fortran66 https://blog.hatena.ne.jp/fortran66/ fortran66のブログ https://fortran66.hatenablog.com/ Fortran2003 半径１の１／４円中に入るか入らないかで円周率を求める古典的な例です。一千万回の試行をしているけれども、組み込み乱数がちゃんとしているかは不明ｗ 試行回数Nに対して、おおむね 1/sqrt(N) に比例して誤差が減るはずなので目安の数字も出力しています。 ソース・プログラム program MonteCarlo implicit none real, parameter :: pi = 4.0 * atan(1.0) integer, parameter :: nn = 10**7 real :: x(nn), y(nn), pp call random_seed() call random_n… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Ffortran66.hatenablog.com%2Fentry%2F20111122%2F1322023222" title="Monte Carlo で円周率 - fortran66のブログ" class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> http://cdn-ak.f.st-hatena.com/images/fotolife/f/fortran66/20111123/20111123133543.gif Hatena Blog https://hatena.blog 2011-11-22 13:40:22 rich https://fortran66.hatenablog.com/entry/20111122/1322023222 1.0 100%