I read that for an open, bounded domain $\Omega \in \mathbb R^n$ the Poisson equation $$-\Delta u = f$$ with given boundary data in fact does always exist and is unique for Dirichlet boundary data or unique up to a constant for Neumann boundary data.
Why do we even bother a variational formulation to obtain a weak solution?