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

wll

2 r9 g9 U# n& b; ^1 d
f(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)
9 U8 q1 ?9 t" ]( M6 f如果一项一项的按顺序写下去，可以实现这个式子，但是很繁琐。
大家有什么办法，用什么函数或者语句可以实现拉格朗日插值多项式？
各位谁有思路给一下提示，先谢过啦^_^

6 V: ]: W0 C) T  L/ X8 _另：
( T- }" ^  q" V  [就比如说我想实现这样一个循环
DO/trm:,i,1,5
. t# }0 m5 f7 K& R$ s1 ?6 ax(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)))
7 R5 O7 f! I  Rtrm:

  h+ q6 @. A4 ]+ K# w: l不论i等于几时，分母都有0项，改怎么办呢？
! {/ ~+ T& y& N2 `: y! q! U