In this paper, we study the two-dimensional Bose-Hubbard model with short- and longrange interactions in the canonical ensemble. Using a Variational Monte-Carlo method, we obtain the phase diagrams containing four different phases: Mott insulator, density wave, superfluid and supersolid. The transition lines are determined using the structure factor and the superfluid density. We observe that the phase diagrams for short and long-range interactions are very similar quantitatively and qualitatively, but show also some significant differences. To evaluate the quality of our Variational approach, we compare our results with results published by other groups, obtained with various methods, such as an exact quantum Monte-Carlo algorithm and Gutzwiller Ansatz.