clear all clf X=-5:0.05:5; [x y]=meshgrid(X,X); var=0.5; miu=rand(3,2); figure(1) hold on plot(0,0,'r*'); for i=1:3 rb=gauss2d(x,y,miu(i,:),var); d(i)=gauss2d(0,0,miu(i,:),var); contour(x,y,rb,5) line([0 miu(i,1)],[0 miu(i,2)]) s=sprintf('%.2f',d(i)); text(miu(i,1), miu(i,2),s) end hold off