function A = simRetina(samples, noise)

% 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;

A = zeros((length(x3)*length(x4)),samples);


noise

for i = 1:samples
    mx3= mx3 + (rand(size(mx3))/10)*noise(1);
    mx= mx + (rand(size(mx))/10)*0;
    Cx= Cx + (rand(size(Cx))/1000)*noise(2);
    Cx= Cx.*eye(size(Cx));
    Cx5= Cx5 + (rand(size(Cx5))/1000)*noise(3);
    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