Precision of find_root

asked 2012-09-17 15:28:51 -0600

TumericTJ gravatar image

Here is a straightforward question:

I am wondering about the precision of find_root. Looking into the documentation (I'm using Sage v. 5.0), there is a comment along with one of the examples that the "precision isn't very good on some machines." At the same time, it says that the routine converges unless it throws an error.

A post response from William Stein clarifies that the function is passed to SciPy which has double precision. So I'm assuming that there may be a way for more precision on some machines, but that one can always count on having at least this much.

edit retag flag offensive close merge delete

Comments

Maybe you should consider mpmath: sage: from mpmath import * sage: mp.dps = 100; mp.pretty = True sage: findroot(sin, 3) 3.141592653589793238462643383279502884197169399375105820974944592307816406286208998628034825342117068

achrzesz gravatar imageachrzesz ( 2012-09-17 19:16:35 -0600 )edit