;
PRO hnmf2, X, Y, XOUX, YOUX
;
; Find the NmF2 and hmF2 based on the electron profile
;
; input
;     X: height
;     Y: Ne profile
;
; output
;     XOUX: hmF2
;     YOUX: NmF2
;
nz=n_elements(x)
for k=nz-2,1,-1 do begin
   if ((y[k] ge y[k-1]) and (y[k] ge y[k+1])) then break
endfor
h=[x[k-1],x[k],x[k+1]]
p=[y[k-1],y[k],y[k+1]]

h12=h(0)*h(0)
h22=h(1)*h(1)
h32=h(2)*h(2)
DELTX=h12*h(1)+h22*h(2)+h32*h(0)-h32*h(1)-h12*h(2)-h22*h(0)
ATX=p(0)*h(1)+p(1)*h(2)+p(2)*h(0)-h(1)*p(2)-h(2)*p(0)-h(0)*p(1)
AX=ATX/DELTX
BTX=h12*p(1)+h22*p(2)+h32*p(0)-h32*p(1)-h12*p(2)-h22*p(0)
BX=BTX/DELTX
CTX=h12*h(1)*p(2)+h22*h(2)*p(0)+h32*h(0)*p(1)-h32*h(1)*p(0)-h12*h(2)*p(1)-h22*h(0)*p(2)
CX=CTX/DELTX
XOUX=-(BX/(2.*AX))
YOUX=-((BX*BX-4.*AX*CX)/(4.*AX))

RETURN
END