> restart;

> eq1:=v2(t)=i2(t)*2.5e3;

> eq2:=v3(t)=i3(t)*1.5e3;

> eq3:=v4(t)=iL2(t)*1.5e3;

> eq4:=vL1(t)=25e-6*diff(iL1(t),t);

> eq5:=vL2(t)=25e-6*diff(iL2(t),t);

> eq6:=i2(t)=-iL1(t);

> eq7:=iL1(t)=i3(t)+iL2(t);

> eq8:=v2(t)=vL1(t)+v3(t);

> eq9:=v3(t)=vL2(t)+v4(t);

> soln:=dsolve({eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9,iL1(0)=3.3e-3,iL2(0)=6.6e-3},{vL1(t),vL2(t),v2(t),v3(t),v4(t),iL1(t),iL2(t),i2(t),i3(t)},method=laplace);

> evalf(soln);

>