function ShowVolume()
global count

% Integrate
% f(x,y) = exp(-15*(x-3)^2 -12*(y-4)^2)+exp(-10*(x-4)^2 -13*(y-3)^2);
% over the triangle with vertices (1,10), (2,-4), and (6,6)...
u = [1 2 6];
v = [10 -4 6];
fVals = [MyF(u(1),v(1)) MyF(u(2),v(2)) MyF(u(3),v(3))];
clc
% Set maximum number of levels at
L = 10;

fprintf('You:\n\n')
disp(' tol Volume Number f-evals ')
disp('-------------------------------------------')
for tol = [.01 .001 .0001 .00001 .000001]
count = 0;
V = Volume(u,v,@MyF,fVals,tol,L);
fprintf('%10.1e %12.10f %5d \n',tol,V,count)
end

fprintf('\n\nCVL:\n\n')
disp(' tol Volume Number f-evals ')
disp('-------------------------------------------')
for tol = [.01 .001 .0001 .00001 .000001]
count = 0;
V = VolumeCVL(u,v,@MyF,fVals,tol,L);
fprintf('%10.1e %12.10f %5d \n',tol,V,count)
end
fprintf('\n\nMaximum number of prismoids = 3^%1d = %1d\n',L,3^L)
fprintf('This is approximately the maximum number of f-evaluations\n')

function z = MyF(x,y)
global count
count=count+1;
z = exp(-15*(x-3)^2 -12*(y-4)^2)+exp(-10*(x-4)^2 -13*(y-3)^2);