z-(1000*(1-z)/(8-4*z+2*z^2-z^3)+C的N {
init:
z =pixel
loop:
z1=z
if |z|<1e20&&|z-1|>1e-12
z = z-(1000*(1-z)/(8-4*z+2*z^2-z^3)+@c)/ (-1000/(8-4*z+2*z^2-z^3)-(1000*(1-z))*(-4+4*z-3*z^2)/(8-4*z+2*z^2-z^3)^2)
else
z=z1
endif
bailout:
|z-z1| >= @bailout
default:
title = "z-(1000*(1-z)/(8-4*z+2*z^2-z^3)+cN"
helpfile = "Uf*.chm"
helptopic = "Html\formulas\standard\julia.html"
$IFDEF VER50
rating = recommended
$ENDIF
param c
caption = "c"
default = (-.96164, .13619)
hint = "Use this parameter to create many different Julia sets. A good \
way to set this parameter is with the Switch, Eyedropper, or \
Explore features."
endparam
param bailout
caption = "Bailout value"
default = 0.0001
min = 0.0
$IFDEF VER40
exponential = true
$ENDIF
hint = "This parameter defines how soon an orbit bails out while \
iterating. Larger values give smoother outlines; values around 4 \
give more interesting shapes around the set. Values less than 4 \
will distort the fractal."
endparam
switch:
type = "Mandelbrot"
power = power
bailout = bailout
}
在画板中搞，又按胡老师与榕坚老师的提示，编写上面代码，不用et改用em作色，结果所得图皆如下图：
