程序如下:%输入系数值
h1=300;h2=10;h3=10;
q=50;l=3000;l0=2400;
syms A1 A2 B1 B2 C1 C2 C3 D1 D2
syms x y m a1 a2
a1=m*pi/l;a2=m*pi/l0;
Rf1=A1*(-sinh(a1*h1))+B1*cosh(a1*h1)+C1*h1*sinh(a1*h1)+D1*(-h1)*cosh(a1*h1)-(4*q*(sin(m*pi/2))^2)/(m*pi*a1^2);
Rf2=A1*cosh(a1*h1)-B1*sinh(a1*h1)-C1*(sinh(a1*h1)/a1+h1*cosh(a1*h1))+D1*(cosh(a1*h1)/a1+h1*sinh(a1*h1));
Rf3=-a1^2*sin(a1*x)*B1-B2*(2*sin(a2*(x-(l-l0)/2))-a2^2*sin(a2*(x-(l-l0)/2))*(x-(l-l0)/2)*(x-(l+l0)/2)+2*a2*cos(a2*(x-(l-l0)/2))*(2*x-l));
Rf4=-a1^2*cos(a1*x)*(A1+D1/a1)-(sin(a2*(x-(l-l0)/2))*(l-2*x)-a2*cos(a2*(x-(l-l0)/2))*(x-(l-l0)/2)*(x-(l+l0)/2))*(A2*a2+D2);
%关于Aj积分
RF1=int(Rf1*diff(Rf1,A1),x,0,l)+int(Rf2*diff(Rf2,A1),x,0,l)+int(Rf3*diff(Rf3,A1),x,0,l)+int(Rf4*diff(Rf4,A1),x,0,l)
我的版本是64位R2011a