Local error control in codes for ordinary differential equations

There are several schemes for the control of local error which are seen in differential equation solvers. The analysis attempts to explain how the selection of a scheme influences the behavior of global error seen in high quality production codes. Two rules of thumb for estimating global errors are given theoretical support when used in conjunction with suitable codes. Substantial numerical experiments support the analysis and conclusions.