We've been having some fun with cubic equations. The other day we figured out that the relations A^2-4=B, when applied transitively three times in a row, generates the roots of a cubic equation. The trick is that you have to find a starting point whereby, after three iterations, you get back to where you've started. You can see what is happening on this graph:
You can see that starting with A=1.6, we generate B= -1.4 and C= -2.4. Let's look at the elementary symmetric polynomials in A, B and C to see what our cubic equation should look like:
A + B + C = -2.2
AB + BC + CA = -2.72
ABC = 5.376
Now, if we remember what Balarka did the other day, he generated an octic equation by substituing the recursive relationship into itself three times, factored out the trivial quadratic, and then factored the resulting sextic into two terms, one of which was:
x^3 + 2x^2 -3x +5
This is really close to the equations you get from my symmetric polynomials in A, B and C....considering I just drew those lines on the graph by eyeball (no cheating), I'm a bit amazed how well it worked. (Remember the signs alternate when you sub in the coefficients...I told Balarka he had a sign error yesterday, but I was wrong.)
It was actually quite a surprise when Balarka's octic factored cleanly. We had started with the recursive relationship A^2 - 2 = B....going through the same steps, generating the octic, factoring out the quadratic, and factoring the resulting sextic...well, we knew that one had to factor because there was already a known factor that was generated from combining the seventh roots of unity. But we didn't actually expect the arbitrary relationship A^2-4 to also factor.
So what about A^2-6? Well, I could have tried graphing it, but then I had a better idea...why not program the relationship into Excel? It's a super-easy spreadsheet calculation...you literally need only five or six cells. You just keep iterating your seed value until you get the recursion working. And of course, with instant correction, you can easily work it out to four or five decimal places. It's kind of fun.
I tried it first on A^2-4, and verified that the cycle did in fact exist...in fact, there were (of course) two cycles, corresponding to the two cubics that we factored out of the sextic. But how could we tell that these were the roots of a cubic equation in integer coefficients? That was easy! I just got Excel to display the values of the sums and products of the roots. At the sime time as the cyclic relationship converged back on itself, you could see the sums and products zeroing in on integer values. We had not only found the roots but also the cubic equation which they came from?
So what about A^2 - 6? (For some reason it seemed reasonable to stick to even values.) This one was different. I found the following triplets with cyclic relationships:
(1.7425, -2.9636, 2.7827)
(2.0547, -1.7782, -2.8380)
but checking the elementary symmetric combinations of these triplets, it was clear that they were not the roots of cubic equations with integer coefficients:
(1.56, -8.56, -14.37)
(-2.56, -4.44, 10.37)
(These triplets came up as integers when we worked with A^2-4). But were the coefficients of the two cubic equations at least something like the roots of a quadratic? Indeed they were. Examining the pair of triplets listed above, you can see by taking sums and products that the coefficients are related as follows:
The x^2 coefficients are given by the relationship x^2 + x - 4 = 0
The x coefficients are given by the relationship x^2 + 13x +38 = 0
The constant coefficients are given by x^2 +4x - 149 = 0
And that's about all I've got to say for now...