In this paper, a new numerical method for solving the fractional Riccati differential equation is presented. The fractional derivatives are described in the Caputo sense. The method is based upon fractional-order Bernoulli functions approximations. First, the fractional-order Bernoulli functions and their properties are presented. Then, an operational matrix of fractional order integration is derived and is utilized to reduce the under study problem to a system of algebraic equations. Error analysis included the residual error estimation and the upper bound of the absolute errors are introduced for this method. The technique and the error analysis are applied to some problems to demonstrate the validity and applicability of our method.