function ydot=f1(t,y)
Um=22e3/sqrt(3)*sqrt(2);
Omega=314.16;
R1=7.260000000000000;
Lh=3.144777098309160e+03;
ydot=(Um*sin(Omega*t)-y*R1)/Lh;