1real function h1(x)
2
3!  sigma=1.0/sqrt(2.0)
4  sigma=1.0
5  xlim=sigma/sqrt(6.0)
6  ax=abs(x)
7  sgnx=1.0
8  if(x.lt.0) sgnx=-1.0
9  if(ax.le.xlim) then
10     h1=x
11  else
12     z=exp(1.0/6.0 - (ax/sigma)**2)
13     h1=sgnx*sqrt(6.0)*sigma*(2.0/3.0 - 0.5*z)
14  endif
15
16  return
17end function h1
18