xtmp(1:n) = x(1:n) + qHlf * ck1(1:n) * ds ck2(1:n) = f(n, s + qHlf * ds, xtmp(1:n)) というのはやりたくなかったので, ck2(1:n) = f(n, s + qHlf * ds, x(1:n) + qHlf * ck1(1:n) * ds) と書いたら, forrtl: warning (402): fort: (1): In call to F, an array temporary was created for argument #3 forrtl: 警告 (402): fort: (1): F への呼び出しで、一時的配列が引数 #3 用に作成されました。 って警告がでるので,ここは仕方なし. ! --+*---1----+----2----+----3----+----4----+----5----