e751.m

matlab算法集 matlab算法集

%------------------------------------------------------------------
% Example 7.5.1: Romberg Integration 
%------------------------------------------------------------------
   clc
   clear
   f = inline ('exp(x)/(exp(1) - 1)','x');
   m = 5;
   a = 0;
   b = 1;
   tol = 0;
   fprintf ('Example 7.5.1: Romberg Integration\n'); 
   [y,j] = romberg (a,b,tol,m,f,1);
   fprintf ('\n');
   show ('y',y)
   show ('Extrapolation level',j)
%------------------------------------------------------------------