Thanks John and David. It looks as if I will try to use python/scipy, unless someone has written wrappers for Sundials for python. Does anybody know if this has been done? John