Численные методы для систем дифференциальных уравнений

1) найти численное решение системы уравнений

2)записать значения времени и фазовых координат для каждого шага интегрирования

3)Построить график изменения входной величины во времени на основании записанных данных

Нужно очень подробное решение на бумаге, которое потом можно было бы обратить в код или сразу код на языке Си. Вариант 14 (строка с номером 14)