In this paper, we discussed the transmission of malaria disease in the work of an SIRS-SI model with treatments are given to humans and mosquitoes. Furthermore, we investigated the Homotopy Perturbation Method (HPM) to construct the approximate analytical solution of the simultaneous non linear ordinary differential equations arise in this model. A stability analysis was then performed and MATLAB numerical simulation was provided to clarify the result. These analytical solutions represent a significant simplification of the system’s description allowing easy curve fitting to experiment. Analytical results are compared with simulation results and satisfactory agreement is noted