FnSF              {
                     ; Version 2002 Patchlevel 4
reset=2002 type=formula formulafile=*.frm formulaname=multifractal
ismand=y function=exp/exp/exp passes=t
center-mag=-1.72529/-0.0175526/4.196079
params=-0.1738944669942329/-0.3000579851680044/384.00001/256.0076/256038\
4.00008/1536.0256/768.01536/5120.01024/38888/2 float=y maxiter=3072
inside=maxiter outside=tdis logmap=yes periodicity=0 rseed=-2436
colors=200U00<9>A00<4>000<23>eeegggiiijjjlllnnn<2>sss<25>BBB999777<2>222\
000012<24>0bo0cq0es<3>0kz<25>09C08A068<2>022000211<22>keOmgPoiQqjRslSunT\
<2>zsX<25>CB7A96875<3>000400<23>o00q00s00<3>z00<15>W00
}
FRM:Multifractal     { ; Albrecht Niekamp Apr, 03
;only integers as first input, second input 5 digits to the right
;real(p2) factor1,border1
;imag(p2) (-)maxiter1,input bas.frm1_2digits input bas.frm2_2digits
;real(p3) bailout,hyb.mandel: number of sides (Many_mods L.Allison)
;imag(p3) factor2,border2
;real(p4) maxiter2,factor3
;imag(p4) border3,maxiter3
;real(p5) 1digit_shape only 5digits: shape,out1,out2,in1,in2
; 1_mand 2_jul 3_hyb.mand 4_hyb.jul +5_reset (not shape) <0_no in/out1
; optional basic frm2: out-nbr 1digit in-nbr 1digit 3_both +5_shape
;imag(p5) 1digit : 1_iter1reset 2_iter2reset 3_both +5_shape effect
; 2digits : many_mods multipl. shutoff +10_in/out1 +20_inout2 +30_both
; optional input2: frequency_2digits level_2digits amplitude_2digits
;
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)
if (dd<5)
start=0
endif
vb=(dd==3)+(dd==4)
da=round((da-dd)*100)
d=trunc(dd/10000)
vb=(d==3)+(d==4)+vb
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-d*10
d=trunc(dd)
inr2=(d>5)
d=d-5*inr2
v3=(d==3)
v4=(d==4)
d2=(d==2)+(d==4)
d=trunc(da/10)
ab=(d>=5)
d=d-5*ab
da=da-10*d
bb1=(d==1)
bb2=(d==2)
d=da
ab=(d>=5)
d=d-5*ab
ab1=(d==1)
ab2=(d==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 by S.Gallet
d=d-freq*1000
level=(trunc(d/100))/100            ;water level
d=d-level*10000
ampl=d/100                          ;amplitude of the wave
angle=real(rotskew*pi/180)
t=exp(-flip(angle))
h=1/real(magxmag)
q=h/0.75*imag(magxmag)
tanskew=tan(imag(rotskew*pi/180))
w=2*q*t
v=2*h*(tanskew+flip(1))*t
z3rd=center+(-q-h*tanskew-flip(h))*t
z=pixel-z3rd
s=imag(conj(w)*z)/imag(conj(w)*v)
if (s<=level)
   da=level-s
   z=z+2*da*(1+ampl*sin(freq*da^0.2))*v
endif
z=z+z3rd
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))       ;frm: Albrecht Niekamp
else
z2=fn1(z)+c                         ;frm: Many_mods Linda Allison
sqz2=cos(z2)
z=c*(1-sqz2)/(1+sqz2)
endif
else if (ab)
z=z*z+c+c*c-p7                     ;frm: Pusk s Istv n
else
z2=z*z                             ;frm: Pusk s Istv n modified
z=z2*z2+(p6+u)*z2+c-p1*(d3==0)
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 (start==1)
if (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
endif
bo<=ba
}