fortran66のブログ

fortran について書きます。

2010-03-03から1日間の記事一覧

Chudonovskyの公式

14桁ずつ数値が決まっていきます。 ■実行結果 収束の様子と真値との差を表示しています。最後に4*ATAN(1)で真値を書き出しています。 ■ソースコード MODULE m_pi IMPLICIT NONE INTEGER, PARAMETER :: kq = KIND(1.0q0) CONTAINS REAL(kq) ELEMENTAL FUNCTIO…