source: trunk/SRC/Interpolation/imoms3.pro

;+
;
; @param xin {in}{required}
;
; @version
; $Id$
;
;-
FUNCTION imoms3, xin
;
  compile_opt idl2, strictarrsubs
;

  x = abs(xin)
  y = fltarr(n_elements(x))

  test1 = where(x LT 1)
  IF test1[0] NE -1 THEN BEGIN
    xtmp = x[test1]
    y[test1] = 0.5*xtmp*xtmp*xtmp -xtmp*xtmp -0.5*xtmp + 1
  ENDIF
  test1 = where(x LT 2 AND x GE 1)
  IF test1[0] NE -1 THEN BEGIN
    xtmp = x[test1]
    y[test1] = (-1./6.)*xtmp*xtmp*xtmp + xtmp*xtmp - (11./6.)*xtmp + 1
  ENDIF

  RETURN, y
END