A function to fit fair (bending energy minimizing) cubic bezier curves through a set of given ordered points in the plane.