this paper presents an operational formulation of the tau method based upon orthogonal polynomials by using a reduced set of matrix operations for the numerical solution of nonlinear multi-order fractional differential equations(fdes). the main characteristic behind the approach using this technique is that it reduces such problems to those of solving a system of non-linear algebraic equations....