New Runge-Kutta algorithms are developed which determine the solution of a system of ordinary differential equations at any point within a given integration step, as well as at the end of each step.

A new class of stochastic Runge–Kutta methods for the weak approximation of the solution of Itô stochastic differential equation systems with a multidimensional Wiener process is introduced. As the ...