z^z^z^z^...

面白いネタ発見!
http://d.hatena.ne.jp/rikunora/20100504/p1
ちょっと計算してみました。 exp(a+ib)=z となるa,bを選び、初期値z, 漸化式 Z_{n+1}=\exp( (a+ib)Z_n) で変換していきます。かなり遠くへ飛んだら発散とみなして計算打ち切り。同じ点に留まるようになったら収束として計算打ち切り。ただし周期 P でP個の点をぐるぐる回る収束もあるのでそれも調べる。絵は発散が黒、周期で収束の場合は周期によって色分け。制限内にどっちにもならない場合はグレー。中心0, 範囲±10を図示してあります。

また初期値zをいろいろ変えて、漸化式はZ_{n+1}=i^ {Z_n}固定でやった場合に収束までの回数で色をつけたジュリア集合はこちら

a=0.3, b=2 のジュリア集合。周期27。白い点がアトラクター。

画像生成コード
http://d.hatena.ne.jp/ita/00010417/p1
http://d.hatena.ne.jp/ita/00010417/p2