Mathematica 9 is now available
Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
2006
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*December
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 2006

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: a Quaternion quadratic level Pisot polynomial

  • To: mathgroup at smc.vnet.net
  • Subject: [mg68637] Re: a Quaternion quadratic level Pisot polynomial
  • From: Roger Bagula <rlbagula at sbcglobal.net>
  • Date: Mon, 14 Aug 2006 06:44:08 -0400 (EDT)
  • References: <ebmtoj$6k1$1@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com

jdturner wrote:

>Roger Bagula wrote:
>  
>
>>I tried all kinds of combinations.
>>q^4-4*q^2-1
>>is the opposite of a Pisot ( one root inside the Uniot Quaternion and
>>the rest outside).
>>
>>(-1+I+J+K)*q^2-2*q-1==0
>>solves to give a set of Pisot like roots.
>>{{-0.34833, 0.385615, -0.34833, -0.15167},
>> {-0.15167, -0.885615, -0.15167, -0.34833},
>>{-0.420313, 0.116972,  0.116972, 0.116972},
>> {-0.0796875, -0.616972, -0.616972, -0.616972}}
>>
>>    
>>
>I have a question for you.  Why four roots, since the equation is
>easily solved by completing the square.  First rewrite the equation as
>
>q^2+c1*q+c0=0, c1=-2/(-1,I+J+k), c0 = -1/(-1+I+J+K)
>
>yielding the roots
>qp = -c1/2+qsqrt( -c0+c1~c1/4 )
>    = [ - 0.0796874903417d0, [ - 0.61697245625346d0,
>         - 0.61697245625346d0,  - 0.61697245625346d0]]
>qm= -c1/2-qsqrt( -c0+c1~c1/4 )
>    = [ - 0.4203125096583d0, [0.11697245625346d0, 0.11697245625346d0,
>   0.11697245625346d0]]
>
>which coorespond to your last two rts.
>
>Jim Turner
>
>  
>

Jim Turner,
That's a good question.
It appears Mathematica solves four equations in four unknowns and gives at least four answers
even with {t_Real, x_Real, y_Real, z_Real}
Peter Pleasants finds :
4x^4 + 4x^3 + 6x^2 + 4x + 1
Has Pisot roots too. Mathematica gives 10 roots... two of which are the Real root
that are the same as this polynomial.
You are welcome to try your method on that.
I was trying to get a method in Mathematica that gave
correct roots for any polynomial
and the specific polynomial that behaved as a Pisot:
1) Polynomial with integer coefficients
2) one real modulus/absolute values greater than one
3) other modulus/absolutue values less than one
Complex limited , Real limited ones ones are already known. 

I can't say that I understand Quaternion roots all the well
even with my work.
Since you are the only one to respond,
I assume others are even less brave.
Roger Bagula


  • Prev by Date: RE: ContourPlot and finite approximation to level sets
  • Next by Date: Re: Postscipt fonts in graphics
  • Previous by thread: Re: a Quaternion quadratic level Pisot polynomial
  • Next by thread: "No more memory available" -- a recurring problem