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