Casio50 » fx-3x50P » Programs 程式集 » Numerical Root-Finding 數值法求根 »

Newton’s Method
牛頓法

本式可以牛頓法求一方程的根

<COMP>
Goto 3 : Lbl 0 : f(X) → Y : Goto 2 :
Lbl 3 : ? → X : ? → D : D + X 10 - 6 ( D = 0 → D : 0 → A :
Lbl 1 : Goto 0 : Lbl 2 : A = 0 → C : X + D → X : A B + C Y → B : C → A : A => Goto 1 :
X – D – D Y ÷ ( Y – B → X ◢ Goto 1
>94 Bytes

操作方法

一般操作 例子
求
f(x) = 0
並以a為初始值 		求
cos^2 (x^2) = x^2 (in radian mode)
並以 0.5 為初始值
整體程式為:
<COMP>
Goto 3 : Lbl 0 : (cos X ² ) ² – X ² → Y : Goto 2 :
Lbl 3 : ? → X : ? → D : D + X 10 - 6 ( D = 0 → D : 0 → A :
Lbl 1 : Goto 0 : Lbl 2 : A = 0 → C : X + D → X : A B + C Y → B : C → A : A => Goto 1 :
X – D – D Y ÷ ( Y – B → X ◢ Goto 1
102 Bytes
(cos^2 (x^2) = x^2可寫成cos^2 (x^2) - x^2 = 0)
啟動程式 按 Prog 1
輸入初始值
a EXE		 0.5 EXE
輸入微分之差值. 不賦值則以 a/1000000 為差值
[d] EXE		 EXE
顯示近似值
x1 EXE
x2 EXE
x3 EXE ...		 0.965580083 EXE
0.812866149 EXE
0.801185672 EXE
0.801070776 EXE
0.801070765 EXE
0.801070765 EXE
0.801070765 EXE
0.801070765 EXE ...
(x = 0.8011 corr. to 4 d.p.)
(正確答案 = 0.801070765425)

記憶體

A 0
B fn
C 0
D d
X xn
Y ?
M  

相關程式

  1. Secant Method / 正割法
  2. Newton’s Method (II) / 牛頓法 (二)
  3. Bisection Method / 二分法
  4. Regular Falsi / 試位法