13/09/14 21:44:06.31 4CwOwE5K
前掲のプログラムをちょっといじって、
REAL CATCH(0:*), POP(0:*)
RK = 100000 ← 初期資源量を強制的に10万に
R = 1.4184 * 0.04 ← 繁殖力を強制的に0.04に
POP(ISTART) = RK
DO 10 I = ISTART + 1, IYEAR
D = POP(I - 1) / RK
POP(I) = POP (I - 1) * (1. + R * (1. - D * D)) - 5000 ← 毎年の捕獲量を強制的に5千に
WRITE (*,*), I, POP (I) ← 画面出力
IF (POP(I) .LE. 0) GOTO 20
10 CONTINUE
STOP ← プログラムを強制終了
STKSIM = POP(IYEAR)
RETURN
出力結果
1 95000.
2 90525.5156
3 86452.6406
4 82691.6094
5 79175.1328
6 75851.2656
というか見れば分かるけど、Rはサブルーチンに入ってきたときから一切変更されてないね。