I'm looking for an NL ODE solver to run in MS

fortran Version 5.1 or

for Fortran Powerstation, if available. I don't have an IMSL library, so

that's why I am looking. Thanks to anyone who can help. Oh, I also need

any directions associated with using this subroutine.