#include <feel/feel.hpp>
Bratu Problem
solve \( -\Delta u + \lambda \exp(u) = 0, \quad u_\Gamma = 0\) on \(\Omega\)