>Magister Scientiae - MSc / This thesis deals with analysis and implementation of a positivity preserving numerical method for a vaccination model for the transmission dynamics of two HIVsubtypes in a given community. The continuous model is analyzed for stability and equilibria. The qualitative information thus obtained is used while designing numerical method(s). Three numerical methods, namely, Implicit Finite Difference Method (IFDM), Non-standard Finite Difference Method (NSFDM) and the Runge-Kutta method of order four (RK4), are designed and implemented. Extensive numerical simulation are carried out to justify theoretical outcomes.