Ask Your Question

Numerical approximations

asked 2015-06-04 01:33:02 +0200

Ducky47 gravatar image

updated 2015-06-04 04:08:18 +0200

kcrisman gravatar image

[x == -2^(1/4), x == 2^(1/4)]

What's the best/quickest way to get numerical approximations of these values of x? Thanks.

edit retag flag offensive close merge delete

3 Answers

Sort by ยป oldest newest most voted

answered 2015-06-04 04:10:46 +0200

kcrisman gravatar image

Try taking the right hand side of each solution.

sage: A =    solve(g(x)==0,x)
sage: [a.rhs().n() for a in A]
[-1.18920711500272, 1.18920711500272]

where n() is short for numerical_approx(). There are other ways of approaching this as well, such as via the solution_dict=True keyword to solve().

edit flag offensive delete link more

answered 2015-06-05 14:41:22 +0200

Ducky47 gravatar image

With my background neither is particularly straight-forward but I am happy to learn about both methods. Thanks.

edit flag offensive delete link more

answered 2015-06-05 13:39:02 +0200

tmonteil gravatar image

In your case, since g can be seen as a polynomial, you can also look at its roots in the real double field (floating-point approximations with 53 bits of precision):

sage: g.polynomial(RDF).roots()
[(-1.189207115002721, 1), (1.189207115002721, 1)]
edit flag offensive delete link more


I see two ordered pairs, I think, but if that's so I don't understand why the y-values are 1 and not 0.

Ducky47 gravatar imageDucky47 ( 2015-06-05 14:49:09 +0200 )edit

OK. "The first of the pair is the root, the second of the pair is its multiplicity."

Ducky47 gravatar imageDucky47 ( 2015-06-05 17:07:38 +0200 )edit

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools


Asked: 2015-06-04 01:33:02 +0200

Seen: 705 times

Last updated: Jun 05 '15