A set of MapleV R.4/5 software routines for calculating the numerical evolu
tion of dynamical systems and flexible plotting the results is presented. T
he package consists of an initial condition generator (on which the user ca
n impose quite general constraints), a numerical solving manager, plotting
commands that allow the user to locate and focus in on regions of possible
interest and, finally, a set of routines that calculate the fractal dimensi
on of the boundaries of those regions. A special feature of the software ro
utines presented here is an optional interface in C, permitting fast numeri
cal integration using standard Runge-Kutta methods, or variations, for high
precision numerical integration. (C) 1999 Elsevier Science B.V. All rights
reserved.