fortran66 https://blog.hatena.ne.jp/fortran66/ fortran66のブログ https://fortran66.hatenablog.com/ Fortran2003 ローレンツ型の関数を‐1から1の範囲で、いくつかの分割数に対して台形公式で数値積分します。出てきた積分値をリチャードソン補外をして改善します。その結果は理屈上シンプソン法によって求めたものと一致します。 実行結果 まずもっとも細かい分割数での座標値を配列生成子で生成し、被積分関数に入れて値を配列として求めます。ここで被積分関数を要素型で宣言することで、引数と関数値のDO LOOPの交換をします。 『y = f( [(x0 + i * h, i = 0, n)] )』 円周率πが求まるはずです。 『integer, parameter :: kd = kind(0.0e0)』の精度を変えることで… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Ffortran66.hatenablog.com%2Fentry%2F20111123%2F1322043773" title="台形公式をリチャードソン補外でシンプソン法に。 - 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/20111123191637.gif Hatena Blog https://hatena.blog 2011-11-23 19:22:53 rich https://fortran66.hatenablog.com/entry/20111123/1322043773 1.0 100%