x=-3:.1:3;
y=x;
x=repmat(x',1,61);
y=repmat(y,61,1);
Z=normpdf(x).*normpdf(y);
surf(x,y,Z)