In this paper, the problem of the spread of a non-fatal disease in a population is solved by using the Hermite collocation method. Mathematical modeling of the problem corresponds to a three-dimensional system of nonlinear ODEs. The presented scheme reduces the problem to a nonlinear algebraic equation system by expanding the approximate solutions by using Hermite polynomials with unknown coefficients. These coefficients of the Hermite polynomials are computed by using the matrix operations of derivatives together with the collocation method. Maple software is used to carry out the computations. In addition, comparison of our method with the Homotopy perturbation method (HPM) and Laplece-Adomian decomposition method (LADM) proves accuracy of solution.