OdeAuto(Real t0, Real t1, Matrix x, Matrix dx(), Matrix input(), Real eps, Real dtsav)epsは許容誤差を示します。省略するとeps = 1.0E-6となりま す。Autoがつくと,ルンゲクッタ法のきざみ幅を自動調節してくれ ます。Ode()の説明で積分のきざみ幅の重要性について述べました。き ざみ幅hの値を小さくすれば,解の精度は上がります。逆に,hを大き くすれば概して精度は下がりますが,関数が急激に変化しないところでは hを大きめに取っても解の精度が落ちることは有りません。そこで, 合理的に積分のきざみ幅を選ぶためには,1ステップを積分した時 に生ずる誤差の評価を行なうことです。
きざみ幅の自動調節では,xをh進める前に,次の1, 2を行なってhの 値を決定します。
積分結果が誤差以内でなかった場合はそのきざみを半分ずつ小さくするた
め時間がかかりますが,逆に関数の変化がゆるやかな場合は,きざみを大
きくするので,全体の計算時間は短くなります。また,積分の誤差を完全
に満たしているため信頼度は上がります。