Hi, I love this one because it keeps you watching. Albrecht -------------------------start par------------------------- Jrula2 { ; Jrula2 ; structures+overlays+partial symmetry ; frm+image Albrecht Niekamp ; source+colours Lee Skinner ; time 11.29.00 ; Fractint Version 2004 Patchlevel 11 reset=2004 type=formula formulafile=mfr.frm formulaname=multifractal ismand=n function=exp/exp/exp passes=t center-mag=+1.09441249371272300/+0.21759644730115380/76.\ 36852/1/-102.499999999999915/-3.17669501814776822e-013 params=-0.1740997192297128/-0.2996141087069308/384.00004\ /256.1106/2048.00016/768.00009/768.01536/256.01025/28424\ .5/12 float=y maxiter=2048 inside=maxiter outside=tdis logmap=10 symmetry=none periodicity=0 rseed=-2436 cyclerange=2/255 colors=000000t2tv2u<3>c5eZ6`U6X<2>F9LAAG7FO4KW0Qc<3>5Xc6\ Zc7`c7aa7c_7eY<3>FcgHciJbl<3>S_v<3>civflvinvlqvotvsww<2>\ zzz<3>rwzovzhsw<3>Gch8_d8Xa<3>8JP8FL8CI<4>XCAaB9fB7<3>zA\ 0<13>SA1QA1NA1<3>D92<6>UXAW`BYcC<3>hrI<3>UnXRm_NlcJkgFfp\ Aaz<3>PJoSFlWAi_5fc0ce7_jFRoNIzb0<2>zz1zzC<2>zzz<3>jqdfo\ _blU<3>Mb8<3>LQ8KM8KJ8<2>I88H44G00<3>TA5WD6ZF7aI8eLAhMBi\ ODjQF<3>qeBriAtm9vq8xu7zz5zz4zz3<8>zY4zV5zS5<3>zE6wE6tD6\ <10>w8Px8Rx7S<3>z5_<3>f5Ma5IW4E<2>_4H<3>P3DN2CK2A<3>805G\ 0BO0G<3>g0Ul0Xq0`v0d<2>y0fz0gy0h<5>s1or2qs2r } ----------------------------start frm---------------------- Frm:multifractal { ; Albrecht Niekamp 24.3.03 ; d=real(p2) fac1=trunc(d) b1=(d-fac1)*100000 ; d=imag(p2) da=d<0 if (da) d=-d endif mi1=trunc(d) d=round((d-mi1)*10000) p6=trunc(d/100)/10 d=d-1000*p6 if (da) p6=-p6 endif p7=d/10 ; d=real(p3) ba=trunc(d) mm=round((d-ba)*100000) ; d=imag(p3) fac2=trunc(d) b2=(d-fac2)*100000 ; d=real(p4) mi2=trunc(d) fac3=round((d-mi2)*100000) ; d=imag(p4) b3=trunc(d) mi3=round((d-b3)*100000) ; da=real(p5) t=da<0 if (t) da=-da endif start=2-t dd=trunc(da) da=round((da-dd)*100) if (dd<5) start=0 vb=(dd==3)+(dd==4) else d=trunc(dd/10000) vb=(d==3)+(d==4) dd=dd-d*10000 d=trunc(dd/1000) or1=d>5 d=d-5*or1 vv1=(d==3) vv2=(d==4) dd1=(d==2)+(d==4) dd=dd-d*1000 d=trunc(dd/100) or2=d>5 d=d-5*or2 vv3=(d==3) vv4=(d==4) dd2=(d==2)+(d==4) dd=dd-d*100 d=trunc(dd/10) inr1=d>5 d=d-5*inr1 v1=(d==3) v2=(d==4) d1=(d==2)+(d==4) dd=dd-10*d d=trunc(dd) inr2=d>5 d=d-5*inr2 v3=(d==3) v4=(d==4) d2=(d==2)+(d==4) endif ; d=trunc(da/10) da=da-10*d ab=d>=5 d=d-5*ab bb1=(d==1) bb2=(d==2) dd=da>=5 da=da-5*dd ab=ab+dd ab1=(da==1) ab2=(da==2) ; d=imag(p5) dd=trunc(d) d=round((d-dd)*1000000) le=d>10000 da=trunc(dd/10) so1=((da==1)+(da==3)==0) so2=((da==2)+(da==3)==0) dd=dd-10*da wo=dd>=5 dd=dd-5*wo ir1=(dd==1)+(dd==3) ir2=(dd==2)+(dd==3) if (le) freq=round((trunc(d/10000))*10) ;frequency_lake effect d=d-freq*1000 level=(trunc(d/100))/100 ;water level d=d-level*10000 ampl=d/100 ;amplitude of the wave u=real(rotskew*pi/180) t=exp(-flip(u)) bo=1/real(magxmag) q=bo/0.75*imag(magxmag) dd=tan(imag(rotskew*pi/180)) d3=2*q*t rs=2*bo*(dd+flip(1))*t z2=center+(-q-bo*dd-flip(bo))*t z=pixel-z2 d=imag(conj(d3)*z)/imag(conj(d3)*rs) if (d<=level) da=level-d z=z+2*da*(1+ampl*sin(freq*da^0.2))*rs endif z=z+z2 endif if (ismand) if (vb) c=0.4*log(sqr(pixel^mm)) z=z*le else c=z*le+pixel*(le==0) z=0 endif else c=p1 z=z*le+pixel*(le==0) endif t=0 d3=(ismand==0) u=d3 : if (vb) if (d3) z=(z*z-real(c))^(z*u+imag(c)) else z2=fn1(z)+c ;Linda Allison q=cos(z2) z=c*(1-q)/(1+q) endif elseif (ab) z=z*z+c+c*c-p7 ;Pusk s Istv n elseif (d3) z2=z*z z=z2*z2+(p6+u)*z2+c else z2=z*z ;Pusk s Istv n z=z2*z2+(p6+u)*z2+c-p1 endif bo=|z| If (wo) t=t+1 if ((t<mi1)+(bo>b1)) u=2*(fn1(t/fac1)) if (ismand) z=z*u if (vb) c=0.4*log(sqr(pixel^mm)) else c=z*le+pixel*(le==0) endif else z=z*le+pixel*(le==0) c=p1*u endif wo=0 if (ir1) t=0 endif endif elseif (start) t=t+1 if (start==2) if (bo>b2) u=2*(fn2(t/fac2)) rs=or1 ab=bb1 if (dd1) if (rs) z=pixel endif c=p1*u vb=vv2 d3=vb else d3=0 z=z*u c=z if (rs) z=0 endif vb=vv1 if (vb) if (so1) mm=mm+mm endif c=0.4*log(sqr(pixel^mm)) endif start=1 endif elseif (t>mi2) u=2*(fn2(t/fac2)) rs=inr1 ab=ab1 if (d1) if (rs) z=pixel endif c=p1*u vb=v2 d3=vb else d3=0 z=z*u c=z if (rs) z=0 endif vb=v1 if (vb) if (so1) mm=mm+mm endif c=0.4*log(sqr(pixel^mm)) endif endif start=1 if (ir2) t=0 endif endif elseif (bo>b3) u=2*(fn3(t/fac3)) ab=bb2 rs=or2 if (dd2) if (rs) z=pixel endif c=p1*u vb=vv4 d3=vb else d3=0 z=z*u c=z if (rs) z=0 endif vb=vv3 if (vb) if (so2) mm=mm+mm endif c=0.4*log(sqr(pixel^mm)) endif start=0 endif elseif (t>mi3) u=2*(fn3(t/fac3)) ab=ab2 rs=inr2 if (d2) if (rs) z=pixel endif c=p1*u vb=v4 d3=vb else d3=0 z=z*u c=z if (rs) z=0 endif vb=v3 if (vb) if (so2) mm=mm+mm endif c=0.4*log(sqr(pixel^mm)) endif endif start=0 endif endif bo<=ba }