分形小记

前阵子被何童靴安利了一个知乎帖子：
有没有一段代码，让你觉得人类的智慧也可以璀璨无比？
图形那块美的不要不要~所以基本上一个一个瞅了瞅，所以本文算是个记录吧~

---

一、Mandelbrot分形图案

（图片来自知乎）
右侧那浑圆的大屁股，成功的引起的我的注意~
原文代码主要如下：

float x=0,y=0;
int k;
for(k=0;k++<256;){
    float a=x*x-y*y+(i-768.0)/512;
    x=a;
    if(x*x+y*y>4)break;
}
return log(k)*47;

/*作者：烧茄子
链接：http://www.zhihu.com/question/30262900/answer/48741026
来源：知乎
著作权归作者所有，转载请联系作者获得授权。*/

代码大致是做了这件事：

用 C(Cx,Cy)代表图像上某一点，对参数x=0,y=0作出如下迭代：
x = x*x - y*y + Cx
y = 2 * x * y + Cy
直至x*x+y*y > 4 为止
然后输出迭代所需的次数作为颜色值

代码对三个通道分别作了一些不同的处理，使之更加好看了一丢丢。

刚刚看到这个迭代式其实我是懵逼的……搜了搜终于发现原来就是复数乘法 =_=……

简单而言，Mandelbrot集合是这样的：
对于复平面上的一个点c，和初始点x，对它做如下迭代：
f(x) = x^2 + c
可以得到一个迭代数列。
使这个数列收敛的所有的点，就组成了 Mandelbrot 集合。

所以终于想通了那段代码的含义了~

然后我试图轻微修改下 f(x)，看能不能得到更有趣的结果~

1. f(x) = t * x^2 + c, t为float

结果是trivial的……因为这个参数t其实就相当于这样一个变换：
t*f(x) = (t*x)^2 + t*c
也就是说，只是影响了收敛判断值，以及把c缩小了t倍。
……不贴效果图了……

2. f(x) = x^3 + c

// 迭代函数：x^3 + c
void