原子动力工学特论.pptVIP

原子動力工学特論 レポート1 交通電子機械工学専攻 齋藤　泰治 目次 ニュートン法について * ニュートン法 方程式f(x)=0の解を求める 処理手順 1. 初期値 x の決定 x = x0　を決める 2. 接線と x 軸の交点の計算 x = x0 における y = f(x) の接線を引き、今度はこの接線と y=0 (x軸)の交点を x1 とする すなわち、 xn+1 = xn - (f(x)/f ’(x)) を計算する 3. 2. を繰り返して値が収束したらそれを解とする ニュートン法 ?初期値ｘ=ｘ0をとる グラフに接線１を引く 接線１がX軸をきる点を次の点、ｘ1とする ｘ＝ｘ1のところで次の接線２を引き どうようにX軸をきる点をｘ2とする 同じことをn回繰り返すと 真の値xに極めて近い値 を得る ｘ0 ｘ1 ｘ2 ｘ≈xn 問1　 を解けばよい より、 従って、解は にある x0=1.000000 i=1 x=1.000000 f(x)= 0.158529015192 df(x)=1.459699 i=2 x=0.891396 f(x)= 0.016637296200 df(x)=1.154467 i=3 x=0.876985 f(x)= 0.000288172633 df(x)=1.114499 i=4 x=0.876726 f(x)= 0.000000092913 df(x)=1.113783 i=5 x=0.876726 f(x)= 0.000000000000 df(x)=1.113783 ans=0.876726 プログラム実行結果1-1 x0=-1.000000 i=1 x=-1.000000 f(x)= 1.841470984808 df(x)=-2.540302 i=2 x=-0.275098 f(x)= 0.347319458830 df(x)=-1.512593 i=3 x=-0.045479 f(x)= 0.047531597880 df(x)=-1.089923 i=4 x=-0.001869 f(x)= 0.001872376843 df(x)=-1.003735 i=5 x=-0.000003 f(x)= 0.000003475719 df(x)=-1.000006 i=6 x=0.000000 f(x)= 0.000000000009 df(x)=-0.999999 ans=0.000000 プログラム実行結果1-2 問２　 を解く また、ｘ＝0　で かつ 従って、解は x0=0.000000 i=1 x=0.000000 f(x)= 0.000000000000 df(x)=-0.999999 ans=0.000000 x0=2.000000 i=1 x=2.000000 f(x)= 0.551482119784 df(x)=0.773962 i=2 x=1.287456 f(x)= -0.003592498639 df(x)=0.684155 i=3 x=1.292707 f(x)= 0.000007627617 df(x)=0.687055 i=4 x=1.292696 f(x)= 0.000000000037 df(x)=0.687049 ans=1.292696 プログラム実行結果2-1,2 x0=5.000000 i=1 x=5.000000 f(x)= -0.276924238464 df(x)=-0.965662 i=2 x=4.713229 f(x)= 0.008136103732 df(x)=-1.008975 i=3 x=4.721292 f(x)= 0.000000405684 df(x)=-1.008864 i=4 x=4.721293 f(x)= 0.000000000000 df(x)=-1.008864 ans=4.721293 Ｘが十分に大きくなれば、exp(-x)は限りなく0に近づくので 解は、 （xが大きいとき） プログラム実行結果2-3 問3 を解く x0=1.000000 i=1 x=1.000000 f(x)= 0.864664716763 df(x)=2.270671 i=2 x=0.619203 f(x)= 0.093566571545 df(x)=1.818098 i=3 x=0.567739 f(x)= 0.001059165386 df(x)=1.778015 i=4 x=0.567143 f(x)= 0.000000126972 df(x)=1.777590 i=5 x=0.56