I am trying to solve a system of first order differential algebraic equations using the subroutine DASPG/DDASPG in IMSL library. To understand how this subroutine works, I copied the example code in the IMSL.PDF file and run the code in the window develop studio. The example is a simple Van der Pol equation and I was hoping to get the exact same results as in the example. But instead, I was prompted with error message: TERMINAL ERROR 7 from DD3SPG, which means the iteration matrix is singular. While the example code produced very smooth results. I wonder why the same exact code cannot reproduce the same results as the example did? Any insight would be greatly appreciated. Thanks
For more complete information about compiler optimizations, see our Optimization Notice.