Symbolsk løsning av nodeligninger med Matlab

% Using Kirchoff's current law (KCL) on a source follower configuration
% to find Vout as a function of Vin
% Kjetil Ullaland, 10. september 1996

ligning1='(Vout-Vgs)/Zc+gm*Vgs+Vout/Rl=0';
ligning2='(Vgs-Vout)/Zc+(Vgs-Vin)/Rs=0';
ligning1=subs(ligning1,'1/(j*w*C)','Zc');
ligning2=subs(ligning2,'1/(j*w*C)','Zc');
pretty(ligning1);
pretty(ligning2);


Vgs_solved=simplify(solve(ligning2,'Vgs'))
pretty(Vgs_solved)
ligning3=subs(ligning1,Vgs_solved,'Vgs');
Vout_solved=simplify(solve(ligning3,'Vout'))
pretty(Vout_solved)