FUNCTION it_gausshist,image

  n_0=127.5
  s=256/6

  hist=HIST_EQUAL(image)
  g_x=EXP(-(hist-n_0)^2/(2*s^2))
  print,MIN(hist),MAX(g_x)
  g2_x=1/(ERRORF(g_x))
  print,MIN(g2_x),MAX(g2_x)
  img=hist(g2_x)
  print,MIN(img),MAX(img)

RETURN,img
END