;+
; 
;-
FUNCTION imoms3, xin

  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