A few questions about solve and exact polynomial roots
Subject: A few questions about solve and exact polynomial roots
From: Mike Valenzuela
Date: Wed, 4 Dec 2013 14:27:04 -0700
Wow...
Thats so cool. I'm looking forward to the next version (or whichever
version your function appears)!
I'm not sure if there are other causes of casus irreducibilis where trig
functions work better. I know that *in general* 5th order and higher order
polynomials do not have closed form solutions in the field of complex
radicals. I am truly clueless whether there are closed form solutions in
terms of trigonometric functions. I know for a fact, that many polynomials
without closed form rational solutions do have trigonometric solutions:
x^5 + x^4 - 4*x^3 - 3*x^2 + 3*x + 1, has roots at 2 * cos(2*k*%pi/11), for
suitable values of k
I've also heard of that ultra-radicals one extention to help solve these.
I really do thank you for your work as I don't have much time as a
grad-student. Unless there is a problem so ugly that I cannot do it by hand
in less time than to code a functional set of rules (and I'm not very
familiar with the terms I would need to make those rules in the first
place), then I usually just do the ugly stuff by hand. So I want to thank
you for your contribution!
On Mon, Dec 2, 2013 at 3:03 PM, Aleksas Domarkas <aleksasd873 at gmail.com>wrote:
> Mike Valenzuela on Dec 2,14:24:50 wrote:
>
> >Hello all,
> >I was messing around with some eigenvectors and values from a simple
> system
> >and came across the *casus irreduciblis* problem: sometimes complex
> >radicals are needed to express real solutions. However, I've been trying
> to
> >reduce the solution to real non-algebraic solutions.
> >(specified here)
> >
> http://en.wikipedia.org/wiki/Casus_irreducibilis#Non-algebraic_solution_in_terms_of_real_quantities
> >(and here)
> >
> http://en.wikipedia.org/wiki/Cubic_function#Trigonometric_.28and_hyperbolic.29_method
>
> >Anyways, I am wondering if maxima natively supports
> >(1) Returning non-algebraic solutions from its solving system, especially
> >if they are simplier.
> >(2) Some combination of trig simplification rules to recover the simple
> >form given in the above link.
>
> >I suppose users could manually use the above identities when applicable,
> >but they would have to be aware of it.
>
>
>
> From my "odes" package (to appear in next Maxima relase)
> function solvet(eq,x)
> returns rectform solution of polynomial equation.
> In "casus irreducibilis" give real solutions expressed
> in trigonometric functions.
>
> Examples:
> (%i2) solvet(x^3-3*x^2+1,x);
> (%o2) [x=2*cos(%pi/9)+1,x=2*cos((5*%pi)/9)+1,x=2*cos((7*%pi)/9)+1]
>
> (%i3) solvet(x^6-3*x^5-3*x^4+12*x^3-3*x^2-6*x+2=0,x);
> (%o3)
> [x=1,x=1-sqrt(3),x=sqrt(3)+1,x=2*cos((2*%pi)/9),x=2*cos((4*%pi)/9),x=2*cos((8*%pi)/9)]
>
> (%i4) solvet(x^3-15*x-5,x);
> (%o4)
> [x=2*sqrt(5)*cos(atan(sqrt(19))/3-(2*%pi)/3),x=2*sqrt(5)*cos(atan(sqrt(19))/3+(2*%pi)/3),x=2*sqrt(5)*cos(atan(sqrt(19))/3)]
>
> (%i5) solve(x^3-15*x-5,x);
> (%o5)
> [x=(-(sqrt(3)*%i)/2-1/2)*((5*sqrt(19)*%i)/2+5/2)^(1/3)+(5*((sqrt(3)*%i)/2-1/2))/((5*sqrt(19)*%i)/2+5/2)^(1/3),x=((sqrt(3)*%i)/2-1/2)*((5*sqrt(19)*%i)/2+5/2)^(1/3)+(5*(-(sqrt(3)*%i)/2-1/2))/((5*sqrt(19)*%i)/2+5/2)^(1/3),x=((5*sqrt(19)*%i)/2+5/2)^(1/3)+5/((5*sqrt(19)*%i)/2+5/2)^(1/3)]
>
> best
>
> Aleksas D
>
> _______________________________________________
> Maxima mailing list
> Maxima at math.utexas.edu
> http://www.math.utexas.edu/mailman/listinfo/maxima
>
>