Another quickie today, a Kerry Mitchell formula with gif and jpg
at the usual http://maxitersfractalfollies.blogspot.com.
Also I finally got around to re-sampling the Aug. 16 post.
fract424.gif { ; green tendrils
; blank
; calctime 0:12:31.88
; created Sep 17, 2010
; Fractint Version 2004 Patchlevel 9
reset=2004 type=formula formulafile=kerrym.frm
formulaname=field3_man passes=1
center-mag=-0.12869856067348150/+0.98735067767556690/237.828
params=10000000/7 float=y maxiter=1500 inside=0 decomp=256
periodicity=0
colors=0100200300400500700800900A00B00D00E00F00G00H00J00K00L00M00N00P00Q\
00R00S00T00V00W00X00Y00Z00`00a00b00c00d00f00g00h00i00k00k01k02k03l04l05l\
05m06m07m08n09n0An0Bo0Bo0Co0Dp0Ep0Fp0Gp0Hq0Hq0Iq0Jr0Kr0Lr0Ms0Ns0Ns0Ot0Pt\
0Qt0Ru0Su0Tu0Tu0Uv0Vv0Wv0Xw0Yw0Zw0Zx0_x0`x0ay0by0cy0dz1dz1ez1ez1fz2fz2gz\
2gz2hz3hz3iz3iz3jz4jz4kz4kz4lz5mz5mz5nz5nz6oz6oz6pz6pz7qz7qz7rz7rz8sz8sz\
8tz8uz9uz9uz9tz9tz9tz9sz9sz8sz8rz8rz8rz8qz8qz7qz7pz7pz7oz7oz7oz6nz6nz6nz\
6mz6mz6mz5lz5lz5lz5kz5kz5kz4jz4jz4iz4iz4iz4hz3hz3hz3gz3gz3gz3fz2fz2fz2ez\
2ez2dz1dz1cz1bz1ay1`y1_y1_x1Zx1Yx1Xw1Ww1Vw1Uv1Uv1Tv1Su1Ru1Qu1Pt1Ot1Ot1Nt\
1Ms1Ls1Ks1Jr1Ir1Ir1Hq1Gq1Fq1Ep1Dp1Cp1Co1Bo1Ao19n18n17n16m06l06k06j06i06h\
06g06f05e05d05c05b05a05`05_05Z04Y04X04W04V04U04T04S04R03P03O03N03M03L03K\
03J03I02H02G02F02E02D02C02B02A0190180170160150140130010
}
frm:field3_man { ; Kerry Mitchell 20jan98
;
; Colors Mandelbrot set by nearest approach to
; period 3 field lines
;
; use decomp=256
; real(p1) = bailout (try 1e12)
; imag(p1) = number of field line to use:
; 0, 1, 2, 3, 4, 5, 6, 7
; 3 iterations per pixel
; variable zc used for calculation, z for coloring
;
zc=0, c=pixel, maxr=real(p1), minr=maxr, iter=1
fieldangle=tan(imag(p1)*2*pi/7):
;
; iteration
; compare tangent of polar angle with desired
; field line angle, update minimum if needed
;
iter=iter+3, zc=sqr(zc)+c, zc=sqr(zc)+c, zc=sqr(zc)+c
rzc=|zc|, tanangle=imag(zc)/real(zc),
r=cabs(fieldangle-tanangle)
if (r<minr)
minr=r
end if
;
; bailout
; set "iteration done" flag (iter=-1)
; use log of minimum difference of angles as
; decomposition angle
;
if ((rzc>maxr)||(iter>=maxit))
iter=-1
angle=log(minr)
z=cos(angle)+flip(sin(angle))
end if
iter>0
}
Roger Alexander