拉格朗日多项式插值遇到的问题

快乐枫

根据彭国伦fortran书中提到的拉格朗多项式日插值（P465），写了如下function。但是有时候可以插值成功，有时候却得到结果NaN，不知什么原因？
是函数写的有问题呢？还是拉格朗日多项式方法自身的不稳定造成的呢？求指点

已知点是X1(N)，Y1(N)，N是已知点个数，求x点处的插值

function lagrange(x,N,X1,Y1)
implicit none

real x
integer N
real X1(N),Y1(N)
real lagrange

real coeff
integer i,j
lagrange=0
do i=1,N
   coeff=1
   do j=1,N
    if(i/=j) coeff=coeff*(x-X1(j))/(X1(i)-X1(j))
   end do
   lagrange=lagrange+coeff*Y1(i)
end do

return
end function

MENGMENGDA

感谢！！！！！！

快乐枫

MENGMENGDA 发表于 2014-7-20 18:57
感谢！！！！！！

知道问题所在吗？？