%----------------------------------------------------------
% Function F0 defined for convenience in computing nuclear
% and two-electron integrals. We define F(0) = 1.
%----------------------------------------------------------
function F = F0(t)
if (t == 0)
    F = 1;
else
    F = 1/2*sqrt(pi/t)*erf(sqrt(t));
end