誤差函數
更新日期: 2011年12月29日
程式速度很快,但準確度較低(誤差少於4×10-4相對準確度)。
注意: e^( 是按 SHIFT ex
程式 (38 bytes,使用記憶A)
8887 ÷ 63473: ?→A:
√( 1 - e^( -A2 (4÷π + Ans A2) ÷ (Ans A2 + 1
例題: 計算誤差函數 erf(0.2)的值。
按 Prog 1 再按 0.2 EXE (顯示答案0.2227)
附錄程式 (由網友 Skwai 提供及經網友roviury 改良)
附錄程式準確度很高,但計算時間較長。
程式一(54 bytes,比程式二快很多)
?→X: MM-: X→Y: While Ans-Y≠AnsM+:
Ans→D:-YX2( 2M - 1) ÷ (2M2 + M→Y:
D+Y: WhileEnd: 2Ans÷√(π
程式二需要在 SD 模式下執行,因此在選擇新程式位置後,按 4 選用SD Lin模式。
注意: 藍色的英文字為統計模式中的變數(Σx 按 Shift 1 2,n 按 Shift 1 3)。
程式二(42 bytes)
?→X: ClrStat: X DT: While Σx-Ans≠Σx:
- AnsX2( 2n - 1) ÷ ( 2n2+n DT: WhileEnd: 2Σx÷√(π
ClrStat: 0 DT: ?→A: 1: While
Σx - Ans≠Σx:
2Σy+1: 2(-1)^(Σy)A^(Ans)÷(AnsΣy!
√(π , 1 DT:
WhileEnd: Σx→A