fortran - Fortran 90/95中ODE IVP系统的有限差分/欧拉方法？

Question

我认为可能会问类似的问题（虽然我没有找到）。是否有任何 Fortran 90/95 代码可以使用有限差分或欧拉方法解决耦合 ODE IVP（我相信它是可用的，但我没有找到）？我有大约 150 个带有 IVP 的耦合 ODE，形式如下，

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

标签： fortranode