Automatic Grid Control in Adaptive BVP Solvers
Grid adaptation in two-point boundary value problems is usually based on mapping a uniform auxiliary grid to the desired nonuniform grid. Here we combine this approach with a new control system for constructing a grid density function I center dot(x). The local mesh width Delta x (j + 1/2) = x (j + 1) -aEuro parts per thousand x (j) with 0 = x (0) < x (1) < ... < x (N) = 1 is computed as Delta x (