Solving a system of first order differential equations (predictor corrector)

[C++]

Hello, I currently have working code to solve a single first order differential equation using a predictor-corrector. I need to modify this to solve a system of first order differential equations. Here is the system I need to solve: y'' = -y' + 6y; y(0)=1; y'(0)=-2, on [0,4] And here...