榕坚

奇怪了这UF，虽然面板提示该语句出错，但最终还是能做出图象来。我就试着用它两端取模来替代它，得到的结果好象是一样的：

柳烟

152# 榕坚
此UF在那个地方？

榕坚

153# 柳烟


就是150#那个。

榕坚

Carr2280(YAXIS) {
; Updated for UF by Erik Reckase, March 2000
                  ; Modified Sylvie Gallet frm. [101324,3444],1996
init:
  pixel2=-abs(real(pixel))+flip(imag(pixel)), p0=.2*pixel2
  m=conj(tanh(conj(sinh(pixel2+flip(.0035/pixel2)))))
  c=6*pixel2+flip(.0055/pixel2) \
    -log(pixel2-1/tan(.001/pixel2))^-4.8
  z=2*m+3*pixel2-flip(.0035/pixel2), p5=.225*pixel2
  d1=imag(.0002515/pixel2), d4 = 4*d1
  z1 = 1.5*z+d1, z2 = 2.25*z+d1, z3 = 3.375*z+d1
  z4 = 5.0625*z+d1, z6 = z4^6
  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)
    c=z=z1
  ELSEIF (iter==l2)
    c=z=z2
  ELSEIF (iter==l3)
    c=z=z3
  ELSEIF (iter==l4)
    z=0, c=0
  ELSE
    z = z - p0
  ENDIF
  IF (iter==l4)
    z = z + z4+d4, c = c + z6+d4
  ELSE
    z = z + d4, c = c + d4
  ENDIF
  z = z*z + .975*c - p5
  iter=iter+1
bailout:
  |real(z)| <= bailout
default:
  title = "Carr 2280"
  periodicity = 0
  maxiter = 500
  magn = 6
  center = (0,0)
  method = multipass
  param p1
    caption = "Iter Limits 1,2"
    default = (100,150)
    hint = "0< Re(P1) < Im(P1) < Re(P2) < Im(P2) <maxiter \
            for best results"
  endparam
  param p2
    caption = "Iter Limits 3,4"
    default = (200,250)
    hint = "0< Re(P1) < Im(P1) < Re(P2) < Im(P2) <maxiter \
            for best results"
  endparam
}

柳烟

此分形计算太复杂，我作成后，扫出的图不对劲，差十万八千里。已经疲惫了，空了再检查一下算式，再看看代码。

下载 (10.77 KB)

2011-4-8 20:37

柳烟

我仔细检查，找出一处计算错误，原来我将代码中的乘看成乘方去算了。改后再看看正常不正常，实在没把握，搞了这么久的分形，这个计算的复杂性数第一。

榕坚

156# 柳烟


这个错误的结果也挺有趣的，我一直弄了快两天了，那个IF inter==l4太奇怪了。今天再仔细琢磨好象有点明白了。现在实在没干劲再弄，等一、两天再看看。

柳烟

又纠正了两处计算错误，还是不对劲。

榕坚

这是我第一次弄出来的结果，可与UF一比较细节处差太多了：

榕坚

还有这个，短短的几行代码弄得我晕头转向：
rsp_0711h1 {
  if (cabs(p1)==0), n=8, else, n=(round(cabs(p1))), endif ;points on star
  nm1=n-1, nm2=n-2, a=n*nm1,
  if (cabs(p2)==0), c=pixel, else, c=p2, endif ;slices with high p1
  if (cabs(p3)==0), r=1.5, else, r=p3, endif  ; looked good with 1.0
  z=pixel, oz=f=f1=f2=0
  :
  oz=z, f=z^n-z-c, f1=n*z^nm1-1, f2=a*z^nm2
  z=z-r*(2*f/f1)/(1+sqrt(1-2*f*f2/f1/f1))
  CABS(z-oz)>.0000001
}