% Prechodne deje relizovane na transformatoru
% a vedeni

% Reseni prubehu proudu nakratko numerickou metodou
pom1=Ls/2+Lh;
pom2=Ls/2+Lh-Lh^2/pom1;
Ukm=Un1/sqrt(3)*sqrt(2)*Ukp/100;
f1kc = @(t,y) [(Ukm*sin(omega*t)-y(1)*Rk/2-Lh*y(2)*Rk/2/pom1)/pom2
               (Lh*Ukm*sin(omega*t)/pom1-Lh*y(1)*Rk/2/pom1-y(2)*Rk/2)/pom2];
[t,y] = ode23(f1kc, 0 : 0.0001 : 0.1, [0 0]);

plot(t,y);
grid on;
