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

Bisection Method (III)
二分法 (三)

本式可以二分法求一方程的根.

<COMP>
Goto 3 : Lbl 0 : f(X) → C : X = A => Goto 4 : Goto 2 :
Lbl 3 : ? → A : ? → B :
Lbl 1 : . 5 ( A + B → D ◢ D → X : Goto 0 :
Lbl 2 : C → Y ◢ A → X ◢ B ◢ Goto 0 :
Lbl 4 : 0 > C Y => D → B : C Y > 0 => D → A : Goto 1
>88 Bytes

操作方法

一般操作 例子
求
f(x) = 0
並以 [a,b] 為初始區間 		求
cos^2 (x^2) = x^2 (in radian mode)
並以 [0.5, 1] 為初始區間
<COMP>
Goto 3 : Lbl 0 : ( cos X ² ) ² − X ² → C : X = A => Goto 4 : Goto 2 :
Lbl 3 : ? → A : ? → B :
Lbl 1 : . 5 ( A + B → D ◢ D → X : Goto 0 :
Lbl 2 : C → Y ◢ A → X ◢ B ◢ Goto 0 :
Lbl 4 : 0 > C Y => D → B : C Y > 0 => D → A : Goto 1
97 Bytes
(cos^2 (x^2) = x^2可寫成cos^2 (x^2) - x^2 = 0)
啟動程式 按 Prog 1
輸入初始值
a EXE b EXE		 0.5 EXE 1 EXE
顯示近似值、其對應函數結果及該近似值之所在區間
x1 EXE f(x1) EXE a1 EXE b1 EXE
x2 EXE f(x2) EXE a2 EXE b2 EXE
x3 EXE f(x3) EXE a3 EXE b3 EXE ...		 0.75 EXE 0.153088258 EXE 0.5 EXE 1 EXE
0.875 EXE -0.24585699 EXE 0.75 EXE 1 EXE
0.8125 EXE -0.036219891 EXE 0.75 EXE 0.875 EXE
0.78125 EXE 0.061141116 EXE 0.75 EXE 0.8125 EXE
0.796875 EXE 0.013121343 EXE 0.78125 EXE 0.8125 EXE
0.8046875 EXE -0.011386487 EXE 0.796875 EXE 0.8125 EXE
...
(正確答案 = 0.801070765425)

記憶體

A an-1
B bn-1
C fn-1
D xn
X xn-1
Y f(an-1)
M  

其他資料

區間 [a,b] 是指變數 x 是在 a 和 b 之間.

相關程式

  1. Bisection Method / 二分法
  2. Bisection Method (II) / 二分法 (二)