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