program c
	double precision a,b,c,d
	a=10.788663115743002
	b=0.37935678234628661
	c = derfc(a*b)
	d = erfc(a*b)
	print *, a,b
	print *, c,d
	stop
	end