% Blatt 6 Aufgabe3

%Gegeben:
tau1=0.06;
delta_tau=25;
alpha=-2.5;
k=360/delta_tau;
rho=999.975;
RE=6.370e6;
omega=7.29e-5;

x_1=-1.0e6;
x_2=-4.0e6;
H=200;

display('Winkel:')
theta=[-10:5:25]'

Dtau_x=(tau1*k*sind(k*(theta-alpha)))/RE;
D2tau_x=(tau1*k^2*cosd(k*(theta-alpha)))/RE^2;


M_x1=(x_1*(2*omega*cosd(theta)).^-1).*(RE*D2tau_x+tand(theta).*Dtau_x);
M_x2=(x_2*(2*omega*cosd(theta)).^-1).*(RE*D2tau_x+tand(theta).*Dtau_x);

M_y=-RE*(2*omega*cosd(theta)).^(-1).*Dtau_x;

display('Mittlere Geschwindigkeiten x_1=1000km und x_2=4000km:')
u_x1=M_x1/(H*rho)
u_x2=M_x2/(H*rho)
v=M_y/(H*rho)

M_y=0;
