matlab求二维概率密度二维概率分布

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: