关于fortran的rk4

你的问题是k是未定义的,所以认为是整数,diff_equation中的dydx是未定义的,所以认为是浮点数。类型不一致,导致号码奇怪。

我强烈建议您添加:

隐式无

然后对于每个程序的每个数组,变量声明它的类型。

不要用维度,永远忘记这个词。

整数数组的定义如下:

整数::i(n)

浮点数组的定义如下:

实数::r(n)