Charles Schwartz
A new computational procedure is offered to provide simple, accurate and flexible methods for using modern computers to give numerical evaluations of the various Bessel functions. The Trapezoidal Rule, applied to suitable integral representations, may become the method of choice for evaluation of the many Special Functions of mathematical physics.
Charles Schwartz
We extend the technique of using the Trapezoidal Rule for efficient evaluation of the Special Functions of Mathematical Physics given by integral representations. This technique was recently used for Bessel functions, and here we treat Incomplete Gamma functions and the general Confluent Hypergeometric Function.
Charles Schwartz
A set of algorithms is presented for efficient numerical calculation of the time evolution of classical dynamical systems. Starting with a first approximation for solving the differential equations that has a "reversible" character, we show how to bootstrap easily to higher order accuracy.The method, first shown for a single particle in one dimension, is then neatly extended to many dimensions and many particles.