function ydot=f1(t,y)
Uk=10;
Um=Uk/100*400e3/sqrt(3)*sqrt(2);
Omega=314.16;
Rk=4.57174;
Ls=0.14478;
ydot=(Um*sin(Omega*t)-y*Rk)/Ls;