Subsections

• 1 問題

2 境界値問題としてのLane-Emden 方程式

Lane-Emden方程式(7.35)を境界値問題として解いてみよう。 適当なフィッティングポイント$\xi_{\rm f}$を設け、 中心からフィッティングポイントまでの解は前節で解説した通りの方法で求め、 その解に、表面、即ち$\theta(\xi)=0$となる $\xi=\xi_{\rm s}$から 積分してきた解が、このフィッティングポイントで一致するようにすれば良い。

表面からの積分は$\xi_{\rm s}$の値を決めた上で行わなければならないのは言うまでもない が、微分方程式は２階だから、 $\xi=\xi_{\rm s}$での初期条件が２つ必要である。 その１つは $\theta(\xi)=0$だが、もう１つの$d\theta/d\xi$は最初はわからない。 そこで、$\xi_{\rm s}$とそこでの$d\theta/d\xi$の値をパラメーターとして与えてやれば、 外側の解をこの２つのパラメーターの関数として扱う事が出来る。 $y_1 = \theta$， $y_2 = \xi^2 d\theta/d\xi$とし、 $y_1^{{\rm (int)}}$と $y_2^{{\rm (int)}}$を中心からの積分で求めた $\xi=\xi_{\rm f}$での値、 $y_1^{{\rm (ext)}}$と $y_2^{{\rm (ext)}}$を表面からの積分で求めた $\xi=\xi_{\rm f}$での値としよう。内側の解はユニークに決まっている。 一方、外側の解は $x \equiv \xi_{\rm s}$と そこでの $z \equiv d\theta/d\xi$の値に依存している。 全体の解となるには、 $\xi=\xi_{\rm f}$で

$$Y(x,z) \equiv y_1^{{\rm (int)}}-y_1^{{\rm (ext)}}(x,z) =0$$ (96)

$$Z(x,z) \equiv y_2^{{\rm (int)}}-y_2^{{\rm (ext)}}(x,z) =0$$ (97)

が満たされなければならない。

手順としては、 先ず、適当に選んだ$x_{\rm t}$と$z_{\rm t}$に対する$Y$と$Z$を 求める。 次に、$z$の値を変えずに、$x$を$x_{\rm t}$から $x_{\rm t}+\delta x$に変えて、 $Y$及び$Z$の変化量を求める。 同様に、$x$の値を変えずに、$z$を$z_{\rm t}$から $z_{\rm t}+\delta z$に変えて、 $Y$及び$Z$の変化量を求める。 これから

$$\partial Y/\partial x = [Y(x_{\rm t}+\delta x,z_{\rm t})-Y(x_{\rm t},z_{\rm t})]/\delta x$$ (98)

$$\partial Y/\partial z = [Y(x_{\rm t},z_{\rm t}+\delta z)-Y(x_{\rm t},z_{\rm t})]/\delta z$$ (99)

を数値的に求める。 今、$x_{\rm t}$、$z_{\rm t}$の値が、真の解$x_0$，$z_0$から $\Delta x$，$\Delta z$ずれているとすると、１次の精度で、

$$Y(x_{\rm t},z_{\rm t}) + (\partial Y/\partial x)\Delta x + (\partial Y/\partial z)\Delta z =0$$ (100)

$$Z(x_{\rm t},z_{\rm t}) + (\partial Z/\partial x)\Delta x + (\partial Z/\partial z)\Delta z =0$$ (101)

となる。そこで、上の２式から$\Delta x$と$\Delta z$を求め、 $x=x_{\rm t}+\Delta x$を改めて$x_{\rm t}$， $z=z_{\rm t}+\Delta z$を改めて$z_{\rm t}$として、 $\Delta x$，$\Delta z$が充分小さくなるまで、逐次計算を繰り返せば良い。

1 問題

境界値問題として、Lane-Emden方程式を解け。