% Reseni prubehu proudu naprazdno numerickou metodou
f1cc = @(t,y) [(Un1/sqrt(3)*sqrt(2)*sin(omega*t)-y*Rk1)/Lh];
[t,y] = ode23(f1cc, 0 : 0.0001 : 0.1, 0, []);

plot(t,y);
grid on;
