% Written by Min-Hua Huang
% June 9, 2011
% Table 1 of "Indefinite Integration of the Gamma Integral".
% Total Number of Significant Digits by Varying $C$.
clear; clc; format long; digits(700);
fprintf('Total number of significant digits and the considerably larger number C.\n');
fprintf('Target of analysis: Gamma(s), s=0.1,0.2,...,1.\n');
record=ones(10,10); kk=1; kkkk=1; precn=700;
for i=0.1:0.1:1
   answer=vpa(gamma(sym(i)),precn); 
   kkk=1;
   fprintf('Gamma(%2.1f), C=10,20,...,100.\nError=',i);
   for j=10:10:100
       result=vpa(sym(j)^(i)*exp(sym(-j))*h1(i-1,-j,precn));
       record(kk,kkk)=vpa(result-answer,precn);
       sdp=floor(-log10(abs(record(kk,kkk))));
       fprintf(' %2u',sdp);
       kkk=kkk+1;
       kkkk=kkkk+1;
   end
   fprintf('.\n');   
   kk=kk+1; 
end    
fprintf('\n');