Introduction - If you have any usage issues, please Google them yourself
Double precision COLSYS subroutine, used to solve ordinary differential linear, nonlinear equations. After expansion version can solve the 1000 highest order for 4 order equation. And achieve satisfactory precision.