Could anyone help me with the following issue.
Using the GSL library in R I define the following code:

#########
library(gsl);
S<-function(n)
{ r<-0:n;
ans<-sum(gsl_sf_choose(n,r)*(-1)^r*2^(2*r)*gamma_inc(6-2*r,2))
ans }
#########
>SS(10)   yields 34.91868
>SS(40)   yields 5.340422
>SS(60)   yields 180.3162
Doing the same computations in maple I get
 
34.918679360927169740821310620770402166885975646756
5.340473872869891061658721253647686930049562214921
2.2269219454888559341895277572725106
 
respectively.  Is this a precision problem?


		
---------------------------------

	[[alternative HTML version deleted]]


