The solutions to the two coupled 1^{st} order differential
equations:

can be obtained by the Runge-Kutta method, correct to fourth order
terms in x, using

The incremental results are:

These equations can easily be extended to systems of N differential
equations. They can also be used in solving 2^{nd} order
differential equations if the functions F and G are taken as follows: