3体問題

日常の数学・物理    ーー それでも地球はまわっている



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


Xn+1 = a・Xn(1-Xn)  0<Xn≦1、0<a≦4 , n→∞

Sub Logistic()     'ロジスティック写像

For a = 2.9 To 4 Step 0.0001

  Xn = Rnd()       '初期値を乱数で実施 0~1

For n世代 = 0 To 10000
  Xn1 = a * Xn * (1 - Xn)
  Xn = Xn1
Next n世代
   
  Cells(6 + cnt, 3) = a
  Cells(6 + cnt, 4) = Xn1

  cnt = cnt + 1
Next a

End Sub

※但し、旧エクセルだと取込める散布図のデータ数が非常に低く、アドインなどで拡張する必要がある。