Расчетные формулы метода Рунге-Кутта 4-го порядка для системы ОДУ 1-го порядка:
где 
;
;
;
m – количество узлов;
– номер функции;
– номер узла;
;
;
;
.
На рис. 8.1 – 8.3 представлен алгоритм решения системы обыкновенных дифференциальных уравнений.
Рис. 8.1. Алгоритм решения системы ОДУ.
Рис. 8.2. Алгоритм расчета новой точки методом Рунге-Кутта 4-го порядка:
Рис. 8.3. Схема процедуры вычисления правой части системы ОДУ 1-го порядка:
