45 real(dp),
intent(in) :: x
50 t = 1.0_dp/(1.0_dp+0.5_dp*z)
52 erfcc = t * exp( -z*z -1.26551223_dp &
53 + t * ( 1.00002368_dp &
54 + t * ( 0.37409196_dp &
55 + t * ( 0.09678418_dp &
56 + t * ( -0.18628806_dp &
57 + t * ( 0.27886807_dp &
58 + t * ( -1.13520398_dp &
59 + t * ( 1.48851587_dp &
60 + t * ( -0.82215223_dp &
61 + t * 0.17087277_dp ) ) ) ) ) ) ) ) )
Declarations of kind types for SICOPOLIS.
real(dp) function erfcc(x)
Computation of the complementary error function erfc(x) = 1-erf(x) with a fractional error everywhere...