% GEP: elektricno vezje

A = [6 -4; -4 5];

B = [-9 5; 5 -6];

% lastne vrednosti in vektorji za GEP
[V,D] = eig(A,B)

% partikularna resitev, ko t --> inf oz I'=0 (0 = B Ip + [0; 6])

Ip = -B\[0;6]

% zacetni pogoj: I(0) = V c + Ip = 0
c = -V\(Ip)

% koncna resitev: I = V [e^(D*t)] c + Ip

I = @(t) c(1)*V(:,1)*exp(D(1,1)*t) + c(2)*V(:,2)*exp(D(2,2)*t) + Ip;

figure
t = linspace(0,10);
plot(t, I(t))