Matlab - Parámetros de 2 puertos: Resolución de sistema de ecuaciones.

clc
clear all;
syms  U3 U4 Ii Io Vi Vo Ib Vce
RB1=22e3;
RB2=68e3;
RC=3.3e3;
RE1=0.47e3;
hie=5e3;
hoe=11e-6;
hfe=200;
hre=2.4e-4;
EC1=solve(((1/RB1)+(1/RB2)+(1/hie))*Vi-(1/hie)*U3==Ii,...
          (hoe+(1/RC))*Vo-hoe*U4==Io-hfe*Ib,...
          (-1/hie)*Vi+(1/hie)*U3==-Ib,...
          -hoe*Vo+((1/RE1)+hoe)*U4==Ib+hfe*Ib,...
          U3-U4==hre*Vce,...
          Vce==Vo-U4,...
	U3,U4,Ib,Vce,Ii,Vo);
gi=subs(EC1.Ii,Io,0)/Vi
gf=subs(EC1.Vo,Io,0)/Vi
gr=subs(EC1.Ii,Vi,0)/Io
go=subs(EC1.Vo,Vi,0)/Io