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

Newton’s Method (II)
牛頓法 (二)

本式可以牛頓法求複函數之根

<CMPLX>
Goto 3: Lbl 0: f(M) → A: Goto 2:
Lbl 3: ? → M: ? → B: Lbl 4: 0 → C:
Lbl 1: Goto 0: Lbl 2:
C = 0 => B M– => A → C => Goto 1:
A = C => B M+ => Goto 5:
B + BC ÷ ( A – C M+:
Lbl 5: M ◢ Goto 4
>81 Bytes

操作方法

一般操作 例子
求根:
f(x)=0


x=k
作初始值
求根:
x^10 + 1/x = 0


x = i
作初始值
整體程式為
<CMPLX>
Goto 3: Lbl 0: M M ³ ³ + M -1 → A: Goto 2:
Lbl 3: ? → M: ? → B: Lbl 4: 0 → C:
Lbl 1: Goto 0: Lbl 2:
C = 0 => B M– => A → C => Goto 1:
A = C => B M+ => Goto 5:
B + BC ÷ ( A – C M+:
Lbl 5: M ◢ Goto 4
86 Bytes

註:
冪運算子 "^" 並不適用於複數
啟動程式 Prog 1
輸入資料
k EXE dx EXE
i EXE 0.1 EXE
顯示近似值
x1 EXE
x2 EXE
x3 EXE ...
0.062311667 + 0.861821773i EXE
0.354295012 + 1.01790267i EXE
0.212871259 + 0.940495832i EXE
0.129042893 + 1.00798382i EXE
0.133029109 + 0.9853219i EXE
0.144559063 + 0.985614051i EXE
0.143970985 + 0.990866787i EXE
0.141878954 + 0.990489856i EXE
0.142038603 + 0.98963589i EXE
...
(正確答案 = 0.14231481085+0.989821431509i)

記憶體

A ?
B dx
C f(xn-1)
M xn

其他資料

xn 非常接近正確值時, 可能會出現「周期性偏差」

相關程式

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