Here's what I get:
(%i1) :lisp (expt 2d0 (/ #c(10d0 0d0) 2));
#C(32.0 0.0)
(%i1) :lisp (gamma-lanczos (+ 1 (/ #c(10d0 0d0) 2)));
#C(119.99999999999996 0.0)
(%i1) :lisp (expt (/ 2 (float pi)) (/ (- 1 (cos (* (float pi) #c(10d0
0d0)))) 4));
#C(1.0 0.0)
Looks the same as yours.
David
BTW, the testsuite now doesn't complain.
On Wed, 2008-11-26 at 10:28 -0500, Raymond Toy wrote:
> David Ronis wrote:
> > Just updated from cvs/trunk on a i686-linux box, running clisp-2.47. I
> > had the following testsuite failure:
> >
> >
> > Running tests in rtest_gamma:
> > ********************** Problem 134 ***************
> > Input:
> > closeto(double_factorial(10.0), 3840.0, 1.0E-12)
> >
> >
> > Result:
> > 1.364242052659392E-12
> >
> > This differed from the expected result:
> > true
>
>
> Can you compute the following values for me?
>
> :lisp (expt 2d0 (/ #c(10d0 0d0) 2))
> :lisp (gamma-lanczos (+ 1 (/ #c(10d0 0d0) 2)))
> :lisp (expt (/ 2 (float pi)) (/ (- 1 (cos (* (float pi) #c(10d0 0d0)))) 4))
>
> This is what I get with clisp:
> #C(32.0 0.0)
> #C(119.99999999999996 0.0)
> #C(1.0 0.0)
>
> So double_factorial should be the product of these, which is
> #C(3839.9999999999986 0.0)
>
> Ray
>
>
> _______________________________________________
> Maxima mailing list
> Maxima at math.utexas.edu
> http://www.math.utexas.edu/mailman/listinfo/maxima
>