% Sphere to stand in for the retina
[X,Y,Z] = sphere(60);

X = mhalf(X,'d')+1;
Y = mhalf(Y,'d')+1;
Z = mhalf(Z,'d');


% Possible gaussian shapes and sizes
mx=[0 0]';
mx1=[2.5 2.5]';
mx2=[1 1]';
mx3=[0.4 0.4]';

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];
Cx5=[0.015 0; 0 0.015];

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;

%x3=1:0.033:2;
%x4=1:0.033:2;

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

noise1 = 2;
noise2 = 2;
noise3 = 2;

mx3= mx3 + (rand(size(mx3))/10)*noise1
mx= mx + (rand(size(mx))/10)*0
Cx= Cx + (rand(size(Cx))/1000)*noise2
Cx= Cx.*eye(size(Cx))
Cx5= Cx5 + (rand(size(Cx5))/1000)*noise3
Cx5= Cx5.*eye(size(Cx5))

g = multiGauss(x3,x4,mx3,Cx5);
g2 = multiGauss(x3,x4,mx,Cx);
%g = g+1;

%x3 = x3+1;
%x4 = x4+1;
%g = (g+1)/100;

g = g/1500;
g2 = g2*15;

g = g+g2;




%sph = surf(X,Y,Z,'EdgeColor','none')
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