erfc.f   [plain text]

c============================================== test.f
                real x, y
		real*8 x1, y1
		x=0.
		y = erfc(x)
		if (y .ne. 1.) call abort

		x=1.1
		y = erfc(x)
		if (abs(y - .1197949) .ge. 1.e-6) call abort

* modified from x=10, y .gt. 1.5e-44 to avoid lack of -mieee on Alphas.
		x=8
		y = erfc(x)
		if (y .gt. 1.2e-28) call abort

		x1=0.
		y1 = erfc(x1)
		if (y1 .ne. 1.) call abort

		x1=1.1d0
		y1 = erfc(x1)
		if (abs(y1 - .1197949d0) .ge. 1.d-6) call abort

		x1=10
		y1 = erfc(x1)
		if (y1 .gt. 1.5d-44) call abort
		end
c=================================================
!output:
!  0.  1.875
!  1.10000002  1.48958981
!  10.  5.00220949E-06
!
!The values should be:
!erfc(0)=1
!erfc(1.1)= 0.1197949
!erfc(10)<1.543115467311259E-044