On 5-May-09, at 5:00 PM, Henryk Trappmann wrote:
> > sage: a = 1/(48*sqrt(1)) - 7/(96*1**(3/2)) + 3/(32*1**(5/2)) - 5/ > (128*1**(7/2)) > sage: a > 1/384 > sage: for k in range(5): print a.n(digits=10-k) > ....: > 0.002604166668 > 0.00260416669 > 0.0026041665 > 0.002604164 > 0.00260419 I cannot guarantee that this is the source of your problem, but: computing with 5 digits of precision usually means that all *intermediate* results are computed to 5 digits, and not that the *final* result is correct to 5 digits. It is possible that you are hitting such rounding errors. Nick --~--~---------~--~----~------------~-------~--~----~ To post to this group, send email to sage-devel@googlegroups.com To unsubscribe from this group, send email to sage-devel-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-devel URLs: http://www.sagemath.org -~----------~----~----~----~------~----~------~--~---
