In this paper, a modern method is presented to solve a class of fractional optimal control problems (FOCPs) indirectly. First, the necessary optimality conditions for the FOCP are obtained in the form of two fractional differential equations (FDEs). Then, the unknown functions are approximated by the hybrid functions, including Bernoulli polynomials and Block-pulse functions based o...