Generation of higher order pseudospectral integration matrices

A new explicit expression of the higher order pseudospectral integration matrices is presented using an explicit formula for computing iterated integrals of Chebyshev polynomials. Applications to initial value problems, boundary value problems, linear integral and integro-differential equations are presented. The present numerical results are in satisfactory agreement with the exact solutions and show the advantages of this method to some other known methods.