Abstract
In this paper, a new HIV CD4+T cells reaction-diffusion model in two dimensions has been introduced. Two novel and efficient positivity preserving finite difference schemes for the numerical solution has been used. The positivity property is of great importance in epidemic models because negative values have no meaning. The stability and consistency of the proposed positivity preserving schemes have been discussed briefly. A comparison of the proposed schemes with an extensively used Euler scheme has been provided. The numerical simulations of all the schemes have been presented with the help of a numerical test and found that the Euler method shows false behavior which is not a part of the continuous system. Moreover, both proposed positivity preserving schemes illustrate the behavior which is consistent with the continuous system.