function F = gaussian_f0(t)
if(t == 0)
    F = 1;
else
    F = 1/2*sqrt(pi/t)*erf(sqrt(t));
end