XLPack ソルバー 2階常微分方程式

2 階常微分方程式の初期値問題の解を求めます.

d2y1/dt2 = -y1/r, d2y2/dt2 = -y2/r (ただし, r = (y1^2 + y2^2)^(3/2))
初期値は, t0 = 0 において y1 = 1, y2 = 0, dy1/dt = 0, dy2/dt = sqrt(3)

を解く例を示します. 従属変数が2つで方程式が2本の2階連立常微分方程式の例になります.

独立変数(t)セル(この場合B8)および従属変数(yi)セル(この場合B9～B10)にはプログラムがそれぞれの値を書き込みます. 微分式(dyi/dt)セル(この場合B13～B14)にはtおよびyiよりd2y1/dt2およびd2y2/dt2を求める数式を入力しておきます. この場合, B13は =-B9*(B9^2+B10^2)^(-3/2), B14は =-B10*(B9^2+B10^2)^(-3/2) とします.

初期値および出力セル範囲(この場合D8:H28)の各行の1列目にtを入力します. 1行目には初期値(この場合 t0, t0におけるy1, y2, dy1/dt, dy2/dt)を入力しておきます.「計算開始」をクリックすると, 各行のtにおけるy1, y2, dy1/dt, dy2/dtを求め2列目以降に出力します. 指定した範囲内であってもtが空白の行に出会うと計算を終了します.

使用するプログラムは Dopn43, Dopn64, Dopn86 または Dopn1210 から選ぶことができます.

計算精度の標準値は 1.0e-8 です. 各ルーチンのRTolおよびATolにこれを設定します.

それぞれのセル範囲は必要な範囲よりも大きく指定することができ, その場合には左上から必要なだけの大きさの範囲を使います.

「保存／読込」についてはこちらを参照ください.

「ヘルプ」をクリックするとネット接続中であればこのページを表示します.

右上の「？」は正しく動作しません. 「ヘルプ」ボタンを使用してください.