matlab如何计算多次方程组,matlab多元多次方程组求解咨询

%%代码：m文件代码：

function eq=myfun(Vd0,Vd1,Vd2,Vd3,Vd4,Vc2,Vc3,Vc4,Vc5,Vb1,Vb2,Vb3,Vb4,Pd1,Pd2,Pd3,Pd4,Pc1,Pc2,Pc3,Pc4)

syms Vd0 Vd1 Vd2 Vd3 Vd4 Vc2 Vc3 Vc4 Vc5 Vb1 Vb2 Vb3 Vb4 Pd1 Pd2 Pd3 Pd4 Pc1 Pc2 Pc3 Pc4

%%global Vd0 Vd1 Vd2 Vd3 Vd4 Vc2 Vc3 Vc4 Vc5 Vb1 Vb2 Vb3 Vb4 Pd1 Pd2 Pd3 Pd4 Pc1 Pc2 Pc3 Pc4

g=9.8

de=0.25

Cd=0.9

Cc=0.2

P0=100000

Ab=48

Ad=60

Ac=60

H=1.5

Kd=1.5

Kc=2.5

Kb=0.2

N=4

Vc1=0

Vc0=0

Vb0=0

Pd0=101000

Pc5=101000

eq(1)=Vc1+(Ab/Ac)^2*Vb1-Vc2;

eq(2)=Vc2+(Ab/Ac)^2*Vb2-Vc3;

eq(3)=Vc3+(Ab/Ac)^2*Vb3-Vc4;

eq(4)=Vc4+(Ab/Ac)^2*Vb4-Vc5;

eq(5)=Vd0-(Ab/Ad)^2*Vb0-Vd1;

eq(6)=Vd1-(Ab/Ad)^2*Vb1