(元のFORTRANプログラムでは、単精度浮動小数点数を用いていたが、 Cプログラムでは、倍精度浮動小数点数を用いている。 そのため計算結果に違いが生じることがある。 特に反復計算では、EPS の値はかなり小さく取れる。)