In Mathematica documentation one is prompted to use a grid with points at the zeros of the Chebyshev polynomials so that Runge's phenomena arising from
DifferenceOrder->"Pseudospectral" are eliminated.
Although there is an explicit example of how to construct such a grid
CGLGrid[x0_, L_, n_Integer /; n > 1] := x0 + 1/2 L (1 - Cos[π Range[0, n - 1]/(n - 1)]) cgrid = CGLGrid[-5, 10.,16];
I have not yet found how to introduce this grid in
Can anyone post a simple example?