Вам придется прочитать несколько книг по численным методам для ОДУ. Я думаю, что вы не можете сделать общий код даже для этого относительно простого случая. У вас здесь y' + F(y,t) = 0, который является ODE первого порядка. Не зная структуры F, мы не можем говорить об одном способе найти решение. Но если F хотя бы непрерывна по Липшицу, вы можете попробовать любой из методов Эйлера. Они не так сложны.
Я могу представить, что ваш код уже использует один из методов Эйлера, поэтому, если вам не важно, что F - какая-то «плохая» функция, вы можете просто изменить свой код, чтобы использовать вместо него F и y (0) = y_0 у / т и (0,1).