We consider the weakly asymmetric simple exclusion process on the discrete space {1, ..., n − 1} (n ∈ ), in contact with stochastic reservoirs, both with density ρ ∈ (0, 1) at the extremity points, and starting from the invariant state, namely the Bernoulli product measure of parameter ρ. Under time diffusive scaling t n 2 and for ρ = 1 2 , when the asymmetry parameter is taken of order 1/ n, we prove that the density fluctuations at stationarity are macroscopically governed by the energy solution of the stochastic Burgers equation with Dirichlet boundary conditions, which is shown to be unique and different from the Cole-Hopf solution.