%Guion de Octave para reolver un circuito RLC con condiciones iniciales Vc(0)<>0
% e Ic(0)=0;

c=0.1e-6; %capacitor de 1uF
l=10e-3;  %ionductancia de 1mH
r=10e3;   %Resistencia de 1K
vc0=10;   %Voltaje inicial en el capacitor
il0=1;    %Corriente inicial en la inductancia
A=[0 1/c;-1/l -r/l];
[X,Lambda]=eig(A);
cc=X\[vc0;il0];
fprintf("vc(t)=(%f)(%f) exp(%f*t)+(%f)(%f) exp(%f*t)\n",cc(1),X(1,1),Lambda(1,1),...
                                                  cc(2),X(1,2),Lambda(2,2));
fprintf("il(t)=(%f)(%f)exp(%f*t)+(%f)(%f) exp(%f*t)\n",cc(1),X(2,1),Lambda(1,1),...
                                                  cc(2),X(2,2),Lambda(2,2)); 
tau=max(abs(1/Lambda(1,1)),abs(1/Lambda(2,2)));

t=linspace(0,5*tau,100);
vc=cc(1)*X(1,1)*exp(Lambda(1,1)*t)+cc(2)*X(1,2)*exp(Lambda(2,2)*t);
il=cc(1)*X(2,1)*exp(Lambda(1,1)*t)+cc(2)*X(2,2)*exp(Lambda(2,2)*t);
plot(t,vc,t,il);