Any help is really appreciated.
Good days
Andrea
Ps:This is destinate to be a wonderful formula, if i fix the julia.... ;)
SlopeDivideJulibrot{
; draws 4-D slices of DivideBrot Julibrots
; Original code by Jim Muth, Slope version and implementation by Andrea Spinozzi
; October2009
Init:
;Mandelbrot Mode
pix=pixel
u=real(pix)
v=imag(pix)
a=pi*real(p1*0.0055555555555556)
b=pi*imag(p1*0.0055555555555556)
g=pi*real(p2*0.0055555555555556)
d=pi*imag(p2*0.0055555555555556)
ca=cos(a)
cb=cos(b)
sb=sin(b)
cg=cos(g)
sg=sin(g)
cd=cos(d)
sd=sin(d)
aa=-(real(p4)-2)
bb=(imag(p4)+0.00000000000000000000001)
p=u*cg*cd-v*(ca*sb*sg*cd+ca*cb*sd)
q=u*cg*sd+v*(ca*cb*cd-ca*sb*sg*sd)
r=u*sg+v*ca*sb*cg
s=v*sin(a)
c=p+flip(q)+p3
z=r+flip(s)+@startseed
if @mode=="Julia"
pix=pixel
u=real(@startseed)
v=imag(@startseed)
a=pi*real(p1*0.0055555555555556)
b=pi*imag(p1*0.0055555555555556)
g=pi*real(p2*0.0055555555555556)
d=pi*imag(p2*0.0055555555555556)
ca=cos(a)
cb=cos(b)
sb=sin(b)
cg=cos(g)
sg=sin(g)
cd=cos(d)
sd=sin(d)
aa=-(real(p4)-2)
bb=(imag(p4)+0.00000000000000000000001)
p=u*cg*cd-v*(ca*sb*sg*cd+ca*cb*sd)
q=u*cg*sd+v*(ca*cb*cd-ca*sb*sg*sd)
r=u*sg+v*ca*sb*cg
s=v*sin(a)
c=p+flip(q)+p3
z=r+flip(s)+pix
endif
Loop:
z=fn1(z)/(z^(aa)+bb)+c
Bailout:
|z| < @bailout
Default:
title="SlopeDivideJulibrot"
periodicity=0
method=multipass
maxiter=250
Heading
caption="Mandelbrot Mode"
visible=(@mode=="Mandelbrot")
endheading
Heading
caption="Julia Mode"
visible=(@mode=="Julia")
endheading
param mode
caption="Current mode"
default=0
enum="Mandelbrot" "Julia"
visible=false
endparam
param switchmode
caption="Switch to.."
default=1
enum="Mandelbrot" "Julia"
visible=false
endparam
Param bailout
caption="Bailout"
default=10000
min=0.0
Endparam
Param startseed
caption="Start Z/Seed Value"
default=(0.0,0.0)
Endparam
Heading
caption="Parameters"
Endheading
Param p1
caption="Parameter 1"
default=(0.0,0.0)
Endparam
Param p2
caption="Parameter 2"
default=(0.0,0.0)
Endparam
Param p3
caption="Center"
default=(0.0,0.0)
Endparam
Param p4
caption="Exponent"
default=(2.0,0.0)
Endparam
Heading
caption="Functions"
Endheading
Func fn1
caption="Function Z"
default=sqr()
Endfunc
Switch:
type = "SlopeDivideJulibrot"
bailout = @bailout
startseed = pixel
p1 = p1
p2 = p2
p3 = p3
p4 = p4
fn1 = fn1
mode = @switchmode
switchmode = @mode
version = @version
}