/*
*  R. Turco
*  v. 1.3
*  2010
*/

/*
   Serie Inversa dei primoriali (p# = p1*p2*p3*...*pn)
   SerRevPrimorial is convergent at S=1.704761... 
   n=7
*/


SerRevPrimorial(n) = local(S=0,P=1);{
 
 forprime(i=2,n,
     P = ProdPrime(i);
     S = S + (1/P) * 1.0;
 );
 return(S);
}

ProdPrime(v)= local(p=1);{
 
 forprime(i=2,v,
          p = p * i;
 );
 return(p);
}

/*
   Serie Inversa dei fattoriali (n! = n*n-1*n-2...)
   SerRevFactorial is convergent at S=2.7182818284590452353602874712... 
   n=7
*/

SerRevFactorial(n) = local(S=0,P=1);{
 
 for(i=0,n,
     P = i!;
     S = S + (1/P) * 1.0;
 );
 return(S);
}