%------------------------------------------------------------------
% Example 7.8.1: Two-Dimensional Integral                          
%------------------------------------------------------------------
   clc
   clear
   global n 
   N = 6;
   fprintf ('Example 7.8.1: Two-Dimensional Integral\n');
   fprintf ('\n  n       I(D)    evaluations\n');
   fprintf ('-----------------------------\n'); 
   for i = 1 : N
      n = 10*i;
      z = simpson (-1,1,n,'fung781');
      fprintf ('%3i %12.7f %8g\n',n,z,n*n);
   end
   fprintf ('-----------------------------\n'); 
   wait
%------------------------------------------------------------------