ロジスティック写像の分岐図をkeisan.casio.jpにUP!

高精度計算サイトkeisan.casio.jpにロジスティック写像
Xn+1 = a * Xn * (1 - Xn)
の分岐図を描く自作式をUPしましたよ。

ロジスティック写像の分岐図

これは実は久々の自信作。計算は超簡単だけれど、グラフにするのが。

println(...)を使うと引数は8個しかとれなくて、すっかすっかのグラフになってしまう。そこでprint(...)を19個使って、最後にprintln()とすると、20個までの系列のグラフが描ける。本当は100系列くらいまでやるともっとそれっぽいんだけれど、図を見てもらえばわかるようにちゃんと窓も計算できてるし、いいかな。

あとは、精度保証計算をすると、何回も繰り返すと有効数字がめちゃくちゃ減るので最初から桁数は多く取っておくのがよろしいかと。