- UID
- 263
- 帖子
- 942
- 精华
- 6
- 积分
- 1117
|
在画板中实现这段代码:
init:
pixel2=-abs(real(pixel))+flip(imag(pixel))
d=log(1/pixel2-conj(cos(1/pixel2)))
c=pixel2^3-conj(.1/pixel2-flip(.01/pixel2))
z=abs(pixel2^12)-flip(.01/pixel2)-conj(.1/pixel2)+.3
zmod = (-.7456,-.132) + flip(|d/80*pixel2|)
cmod = |c/7.5| + zmod
c1=1.5*z^1.2,c2=2.25*conj(z)
c3=3.375*flip(z),c4=5.0625*flip(z)
l1=round(real(p1)),l2=round(imag(p1))
l3=round(real(p2)),l4=round(imag(p2))
float bailout=16,int iter=0
loop:
IF (iter==l1)
cmod = |c1/7.5| + zmod, z = 0
ELSEIF (iter==l2)
cmod = |c2/7.5| + zmod, z = 0
ELSEIF (iter==l3)
cmod = |c3/7.5| + zmod, z = 0
ELSEIF (iter==l4)
cmod = |c4/7.5| + zmod, z = 0
ENDIF
z = z^2 + cmod
iter=iter+1
bailout:
|z|<=bailout
default:
title = "Carr 3202"
periodicity = 0
maxiter = 700
magn = 1.4
center = (0,0)
method = multipass
param p1
caption = "Iter Limits 1,2"
default = (200,300)
hint = "0< Re(P1) < Im(P1) < Re(P2) < Im(P2) <maxiter \
for best results"
endparam
param p2
caption = "Iter Limits 3,4"
default = (400,550)
hint = "0< Re(P1) < Im(P1) < Re(P2) < Im(P2) <maxiter \
for best results"
endparam
} |
|