Status: Accepted
Owner: ----
Labels: Type-Defect Priority-Medium
New issue 2878 by [email protected]: numerical evaluation problem
http://code.google.com/p/sympy/issues/detail?id=2878
This is confusing to me:
a = integrate(sin(x)**10, x).diff(x).simplify() - sin(x)**10
a
4 2 4
63*sin (x) 63*sin (x) 63*cos (x) 63
- ---------- + ---------- + ---------- - ---
128 64 128 128
a.subs(x, S("0.6")).n(100)
0.0000000000000000555111512312578270211815834045410156250000000000000000000000
0000000000000000000000000000000000000000
a.subs(x, S(6)/10).n(100)
-.0e-213
Why does the 0.6, which should be exactly 6/10 return a different result?
Could it be, that 0.6 is interpreted as a floating point number with just
16 digits of precision?
--
You received this message because you are subscribed to the Google Groups
"sympy-issues" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sympy-issues?hl=en.
