[math-fun] Mathematica nth root of unity question
I learned (or was reminded) only recently that the real and imaginary parts of any nth root of unity can each be expressed in terms of (an iterated rational expression in) radicals of integers (using only roots lower than the nth). 1) Is there an easy proof of this? 2) How do I get Mathematica to output these expressions? (Everything I've tried so far gets only the dumbest error messages.) For instance, Wikipedia lists the real part of (exp(2pi*i/7)) as follows: cos(2pi/7) = (-1 + ((7 + 21*sqrt(-3)/2)^(1/3) + ((7 - 21*sqrt(-3)/2)^(1/3))/6 . I'm confident that several of the 2^2 3^2 = 36 choices of roots give the correct answer, but: which ones? This may be easy to guess, but how about cos(6pi/23) ? 3) There ought to be a notation that doesn't leave these choices unspecified. —Dan
On 2015-11-24 19:42, Dan Asimov wrote:
I learned (or was reminded) only recently that the real and imaginary parts of any nth root of unity can each be expressed in terms of (an iterated rational expression in) radicals of integers (using only roots lower than the nth).
1) Is there an easy proof of this?
2) How do I get Mathematica to output these expressions? (Everything I've tried so far gets only the dumbest error messages.)
For instance, Wikipedia lists the real part of (exp(2pi*i/7)) as follows:
cos(2pi/7) = (-1 + ((7 + 21*sqrt(-3)/2)^(1/3) + ((7 - 21*sqrt(-3)/2)^(1/3))/6
In[43]:= ToRadicals[Cos[\[Pi]/7]] Out[43]= -(1/2) (-1)^(6/7) (1 + (-1)^(2/7)) In[44]:= Developer`TrigToRadicals[Cos[\[Pi]/7]] During evaluation of In[44]:= Developer`TrigToRadicals::obs: Developer`TrigToRadicals has been superseded by ToRadicals, [I don't see how. --rwg] and is now obsolete. It will not be included in future versions of the Wolfram Language. >> Out[44]= 1/2 (1/ 3 (1/2 (1 - I Sqrt[7]) - ((-1 + I Sqrt[3]) (1/2 (-1 + I Sqrt[7]) + 1/2 (-1 - I Sqrt[7]) (1/2 (-1 + I Sqrt[3]) + 1/4 (-1 + I Sqrt[3])^2)))/(2 (6 + 3/4 (-1 + I Sqrt[3]) (-1 + I Sqrt[7]) + 1/2 (-1 - I Sqrt[7]) (1 + 3/4 (-1 + I Sqrt[3])^2))^( 1/3)) - 1/ 4 (-1 + I Sqrt[3])^2 (6 + 3/4 (-1 + I Sqrt[3]) (-1 + I Sqrt[7]) + 1/2 (-1 - I Sqrt[7]) (1 + 3/4 (-1 + I Sqrt[3])^2))^(1/3)) + 1/3 (1/2 (1 + I Sqrt[7]) - ((-1 + I Sqrt[3])^2 (1/2 (-1 - I Sqrt[7]) + 1/2 (-1 + I Sqrt[7]) (1/2 (-1 + I Sqrt[3]) + 1/4 (-1 + I Sqrt[3])^2)))/(4 (6 + 3/4 (-1 + I Sqrt[3]) (-1 - I Sqrt[7]) + 1/2 (-1 + I Sqrt[7]) (1 + 3/4 (-1 + I Sqrt[3])^2))^( 1/3)) - 1/ 2 (-1 + I Sqrt[3]) (6 + 3/4 (-1 + I Sqrt[3]) (-1 - I Sqrt[7]) + 1/2 (-1 + I Sqrt[7]) (1 + 3/4 (-1 + I Sqrt[3])^2))^(1/3))) In[47]:= MinimalPolynomial[Cos[2 \[Pi]/7]]@x Out[47]= -1 - 4 x + 4 x^2 + 8 x^3 In[50]:= Solve[0 == %] Out[50]= {{x -> 1/6 (-1 + 7^(2/3)/(1/2 (1 + 3 I Sqrt[3]))^( 1/3) + (7/2 (1 + 3 I Sqrt[3]))^(1/3))}, {x -> -(1/6) - ( 7^(2/3) (1 + I Sqrt[3]))/(6 2^(2/3) (1 + 3 I Sqrt[3])^(1/3)) - 1/12 (1 - I Sqrt[3]) (7/2 (1 + 3 I Sqrt[3]))^(1/3)}, {x -> -(1/ 6) - (7^(2/3) (1 - I Sqrt[3]))/( 6 2^(2/3) (1 + 3 I Sqrt[3])^(1/3)) - 1/12 (1 + I Sqrt[3]) (7/2 (1 + 3 I Sqrt[3]))^(1/3)}}
.
I'm confident that several of the 2^2 3^2 = 36 choices of roots give the correct answer, but: which ones? This may be easy to guess, but how about cos(6pi/23) ?
In[51]:= Developer`TrigToRadicals[Cos[6 \[Pi]/23]] During evaluation of In[51]:= Developer`TrigToRadicals::obs: Developer`TrigToRadicals has been superseded by ToRadicals, and is now obsolete. It will not be included in future versions of the Wolfram Language. >> Out[51]= <Enormous expression> In[52]:= Short[%] Out[52]//Short= 1/2 (1/11 (1/2 (-1+I Sqrt[23])+<<9>>+(<<1>> <<1>>)/(9<<5>>5<<1>><<1>>))+<<1>>) --rwg
3) There ought to be a notation that doesn't leave these choices unspecified.
"The ordering used by Root[f,k] takes real roots to come before complex ones, and takes complex conjugate pairs of roots to be adjacent. "
—Dan
Thanks, WDS and RWG for your answers! RWG: At the bottom, what do the expressions <<4k+1>> mean? (And I hope you'll report to WR that "Developer`TrigToRadicals has been superseded by ToRadicals" does *not* appear to be the case.) —Dan P.S. I asked about cos(6pi/23) in particular because Kummer found that the cyclotomic field Q(exp(2pi/23)) is the first one whose ring of algebraic integers is *not* a unique factorization domain. Thereby showing Lamé's putative proof of FLT invalid, as Liouville suspected (http://sites.mathdoc.fr/JMPA/PDF/JMPA_1847_1_12_A15_0.pdf <http://sites.mathdoc.fr/JMPA/PDF/JMPA_1847_1_12_A15_0.pdf>).
On Nov 25, 2015, at 9:55 AM, rwg <rwg@sdf.org> wrote:
On 2015-11-24 19:42, Dan Asimov wrote:
I learned (or was reminded) only recently that the real and imaginary parts of any nth root of unity can each be expressed in terms of (an iterated rational expression in) radicals of integers (using only roots lower than the nth). 1) Is there an easy proof of this? 2) How do I get Mathematica to output these expressions? (Everything I've tried so far gets only the dumbest error messages.) For instance, Wikipedia lists the real part of (exp(2pi*i/7)) as follows: cos(2pi/7) = (-1 + ((7 + 21*sqrt(-3)/2)^(1/3) + ((7 - 21*sqrt(-3)/2)^(1/3))/6 In[43]:= ToRadicals[Cos[\[Pi]/7]]
Out[43]= -(1/2) (-1)^(6/7) (1 + (-1)^(2/7))
In[44]:= Developer`TrigToRadicals[Cos[\[Pi]/7]]
During evaluation of In[44]:= Developer`TrigToRadicals::obs: Developer`TrigToRadicals has been superseded by ToRadicals, [I don't see how. --rwg] and is now obsolete. It will not be included in future versions of the Wolfram Language. >>
Out[44]= 1/2 (1/ 3 (1/2 (1 - I Sqrt[7]) - ((-1 + I Sqrt[3]) (1/2 (-1 + I Sqrt[7]) + 1/2 (-1 - I Sqrt[7]) (1/2 (-1 + I Sqrt[3]) + 1/4 (-1 + I Sqrt[3])^2)))/(2 (6 + 3/4 (-1 + I Sqrt[3]) (-1 + I Sqrt[7]) + 1/2 (-1 - I Sqrt[7]) (1 + 3/4 (-1 + I Sqrt[3])^2))^( 1/3)) - 1/ 4 (-1 + I Sqrt[3])^2 (6 + 3/4 (-1 + I Sqrt[3]) (-1 + I Sqrt[7]) + 1/2 (-1 - I Sqrt[7]) (1 + 3/4 (-1 + I Sqrt[3])^2))^(1/3)) + 1/3 (1/2 (1 + I Sqrt[7]) - ((-1 + I Sqrt[3])^2 (1/2 (-1 - I Sqrt[7]) + 1/2 (-1 + I Sqrt[7]) (1/2 (-1 + I Sqrt[3]) + 1/4 (-1 + I Sqrt[3])^2)))/(4 (6 + 3/4 (-1 + I Sqrt[3]) (-1 - I Sqrt[7]) + 1/2 (-1 + I Sqrt[7]) (1 + 3/4 (-1 + I Sqrt[3])^2))^( 1/3)) - 1/ 2 (-1 + I Sqrt[3]) (6 + 3/4 (-1 + I Sqrt[3]) (-1 - I Sqrt[7]) + 1/2 (-1 + I Sqrt[7]) (1 + 3/4 (-1 + I Sqrt[3])^2))^(1/3)))
In[47]:= MinimalPolynomial[Cos[2 \[Pi]/7]]@x
Out[47]= -1 - 4 x + 4 x^2 + 8 x^3
In[50]:= Solve[0 == %]
Out[50]= {{x -> 1/6 (-1 + 7^(2/3)/(1/2 (1 + 3 I Sqrt[3]))^( 1/3) + (7/2 (1 + 3 I Sqrt[3]))^(1/3))}, {x -> -(1/6) - ( 7^(2/3) (1 + I Sqrt[3]))/(6 2^(2/3) (1 + 3 I Sqrt[3])^(1/3)) - 1/12 (1 - I Sqrt[3]) (7/2 (1 + 3 I Sqrt[3]))^(1/3)}, {x -> -(1/ 6) - (7^(2/3) (1 - I Sqrt[3]))/( 6 2^(2/3) (1 + 3 I Sqrt[3])^(1/3)) - 1/12 (1 + I Sqrt[3]) (7/2 (1 + 3 I Sqrt[3]))^(1/3)}}
. I'm confident that several of the 2^2 3^2 = 36 choices of roots give the correct answer, but: which ones? This may be easy to guess, but how about cos(6pi/23) ?
In[51]:= Developer`TrigToRadicals[Cos[6 \[Pi]/23]]
During evaluation of In[51]:= Developer`TrigToRadicals::obs: Developer`TrigToRadicals has been superseded by ToRadicals, and is now obsolete. It will not be included in future versions of the Wolfram Language. >>
Out[51]= <Enormous expression> In[52]:= Short[%] Out[52]//Short= 1/2 (1/11 (1/2 (-1+I Sqrt[23])+<<9>>+(<<1>> <<1>>)/(9<<5>>5<<1>><<1>>))+<<1>>) --rwg
3) There ought to be a notation that doesn't leave these choices unspecified. "The ordering used by Root[f,k] takes real roots to come before complex ones, and takes complex conjugate pairs of roots to be adjacent. " —Dan
_______________________________________________ math-fun mailing list math-fun@mailman.xmission.com https://mailman.xmission.com/cgi-bin/mailman/listinfo/math-fun
On 2015-11-25 11:08, Dan Asimov wrote:
Thanks, WDS and RWG for your answers!
RWG: At the bottom, what do the expressions <<4k+1>> mean?
In the doc for Short, Properties & Relations (2) Shortened outputs contain Skeleton objects: In[1]:= Short[Range[100]] Out[1]//Short= {1,2,3,4,5,6,7,8,9,10,11,12,13,<<75>>,89,90,91,92,93,94,95,96,97,98,99,100} In[2]:= Skeleton[75] Out[2]= <<75>> The actual notebook shows a "large output" object with various viewing controls, but clicking on or copying skeletonized subexpressions does not seem to selectively deskeletonize them.
(And I hope you'll report to WR that
"Developer`TrigToRadicals has been superseded by ToRadicals"
does *not* appear to be the case.)
I was going to, but then "realized" (hallucinated?) a simple workaround, which I now forget. I guess I'll have to meekly request it. At least you can get the degree: In[95]:= MinimalPolynomial[Cos[6 \[Pi]/23]] Out[95]= -1 - 12 #1 + 60 #1^2 + 280 #1^3 - 560 #1^4 - 1792 #1^5 + 1792 #1^6 + 4608 #1^7 - 2304 #1^8 - 5120 #1^9 + 1024 #1^10 + 2048 #1^11 & and a handle on the actual root: In[93]:= RootApproximant[Cos[6 \[Pi]/23`69]] Out[93]= Root[-1 - 12 #1 + 60 #1^2 + 280 #1^3 - 560 #1^4 - 1792 #1^5 + 1792 #1^6 + 4608 #1^7 - 2304 #1^8 - 5120 #1^9 + 1024 #1^10 + 2048 #1^11 &, 9]
—Dan
P.S. I asked about cos(6pi/23) in particular because Kummer found that the cyclotomic field Q(exp(2pi/23)) is the first one whose ring of algebraic integers is *not* a unique factorization domain. Thereby showing Lamé's putative proof of FLT invalid, as Liouville suspected (http://sites.mathdoc.fr/JMPA/PDF/JMPA_1847_1_12_A15_0.pdf <http://sites.mathdoc.fr/JMPA/PDF/JMPA_1847_1_12_A15_0.pdf>).
Interesting! You might check out NumberFieldRootsOfUnity, NumberFieldFundamentalUnits, ToNumberField, AlgebraicNumber..., AlgebraicInteger..., but NumberFieldClassNumber[%93] NumberFieldClassNumber::unkn: The class number of the number field generated by Root[-1-12 #1+60 #1^2+280 #1^3-560 #1^4-1792 #1^5+1792 #1^6+4608 #1^7-2304 #1^8-5120 #1^9+1024 #1^10+2048 #1^11&,9,0] is not yet available. >> --rwg
On Nov 25, 2015, at 9:55 AM, rwg <rwg@sdf.org> wrote:
On 2015-11-24 19:42, Dan Asimov wrote:
I learned (or was reminded) only recently that the real and imaginary parts of any nth root of unity can each be expressed in terms of (an iterated rational expression in) radicals of integers (using only roots lower than the nth). 1) Is there an easy proof of this? 2) How do I get Mathematica to output these expressions? (Everything I've tried so far gets only the dumbest error messages.) For instance, Wikipedia lists the real part of (exp(2pi*i/7)) as follows: cos(2pi/7) = (-1 + ((7 + 21*sqrt(-3)/2)^(1/3) + ((7 - 21*sqrt(-3)/2)^(1/3))/6 In[43]:= ToRadicals[Cos[\[Pi]/7]]
Out[43]= -(1/2) (-1)^(6/7) (1 + (-1)^(2/7))
In[44]:= Developer`TrigToRadicals[Cos[\[Pi]/7]]
During evaluation of In[44]:= Developer`TrigToRadicals::obs: Developer`TrigToRadicals has been superseded by ToRadicals, [I don't see how. --rwg] and is now obsolete. It will not be included in future versions of the Wolfram Language. >>
Out[44]= 1/2 (1/ 3 (1/2 (1 - I Sqrt[7]) - ((-1 + I Sqrt[3]) (1/2 (-1 + I Sqrt[7]) + 1/2 (-1 - I Sqrt[7]) (1/2 (-1 + I Sqrt[3]) + 1/4 (-1 + I Sqrt[3])^2)))/(2 (6 + 3/4 (-1 + I Sqrt[3]) (-1 + I Sqrt[7]) + 1/2 (-1 - I Sqrt[7]) (1 + 3/4 (-1 + I Sqrt[3])^2))^( 1/3)) - 1/ 4 (-1 + I Sqrt[3])^2 (6 + 3/4 (-1 + I Sqrt[3]) (-1 + I Sqrt[7]) + 1/2 (-1 - I Sqrt[7]) (1 + 3/4 (-1 + I Sqrt[3])^2))^(1/3)) + 1/3 (1/2 (1 + I Sqrt[7]) - ((-1 + I Sqrt[3])^2 (1/2 (-1 - I Sqrt[7]) + 1/2 (-1 + I Sqrt[7]) (1/2 (-1 + I Sqrt[3]) + 1/4 (-1 + I Sqrt[3])^2)))/(4 (6 + 3/4 (-1 + I Sqrt[3]) (-1 - I Sqrt[7]) + 1/2 (-1 + I Sqrt[7]) (1 + 3/4 (-1 + I Sqrt[3])^2))^( 1/3)) - 1/ 2 (-1 + I Sqrt[3]) (6 + 3/4 (-1 + I Sqrt[3]) (-1 - I Sqrt[7]) + 1/2 (-1 + I Sqrt[7]) (1 + 3/4 (-1 + I Sqrt[3])^2))^(1/3)))
In[47]:= MinimalPolynomial[Cos[2 \[Pi]/7]]@x
Out[47]= -1 - 4 x + 4 x^2 + 8 x^3
In[50]:= Solve[0 == %]
Out[50]= {{x -> 1/6 (-1 + 7^(2/3)/(1/2 (1 + 3 I Sqrt[3]))^( 1/3) + (7/2 (1 + 3 I Sqrt[3]))^(1/3))}, {x -> -(1/6) - ( 7^(2/3) (1 + I Sqrt[3]))/(6 2^(2/3) (1 + 3 I Sqrt[3])^(1/3)) - 1/12 (1 - I Sqrt[3]) (7/2 (1 + 3 I Sqrt[3]))^(1/3)}, {x -> -(1/ 6) - (7^(2/3) (1 - I Sqrt[3]))/( 6 2^(2/3) (1 + 3 I Sqrt[3])^(1/3)) - 1/12 (1 + I Sqrt[3]) (7/2 (1 + 3 I Sqrt[3]))^(1/3)}}
. I'm confident that several of the 2^2 3^2 = 36 choices of roots give the correct answer, but: which ones? This may be easy to guess, but how about cos(6pi/23) ?
In[51]:= Developer`TrigToRadicals[Cos[6 \[Pi]/23]]
During evaluation of In[51]:= Developer`TrigToRadicals::obs: Developer`TrigToRadicals has been superseded by ToRadicals, and is now obsolete. It will not be included in future versions of the Wolfram Language. >>
Out[51]= <Enormous expression> In[52]:= Short[%] Out[52]//Short= 1/2 (1/11 (1/2 (-1+I Sqrt[23])+<<9>>+(<<1>> <<1>>)/(9<<5>>5<<1>><<1>>))+<<1>>) --rwg
3) There ought to be a notation that doesn't leave these choices unspecified. "The ordering used by Root[f,k] takes real roots to come before complex ones, and takes complex conjugate pairs of roots to be adjacent. " —Dan
_______________________________________________ math-fun mailing list math-fun@mailman.xmission.com https://mailman.xmission.com/cgi-bin/mailman/listinfo/math-fun
_______________________________________________ math-fun mailing list math-fun@mailman.xmission.com https://mailman.xmission.com/cgi-bin/mailman/listinfo/math-fun
participants (3)
-
Dan Asimov -
Dan Asimov -
rwg