“…For the approximation of the second derivative, we can either use the first derivative operator twice, or use u xx ≈ P −1 (−M +BD)u, where M + M T ≥ 0,B is given by (40), and D is a first-derivative approximation at the boundaries, i.e., 2, 3, . .…”