Quit[]

G=667384*^-16 (m^3)(kg^(-1))(sek^(-2));
c=299792458m (sek^(-1));
r=6378000m;
h=20000000m;
M=5972*^21kg;
ω=72921*^-9(sek^(-1));
f=1023/100MHz;

f1[R_]:=(Sqrt[1-((2 G M)/(R (c^2)))])
f2[n_]:=(Sqrt[1-(((v[[n]])^2)/(c^2))])

v:={(Sqrt[((G M)/(r+h))]),  ω(r+h), ω r, 0};

Δ1:=1-f1[r]/f1[r+h]; N[Δ1]

Δ2:=1-f2[w]/f2[u]; w=3; u=1; N[Δ2]

Δ:=(1+Δ1)(1+Δ2); N[Δ-1]

F=N[{f, f/Δ}, 13]

plot:=Plot[{Δ-1, Δ1, Δ2}, {h,0,2*^8m}, PlotRange->All, Frame->True, PlotStyle->{Green, Blue, Red}, GridLines->{{sol}, {}}];
sol:=NSolve[Δ==1&&h>0, h][[1,1,2]];
Clear[h,u]; u=2; m=1;
"h (stationär)"->sol"m"
plot
w=4; u=1;
"h (orbital)"->sol"m"
plot
w=3; u=1;
"h (orbital+rotation)"->sol"m"
plot

(* Simon Tyran, Vienna. Yukterez, 2015 *)