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

pixels = W:1024 H:1024

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

zr = x
zi = y

m = 1.0
n = 2.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)
dr = complex_div_r(z4r - m, z4i - m, n*z3r, n*z3i)
di = complex_div_i(z4r - m, z4i - m, n*z3r, n*z3i)
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)
dr = complex_div_r(z4r - m, z4i - m, n*z3r, n*z3i)
di = complex_div_i(z4r - m, z4i - m, n*z3r, n*z3i)
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)
dr = complex_div_r(z4r - m, z4i - m, n*z3r, n*z3i)
di = complex_div_i(z4r - m, z4i - m, n*z3r, n*z3i)
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)
dr = complex_div_r(z4r - m, z4i - m, n*z3r, n*z3i)
di = complex_div_i(z4r - m, z4i - m, n*z3r, n*z3i)
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)
dr = complex_div_r(z4r - m, z4i - m, n*z3r, n*z3i)
di = complex_div_i(z4r - m, z4i - m, n*z3r, n*z3i)
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)
dr = complex_div_r(z4r - m, z4i - m, n*z3r, n*z3i)
di = complex_div_i(z4r - m, z4i - m, n*z3r, n*z3i)
zr = zr - dr
zi = zi - di


r = zr
g = zi
b = -zi
