clear all; close all;
pkg load symbolic

syms s t V V0 deni I0
v0=2.5; i0=-1; L=1; C=1; R1=2; R2=1;

% η Sol είναι η τάση V από κομβική ανάλυση
V0 = (V+L*i0)/(s+1)
eq = (V-5/(s+2))/R1 + s*(V-(v0/s+0.5*V0)) + (V+L*i0)/(s+1)
Sol = solve(eq,V)

% deni είναι ο παρονομαστής και Sdeni οι ρίζες του
deni = 2*s^3+6*s^2+7*s+6
Sdeni = solve(deni,s)

% το ζητούμενο ρεύμα I0 και το i0(t) στη μεταβλητή ii
I0 = (Sol+L*i0)/(s+1)
ii = ilaplace(I0,s,t)
