Я думаю, что может быть задан похожий вопрос (хотя я не нашел ни одного). Существует ли какой-либо код на Fortran 90/95 для решения связанных IVP ODE с использованием метода конечных разностей или метода Эйлера (я верю, что он доступен, хотя я не нашел)? У меня есть около 150 связанных ODE с IVP следующей формы,
dx1(t)/dt=a01+a11.x2(t)+a12.x3(t)-b11.x1(t)
dx2(t)/dt=a02+a22.x3(t)-b22.x2(t)
dx3(t)/dt=a03+a31.x1(t)+a32.x2(t)-b13.x3(t)
где a и b - это константы, t - время,
x1(t=0)=0,x2(t=0)=0,x3(t=0)=0