PUT zy=2*zx*zy+cy; zx=zx2-zy2+cx; rr=ln(ln(zx*zx+zy*zy)/4)/ln(2); t=(n-rr)/(2*pi); r=128*(1+sin(t)); g=128*(1-cos(t)); b=128*(1-sin(t));