数値計算プログラム

補間法は、他の解法の部分処理として使用される事が多いと思っています。
逆にいえば、基本的な手法です。

dim x(n),y(n)
s:=0
f:=1
i:=1
repeat
d:=1
j:=1
repeat
if j=i then d:=d*(x-x(j)):go label1
d:=d*(x(i)-x(j)
label1
until j>=n
s:=s+y(i)/d
f:=f*(x-x(i))
until i>=n
end