Here is an example where mpfr-gamma and pari-gamma give different results:

x=1
for i in range(1,400):
    x=x*1.101
    if gamma(ComplexField(prec=1000)(x))-gamma(x) <> 0:
        print x

outputs

4.71129468628944e16

--
Jori Mäntysalo

--
You received this message because you are subscribed to the Google Groups 
"sage-devel" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to sage-devel+unsubscr...@googlegroups.com.
To post to this group, send email to sage-devel@googlegroups.com.
Visit this group at http://groups.google.com/group/sage-devel.
For more options, visit https://groups.google.com/groups/opt_out.

Reply via email to