title: How to do a polynomial fit with fixed points
body: <p>I have been doing some fitting in python using numpy (which uses least squares).</p>
  
  <p>I was wondering if there was a way of getting it to fit data while forcing it through some fixed points?  If not is there another library in python (or another language i can link to - eg c)?</p>
  
  <p><strong>NOTE</strong> I know it's possible to force through one fixed point by moving it to the origin and forcing the constant term to zero, as is noted here, but was wondering more generally for 2 or more fixed points :</p>
  
  <p><a href="http://www.physicsforums.com/showthread.php?t=523360">http://www.physicsforums.com/showthread.php?t=523360</a></p>