柳烟

好象UF中的此分形没有用逃逸时间算法.从baiout看,那个P始终为1

柳烟

这网站不地道，我也试了个，结果差十万八千里。主要是不知 M集的C是不是f(z)+c中的c呢？还是c在其它位置，另就是z到底定位于何处？这组分形庭诱人。

柳烟

我可以试试，简化代码要花时间，因为误删一句或一个字符，可能导致扫不出图。我整出后，将代码帖在此，等待……

柳烟

缩水成这样了，不知是不是你要的效果：
Deviation-Mandelbrot {

Init:

  rotx=@rotx
  roty=@roty
  rotz=@rotz
  rotw=@rotw

    complex z1=real(@start)
    complex z2=imag(@start)
    complex c1=real(pixel)
    complex c2=imag(pixel)

    z=z1+flip(z2)
    c=c1+flip(c2)

    ry=imag(z)*cos(pi*rotx/180)-real(c)*sin(pi*rotx/180)
    rz=imag(z)*sin(pi*rotx/180)+real(c)*cos(pi*rotx/180)
    z=real(z)+ry*1i
    c=rz+imag(c)*1i

    rx=real(c)*sin(pi*roty/180)+real(z)*cos(pi*roty/180)
    rz=real(c)*cos(pi*roty/180)-real(z)*sin(pi*roty/180)
    z=rx+imag(z)*1i
    c=rz+imag(c)*1i

    ry=real(z)*cos(pi*rotz/180)-imag(c)*sin(pi*rotz/180)
    rw=real(z)*sin(pi*rotz/180)+imag(c)*cos(pi*rotz/180)
    z=ry+imag(z)*1i
    c=real(c)+rw*1i

    rx=imag(c)*sin(pi*rotw/180)+imag(z)*cos(pi*rotw/180)
    rw=imag(c)*cos(pi*rotw/180)-imag(z)*sin(pi*rotw/180)
    z=real(z)+rx*1i
    c=real(c)+rw*1i

    z1=real(z)
    z2=imag(z)
    c1=real(c)
    c2=imag(c)

Loop:

  a=z1*z1
  b=z2*z2
  d=z1*z2

  a=(a+flip(real(p4)))^@p6
  b=(b+flip(imag(p4)))^@p7
  d=(d+flip(real(@p8)))^@p10

  z3=(a*p5-b*@p9)*p2
  z4=d*p1
  z5=b*p3


   z1=z3-z5+c1
   z2=z4-z5+c2

   z=z1-z2

Bailout:

  |z1|+|z2|<=@bailout

Default:

title="Deviation-Mandelbrot"
method=multipass
periodicity=0
maxiter=1000

param bailout
  caption="Bailout Value"
  default=128.0
endparam

param rotx
  caption="Rotation x"
  default=0.0
endparam
param roty
  caption="Rotation y"
  default=0.0
endparam
param rotz
  caption="Rotation z"
  default=0.0
endparam
param rotw
  caption="Rotation w"
  default=0.0
endparam

param start
  caption="Start z"
  default=(0.0,0.0)
endparam

param p4
  caption="Parameter 1"
  default=(0.0,0.0)
endparam
param p8
  caption="Parameter 2"
  default=(0.0,0.0)
endparam
param p1
  caption="Multiply 1"
  default=(2.0,0.0)
endparam
param p2
  caption="Multiply 2"
  default=(1.0,0.0)
endparam
param p3
  caption="Multiply 3"
  default=(1.0,0.0)
endparam
param p5
  caption="Multiply 4"
  default=(1.0,0.0)
endparam
param p9
  caption="Multiply 5"
  default=(1.0,0.0)
endparam
param p6
  caption="Exponent 1"
  default=(1.0,0.0)
endparam
param p7
  caption="Exponent 2"
  default=(1.0,0.0)
endparam
param p10
  caption="Exponent 3"
  default=(1.0,0.0)
endparam
param version
  caption="Version"
  default=1.2
  visible=false
endparam

}

柳烟

删掉后，我看了，确实一样。那些旋转量好象是将放倒的M集起扯裂的作用，价值好象不大。

柳烟

请问一下，这M集是不是在网站所给的解析式后加上一个c的M集？z的初始值定位于何处？

柳烟

221# 榕坚
我算出的导函数是：-4(1-z)^2*[z^3-3*z^2+8*z+6]，令其为0,得到z的零点为1,但定位于此点，我在UF中试了，不行。另外我看到三次方程，解起来觉得麻烦，没有去解。零点是由三次方程解得的吗？

柳烟

223# 榕坚
此分形有点怪，阈值过小过大都不行。好象分式函数的M集都有这个特点，阈值要适中。