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

Regular Falsi
試位法

本式可以試位法求一方程的根

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

操作方法

一般操作 例子
求
f(x) = 0
並以 [a,b] 為初始區間		 求
cos^2 (x^2) = x^2 (in radian mode)
並以 [0.5, 1] 為初始區間
整體程式為
<COMP>
Goto 3 : Lbl 0 : ( cos A ² ) ² − A ² → C : X = B => Goto 4 : X = A => Goto 5 : Goto 2 :
Lbl 3 : ? → A : ? → B : Lbl 1 : B → X : Goto 0 :
Lbl 4 : C → Y : A → X : Goto 0 :
Lbl 5 : ( A Y – B C ) ÷ ( Y – C → D ◢ D → X : Goto 0 :
Lbl 2 : 0 > C Y => D → A : C Y > 0 => D → B : Goto 1
118 Bytes
(cos^2 (x^2) = x^2可寫成cos^2 (x^2) - x^2 = 0)
啟動程式 按 Prog 1
輸入初始值
a EXE b EXE		 1 EXE 0.5 EXE
顯示近似值
x1 EXE
x2 EXE
x3 EXE ...		 0.746549032 EXE
0.793952743 EXE
0.800217528 EXE
0.800969717 EXE
0.801058815 EXE
0.801069352 EXE
0.801070598 EXE
0.801070745 EXE
0.801070762 EXE
0.801070764 EXE
0.801070765 EXE
0.801070765 EXE
0.801070765 EXE ...
(正確答案 = 0.801070765425)

記憶體

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

其他資料

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

相關程式

  1. Secant Method / 正割法
  2. Newton’s Method / 牛頓法
  3. Bisection Method / 二分法