[sage-devel] elliptic_e(0.5, 0.1) differs from Mathematica 7 by about 0.04%.
Subject: [sage-devel] elliptic_e(0.5, 0.1) differs from Mathematica 7 by about 0.04%.
From: Valery Pipin
Date: Mon, 10 Aug 2009 13:38:31 +0900
> Dr. David Kirkby wrote:
> > I've built Sage on a sun4u SPARC machine, using Maxima 5.19.0 and ECL
> > 9.8.1
> >
> > Sage fails about 15 tests, one of which is below. I computed the result
> > in Mathematica too, and find the result on the SPARC does differ from
> > what Sage expects, and the result Sage expects is much closer to the
> > result from Mathematica.
> >
> > sage: elliptic_eu (0.5, 0.1)
> > Expected:
> > 0.496054551287
> > Got:
> > 0.495848403419
>
> sorry, I posted the wrong result - there is this failure, but the one I
> meant to post was:
>
> sage: elliptic_e(0.5, 0.1)
> Expected:
> 0.498011394499
> Got:
> 0.497801100392
>
> Mathematica gives: 0.49801139449883153311546104061744810584963105068054
[va at colombo py-q]$ rmaxima
Maxima 5.18.1 http://maxima.sourceforge.net
Using Lisp SBCL 1.0.29
Distributed under the GNU Public License. See the file COPYING.
Dedicated to the memory of William Schelter.
The function bug_report() provides bug reporting information.
(%i1) elliptic_e(0.5, 0.1);
(%o1) .4980113944988315