function [nn,err] = errMC()

  nn = [10, 25, 50, 75, 100, 250, 500];

n = length(nn); rep = 100;

err = zeros(1,n);

for j = 1:n

    s = 0;
    for k=1:rep
	[v, vAna] = euroMC(nn(j));
	s = s+abs(v-vAna); 
    end
    err(j) = s/rep;
end
