%I #9 Jan 19 2015 03:42:50
%S 7,1,8,8,0,8,2,7,2,8,9,7,6,0,3,2,7,0,2,0,8,2,1,9,4,3,4,5,1,2,4,7,5,8,
%T 7,1,8,5,5,9,3,0,1,7,6,3,9,6,8,4,3,7,1,6,2,4,1,0,0,3,5,6,9,9,4,2,3,4,
%U 3,8,4,0,3,8,7,1,8,6,3,8,1,2,8,5,6,9,8,4,1,9,0,6,2,9,9,0,7,8,8,2,4,7,8,6,0
%N Decimal expansion of Pi! (or Gamma(Pi+1)).
%C Some curiosities:
%C Gamma(Pi) = 2.2880377953400324179595889...
%C c = Pi! - Gamma(Pi)= (Pi-1)*Gamma(Pi) = 4.900044933636000284122605...
%C c - 4.9 = 0.0000449336360002841226054360645...
%C d = 20+Pi-exp(Pi) = 0.00090002081052423273355701533095550...
%C d/(c-4.9) = 20.030001812418246261112320495414095...
%C d/(c-4.9) - 20 = 0.0300018124182462611123204954140952...
%C (d/(c-4.9) - 20)^2 = 0.000900108748379635566119435708756...
%C (Pi!)^2 = 51.66853331860352960055302...
%C 3*(Pi!)^2 = 155.005599955810588801659...
%C 3*(Pi!)^2/Pi^3 = 4.999168444313570423878747...
%e 7.188082728976032702082194345124758718559301763968437162410035699423438403...
%t RealDigits[N[Pi!,200]] (* _Vladimir Joseph Stephan Orlovsky_, May 27 2010 *)
%o (PARI) \p 110; x=gamma(Pi+1);y=x/10^ceil(log(x)/log(10)); for(n=1,110,z=y*10;w=floor(z);print1(w,",");y=z-w)
%Y Cf. A000796, A111197.
%K cons,nonn
%O 1,1
%A _Gerald McGarvey_, Nov 01 2005