function ydot=f1(t,y)
Um=13/100*400e3/sqrt(3)*sqrt(2);
Omega=314.16;
Rk=3.2;
Ls=0.26464;
ydot=(Um*sin(Omega*t)-y(1)*Rk)/Ls;