R/erfc.R

In NORMT3: Evaluates Complex Erf, Erfc, Faddeeva, and Density of Sum of Gaussian and Student's t

"erfc" <-
function(z){

ans <- .C("IPerfcvec",
		x=as.double(Re(z)),
		y=as.double(Im(z)),
		ansx=as.double(Re(z)),
		ansy=as.double(Im(z)),
		n = as.integer(length(z)),
		error = as.integer(0),
		PACKAGE="NORMT3")
if (ans$error != 0)
	stop(paste("Error code from TOMS 680 was ", ans$error))

return(complex(real=ans$ansx, imaginary=ans$ansy))

}