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