Collatz問題 #26

#7の散布図と52527からの経路を重ね合わせてみる。
散布図の方のプロットが無い部分の表示はカットする。
散布図のためのデータを作成する処理

$ ./collatz5 > i65536.dat

は#7と同じである。

m = 339
plot [1:100000] 'i65536.dat' every 2 pt 7 ps 0.5, '' every 2::1 lt 3 pt 7 ps 0.5, 'p52527.dat' using 1:(up($2, $3, m)) w l lt 1, '' using 1:(down($2, $3, m)) w l lt 3

初期値97からの経路についても同様の表示を行ってみる。

m = 118
plot [1:10000] 'i100.dat' every 2 pt 7 ps 0.5, '' every 2::1 lt 3 pt 7 ps 0.5, 'p97.dat' using 1:(up($2, $3, m)) w l lt 1, '' using 1:(down($2, $3, m)) w l lt 3

i100.datは#5で作成した初期値100以下の散布図用のデータである。
もしくは、i65536.datの先頭100行に相当するので、そこから切り出してもよい。

さらに表示範囲を狭める。

m = 118
plot [1:100] 'i100.dat' every 2 pt 7, '' every 2::1 lt 3 pt 7, 'p97.dat' using 1:(up($2, $3, m)) w l lt 1, '' using 1:(down($2, $3, m)) w l lt 3