3.4.3 Advanced Topic: Algebraic Numbers
When you enter a Root object, the polynomial that appears in it is automatically reduced to a minimal form.
The representation of algebraic numbers.
In:= Root[24 - 2 # + 4 #^5 &, 1]
This extracts the pure function which represents the polynomial, and applies it to x.
Root objects are the way that Mathematica represents algebraic numbers. Algebraic numbers have the property that when you perform algebraic operations on them, you always get a single algebraic number as the result.
Here is the square root of an algebraic number.
In:= Sqrt[Root[2 - # + #^5 &, 1]]
RootReduce reduces this to a single Root object.
Here is a more complicated expression involving an algebraic number.
In:= Sqrt + Root[2 - # + #^5 &, 1]^2
Again this can be reduced to a single Root object, albeit a fairly complicated one.
In this simple case the Root object is automatically expressed in terms of radicals.
Operations on algebraic numbers.
In:= Root[#^2 - # - 1 &, 1]
When cubic polynomials are involved, Root objects are not automatically expressed in terms of radicals.
In:= Root[#^3 - 2 &, 1]
ToRadicals attempts to express all Root objects in terms of radicals.
If Solve and ToRadicals do not succeed in expressing the solution to a particular polynomial equation in terms of radicals, then it is a good guess that this fundamentally cannot be done. However, you should realize that there are some special cases in which a reduction to radicals is in principle possible, but Mathematica cannot find it. The simplest example is the equation , but here the solution in terms of radicals is very complicated. The equation is another example, where now
This gives a Root object involving a degree six polynomial.
is a solution.
In:= RootReduce[2^(1/3) + Sqrt]
Even though a simple form in terms of radicals does exist, ToRadicals does not find it.
Beyond degree four, most polynomials do not have roots that can be expressed at all in terms of radicals. However, for degree five it turns out that the roots can always be expressed in terms of elliptic or hypergeometric functions. The results, however, are typically much too complicated to be useful in practice.
This computes the sum of the reciprocals of the roots of
Sums of roots.
In:= RootSum[(1 + 2 # + #^5)&, (1/#)&]
Now no explicit result can be given in terms of radicals.
In:= RootSum[(1 + 2 # + #^5)&, (# Log[1 + #])&]
This expands the RootSum into a explicit sum involving Root objects.