四點拉格朗日插值多項式 ( Lagrange Interpolating Polynomial for 4 points )

( 程式版本：1.0，最後更新日期 14 AUG 2003。)

LBL 0	R/S	STO 1	STO 0	R/S	STO 8
STO 7	R/S	STO 2	STO –0	STO 9	R/S
STO 6	STO –7	RCL 0	STO ÷7	R/S	STO 3
STO –9	STO 4	R/S	STO 5	STO –6	RCL 9
STO ÷6	STO +0	R/S	STO –4	R/S	STO –5
RCL 4	STO ÷5	STO +9	RCL 6	STO –5	RCL 7
STO –6	RCL 9	+/-	STO ÷5	RCL 0	+/-
STO ÷6	RCL 4	STO +0	RCL 6	STO –5	RCL 0
+/-	STO ÷5	RCL 7	×	RCL 1	STO 9
－	RCL 1	×	RCL 2	STO +9	×
STO 0	RCL 6	+	RCL 5	×	RCL 0
×	RCL 3	=	STO –8	RCL 6	×
RCL 9	－	RCL 5	×	(	RCL 0
+	RCL 9	×	RCL 3	=	STO –7
RCL 5	R/S	×	(	RCL 9	+
RCL 3	=	STO –6	RCL 6	R/S	RCL 7
R/S	RCL 8	PRGM	98 steps

例：計算經過 (-1, 0), (0, 0), (1, 2) 及 (2, 18) 四點的拉格朗日插值多項式。

按 XEQ 0 ，再按  1 +/- R/S  0 R/S  0 R/S  0 R/S  1 R/S  2 R/S  2 R/S  18 R/S

顯示    2   ( x³ 的係數 ( coefficient ) )
再按 R/S  顯示    1  ( x² 的係數 )
再按 R/S  顯示    –1  ( x 的係數 )
再按 R/S  顯示    0  ( 常數項 )

亦即是說所求的多項式是  2x³ + x² – x

程式執行完成後，按 RCL 5 會顯示多項式 x³ 的係數，按 RCL 6 會顯示多項式 x² 的係數，按 RCL 7 會顯示 x 的係數，而按 RCL 8 會顯示常數項。

如果輸入的四點中有任何兩點的 x 座標 ( x-coordinate ) 相同，則程式會顯示錯誤 Error-Func。