On the backward bifurcation of a vaccination model with non-linear incidence