u=[95.800089.600080.0000129.4000]fangcha=[3.2889-0.9778-0.3333-0.0222;-0.97780.48890.2222-0.2667;-0.33330.22221.1111-0.3333;-0.0222-0.2667-0.33330.7111]c=inv(fangcha)%求方差的逆矩阵symsx1x2x3x4d;a=[(x1-u(1,1))(x2-u(1,2))

　　u=[95.800089.600080.0000129.4000]

　　fangcha=[3.2889-0.9778-0.3333-0.0222;-0.97780.48890.2222-0.2667;-0.33330.22221.1111-0.3333;-0.0222-0.2667-0.33330.7111]

　　c=inv(fangcha)%求方差的逆矩阵

　　symsx1x2x3x4d;

　　a=[(x1-u(1,1))(x2-u(1,2))(x3-u(1,3))(x4-u(1,3))];

　　b=[(x1-u(1,1));(x2-u(1,2));(x3-u(1,3));(x4-u(1,3))];

　　d=(2*pi)^-4/2*(abs(fangcha)^-12)*exp{-12*a*c*b}

　　Errorusing==>mldivide

　　Matrixdimensionsmustagree.

　　Errorin==>Untitled00at10

　　d=(2*pi)^-4/2*(abs(fangcha)^-12)*exp{-12*a*c*b}