一元五次方程

程式可以計算一元五次方程式(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  3Abs Ans cos (arg iAns÷3) => Goto 2:

3iAns + 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)

 

註: 程式使用牛頓法計算第一個根的值,可能需要較長的時間,另外若果第一個根為重根,誤差可能會較大。

 

返回 fx-3650P 程式集

Free Web Hosting