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

ВХОДНЫЕ ДАННЫЕ:
N-число уравнений; h-шаг интегрирования; t-начальное
время; y[i]-массив начальных значений переменных
(i=0,1,2,...,N); функция вычисления правых частей
дифференциальных уравнений-f;
ВЫХОДНЫЕ ДАННЫЕ:
Результат после одного шага решения в массиве начальных
значений переменных y[i];текущее время-t.
ВНУТРЕННИЕ ПЕРЕМЕННЫЕ:
py[i]-массив значений правых частей;
y0[i],y1[i]-масивы для хранения промежуточных
значений переменных; t0-время в начале шага;
i,N-целочисленные переменные.
При использовании данной программы в качестве
вызываемой функции необходимо описать одномерные
массивы y,y0,y1,py размером N (все описания даются
с учётом нулевых значений индексов). Объявить типы
переменных t,h,y,i,N.

Исходники+исполняемый файл