[SciPy-user] optimize, line fit of convex functions

I'd like to extrapolate some points from a fitted line. Given that true line is guaranteed to be convex, are there preferred methods for fitting the data?