Ascher, U., R. Mattheij, and R. Russell, Numerical Solution of Boundary Value Problems for Ordinary Differential Equations, SIAM, Philadelphia, PA, 1995, p. 372.
 Cebeci, T. and H. B. Keller, "Shooting and Parallel Shooting Methods for Solving the Falkner-Skan Boundary-layer Equation," J. Comp. Phys., Vol. 7, 1971, pp. 289-300.
 Hairer, E., and G. Wanner, Solving Ordinary Differential Equations II, Stiff and Differential-Algebraic Problems, Springer-Verlag, Berlin, 1991, pp. 5-8.
 Hindmarsh, A. C., "LSODE and LSODI, Two New Initial Value Ordinary Differential Equation Solvers," SIGNUM Newsletter, Vol. 15, 1980, pp. 10-11.
 Hindmarsh, A. C., and G. D. Byrne, "Applications of EPISODE: An Experimental Package for the Integration of Ordinary Differential Equations," Numerical Methods for Differential Systems, L. Lapidus and W. E. Schiesser eds., Academic Press, Orlando, FL, 1976, pp 147-166.
 Ottesen, J. T., "Modelling of the Baroflex-Feedback Mechanism with Time-Delay," J. Math. Biol., Vol. 36, 1997.
 Shampine, L. F., Numerical Solution of Ordinary Differential Equations, Chapman & Hall Mathematics, 1994.
 Shampine, L. F., and M. K. Gordon, Computer Solution of Ordinary Differential Equations, W.H. Freeman & Co., 1975.
 Skeel, R. D. and M. Berzins, "A Method for the Spatial Discretization of Parabolic Equations in One Space Variable," SIAM Journal on Scientific and Statistical Computing, Vol. 11, 1990, pp. 1-32.
 W.H. Enright and H. Hayashi, "The Evaluation of Numerical Software for Delay Differential Equations," R. Boisvert (Ed.), The Quality of Numerical Software: Assessment and Enhancement, Chapman & Hall, London, 1997, pp. 179-192.