In this paper, a new definition for the fractional order operator called the Caputo-Fabrizio (CF) fractional derivative operator without singular kernel has been numerically approximated using the two-point finite forward difference formula for the classical first-order derivative of the function f (t) appearing inside the integral sign of the definition of the CF operator. Thus, a numerical differentiation formula has been proposed in the present study. The obtained numerical approximation was found to be of first-order convergence, having decreasing absolute errors with respect to a decrease in the time step size h used in the approximations. Such absolute errors are computed as the absolute difference between the results obtained through the proposed numerical approximation and the exact solution. With the aim of improved accuracy, the two-point finite forward difference formula has also been utilized for the continuous temporal mesh. Some mathematical models of varying nature, including a diffusion-wave equation, are numerically solved, whereas the first-order accuracy is not only verified by the error analysis but also experimentally tested by decreasing the time-step size by one order of magnitude, whereupon the proposed numerical approximation also shows a one-order decrease in the magnitude of its absolute errors computed at the final mesh point of the integration interval under consideration.