Automatic Grid Control in Adaptive BVP Solvers
Modern adaptive techniques in two-point boundary value problems generate the mesh by constructing a function that maps a uniform grid to the desired nonuniform grid. This paper describes a new control algorithm for constructing a grid density function $phi(x)$, such that the local mesh width $Delta x_{j+1/2}=x_{j+1}-x_j$ is computed as $Delta x_{j+1/2} = varepsilon_N / varphi_{j+1/2}$. Here $varep
