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

0

program-id. logistic.
data division.
working-storage section.
01 i pic 9999.
01 a pic 9v999999.
01 x pic 9v999999 value 0.8.
procedure division.
perform varying a from 2.0 by 0.002 until a > 4
perform varying i from 0 by 1 until i > 100
compute x = a * x * (1 - x)
display a " " x
end-perform
end-perform
stop run.

gnuplot -p -e "plot 'tmp' pt 5 lt 2 ps 0"

gnuplot -p -e "plot 'tmp' pt 5 lt 1 ps 0"

gnuplot -p -e "set size ratio -1; plot pt 5 lt rgb 'blue' ps 0"

Default