e771.m

matlab算法集 matlab算法集

%------------------------------------------------------------------
% Example 7.7.1: End Point Singularity                             
%------------------------------------------------------------------
   clc
   clear
   N = 10;
   a = 0;
   b = 1;
   f = inline ('1/sqrt(x)','x');
   fprintf ('Example 7.7.1: End Point Singularity\n');		
   fprintf ('\n   n      I(0,1)    percent error\n');
   fprintf ('---------------------------------\n');
   for i = 0 : N
      n = 2^i;
      y = midpoint (a,b,n,f);
      fprintf ('%4i %12.7f   %12.7f\n',n,y,100*(y-2)/2);
   end
   fprintf ('---------------------------------\n');
   wait
%------------------------------------------------------------------