mx=[0 0]';
mx1=[2.5 2.5]';
mx2=[1 1]';

Cx=[1 0; 0 1];
Cx1=[0.1 0; 0 0.1];
Cx2=[2 0; 0 2];
Cx3=[3 0; 0 3];
Cx4=[0.001 0; 0 0.001];

x1=-3:0.1:3;
x2=-3:0.1:3;

%x3 = -10:0.1:10;
%x4 = -10:0.1:10;
x3=-3:0.1:3;
x4=-3:0.1:3;

g = multiGauss(x3,x4,mx,Cx);
g1 = multiGauss(x3,x4,mx,Cx1);
g2 = multiGauss(x3,x4,mx,Cx2);
g3 = multiGauss(x3,x4,mx,Cx3);
g4 = multiGauss(x3,x4,mx,Cx4);
g5 = multiGauss(x3,x4,mx1,Cx4);
g6 = multiGauss(x3,x4,mx2,Cx1);

% Adjust the relative magnitude of the Gaussians
g7 = g6+g3*1000;

img = surf(x3,x4,g7,'EdgeColor','none')
%colormap(gray);
view(0,-90)
%rotate(img,[0,90],90)
%viewmtx(0,90,60)
axis off
%axis equal
axis square
%axis equal
axis vis3d
%axis image
axis fill