Hi, this is a major update for the Multifractal formula. The inputs are a lot easier and a problem is solved that in some cases caused that the inside options were not executed. This means, this update is not completely compatible, so dont remove the previous version. To avoid confusion it has a different name. The images all use the same shape, but with different insides. A lot of interesting border zooms are guaranteed. Janusbird3 uses passes=1 as passes=t causes coloring errors. Fractal Greetings, Al. --------------------------------start-par---------------------------------- Janusbird { ; Albrecht Niekamp May, 03 ; Time 0.39.46.32 ; Version 2002 Patchlevel 5 reset=2002 type=formula formulafile=mult239.frm formulaname=multifractalx ismand=y function=exp/sin/exp passes=t center-mag=-0.488627/4.44089e-016/1.086909 params=-0.2724234748374889/0.09260841170049135/384.00001/255.4505/2048.0\ 0008/768.0256000000001/768.01536/2560.01025/12297.72/0 float=y maxiter=2048 inside=maxiter outside=tdis logmap=4 periodicity=0 rseed=-2436 colors=200pppqqqsss<25>BBB999777<2>222000012<24>0bo0cq0es<3>0kz<25>09C08\ A068<2>022000211<22>keOmgPoiQqjRslSunT<2>zsX<25>CB7A96875<3>000400<23>o0\ 0q00s00<3>z00<25>C00A00800<3>000<23>eeegggiiijjjlllnnn } Janusbird2 { ; Julia as inside ; Albrecht Niekamp May, 03 ; Version 2002 Patchlevel 5 reset=2002 type=formula formulafile=mult239.frm formulaname=multifractalx ismand=y function=exp/exp/exp passes=t center-mag=-0.488627/4.44089e-016/1.086909 params=-0.2724234748374889/0.09260841170049135/384.00001/255.4505/2048.0\ 0008/768.0256000000001/768.01536/2560.01025/12224.71/0 float=y maxiter=2048 inside=maxiter outside=tdis logmap=4 periodicity=0 rseed=-2436 colors=200pppqqqsss<25>BBB999777<2>222000012<24>0bo0cq0es<3>0kz<25>09C08\ A068<2>022000211<22>keOmgPoiQqjRslSunT<2>zsX<25>CB7A96875<3>000400<23>o0\ 0q00s00<3>z00<25>C00A00800<3>000<23>eeegggiiijjjlllnnn } Janusbird3 { ; Inside many_mods and julia mix ; Albrecht Niekamp May, 03 ; Version 2002 Patchlevel 5 reset=2002 type=formula formulafile=mult239.frm formulaname=multifractalx ismand=y function=exp/exp/exp passes=1 center-mag=-0.488627/4.44089e-016/1.086909 params=-0.2724234748374889/0.09260841170049135/384.00001/255.4505/2048.0\ 0016/768.0256000000001/768.01536/2560.01025/12284.71/11 float=y maxiter=2048 inside=maxiter outside=tdis logmap=4 periodicity=0 rseed=-2436 colors=200pppqqqsss<25>BBB999777<2>222000012<24>0bo0cq0es<3>0kz<25>09C08\ A068<2>022000211<22>keOmgPoiQqjRslSunT<2>zsX<25>CB7A96875<3>000400<23>o0\ 0q00s00<3>z00<25>C00A00800<3>000<23>eeegggiiijjjlllnnn } ----------------------------------------start-frm-------------------------- FRM:Multifractalx { ; Albrecht Niekamp Update 260503 ;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 digit1: 1_out1 2_out2 3_both +5_shape ; digit2: 1_in1 2_in2 3_both ;imag(p5) digit1 : many_mods multipl.shutoff 1_in/out1 2_in/out2 3_both ; digit2 :1_iter1reset 2_iter2reset 3_both +5_shape effect ; 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) 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) dd=dd-d*1000 or1=d>5 d=d-5*or1 vv1=(d==3) vv2=(d==4) dd1=(d==2)+(d==4) d=trunc(dd/100) dd=dd-d*100 or2=d>5 d=d-5*or2 vv3=(d==3) vv4=(d==4) dd2=(d==2)+(d==4) d=trunc(dd/10) dd=dd-d*10 inr1=d>5 d=d-5*inr1 v1=(d==3) v2=(d==4) d1=(d==2)+(d==4) inr2=dd>5 dd=dd-inr2*5 v3=(dd==3) v4=(dd==4) d2=(dd==2)+(dd==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 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 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 } -------------------------------------end-------------------------------