function ydot=f1c(t,y, Omega, Um, Rk1, Lh)
ydot=(Um*sin(Omega*t)-y*Rk1)/Lh;
end
