function ydot=f1(t,y)
Un1=110;
Um=Un1/sqrt(3)*sqrt(2);
omega=2*pi*50;
Rk1 =  0.60500;
Lh =  38.709;
ydot=(Um*sin(omega*t)-y*Rk1)/Lh;
end