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