function Disparity = construct_volumes(General, Volumes)

% Possible gaussian shapes and sizes

noise = [2 2 2];
samples = 20;

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

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

x3=-1:0.020:1;
x4=-1:0.020: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

Disparity = A.^(-1);


% Verification statements. Remove when working.
figure(1)
surf(reshape(Disparity(:,1),sqrt(size(Disparity,1)), sqrt(size(Disparity,1))), '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