This paper presents a two-step Bernstein induced scheme for the numerical solution of some COVID-19 models. The scheme is developed via collocation and interpolation techniques invoked on Bernstein polynomials; the proposed scheme is consistent, between orders four and three. This method can estimate the approximate solution at step points simultaneously by using variable step size. A numerical implementation of the scheme was used on the COVID-19 model, and this showed that the scheme can be conveniently applied to some mathematical models of COVID-19.