In this paper, biochemical reaction problem is given in the form of a system of non-linear differential equations involving Caputo fractional derivative. The aim is to suggest an instrumental scheme to approximate the solution of this problem. To achieve this goal, the fractional derivation terms are expanded as the elements of shifted Legendre scaling functions. Then, applying operational matrix of fractional integration and collocation technique, the main problem is transformed to a set of non-linear algebraic equations. This obtained algebraic system can be solved by available standard iterative procedures. Numerical results of applying the proposed method are investigated in details
)