e772.m

matlab算法集 matlab算法集

%------------------------------------------------------------------
% Example 7.7.2: Infinite Limit of Integration                     
%------------------------------------------------------------------
   clc
   clear
   N = 9;
   a = 0;
   b = 1;
   f = inline ('4/(pi*(4 + x.^2))','x');
   g = inline ('4/(pi*(1 + 4*z.^2))','z');
   fprintf ('Example 7.7.2: Infinite Limit of Integration\n');
   fprintf ('\n n      I(0,1)     percent error\n');
   fprintf ('--------------------------------\n'); 
   for n = 1 : N
      y = gausquad(a,b,n,f) + gausquad(a,b,n,g); 
      fprintf ('%2i %12.7f   %12.7f\n',n,y,100*(y-1));
   end
   fprintf ('--------------------------------\n'); 
   wait
%------------------------------------------------------------------