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----