Конечная разность / метод Эйлера для системы ODE IVP в Fortran 90/95?
Думаю, может быть задан аналогичный вопрос (хотя я его не нашел). Есть ли какой-либо код Fortran 90/95 для решения связанных ODE IVP с использованием метода конечных разностей или метода Эйлера (я действительно считаю, что он доступен, но я не нашел)? У меня около 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
Вам нужен Эйлер или какой-нибудь решатель вроде rk4, rk45 / dopri, rkf45, ...?
Как упоминает @HighPerformanceMark, политика stackoverflow заключается в том, чтобы не отвечать на такие вопросы, как «советовать, какое программное обеспечение выбрать». Существует множество библиотек Fortran для научных вычислений, в том числе для решения ODE, вы, вероятно, не спрашивали Google "обыкновенное дифференциальное уравнение Fortran".
Мне нужен метод Эйлера или метод конечных разностей. Спасибо. Я сделаю это с помощью библиотеки ..
Метод Эйлера - это метод решения ОДУ. Там настолько банально, что программировать почти нечего. Вместо этого вы найдете библиотеки в основном для методов Рунге-Кутта. Термин «конечная разность» является широким термином для многих различных методов, здесь нет метод конечных разностей. Кроме того, этот термин обычно используется в контексте PDE, а не ODE. Предлагаю узнать больше о математике этого. Тогда вы увидите, что явный метод Эйлера совершенно тривиален для программирования.
@VladimirF Спасибо. Я изучаю математику ...
Обратите внимание, что вопросы, касающиеся сторонних ресурсов, здесь явно не по теме, если вы не уверены в этом, обратитесь в справочный центр, где все это объясняется. И спрашивать, где вы должны задать такой вопрос, тоже было бы не по теме.