GRIP中如何实现拉格朗日插值函数？？

wll

, F% n* z8 n# b- g0 X1 Y. j
1 h* }6 w* T  a' d# Q2 |5 m% F( Wf(r)=x1*L1(r)+x2*L2(r)+...+x8*L8(r)
其中Li(r)=(r-r1)...(r-r(i-1))(r-r(i+1))...(r-rn)/(ri-r1)...(ri-r(i-1))(ri-r(i+1))...(ri-rn)
如果一项一项的按顺序写下去，可以实现这个式子，但是很繁琐。
大家有什么办法，用什么函数或者语句可以实现拉格朗日插值多项式？
# k; r% C4 ?9 [$ X* q/ D$ `各位谁有思路给一下提示，先谢过啦^_^

另：
# t  L! F) ?' Z3 O+ I' E) r( B9 }! d就比如说我想实现这样一个循环
& w4 X6 o9 D, s& ~: DDO/trm:,i,1,5
x(i)=n(1)*(e(i)-m(1))*(e(i)-m(2))*(e(i)-m(3))/((m(i)-m(1))*(m(i)-m(2))*(m(i)-m(3)))
trm:
( W8 F% n0 g6 t! p, |! W- _$ l
不论i等于几时，分母都有0项，改怎么办呢？