Group Abstract Group Abstract

Message Boards Message Boards

0
|
12.4K Views
|
9 Replies
|
0 Total Likes
View groups...
Share
Share this post:
GROUPS:

Solve for third degree polynomial

Posted 10 years ago

Hello,

The following 3rd degree polynomial has 3 real roots. All roots can be written as multiple of Cosines. Using Solve, I cannot obtain all three roots expressed in a simple manner. Any idea ? Thanks for your help.

Solve[x^3-3x-1==0,x]

Pierre

POSTED BY: Pierre Henrotay
9 Replies
In[64]:= Solve[x^3 - 3 x - 1 == 0, x, Reals]

Out[64]= {{x -> Root[-1 - 3 #1 + #1^3 &, 1]}, {x ->  Root[-1 - 3 #1 + #1^3 &, 2]}, {x -> Root[-1 - 3 #1 + #1^3 &, 3]}}

In[67]:= N[{{x -> Root[-1 - 3 #1 + #1^3 &, 1]}, {x -> Root[-1 - 3 #1 + #1^3 &, 2]}, {x -> Root[-1 - 3 #1 + #1^3 &, 3]}}]

Out[67]= {{x -> -1.53209}, {x -> -0.347296}, {x -> 1.87939}}
In[68]:= {{x -> -1.532088886237956`}, {x -> -0.3472963553338607`}, {x \-> 1.8793852415718169`}}
POSTED BY: Simon Cadrin
Posted 10 years ago

@Pierre If you are ok with numerical solution, here it is..

NSolve[x^3 - 3 x - 1 == 0, x,10]
{{x -> -1.532088886}, {x -> -0.3472963553}, {x -> 1.879385242}}
POSTED BY: Okkes Dulgerci
Posted 10 years ago
POSTED BY: Pierre Henrotay
POSTED BY: John Doty
Posted 10 years ago
ExpToTrig[ToRadicals[Root[-1 - 3 #1 + #1^3 &, 1]]] // Re
POSTED BY: Okkes Dulgerci
POSTED BY: Udo Krause
Posted 10 years ago

Thank you for pointing out my misunderstanding. I saw Cubics false in the details and assumed that was the default value.

In[1]:= FullSimplify[x /. Solve[x^3 - 3 x - 1 == 0, x]]

Out[1]= {2 Cos[\[Pi]/9], Root[-1 - 3 #1 + #1^3 &, 1], -Cos[\[Pi]/9] + Sqrt[3] Sin[\[Pi]/9]}

In[2]:= ExpToTrig[ToRadicals[Root[-1 - 3 #1 + #1^3 &, 1]]]

Out[2]= -(1/2) Cos[\[Pi]/9] + 1/2 I Sqrt[3] Cos[\[Pi]/9] - Cos[(2 \[Pi])/9] - 1/2 I Sin[\[Pi]/9] - 1/2 Sqrt[3] Sin[\[Pi]/9] - I Sin[(2 \[Pi])/9]

FullSimplify only the complex terms

In[3]:= -(1/2) Cos[\[Pi]/9] + FullSimplify[1/2 I Sqrt[3] Cos[\[Pi]/9] - 1/2 I Sin[\[Pi]/9] - I Sin[(2 \[Pi])/9]] - Cos[(2 \[Pi])/9] - 1/2 Sqrt[3] Sin[\[Pi]/9]

Out[3]= -(1/2) Cos[\[Pi]/9] - Cos[(2 \[Pi])/9] - 1/2 Sqrt[3] Sin[\[Pi]/9]
POSTED BY: Bill Simpson
Posted 10 years ago

Thanks. But Cubics->True is apparently the default for Solve (but not for Reduce). The first two roots are addressed by: Solve[x^3 - 3 x - 1 == 0, x] // ComplexExpand // FullSimplify

But the third one is the one which causes me trouble indeed. I tried various combinations (ComplexExpand, FullSimplify etc..). No success so far.

POSTED BY: Pierre Henrotay
Posted 10 years ago

In the help page for Solve click on Details and Options and read all the extra documentation for Solve.

In that find that including the option Cubics->True will give you what you want

Solve[x^3 - 3 x - 1 == 0, x, Cubics -> True]

Using FullSimplify on the result will translate two out of the three into trig form.

Using ExpToTrig and FullSimplify in various ways on parts of the third one can get it into trig form.

POSTED BY: Bill Simpson
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
Attachments
Remove
or Discard