fractional differential equations. The model is formulated in terms of the left and right Caputo
fractional derivatives. We determine the necessary conditions for the optimality of the
controlled dynamical system. The forward-backward sweep method with generalized euler
scheme is applied to numerically compute the solutions of the optimality system.