clear
mu=[0?0];
sigma=[1?0;0?1];
x=-5:0.1:5;
y=-5:0.1:5;
for?i=1:length（x（1:））
????for?j=1:length（y（1:））
????????p（ij）=mvnpdf（[x（i）?y（j）]musigma）;
????end?
end
subplot（121）
mesh（xyp）;
title（‘二維正態概率密度函數‘）
u=-5: