mx=[0 0]';
mx3=[0.4 0.4]';

Cx=[1 0; 0 1];
Cx5=[0.015 0; 0 0.015];

x3=-1:0.066:1;
x4=-1:0.066:1;

noise = [2 2 2];

G = mixtureGaussians(x3,x4,mx,mx3,Cx,Cx5,noise);

gss = surf(x3,x4,G,'EdgeColor','none')

%colormap(gray);
%view(0,-90)
view(0,30)
%rotate(img,[0,90],90)
%viewmtx(0,90,60)
%axis off
axis equal
axis square
%axis equal
axis vis3d
%axis image
axis fill