mirror of
https://github.com/saitohirga/WSJT-X.git
synced 2024-11-22 20:28:42 -05:00
7d23849786
git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@8310 ab8295b8-cf94-4d9e-aec4-7959e3be5d79
18 lines
298 B
Fortran
18 lines
298 B
Fortran
real function h1(x)
|
|
|
|
! sigma=1.0/sqrt(2.0)
|
|
sigma=1.0
|
|
xlim=sigma/sqrt(6.0)
|
|
ax=abs(x)
|
|
sgnx=1.0
|
|
if(x.lt.0) sgnx=-1.0
|
|
if(ax.le.xlim) then
|
|
h1=x
|
|
else
|
|
z=exp(1.0/6.0 - (ax/sigma)**2)
|
|
h1=sgnx*sqrt(6.0)*sigma*(2.0/3.0 - 0.5*z)
|
|
endif
|
|
|
|
return
|
|
end function h1
|