function A = simRetina(samples)

% Possible gaussian shapes and sizes

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;

s = size(x3);
s
A = zeros(x3*x4,samples);

size(A)

for i = 1:size
    mx3= mx3 + (rand(size(mx3))/10)*2
    mx= mx + (rand(size(mx))/10)*2
    Cx= Cx + (rand(size(Cx))/1000)*2
    Cx= Cx.*eye(size(Cx))
    Cx5= Cx5 + (rand(size(Cx5))/1000)*2
    Cx5= Cx5.*eye(size(Cx5))

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

    g = g/1500;
    g2 = g2*15;
    g = g+g2;
    
    g = g(:);
    A(:,i) = g;
end