% Reseni prubehu proudu naprazdno numerickou metodou

% Lh*i' + Rk1*i = Um*sin(omega*t)
f1i = @(t,y,yp) [Lh*yp + Rk1*y - Un1/sqrt(3)*sqrt(2)*sin(omega*t)];

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

plot(t,y);
grid on;
