Logistic map & Lorenz attractor 昔のインターフェースを流用して、作図してみます。 fortran66.hatenablog.com fortran66.hatenablog.com 横に縞が入るのが、昔の電送写真のようでいいですね。 プログラム device.f90 module device implicit none type, a…

Monte Carlo で素朴に円周率を計算してみます 1000 点での見積もりなので二桁ほどの精度もありません。 でも点数を増やすと図がつぶれてしまいますｗ 実行結果 Monte Carlo: estimated pi = 3.19600010 ⡟⡛⠛⡛⠛⠛⠛⢿⠭⠭⣩⣉⣙⢉⠉⠙⠉⡩⠉⢉⠋⠉⠍⠉⠍⠝⠉⠉⠉⠉⠉⠉⠉⣙⠝⠉⠉⠍⠉⠉⠉⢙⠉⠩⠉⡉⠉⠉⠉⢹⠀…