[SciPy-user] Interpolation polynomials

[Neilen Marais]

Hi,

I'm looking for functions that evalue to the Lagrangian interpolation
polynomials, i.e.

L_i = \prod_{j=0..p, j!=i}(x-Xj)/(Xi-Xj) 

where Xi are the p+1 interpolation points. I.e., I'd like to pass in p+1
interpolation points and get back p+1 polynomials L_i that are each zero at
all interpolation points except Xi where the value is 1. Is there an easy
way to construct these polynomials?

Thanks
Neilen