% Reseni prubehu proudu nakratko numerickou metodou

% Ls*i' + Rk*i = Uk*sin(omega*t)
f1ki = @(t,y,yp) [Ls*yp + Rk*y - Un1/sqrt(3)*sqrt(2)*Ukp/100*sin(omega*t)];

[t,y] = ode15i(f1ki, 0 : 0.0001 : 0.1, 0, 0);

plot(t,y);
grid on;
