We present a numerical approach for the solution of the dissipative Gross-Pitaevskii equation coupled to the reservoir equation governing the exciton-polaritons Bose-Einstein condensation. It is based on the finite difference method applied to space variables and on the fourth order Range-Kutta algorithm applied to the time variable. Numerical tests illustrate the stability and accuracy of the proposed scheme. Then results on the behavior of the condensate under large Gaussian pumping and around the threshold are presented. We determine the threshold through the particular behavior of the self-energy and characterize it by tracking the establishment time of the steady state.