function?r=houyan（aausigma）
pi=3.1415926;
r=double（（1./（sqrt（2*pi）*det（sigma）^0.5））*exp（-0.5*（aa-u）‘*inv（sigma）*（aa-u