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

Secant Method
正割法

本式可以正割法求一方程的根

<COMP>
Goto 3 : Lbl 0 : f(A) → C : Goto 2 :
Lbl 3 : ? → A : ? → B : Lbl 1 : Goto 0:
Lbl 2 : ( B C – A D ) ÷ ( C – D → X :
C → D : A → B : X → A ◢ Goto 1
>61 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 : Goto 2 :
Lbl 3 : ? → A : ? → B : Lbl 1 : Goto 0:
Lbl 2 : ( B C – A D ) ÷ ( C – D → X :
C → D : A → B : X → A ◢ Goto 1
70 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.5 EXE
0.746549032 EXE
0.822920645 EXE
0.800003098 EXE
0.801051419 EXE
0.801070782 EXE
0.801070765 EXE
0.801070765 EXE ...
(x = 0.8011 corr. to 4 d.p.)
(正確數值 = 0.801070765425)

記憶體

A xn
B xn-1
C fn
D fn
X xn
Y  
M

相關程式

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