sin(2*pi*b)/b*exp(-0.5)*Hypergeometric1F1[-b/2, 1-b/2, 1/2] where b=0.346