I've done some test runs & compared output w/ C version, and it seems to be working in SBCL. I've tested in CLISP also, but it errors out. In other Lisp environments, YMMV. I'm looking into other integration algorithms as well, & I'd like to get together a simple visualization package to animate the results. Comments/patches welcome from any interested parties.