(z^2-Q^2)(z^2-c^2)系列分形

lnszdzg

前几天，老师们完了(z-1)(z+1)(z-b+a)(z-b-a) 的牛迭M集征解(越玩越有趣)
今天，在Mathcad中也实现了这个，而且发现：这一系列的结构是相似的，可以统一在一个分形之中。

下载 (14.48 KB)

2015-1-28 21:49

下载 (12.52 KB)

2015-1-28 21:49

下载 (17.06 KB)

2015-1-28 21:49

下载 (19.34 KB)

2015-1-28 21:49

下载 (17.87 KB)

2015-1-28 21:49

下载 (18.75 KB)

2015-1-28 21:49

下载 (17.72 KB)

2015-1-28 21:49

下载 (17.75 KB)

2015-1-28 21:50

下载 (18.83 KB)

2015-1-28 21:49

下载 (17.78 KB)

2015-1-28 21:49

lnszdzg

而且，还有一些很有意思的变种

下载 (15.04 KB)

2015-1-28 22:01

下载 (18.69 KB)

2015-1-28 22:01

下载 (20.18 KB)

2015-1-28 22:01

应剑客之约，让我邀请大家到Mathcad吧做客。
希望几位老师能到Mathcad吧看看
http://tieba.baidu.com/f?kw=mathcad&fr=ala0

lnszdzg

下载 (14.9 KB)

2015-1-28 22:03

UF变种

lnszdzg

下载 (18.12 KB)

2015-1-28 22:04

下载 (17.3 KB)

2015-1-28 22:04

下载 (25.46 KB)

2015-1-28 22:04

lnszdzg

下载 (30.65 KB)

2015-1-28 22:44

柳烟

杜老兄越玩越高超，快入化境了，祝贺。

lnszdzg

下载 (51.17 KB)

2015-1-29 15:52

也是变种

lnszdzg

UF代码，很简单，却很有效哈
Newton_Mandelbrot1 {
init:
  complex Q = @Q
  int p = @power
  c = #pixel
  z = sqrt((Q^2+c^p)/6 )
oop:
  err = (z^4-(c^p+Q^2)*z^2+c^p)/(4*z^3-2*(c^p+Q^2)*z)
  z = z - err
bailout:
  |err| > @bailout
default:
  title = "Newton_Mandelbrot1"
  maxiter = 2000
  periodicity = 0
  center = (0,0)
  magn = 1.0
param power
    caption = "power"
    default = 1
  endparam
  param Q
    caption = "Q"
    default = (1.0,0.0)
  endparam
param bailout
    caption = "Bailout"
    default = 0.000001
endparam
}

lnszdzg

“蝴蝶结”的个数有p确定，如果增加变换，得到的图形也很美。

lnszdzg

邀请大家到Mathcad吧做客。
http://tieba.baidu.com/f?kw=mathcad&fr=ala0