This thesis deals with analysis and implementation of a positivity preserving numerical method for a vaccination model for the transmission dynamics of two HIV subtypes in a given community. The continuous model is analyzed ...