function res = iter_f_psi(u)

  global N L ep LAM

  res = L*u - LAM*u - ep*u.^3;

