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

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

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