Jim Muth wrote:
FOTD -- July 20, 2009 (Rating 8) Possible new formula alert! ....created by the DivideBrot series of formulas, with a mixture of Z^101 and Z^2. The DivideJulibrot formula included in the parameter file might be posted for the first time in this FOTD. I'm not sure.
Actually, this is the third version by that same name. Each one has very slight differences. Here are all three versions: frm:DivideJulibrot { ; Jim Muth (first version) ; draws slices of DivideBrot Julibrots pix=pixel, u=real(pix), v=imag(pix), a=pi*real(p1*0.0055555555555556), b=pi*imag(p1*0.0055555555555556), g=pi*real(p2*0.0055555555555556), d=pi*imag(p2*0.0055555555555556), ca=cos(a), cb=cos(b), sb=sin(b), cg=cos(g), sg=sin(g), cd=cos(d), sd=sin(d), aa=real(p5)-2, bb=imag(p5)+0.00000000000000000001, p=u*cg*cd-v*(ca*sb*sg*cd+ca*cb*sd), q=u*cg*sd+v*(ca*cb*cd-ca*sb*sg*sd), r=u*sg+v*ca*sb*cg, s=v*sin(a), c=p+flip(q)+p3, z=r+flip(s)+p4: z=sqr(z)/(z^(-aa)+bb)+c |z| < 1000000 } frm:DivideJulibrot { ; Jim Muth (second version) ; draws 4-D slices of DivideBrot Julibrots pix=pixel, u=real(pix), v=imag(pix), a=pi*real(p1*0.0055555555555556), b=pi*imag(p1*0.0055555555555556), g=pi*real(p2*0.0055555555555556), d=pi*imag(p2*0.0055555555555556), ca=cos(a), cb=cos(b), sb=sin(b), cg=cos(g), sg=sin(g), cd=cos(d), sd=sin(d), aa=-(real(p5)-2), bb=imag(p5+0.00000000000000000001), p=u*cg*cd-v*(ca*sb*sg*cd+ca*cb*sd), q=u*cg*sd+v*(ca*cb*cd-ca*sb*sg*sd), r=u*sg+v*ca*sb*cg, s=v*sin(a), c=p+flip(q)+p3, z=r+flip(s)+p4: z=sqr(z)/(z^(aa)+bb)+c |z| < 1000000 } frm:DivideJulibrot { ; Jim Muth (third version) ; draws 4-D slices of DivideBrot Julibrots pix=pixel, u=real(pix), v=imag(pix), a=pi*real(p1*0.0055555555555556), b=pi*imag(p1*0.0055555555555556), g=pi*real(p2*0.0055555555555556), d=pi*imag(p2*0.0055555555555556), ca=cos(a), cb=cos(b), sb=sin(b), cg=cos(g), sg=sin(g), cd=cos(d), sd=sin(d), aa=-(real(p5)-2), bb=(imag(p5)+0.00000000000000000000001), p=u*cg*cd-v*(ca*sb*sg*cd+ca*cb*sd), q=u*cg*sd+v*(ca*cb*cd-ca*sb*sg*sd), r=u*sg+v*ca*sb*cg, s=v*sin(a), c=p+flip(q)+p3, z=r+flip(s)+p4: z=sqr(z)/(z^(aa)+bb)+c |z| < 1000000 } Sincerely, P.N.L. ------------------------------------------------- http://home.att.net/~Paul.N.Lee/PNL_Fractals.html http://www.Nahee.com/Fractals/