The domain of attraction of a nonlinear differential equations is the region of initial points of solution tending to the equilibrium points of the systems as the time going. Determining the domain of attraction is one of the most important problems to investigate nonlinear dynamical systems. In this article, we first present two algorithms to determine the domain of attraction by using the moment matrices. In addition, as an application we consider a class of SIRS infection model and discuss asymptotical stability by Lyapunov method, and also estimate the domain of attraction by using the algorithms.