Abstract
An explicit unconditionally stable scheme is proposed for solving time-dependent partial differential equations. The application of the proposed scheme is given to solve the COVID-19 epidemic model. This scheme is first-order accurate in time and second-order accurate in space and provides the conditions to get a positive solution for the considered type of epidemic model. Furthermore, the scheme’s stability for the general type of parabolic equation with source term is proved by employing von Neumann stability analysis. Furthermore, the consistency of the scheme is verified for the category of susceptible individuals. In addition to this, the convergence of the proposed scheme is discussed for the considered mathematical model.