% Reseni prubehu proudu naprazdno numerickou metodou
tspan = 0 : 0.0001 : 0.1;
y0 = 0;
[t,y] = ode23(@(t,y)f1c(t,y, Omega, Um, Rk1, Lh), tspan, y0,[]);

plot(t,y);
grid on;
