% script: hw2p1test.m
% test file for function DoubleIntEst.m

clc

n=[100 1000 5000 10000];
est=zeros(1:length(n));

disp('   n    Integral Estimate');
disp('-------------------------');

for i=1:length(n)
  est(i)=DoubleIntEst(n(i));
  disp(sprintf('%5.0f     %6.5f',n(i),est(i)))
end

  