06/02/09 02:18:09
>>141 ほれ!ティンティンうpねw
微分方程式の積分を最高簡単にしておいた。誤差が溜まって結果怪しい。
ルンゲ=クッタとか真面目なのに自分で治しとけw
PROGRAM oppai
PARAMETER (h = 0.2, a = 1.0, b = -1.0, c = -1.0, d = 1.0)
t = 0.0
x = 3.0
y = 1.0
DO 10 i = 1, 100
t = i * h
dx = (a + b * y) * x * h
dy = (c + d * x) * y * h
x = MAX(0.0, x + dx)
y = MAX(0.0, y + dy)
WRITE(6, *) x, y, t, i
10 CONTINUE
END
>>142 インデント
結果はというと、t=10くらいで、生命体xが滅亡し、その後餌のなくなった
生命体yもすぐに滅亡したw