Hi,<br><br>You can can create p+1 polynomials with the poly() function, and then scale them by evaluating each of them at the remaining point.<br><br>Matthieu<br><br><div><span class="gmail_quote">2007/8/21, Neilen Marais &lt;
<a href="mailto:nmarais@sun.ac.za">nmarais@sun.ac.za</a>&gt;:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hi,<br><br>I&#39;m looking for functions that evalue to the Lagrangian interpolation
<br>polynomials, i.e.<br><br>L_i = \prod_{j=0..p, j!=i}(x-Xj)/(Xi-Xj)<br><br>where Xi are the p+1 interpolation points. I.e., I&#39;d like to pass in p+1<br>interpolation points and get back p+1 polynomials L_i that are each zero at
<br>all interpolation points except Xi where the value is 1. Is there an easy<br>way to construct these polynomials?<br><br>Thanks<br>Neilen<br><br>_______________________________________________<br>SciPy-user mailing list
<br><a href="mailto:SciPy-user@scipy.org">SciPy-user@scipy.org</a><br><a href="http://projects.scipy.org/mailman/listinfo/scipy-user">http://projects.scipy.org/mailman/listinfo/scipy-user</a><br></blockquote></div><br>