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

注意：此程式只能在 HP-32SII 上使用。

LBL A	R/S	STO A	STO J	R/S	STO G
STO H	R/S	STO B	STO –J	STO I	R/S
STO F	STO –G	RCL J	STO ÷G	R/S	STO C
STO –I	STO D	R/S	STO E	STO –F	RCL I
STO ÷F	STO +J	R/S	STO –D	R/S	STO –E
RCL D	STO ÷E	STO +I	RCL F	STO –E	RCL G
STO –F	RCL I	+/-	STO ÷E	RCL J	+/-
STO ÷F	RCL D	STO +J	RCL F	STO –E	RCL J
+/-	STO ÷E	RCL A	STO I	RCL ×G	RCL B
STO +I	RCL ×A	STO J	RCL ×F	－	RCL E
RCL ×J	RCL ×C	+	STO –H	RCL F	RCL ×I
RCL I	RCL ×C	RCL +J	RCL ×E	－	STO –G
RCL I	RCL +C	RCL E	R/S	×	STO –F
RCL F	R/S	RCL G	R/S	RCL H	PRGM
124.5 bytes

( 程式連同數字記憶共使用 204.5 bytes。)
( 如果日後想繼續輸入其他程式，則請在離開 PRGM 模式前按 RTN。)

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

按 XEQ A，再按  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 E 會顯示多項式 x³ 的係數，按 RCL F 會顯示多項式 x² 的係數，按 RCL G 會顯示 x 的係數，而按 RCL H 會顯示常數項。

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