Legendre polynomials' first derivatives have been used as the basis function via the pseudo-Galerkin spectral method. Operational matrices for derivatives have been used and extended to deal with the system of ordinary ...