function ydot=f1c(t,y, omega, Un1, Rk1, Lh)
Um=Un1/sqrt(3)*sqrt(2);
ydot=(Um*sin(omega*t)-y*Rk1)/Lh;
end
