Generation of pseudospectral differentiation matrices I

Research output: Contribution to journalArticlepeer-review

84 Scopus citations


We present a simple method for computing n × n pseudospectral differentiation matrices of order p in script type O:(pn2) operations for the case of quasi-polynomial approximation. The algorithm is based on Fornberg's finite difference algorithm and is numerically stable. A Fortran implementation is included. A necessary and sufficient condition for Dp = D1p is also given.

Original languageEnglish (US)
Pages (from-to)1640-1657
Number of pages18
JournalSIAM Journal on Numerical Analysis
Issue number4
StatePublished - 1997


  • Collocation polynomials
  • Pseudospectral methods

ASJC Scopus subject areas

  • Numerical Analysis
  • Computational Mathematics
  • Applied Mathematics


Dive into the research topics of 'Generation of pseudospectral differentiation matrices I'. Together they form a unique fingerprint.

Cite this