Numerical approximations
g(x)=x^2-sqrt(2)
solve(g(x)==0,x)
[x == -2^(1/4), x == 2^(1/4)]
What's the best/quickest way to get numerical approximations of these values of x? Thanks.
g(x)=x^2-sqrt(2)
solve(g(x)==0,x)
[x == -2^(1/4), x == 2^(1/4)]
What's the best/quickest way to get numerical approximations of these values of x? Thanks.
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()
.
With my background neither is particularly straight-forward but I am happy to learn about both methods. Thanks.
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)]
Please start posting anonymously - your entry will be published after you log in or create a new account.
Asked: 2015-06-04 01:33:02 +0100
Seen: 977 times
Last updated: Jun 05 '15