初めまして。
Fortranというよりは数値計算の分野の質問かも知れませんが、
もし経験者がいらっしゃいましたらご教授いただきたく存じます。
ほぼスレ違いですが失礼します。

現在Fortran90で6段5次ルンゲ・クッタ-フェールベルグ公式による常微分方程式計算をしているのですが、
誤差が小さいとされるこの公式がを自前でプログラム作成しましたところ、
それまでに使用していた4段4次ルンゲ・クッタ法より遥かに大きな誤差が出てしまいました。
係数の誤りは何度も確認しましたがそのようには思えません。
他に考えられる原因はありますでしょうか。