Abstract
This paper is devoted to present the approximate solutions with helping of an efficient numerical method for the non-linear coupled system of dynamical marriage model in the fractional of Riemann-Liouville sense (FDMM). The proposed system describes the dynamics of love affair between a couple. The proposed method is dependent on the use of useful properties of the operational matrices of Bernstein polynomials. The operational matrices for the fractional integration in the Riemann-Liouville sense and the product are used to reduce FDMM to the solution of non-linear system of algebraic equations using Newton iteration method. Numerical simulation is given to show the validity and the accuracy of the proposed algorithm. We introduce a comparison with the obtained solution using Runge-Kutta method. (C) 2017 all rights reserved.