边界条件微分方程MATLAB

clc
clear
syms X1 X2 w t Y1 Y2 Z1 Z2 M m0 Cx Kx ly r Cy Ky J J0 lx lox loy

[X1X2Y1Y2Z1Z2]=solve（‘（M+m0）*（-sin（w*t）*w^2*X1-cos（w*t）*w^2*X2）+Cx*（cos（w*t）*w*X1-sin（w*t）*w*X2）+Kx*x+ly*Cx*（cos（w*t）*w*Z1-sin（w*t）*w*Z2）+ly*Kx*z=m0*r*w^2*cos（w*t）‘...
    ‘（M+m0）*（-sin（w*t）*w^2*Y1-cos（w*t）*w^2*Y2）+Cy*（cos（w*t）*w*Y1-sin（w*t）*w*Y2）+Ky*y-ly*Cy*（cos（w*t）*w*Z1-sin（w*t）*w*Z2）-ly*Ky*z=m0*r*w^2*sin（w*t）‘...
    ‘（J+J0）*（-sin（w*t）*w^2*Z1-cos（w*t）*w^2*