indexing must appear last in an index expression出错求助

ssaw0920

syms Phy  phy1 phy2 phy3  xi phyi xi1 xi2 xi3 phyi1 phyi2 phyi3 L a b pho1 xi1t xi2t xi3t phyi1t phyi2t phyi3t ...
    xi1tt xi2tt xi3tt phyi1tt phyi2tt phyi3tt phy1t phy2t phy3t phy1tt phy2tt phy3tt...
    x11 x12 x13 phy11 phy12 phy13 x11t x12t x13t phy11t phy12t phy13t x11tt x12tt x13tt phy11tt phy12tt phy13tt...
    x21 x22 x23 phy21 phy22 phy23 x21t x22t x23t phy21t phy22t phy23t x21tt x22tt x23tt phy21tt phy22tt phy23tt...
    x31 x32 x33 phy31 phy32 phy33 x31t x32t x33t phy31t phy32t phy33t x31tt x32tt x33tt phy31tt phy32tt phy33tt
H1 =[1/2*L*a*b*pho1*(-2*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x11+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x12+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x13)*sin(conj(phy3))*x11+4*conj(conj(D(conj)(0))*x13)*cos(conj(phy2))*x11-4*conj(conj(D(conj)(0))*x13)*cos(conj(phy2))*x12+4*conj(conj(D(conj)(0))*x13)*sin(conj(phy2))*x13-2*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x21+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x22+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x23)*sin(conj(phy3))*x21-2*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x21+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x22+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x23)*cos(conj(phy3))*x22+2*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x21+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x22+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x23)*cos(conj(phy3))*x23+2*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x22+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x23)*cos(conj(phy3))*x21-2*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x22+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x23)*sin(conj(phy3))*x22-2*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x22+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x23)*cos(conj(phy3))*x23+4*conj(conj(D(conj)(0))*x23)*cos(conj(phy2))*x21-4*conj(conj(D(conj)(0))*x23)*cos(conj(phy2))*x22+4*conj(conj(D(conj)(0))*x23)*sin(conj(phy2))*x23-2*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x31+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x32+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x33)*sin(conj(phy3))*x31-2*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x31+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x32+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x33)*cos(conj(phy3))*x32+2*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x31+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x32+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x33)*cos(conj(phy3))*x33+2*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x32+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x33)*cos(conj(phy3))*x31-2*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x32+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x33)*sin(conj(phy3))*x32-2*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x32+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x33)*cos(conj(phy3))*x33+4*conj(conj(D(conj)(0))*x33)*cos(conj(phy2))*x31-4*conj(conj(D(conj)(0))*x33)*cos(conj(phy2))*x32+4*conj(conj(D(conj)(0))*x33)*sin(conj(phy2))*x33+sin(conj(phy3))*L^2*D(conj)(0)-cos(conj(phy3))*a^2*D(conj)(-1)+cos(conj(phy3))*b^2*D(conj)(1)+sin(conj(phy3))*a^2*D(conj)(0)-cos(conj(phy3))*b^2*D(conj)(-1)+sin(conj(phy2))*b^2*D(conj)(0)-2*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x11+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x12+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x13)*cos(conj(phy3))*x12+2*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x11+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x12+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x13)*cos(conj(phy3))*x13+2*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x12+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x13)*cos(conj(phy3))*x11-2*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x12+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x13)*sin(conj(phy3))*x12-2*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x12+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x13)*cos(conj(phy3))*x13),1/12*L*a*b*pho1*(-12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x11+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x12+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x13)*sin(conj(phy3))*x11+24*conj(conj(D(conj)(0))*x13)*cos(conj(phy2))*x11-24*conj(conj(D(conj)(0))*x13)*cos(conj(phy2))*x12+24*conj(conj(D(conj)(0))*x13)*sin(conj(phy2))*x13-12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x21+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x22+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x23)*sin(conj(phy3))*x21-12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x21+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x22+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x23)*cos(conj(phy3))*x22+12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x21+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x22+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x23)*cos(conj(phy3))*x23+12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x22+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x23)*cos(conj(phy3))*x21-12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x22+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x23)*sin(conj(phy3))*x22-12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x22+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x23)*cos(conj(phy3))*x23+24*conj(conj(D(conj)(0))*x23)*cos(conj(phy2))*x21-24*conj(conj(D(conj)(0))*x23)*cos(conj(phy2))*x22+24*conj(conj(D(conj)(0))*x23)*sin(conj(phy2))*x23-12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x31+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x32+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x33)*sin(conj(phy3))*x31-12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x31+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x32+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x33)*cos(conj(phy3))*x32+12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x31+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x32+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x33)*cos(conj(phy3))*x33+12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x32+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x33)*cos(conj(phy3))*x31-12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x32+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x33)*sin(conj(phy3))*x32-12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x32+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x33)*cos(conj(phy3))*x33+24*conj(conj(D(conj)(0))*x33)*cos(conj(phy2))*x31-24*conj(conj(D(conj)(0))*x33)*cos(conj(phy2))*x32+24*conj(conj(D(conj)(0))*x33)*sin(conj(phy2))*x33+6*sin(conj(phy3))*L^2*D(conj)(0)-6*cos(conj(phy3))*a^2*D(conj)(-1)+6*cos(conj(phy3))*b^2*D(conj)(1)+6*sin(conj(phy3))*a^2*D(conj)(0)-6*cos(conj(phy3))*b^2*D(conj)(-1)+6*sin(conj(phy2))*b^2*D(conj)(0)-L^2*sin(conj(phy2))-12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x11+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x12+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x13)*cos(conj(phy3))*x12+12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x11+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x12+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x13)*cos(conj(phy3))*x13+12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x12+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x13)*cos(conj(phy3))*x11-12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x12+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x13)*sin(conj(phy3))*x12-12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x12+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x13)*cos(conj(phy3))*x13),1/12*L*a*b*pho1*(-a^2*sin(conj(phy3))-L^2*sin(conj(phy3))-12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x11+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x12+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x13)*sin(conj(phy3))*x11+24*conj(conj(D(conj)(0))*x13)*cos(conj(phy2))*x11-24*conj(conj(D(conj)(0))*x13)*cos(conj(phy2))*x12+24*conj(conj(D(conj)(0))*x13)*sin(conj(phy2))*x13-12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x21+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x22+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x23)*sin(conj(phy3))*x21-12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x21+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x22+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x23)*cos(conj(phy3))*x22+12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x21+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x22+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x23)*cos(conj(phy3))*x23+12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x22+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x23)*cos(conj(phy3))*x21-12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x22+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x23)*sin(conj(phy3))*x22-12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x22+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x23)*cos(conj(phy3))*x23+24*conj(conj(D(conj)(0))*x23)*cos(conj(phy2))*x21-24*conj(conj(D(conj)(0))*x23)*cos(conj(phy2))*x22+24*conj(conj(D(conj)(0))*x23)*sin(conj(phy2))*x23-12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x31+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x32+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x33)*sin(conj(phy3))*x31-12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x31+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x32+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x33)*cos(conj(phy3))*x32+12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x31+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x32+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x33)*cos(conj(phy3))*x33+12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x32+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x33)*cos(conj(phy3))*x31-12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x32+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x33)*sin(conj(phy3))*x32-12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x32+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x33)*cos(conj(phy3))*x33+24*conj(conj(D(conj)(0))*x33)*cos(conj(phy2))*x31-24*conj(conj(D(conj)(0))*x33)*cos(conj(phy2))*x32+24*conj(conj(D(conj)(0))*x33)*sin(conj(phy2))*x33+6*sin(conj(phy3))*L^2*D(conj)(0)-6*cos(conj(phy3))*a^2*D(conj)(-1)+6*cos(conj(phy3))*b^2*D(conj)(1)+6*sin(conj(phy3))*a^2*D(conj)(0)-6*cos(conj(phy3))*b^2*D(conj)(-1)+6*sin(conj(phy2))*b^2*D(conj)(0)-12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x11+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x12+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x13)*cos(conj(phy3))*x12+12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x11+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x12+2*conj(Re(D(conj)(1))+i*Im(D(conj)(1)))*x13)*cos(conj(phy3))*x13+12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x12+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x13)*cos(conj(phy3))*x11-12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x12+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x13)*sin(conj(phy3))*x12-12*conj(-2*conj(Re(D(conj)(0))+i*Im(D(conj)(0)))*x12+2*conj(Re(D(conj)(-1))+i*Im(D(conj)(-1)))*x13)*cos(conj(phy3))*x13),                                                                                                                                                                            1/6*L*a*b*pho1*(sin(abs(phy3)^2/phy3)*L^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*a^2*D(conj)(-1)+cos(abs(phy3)^2/phy3)*b^2*D(conj)(1)+sin(abs(phy3)^2/phy3)*a^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*b^2*D(conj)(-1)+sin(abs(phy2)^2/phy2)*b^2*D(conj)(0)),                                                                                                                                                                           1/6*L*a*b*pho1*(sin(abs(phy3)^2/phy3)*L^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*a^2*D(conj)(-1)+cos(abs(phy3)^2/phy3)*b^2*D(conj)(1)+sin(abs(phy3)^2/phy3)*a^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*b^2*D(conj)(-1)+sin(abs(phy2)^2/phy2)*b^2*D(conj)(0)),                                                                                                                                                                            1/6*L*a*b*pho1*(sin(abs(phy3)^2/phy3)*L^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*a^2*D(conj)(-1)+cos(abs(phy3)^2/phy3)*b^2*D(conj)(1)+sin(abs(phy3)^2/phy3)*a^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*b^2*D(conj)(-1)+sin(abs(phy2)^2/phy2)*b^2*D(conj)(0)),                                                                                                                                                                             1/6*L*a*b*pho1*(sin(abs(phy3)^2/phy3)*L^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*a^2*D(conj)(-1)+cos(abs(phy3)^2/phy3)*b^2*D(conj)(1)+sin(abs(phy3)^2/phy3)*a^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*b^2*D(conj)(-1)+sin(abs(phy2)^2/phy2)*b^2*D(conj)(0)),                                                                                                                                                                            1/6*L*a*b*pho1*(sin(abs(phy3)^2/phy3)*L^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*a^2*D(conj)(-1)+cos(abs(phy3)^2/phy3)*b^2*D(conj)(1)+sin(abs(phy3)^2/phy3)*a^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*b^2*D(conj)(-1)+sin(abs(phy2)^2/phy2)*b^2*D(conj)(0)),                                                                                                                                                                             1/6*L*a*b*pho1*(sin(abs(phy3)^2/phy3)*L^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*a^2*D(conj)(-1)+cos(abs(phy3)^2/phy3)*b^2*D(conj)(1)+sin(abs(phy3)^2/phy3)*a^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*b^2*D(conj)(-1)+sin(abs(phy2)^2/phy2)*b^2*D(conj)(0)),                                                                                                                                                                              1/6*L*a*b*pho1*(sin(abs(phy3)^2/phy3)*L^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*a^2*D(conj)(-1)+cos(abs(phy3)^2/phy3)*b^2*D(conj)(1)+sin(abs(phy3)^2/phy3)*a^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*b^2*D(conj)(-1)+sin(abs(phy2)^2/phy2)*b^2*D(conj)(0)),                                                                                                                                                                             1/6*L*a*b*pho1*(sin(abs(phy3)^2/phy3)*L^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*a^2*D(conj)(-1)+cos(abs(phy3)^2/phy3)*b^2*D(conj)(1)+sin(abs(phy3)^2/phy3)*a^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*b^2*D(conj)(-1)+sin(abs(phy2)^2/phy2)*b^2*D(conj)(0)),                                                                                                                                                                              1/6*L*a*b*pho1*(sin(abs(phy3)^2/phy3)*L^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*a^2*D(conj)(-1)+cos(abs(phy3)^2/phy3)*b^2*D(conj)(1)+sin(abs(phy3)^2/phy3)*a^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*b^2*D(conj)(-1)+sin(abs(phy2)^2/phy2)*b^2*D(conj)(0)),                                                                                                                                                                               1/6*L*a*b*pho1*(sin(abs(phy3)^2/phy3)*L^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*a^2*D(conj)(-1)+cos(abs(phy3)^2/phy3)*b^2*D(conj)(1)+sin(abs(phy3)^2/phy3)*a^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*b^2*D(conj)(-1)+sin(abs(phy2)^2/phy2)*b^2*D(conj)(0)),                                                                                                                                                                             1/6*L*a*b*pho1*(sin(abs(phy3)^2/phy3)*L^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*a^2*D(conj)(-1)+cos(abs(phy3)^2/phy3)*b^2*D(conj)(1)+sin(abs(phy3)^2/phy3)*a^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*b^2*D(conj)(-1)+sin(abs(phy2)^2/phy2)*b^2*D(conj)(0)),                                                                                                                                                                               1/6*L*a*b*pho1*(sin(abs(phy3)^2/phy3)*L^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*a^2*D(conj)(-1)+cos(abs(phy3)^2/phy3)*b^2*D(conj)(1)+sin(abs(phy3)^2/phy3)*a^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*b^2*D(conj)(-1)+sin(abs(phy2)^2/phy2)*b^2*D(conj)(0)),                                                                                                                                                                             1/6*L*a*b*pho1*(sin(abs(phy3)^2/phy3)*L^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*a^2*D(conj)(-1)+cos(abs(phy3)^2/phy3)*b^2*D(conj)(1)+sin(abs(phy3)^2/phy3)*a^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*b^2*D(conj)(-1)+sin(abs(phy2)^2/phy2)*b^2*D(conj)(0)),                                                                                                                                                                            1/6*L*a*b*pho1*(sin(abs(phy3)^2/phy3)*L^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*a^2*D(conj)(-1)+cos(abs(phy3)^2/phy3)*b^2*D(conj)(1)+sin(abs(phy3)^2/phy3)*a^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*b^2*D(conj)(-1)+sin(abs(phy2)^2/phy2)*b^2*D(conj)(0)),                                                                                                                                                                             1/6*L*a*b*pho1*(sin(abs(phy3)^2/phy3)*L^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*a^2*D(conj)(-1)+cos(abs(phy3)^2/phy3)*b^2*D(conj)(1)+sin(abs(phy3)^2/phy3)*a^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*b^2*D(conj)(-1)+sin(abs(phy2)^2/phy2)*b^2*D(conj)(0)),                                                                                                                                                                              1/6*L*a*b*pho1*(sin(abs(phy3)^2/phy3)*L^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*a^2*D(conj)(-1)+cos(abs(phy3)^2/phy3)*b^2*D(conj)(1)+sin(abs(phy3)^2/phy3)*a^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*b^2*D(conj)(-1)+sin(abs(phy2)^2/phy2)*b^2*D(conj)(0)),                                                                                                                                                                             1/6*L*a*b*pho1*(sin(abs(phy3)^2/phy3)*L^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*a^2*D(conj)(-1)+cos(abs(phy3)^2/phy3)*b^2*D(conj)(1)+sin(abs(phy3)^2/phy3)*a^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*b^2*D(conj)(-1)+sin(abs(phy2)^2/phy2)*b^2*D(conj)(0)),                                                                                                                                                                             1/6*L*a*b*pho1*(sin(abs(phy3)^2/phy3)*L^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*a^2*D(conj)(-1)+cos(abs(phy3)^2/phy3)*b^2*D(conj)(1)+sin(abs(phy3)^2/phy3)*a^2*D(conj)(0)-cos(abs(phy3)^2/phy3)*b^2*D(conj)(-1)+sin(abs(phy2)^2/phy2)*b^2*D(conj)(0));

为什么出现如下错误信息提示：

??? Error: File: H1H2.m Line: 6 Column: 41
()-indexing must appear last in an index expression.

该式子是程序运行出来的，是一个21*1的矩阵

[ 本帖最后由 eight 于 2008-3-21 23:10 编辑 ]

sigma665

你的表达式还不是一般的复杂

eight

除了你本人以外，恐怕只有上帝才能打救你。

ps：你贴的公式既然是运行某程序后得出来的结果，为何会出现错误？难道你运行的不是程序，而是结果？

rencs1

是不是缺少运算符号。。。