function ydot=f1(t,y)
Un=400;
Um=Un/sqrt(3)*sqrt(2);
Omega=2*pi*50;
Rk1=0.8;
Lh=53.389;
ydot=(Um*sin(Omega*t)-y*Rk1)/Lh;
end
