erfcc.f90

来自「FORTRANvisualfortran常用数值算法集及源码」· F90 代码 · 共 11 行

FUNCTION erfcc(x)
REAL erfcc,x
REAL t,z
z=abs(x)
t=1./(1.+0.5*z)
erfcc=t*exp(-z*z-1.26551223+t*(1.00002368+t*(.37409196+t*&
(.09678418+t*(-.18628806+t*(.27886807+t*(-1.13520398+t*&
(1.48851587+t*(-.82215223+t*.17087277)))))))))
if (x<0.) erfcc=2.-erfcc
END FUNCTION erfcc