% Written by Min-Hua Huang
% June 9, 2011
% Test equation (23) of "Indefinite Integration of the Gamma Integral".
clear; clc; format long; error=ones(1,2); kkk=1;
fprintf('Estimating the upper limit of the error resulted from the \na considerably larger C.\n');
fprintf('\nThe taget of analysis: E1(x), the first-order exponential integral.\n\n');
for x=10:10:100
error(kkk,1)=x;
error(kkk,2)=quadv('f1expo',x,1000000,1e-15);
fprintf('When C=%3i, the error is %5.2e.\n',x,error(kkk,2));
kkk=kkk+1;
end    
fprintf('\n');