Unequal intervals four-point derivatives

Higher-order methods Chap. 9, Sect. 9.2.2 for multipoint discretisations. The four-point variant with unequal intervals is probably optimal the system can be solved using an extended Thomas algorithm without difficulty. Numerov methods (Sect. 9.2.7) can achieve higher orders with only three-point approximations to the spatial second derivative. They are not trivial to program. 

This program is again a Cottrell simulation using second-order extrapolation based on the Bl (Laasonen) method and unequal intervals, but in contrast with the above program C0TT EXTRAP, this one makes use of the four-point spatial derivative approximation, and the GU-function. It performs a little better than the above program, at little extra programming effort. 

Unequal intervals Chap. 7. These are essential for most programs. The second spatial derivative requires four points if second-order is wanted (and is recommended). With four-point discretisation, an efficient extended Thomas algorithm can be used, obviating the need for a sparse solver. Very few points can then be used across the concentration profile. For two-dimensional simulations, direct three-point discretisation on the unequally spaced grid was shown to be comparable with using transformation and discretisation in transformed space. 

The program LSV4IRC is a simulation of a reversible reaction with input values of p (dimensionless uncompensated resistance) and Yc (dimensionless double layer capacity). Unequal intervals are used, with asymmetric four-point second spatial derivatives, and second-order extrapolation in the time direction. The nonlinear set of six equations for the boundary values is solved by Newton-Raphson iteration. Some results are seen in Chap. 11. 

---