In this thesis, a Non-Standard Finite Difference (NSFD) method based on the θ− method is formulated for solving the mathematical model of the Beddington-DeAngelis functional response predator-prey dynamical system. The ...