When reading the current threads about root finding methods and Newton fractals,
   https://fractalforums.org/fractal-mathematics-and-new-theories/28/finding-all-roots-of-a-polynomial-with-mathematical-guarantee/2959 , and
   https://fractalforums.org/fractal-mathematics-and-new-theories/28/newton-hines-fractals/2737
I remembered that I wanted to write something about own experiences on variations an mutations on Newton style fractals.

Since I have now a bit more time than usual I would like to start a loose series about several variations in that area.

My approach contains much less mathematics than for instance the above linked threads. It is more "try it, and look what happens".
In several cases one may even doubt a matematical meaningfulness.
But anyway - as long as the resulting pictures are nice enough it doesn't matter to me ;-)

In fact the root of these variations goes back to times where calculation was  still expensive.
The thought was: 
Take a simple method that shows fast results, and then extend it systematically.
After such extension it may be possible to vary a parameter, or to change a power of the polynomial - without disturbing the complete picture (which is usually the problem with mutating simple formulas).

You will find extensions that contain for instance applications of the Newton methods on functions with multiple variables, some use of n-dimensional numbers, or matrices, or cases in which the Newton method is applied to the result of a Newton method, applied to the result of...

There is one hint which may look a bit mathematically though:
To explain what I have done I will use a specific notation.
When I write
   \( x \leftarrow N[N[f(x,y) , x] - x , y] \)
then I mean:
   1.  Take a function with two variables ( \( f(x,y) \) ) , apply the Newton method with respect to x when it comes to the differentiation ( \( N[... , x]  \) )
   2.  Subtract x from the resulting formula
   3.  Take the function you have received in 2. , and apply the Newton method with respect to y when it comes to the differentiation ( \( N[... , y]  \) )
   4.  Iterate the resulting formula (only x in this case) and look whether you get a fractal.

I will add formulas and parameters created in UF5. 
I will attach those files to each of my entries. Each will contain the params/formulas from the previous entries. This is because I had difficulties to replace attachments in existing entries.

Please do not hesitate to add your own examples, or other variations. - Have fun :-)

Linkback: https://fractalforums.org/fractal-mathematics-and-new-theories/28/slightly-extended-newton-method-fractals/3067/

Today I want to share part 3-5 of this series.

Topic is:  Extension by using higher dimensional numbers.

3D'ish triplex and quaternion variations are not discussed here (For these I had shared code for MB3D on another post).

Here I would here like to go back to the 2D roots.
In fact todays variations have their roots in my very first contact to this topic:
There is an ooold article (sorry, in German only; but the formulas should speak for themselves) that originally brought me to the idea:
Gabriele Buhren, Fraktale aus Polynomen

All of today's variations use a common calculation of root findings for \( z^3+1=0 \) . That was enough effort already ;-)

The shapes are not very spectacular. But I think the code provides a good start point for further mutations.

4-dimensional variation (not quaternion) of Newton calculation

The definition of the numbers has been discussed on other threads already. 
In short:
\( q = (a+bi) + (c+di)*j = z1+z2*j \)
\( j^2 = -i \)
\( ij = ji  \)


   zs1=z1,zs2=z2
   h=3*sqr(z1*z1-i*z2*z2)
   h1=z1-@mutil1*(z1/3+(c1*z1*z1+i*z2*(c1*z2-2*c2*z1))/h)
   z2=z2-@mutil2*(z2/3+(z1*(c2*z1-2*c1*z2)+i*c2*z2*z2)/h)
   z1=h1
   #z = @z1_weight*z1*z1+@z2_weight*z2*z2

The pictures that contain *mut* in their names are from a simple mutation that make a mistake in the calculation of the quotient q.

UF params will be attached to the last of today's posts...

6-dimensional numbers, variant 2

Definition of the numbers:
  \( s = z1 + j*z2 + k*z3 \)
  \( j^2= -(1+i); k^2= (1+i) \)
   limited commutativity:  \( jk=-kj= 1-i \)
   conjugate: \( s' = z1 - j*z2 - k*z2 \)
   definition of the division:  \( c/s = (s`*c)/(s`*s) \)

Calculation.
Use it with care. I think  there is a mistake - I don't get the 9-color basic star as for version 1.
So if anyone is bored, and want to try to re-calculate, please feel free...   

    zs1=z1, zs2=z2, zs3=z3
    z1q=z1*z1, z2q=(1+i)*z2*z2, z3q=(1+i)*z3*z3
    q= z1q+z2q-z3q, q=3*q*q
    s1= z1q-z2q+z3q
    s2= -2*z1*z2
    s3= -2*z1*z3
    z1= z1 - @mutil1*(z1/3+(c1*s1+(1+i)*(c3*s3-c2*s2)-(1-i)*(c2*s3-c3*s2))/q)
    z2= z2 - @mutil2*(z2/3+(c2*s1+c1*s2)/q)
    z3= z3 - @mutil3*(z3/3+(c3*s1+c1*s3)/q)
    #z = @z1_weight*z1*z1 + @z2_weight*z2*z2 + @z3_weight*z3*z3
[\code]


Very nice images. These generalized Newton methods are very fertile.