mu=0:0.05:8;
sigma2=0.05:0.01:10;
clf;
[x,y]=meshgrid(mu,sigma2);
z=(1./sqrt(2*pi*y)).*exp(-((4-x).^2)./(2*y));
z1=(1./sqrt(2*pi*y)).*exp(-((5-x).^2)./(2*y));
z2=(1./sqrt(2*pi*y)).*exp(-((3-x).^2)./(2*y));

zall=(z.*z1).*z2;
clf;
mesh(x,y,z);
xlabel('\mu');
ylabel('\sigma^2');
zlabel('z');

pause;

mesh(x,y,zall);
xlabel('\mu');
ylabel('\sigma^2');
zlabel('z');

