本研究では、自己調整型の多段階ルンゲ・クッタ法の効率的な実装アプローチを提案し、これらの手法の安定性解析を任意の副ステップ数の場合に拡張した。また、物理的に意味のある模型問題を提案し、異なる補間手法や副ステップ数の影響を調べた。数値実験の結果は、提案手法の有効性を示している。