matlab解八元方程组,求解八元方程的问题

程序如下：clcclear

syms m n x y Am Bm Cm Dm En Fn Gn Hn

v=0.167

E=3E10

a=3

b=8.8

h=0.5

A=m*pi/a

B=n*pi/b/2

D=E*h^3/12/(1-v^2)

q=18.7*10^3

Rmn=32/q/pi^3/D*sin(m*pi/2)^2*sin(0.7*n*pi/2)*sin(0.2*n*pi/4)/m/n/(A^2+B^2)^2

a1=2*A*cosh(A*b)*sin(n*pi/2)/b/(A^2+B^2)

a2=2*A*sinh(A*b)*sin(n*pi/2)/b/(A^2+B^2)

a3=-2*(2*A^2*B/b/(A^2+B^2)^2+A*sinh(A*b)*sin(n*pi/2)*(A^2-B^2)/b/(A^2+B^2)^2-A^2*cosh(A*b)*sin(n*pi/2)/(A^2+B^2))

a4=-2*(A*cosh(A*b)*sin(n*pi/2)*(A^2-B^2)/b/(A^2+B^2)^2-A^2*sinh(A*b)*sin(n*pi/2)/(A^2+B^2))

b1=-2*A*sinh(B*a)*cos(m*pi)/a/(A^2+B^2)

b2=-2*A*(cosh(B*a)*cos(m*pi)-1)/a/(A^2+B^2)

b3=-2*A*B*sinh(B*a)*cos(m*pi)/(A^2+B^2)+4*A*B^2*(cosh(B*a)*cos(m*pi)-1)/A/(A^2+B^2)^2

b4=4*A*B^2*sinh(B*a)*cos(m*pi)/A/(A^2+B^2)^2-2*A*B*cosh(B*a)*cos(m*pi)/(A^2+B^2)

Z1=(symsum((B^2*(En*(v-1)*b1+Fn*(v-1)*b2+Gn*(2*v*b2+(v-1)*b3)+Hn*(2*v*b1+(v-1)*b4))-(B^2+v*A^2)*Rmn)*sin(n*pi/2),n,1,10)-symsum((B^2*(En*(v-1)*b1+Fn*(v-1)*b2+Gn*(2*v*b2+(v-1)*b3)+Hn*(2*v*b1+(v-1)*b4))-(B^2+v*A^2)*Rmn)*sin(n*pi/2)*(-1)^n,n,1,10))/2

Z2=symsum((A*(Am*a1+Bm*a2+Cm*a3+Dm*a4)+A*Rmn),m,1,10)

Z3=symsum((A*(Am*a1+Bm*a2+Cm*a3+Dm*a4)+A*Rmn)*cos(m*pi),m,1,10)

Z4=(symsum((B*(Am*a1+Bm*a2+Cm*a3+Dm*a4+Rmn)),n,1,10)-symsum((B*(Am*a1+Bm*a2+Cm*a3+Dm*a4+Rmn))*(-1)^n,n,1,10))/2

s1=solve(Fn,En*sinh(B*a)+Fn*cosh(B*a)+Gn*B*a*sinh(B*a)+Hn*B*a*cosh(B*a),...

Bm,Am*(v-1)*cosh(A*b)+Bm*(v-1)*sinh(A*b)+Cm*((v+1)*sinh(A*b)+(v-1)*A*b*cosh(A*b))+Dm*((v+1)*cosh(A*b)+(v-1)*A*b*sinh(A*b)),...

A^2*(Am*(1-v)*sinh(A*b)+Bm*(1-v)*cosh(A*b)+Cm*(2*cosh(A*b)+(1-v)*A*b*sinh(A*b))+Dm*(2*sinh(A*b)+(1-v)*A*b*cosh(A*b)))+Z1,...

B*(En+Hn)+Z2,...

B*(En*cosh(B*a)+Fn*sinh(B*a)+Gn*(sinh(B*a)+B*a*cosh(B*a))+Hn*(cosh(B*a)+B*a*sinh(B*a)))+Z3,...

A*(Am+Dm)+Z4,...

Am,Bm,Cm,Dm,En,Fn,Gn,Hn)

Am=eval(s1.Am)

Bm=eval(s1.Bm)

Cm=eval(s1.Cm)

Dm=eval(s1.Dm)

En=eval(s1.En)

Fn=eval(s1.Fn)

Gn=eval(s1.Gn)

Hn=eval(s1.Hn)