一元五次方程
程式可以計算一元五次方程式(quintic equations)的五個根(包括複數根),程式需要使用兩個程式位置。
程式編寫日期: 2006年8月5日
程式第一部份 (COMP 模式,141 bytes)
?→X: ?→B: ?→C: ?→A: ?→D: ?→Y: π→M:
Lbl 0: (XM^5 + BM^4 + CM3 + AM2 + DM + Y)
÷ (5XM^4 + 4BM3 + 3CM2 + 2AM + D M- => Goto 0:
- (B┘X + M)┘4→B: C┘X - 4BM→C: A┘X + CM→A:
M→Y: 3B2 2 - CB2 - AB - D┘X - AM→M: Y
程式第二部份 (CMPLX 模式,217 bytes)
注意 : 若果不需要記存方程的答案,程式中綠色部份可以不輸入。
程式需要在 CMPLX 模式下執行,因此在輸入程式前請先按 Mode 2。
CB┘2 - 2B3 + A┘4→A: B2 - C┘6→C:
A => Goto 0: 3C→C: C + √( C2 + M→A: MM-: Goto 1:
Lbl 0: CM + C3 - A2: Ans - √(Ans2 + (M┘3 - C2)3:
iAns => 2 3√Abs Ans cos (arg iAns÷3) => Goto 2:
3√iAns + 3√(2A2 - 2CM - 2C3 - iAns: Lbl 2: √(Ans┘2 + C→M:
3C - M2→C: C + A÷M→A: Lbl 1: 2C - A→C:
Abs A => √Ans∠ . 5arg A: Ans - M + B→A◢
Abs C => √Ans∠ . 5arg C: Ans + M + B→C◢
2B - 2M - A→M◢ 4B - A - Ans - C→B
以下例題假設程式第一部份及第二部份分別儲存於P1及P2。
例題1: 解 x5 - 3x4 - 23x3 + 51x2 + 94x - 120 = 0
按 Prog 1 1 EXE - 3 EXE - 23 EXE 51 EXE 94 EXE - 120 EXE (顯示第一根為3)
再按 Prog 2 (顯示第二個根為 - 2) EXE (顯示第三個根為 5)
EXE (顯示第三個根為 - 4) EXE (顯示第五個根為 1)
例題2: 解 3x5 - 2x4 - 22x3 + 27x2 - 10x + 24 = 0
按 Prog 1 3 EXE -2 EXE -22 EXE 27 EXE -10 EXE 24 EXE (顯示第一根為2)
再按 Prog 2 (顯示第二個根為1.86977)
EXE ( 此時計算機右上角出現R<=>I,表示為複數解)
(顯示第三個根實數部為 - 0.101554)
Shift Re<=>Im (顯示第 三個根虛數部為 0.838323 i)
EXE (顯示第四個根為 - 3)
EXE (顯示第五個根的實數部為 - 0.101554)
Shift Re<=>Im (顯示第 五個根虛數部為 - 0.838323 i)
註: 程式使用牛頓法計算第一個根的值,可能需要較長的時間,另外若果第一個根為重根,誤差可能會較大。