(%i1) | load(odes); |
(%i2) | eq: x^3 - 8*x^2 + x + 9 = 0; |
(%i3) | solvet(eq,x); |
(%i4) | solve(eq,x); |
Attempting to verify, that the absolute of the cube root,
is actually the square root of (61/9)...
(%i5) | expr : ((sqrt(5003)*%i)/6+709/54); |
(%i6) | abs(expr)^(1/3); |
Because this time, the base was real, its cube root
was also treated as real.
(%i7) | (61/9)/(%o6); |
Therefore, two terms are being added, that have the same absolute.
Is Maxima able to give me the numbers from both forms?
(%i8) | float(%o3); |
(%i9) | float(%o4); |
Conclusion:
If Maxima is unable to compute the numeric
version of one solution set, how is the
average user expected to do so? And yet, how
is a set of exact, non-numeric roots to become
useful, if their approximate, numeric values
cannot be verified?
(%i10) | s : map(rhs,%o8); |
(%i11) |
ProductSeries(list) := block ( product : 1, for elem in list do ( product : product * -elem ), product )$ |
(%i12) | ProductSeries(s); |
The constant term was (+9)...
The following Maxima function does, what my own suggested
program also does: Compute numeric approximations only:
(%i13) | allroots(eq); |
(%i14) | s : map(rhs,%o13)$ |
(%i15) | ProductSeries(s); |