#https://www.shadertoy.com/view/ltB3WW

pixels = W:1024 H:1024

x = from -2 to 2 W
y = from -2 to 2 H

zr = x
zi = y

m = 1.0
n = 3.0

z2r = complex_mul_r(zr, zi, zr, zi)
z2i = complex_mul_i(zr, zi, zr, zi)
z3r = complex_mul_r(z2r, z2i, zr, zi)
z3i = complex_mul_i(z2r, z2i, zr, zi)
z4r = complex_mul_r(z2r, z2i, z2r, z2i)
z4i = complex_mul_i(z2r, z2i, z2r, z2i)
z5r = complex_mul_r(z3r, z3i, z2r, z2i)
z5i = complex_mul_i(z3r, z3i, z2r, z2i)
dr = complex_div_r(z5r - m, z5i - m, n*z4r, n*z4i)
di = complex_div_i(z5r - m, z5i - m, n*z4r, n*z4i)
zr = zr - dr
zi = zi - di

z2r = complex_mul_r(zr, zi, zr, zi)
z2i = complex_mul_i(zr, zi, zr, zi)
z3r = complex_mul_r(z2r, z2i, zr, zi)
z3i = complex_mul_i(z2r, z2i, zr, zi)
z4r = complex_mul_r(z2r, z2i, z2r, z2i)
z4i = complex_mul_i(z2r, z2i, z2r, z2i)
z5r = complex_mul_r(z3r, z3i, z2r, z2i)
z5i = complex_mul_i(z3r, z3i, z2r, z2i)
dr = complex_div_r(z5r - m, z5i - m, n*z4r, n*z4i)
di = complex_div_i(z5r - m, z5i - m, n*z4r, n*z4i)
zr = zr - dr
zi = zi - di

z2r = complex_mul_r(zr, zi, zr, zi)
z2i = complex_mul_i(zr, zi, zr, zi)
z3r = complex_mul_r(z2r, z2i, zr, zi)
z3i = complex_mul_i(z2r, z2i, zr, zi)
z4r = complex_mul_r(z2r, z2i, z2r, z2i)
z4i = complex_mul_i(z2r, z2i, z2r, z2i)
z5r = complex_mul_r(z3r, z3i, z2r, z2i)
z5i = complex_mul_i(z3r, z3i, z2r, z2i)
dr = complex_div_r(z5r - m, z5i - m, n*z4r, n*z4i)
di = complex_div_i(z5r - m, z5i - m, n*z4r, n*z4i)
zr = zr - dr
zi = zi - di

z2r = complex_mul_r(zr, zi, zr, zi)
z2i = complex_mul_i(zr, zi, zr, zi)
z3r = complex_mul_r(z2r, z2i, zr, zi)
z3i = complex_mul_i(z2r, z2i, zr, zi)
z4r = complex_mul_r(z2r, z2i, z2r, z2i)
z4i = complex_mul_i(z2r, z2i, z2r, z2i)
z5r = complex_mul_r(z3r, z3i, z2r, z2i)
z5i = complex_mul_i(z3r, z3i, z2r, z2i)
dr = complex_div_r(z5r - m, z5i - m, n*z4r, n*z4i)
di = complex_div_i(z5r - m, z5i - m, n*z4r, n*z4i)
zr = zr - dr
zi = zi - di

z2r = complex_mul_r(zr, zi, zr, zi)
z2i = complex_mul_i(zr, zi, zr, zi)
z3r = complex_mul_r(z2r, z2i, zr, zi)
z3i = complex_mul_i(z2r, z2i, zr, zi)
z4r = complex_mul_r(z2r, z2i, z2r, z2i)
z4i = complex_mul_i(z2r, z2i, z2r, z2i)
z5r = complex_mul_r(z3r, z3i, z2r, z2i)
z5i = complex_mul_i(z3r, z3i, z2r, z2i)
dr = complex_div_r(z5r - m, z5i - m, n*z4r, n*z4i)
di = complex_div_i(z5r - m, z5i - m, n*z4r, n*z4i)
zr = zr - dr
zi = zi - di

z2r = complex_mul_r(zr, zi, zr, zi)
z2i = complex_mul_i(zr, zi, zr, zi)
z3r = complex_mul_r(z2r, z2i, zr, zi)
z3i = complex_mul_i(z2r, z2i, zr, zi)
z4r = complex_mul_r(z2r, z2i, z2r, z2i)
z4i = complex_mul_i(z2r, z2i, z2r, z2i)
z5r = complex_mul_r(z3r, z3i, z2r, z2i)
z5i = complex_mul_i(z3r, z3i, z2r, z2i)
dr = complex_div_r(z5r - m, z5i - m, n*z4r, n*z4i)
di = complex_div_i(z5r - m, z5i - m, n*z4r, n*z4i)
zr = zr - dr
zi = zi - di

r = zr
g = zi
b = -zi
