function ydot=f1(t,y)
Um=220e3/sqrt(3)*sqrt(2);
Omega=314.16;
R1=0.605;
Lh=168.1;
ydot=(Um*sin(Omega*t)-y*R1)/Lh;