function res = iter_f_lam(lam)

global N L ep U

  res1 = L*U - lam*U - ep*U.^3;
  res = norm(res1,2);
