In this paper, a computational approach based on the operational matrices in conjunction with orthogonal shifted Legendre polynomials (OSLPs) is designed to solve numerically the multi-order partial differential equations ...