解常微分方程的fortran90代码。挺好用的。