fortran66 https://blog.hatena.ne.jp/fortran66/ fortran66のブログ https://fortran66.hatenablog.com/ Fortran 統計力学などでは、よく階乗の対数を Stirling の公式で近似しますが、どのくらいの精度で近似になっているのか見てみます。ここで Stirling の公式とは、n が大きいところで、 というものです。n のオーダーに比例する項だけが問題になる時は、ここまでで十分となります。計算では、高次の項を [1/n^7] のオーダーまでとった場合の差分も求めています。大体倍精度の場合は、この程度までとれば十分のはずです。 MODULE m_fact IMPLICIT NONE INTEGER, PARAMETER :: kd = SELECTED_REAL_KIND(8) CONTAINS REAL(… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Ffortran66.hatenablog.com%2Fentry%2F20080302%2F1204387713" 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/20080302/20080302010453.gif Hatena Blog https://hatena.blog 2008-03-02 01:08:33 rich https://fortran66.hatenablog.com/entry/20080302/1204387713 1.0 100%