archive/share/unknown/fact.dft

   1
   2 /*  A neat example of how to exponentiate a series.
   3     Although, TAYLOR now has what follows built in.
   4     So this is now for reference only. */
   5
   6 C[0]:1$
   7 C[N]:=RATSIMP(SUM((-1)^(J+1)*S[J+1]*C[N-1-J],J,0,N-1)/N)$
   8 S[1]:%GAMMA$
   9 S[N]:=ZETA(N)$
  10 DEFTAYLOR(Z!,SUM(C[N]*Z^N,N,0,INF))$
  11 DEFTAYLOR(GAMMA(Z),SUM(C[N]*Z^(N-1),N,0,INF))$
  12 NUMERVAL(%GAMMA,.577215665)$
  13 DECLARE(%GAMMA,CONSTANT)$