Today's fractal is a remix of my Oct. 24 post. Same formula and palette, different parameters. Like the Oct.24 post ready made at http://maxitersfractalfollies.blogspot.com Twilight Remix { ; fract484.gif ; blank ; calctime 0:17:07.88 ; created Oct 26, 2010 ; Fractint Version 2004 Patchlevel 10 reset=2004 type=formula formulafile=kerrym.frm formulaname=hermanm_man-polar center-mag=-1.89610111352240000/+0.00222840388810970/9616.574/1/22.49999\ 99999827232/-2.86882184674652763e-011 params=4/2/1/60/3/0 float=y maxiter=1500 inside=0 proximity=1 outside=fmod decomp=256 colors=000302403503504604705805806907A07B08B08C09D1AE1AF1BF1BG1CH1CI1DI1\ EJ1EK1FL1FL1GM1GN1HO2IN2KM2NK2QJ2TH2VG2YE2`D2cB1fA1i81l71n51q41t21w00z03\ z07z0Bz0Fz0Jz0Nz0Rz0Vz0Zz0bz0fz0jz0nz0rz0vz0zz0xz0vz0tz0rz0pz0nz0lz0jz0h\ z0fz0dz0bz0`z0Zz0Xz0Wz0Uz0Sz0Qz0Oz0Mz0Kz0Iz0Gz0Ez0Cz0Az08z06z04z02z00z00\ y10x20v20u30s40r50p50o61n71l81k91i91hA1gB1eC1dC1bD1aE1_F1ZF1YG1WH1VI1TJ2\ SJ2RK2PL2OM2MM2LN2JO2IM2GJ2EH2CE1AC18916614312000101302503605806A07C19D1\ 9E1AF1AF1BG1BH1CI1DI1DJ1EK1EL1FL1FM1GN1HO2IO2JQ2IS2HV2GX2F_2Da2Cd2Bf2Ah1\ 9k18m17p15r14u13w12z00z20z50z70zA0zC0zF0zH0zK0zM0zP0zR0zU0zX0zZ0za0zc0zf\ 0zh0zk0zm0zp0zr0zu0zw0zz0zw0zu0zr0zo0zl0zj0zg0zd0za0z_0zX0zU0zS0zP0zM0zJ\ 0zH0zE0zB0z80z60z30z00x01v02s03q04n06l07i08g19e1Ab1B`1CY1EW1FT1GR1HO2JN2\ IM2HK2GJ2FH2DG2CE2BD2AB19A18817715514413212000000101201 } frm:hermanm_man-polar { ; Kerry Mitchell 16feb98 ; ; real(p1) = z exponent (use integer >= 2; m=n-1) ; imag(p1) = g exponent (integers) ; real(p2) = alpha magnitude (try 1) ; imag(p2) = alpha polar angle (try integers) ; real(p3) = critical point selector (>0 for positive root) ; imag(p3) = unused (<0 for negative root) ; use decomp=256 ; zero and infinity bailouts hardcoded to 1e-6, 1e6 ; coloring speed hardcoded to 4 ; c=pixel, iter=1, n=real(p1), m=imag(p1), nfac=2*n-1 maxr=1e6, minr=1/maxr, speed=4*pi/128 r=real(p2), t=imag(p2), alpha=r*(cos(t)+flip(sin(t))) oln=1/log(n), fac=log(0.5*log(maxr)) c2=sqr(c), hypnum=sqr(n)+sqr(m), pn=1 hypden=sqr(n-m), hypfac=hypnum/hypden if (real(p3)<0) pn=-1 end if if (real(c2)>hypfac) pn=-pn end if if (imag(c)<0) pn=-pn end if afac=c*n, bfac=c2*(n-m)+(n+m), cfac=c*n d=sqrt(bfac*bfac-4*afac*cfac) z=(bfac+pn*d)/(2*afac) : g=(z-c)/(1-c*z), z=alpha*z^n*g^m iter=iter+1, r=|z| ; ; orbit trap around 0 ; renormalize iteration count via decomp angle ; set "iteration done" flag (iter=-1) ; if (r<minr) angle=(iter+oln*(fac-log(log(cabs(z)))))*speed z=cos(angle)+flip(sin(angle)) iter=-1 end if ; ; orbit trap around infinity ; renormalize iteration count via decomp angle ; add pi to angle to separate from 0 orbit trap ; set "iteration done" flag (iter=-1) ; if (r>maxr) angle=(iter+oln*(fac-log(log(cabs(z)))))*speed angle=angle+pi z=cos(angle)+flip(sin(angle)) iter=-1 end if iter>0 } Roger Alexander