菜鸟用matlab求解超越方程，求大神修改错误程序

tongyaodj

>> clear
k1=inf; k2=inf; kt1= inf; kt2= inf; L=100; E=2.2*10^10; b=0.01; h=0.01; I=b*h^3/12; rou=7860; A=b*h; m0=rou*A*L; j0=m0*L^2;M0=100;J0=100;
K1=k1*L^3/( E*I);
K2=k2*L^3/( E*I);
B1=kt1*L/(E*I);       
B2=kt2*L/(E*I);
syms x a1n a2n b1n b2n Z1 Z2 Z3 Z4 Z5 Z6 f
a1n=K1/(x)^3;
a2n=1/(x)^3*(K2-M0*(x)^4/m0);
b1n=B1/x;
b2n=1/x*(-B2+J0*(x)^4/j0);
Z1=b1n*b2n*(a1n+a2n)+(b1n-b2n);
Z2=a1n*a2n*(b1n-b2n)-(a1n+a2n);
Z3=a1n*a2n+b1n*b2n;
Z4=1-a1n*a2n*b1n*b2n;
Z5=a2n*b2n-a1n*b1n;
Z6=a1n*b2n-a2n*b1n;
f=Z1*(cos(x)*sinh(x)+sin(x)*cosh(x))+Z2*(cos(x)*sinh(x)- sin(x)*cosh(x))-2*Z3*sin(x)*sinh(x)+Z4*(cos(x)*cosh(x)-1)+Z5*(cos(x)*cosh(x)+1)+2*Z6*cos(x)*cosh(x);
ans=solve(f)

Warning: Explicit solution could not be found.
> In solve at 140
  In sym.solve at 49