G=667384*^-16(m^3)(kg^(-1))(sek^(-2));  (* Newtons Konstante *)
c=299792458m (sek^(-1)); (* Einsteins Konstante *)

M=5972*^21kg; (* Erdmasse *)
ω=72921*^-9(sek^(-1)); (* Winkelgeschwindigkeit *)

rÄ=6378km; (* Äquatorialer Radius *)
rM=6371km; (* Mittlerer Erdradius *)
rP=6357km; (* Polradius *)
rΦ=(Sqrt[(((Sin[φ]rP))^2)+(((Cos[φ]rÄ))^2)]);  (* Winkelradius *)

gÄ=97803267714*^-10m/(sek^2); (* g am Äquator, gemessen *)
ex=(Sqrt[(((rÄ)^2)-(rP^2))/((rÄ)^2)]); (* Exzentrizität *)
k=193185138639*^-14; (* Koeffizient *)
χ1=gÄ ((1+k ((Sin[φ])^2))/(Sqrt[1-(ex^2) ((Sin[φ])^2)])); χ2=(ω^2)rΦ ((Cos[φ])^2); 
(* Schwerkraft, Fliehkraft *)

km=1000m; m=1; sek=1; kg=1;Plot[{ χ1+χ2,χ1}, {φ, -π/2, +π/2}, 
PlotStyle->{{Blue,Dashed}, Green, Red}, Frame->True]
 (* Erdanziehung in m/(sek^2) *)

τg=(Sqrt[1-2(χ1+χ2) rΦ (c^(-2))]);(* gravitative ZD *)
τs=(Sqrt[1-(((ω rΦ Cos[φ]))^2) (c^(-2))]); (* geschwindigkeitsabhängige ZD*)
Δτ=τg τs;(* Insgesamte Zeitdilatation *)

Plot[{Δτ-1}, {φ,-π/2, π/2}, 
PlotStyle->{Red}, Frame->True]
(* Zeitdifferenz in sek/sek *)