% Reseni prubehu proudu nakratko numerickou metodou
f1kcc = @(t,y) [(Un1/sqrt(3)*sqrt(2)*Ukp/100*sin(omega*t)-y*Rk)/Ls];
[t,y] = ode23(f1kcc, 0 : 0.0001 : 0.1, 0, []);

plot(t,y);
grid on;
