程序上的问题（求助）

qianen

syms s
mb=[3577,0;0,3577];
cb=[140 1.24;1.24,140]*(10e4);
kb=[270.48 -1.28;1.28 270.48]*(10e8);
k=[196 1.12;1.12 196]*(10e7);
c=[14.5 2.31;2.31 14.5]*(10e5);
a=k+s*c;
b=a+kb+s*cb+s^2*mb;
c=kb+s*cb+s^2*mb;
d=a/b*c;
K=vpa(d,10);
m1=2.94e3;
m27=2.94e3;
mi=5.88e3;
li=0.4;
EJ=1.3514e11;
L=[li,li,li,li,li,li,li,li,li,li,li,li,li,li,li,li,li,li,li,li,li,li,li,li,li,li,0];
M=[m1,mi,mi,mi,mi,mi,mi,mi,mi,mi,mi,mi,mi,mi,mi,mi,mi,mi,mi,mi,mi,mi,mi,mi,mi,mi,m27];
KXX=[K(1,1),0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,K(1,1)];
KXY=[K(1,2),0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,K(1,2)];
KYX=[K(2,1),0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,K(2,1)];
KYY=[K(2,2),0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,K(2,2)];
for i=1:27
    u11=[1,L(i),0,0;
        0,1,0,0;
        0,0,1,L(i);
        0,0,0,1];
    u21=[L(i)^2/(2*EJ),L(i)^3/(6*EJ),0,0;
        L(i)/(2*EJ),L(i)^2/(2*EJ),0,0;
        0,0,L(i)^2/(2*EJ),L(i)^3/(6*EJ);
        0,0,L(i)/(2*EJ),L(i)^2/(2*EJ)];
    u12=[-L(i)*(M(i)*s^2+KXX(i)),0,-L(i)*KXY(i),0;
        -(M(i)*s^2+KXX(i)),0,-KXY(i),0;
        -L(i)*KYX(i),0,-L(i)*(M(i)*s^2+KYY(i)),0;
        -KYX(i),0,-(M(i)*s^2+KYY(i)),0];
    u22=[1-(M(i)*s^2+KXX(i))*L(i)^3/(6*EJ),L(i),-KXY(i)*L(i)^3/(6*EJ),0;
        -(M(i)*s^2+KXX(i))*L(i)^2/(2*EJ),1,-KXY(i)*L(i)^2/(2*EJ),0;
        -KYX(i)*L(i)^3/(6*EJ),0,1-(M(i)*s^2+KYY(i))*L(i)^3/(6*EJ),L(i);
        -KYX(i)*L(i)^2/(2*EJ),0,-(M(i)*s^2+KYY(i))*L(i)^2/(2*EJ),1];
    T(:,:,i)=[u11 u12;u21 u22];
end
H=T(:,:,1);
for n=2:27;
    H=T(:,:,n)*H;
end
H
请问下这个程序为什么错误？

ChaChing

建议给出出错提示

qianen

H=T(:,:,n)*H;
出错就在这个语句

F_Inzaghi

H=T(:,:,n).*H;