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