Используя iterp1d от scipy, вы можете подогнать куб к некоторым данным.Но как получить коэффициенты кубического полинома, который вы только что установили?
В качестве тривиального примера:
from scipy.interpolate import interp1d
x = [1, 2, 3.5, 4]
y = [2, 3, 4, 5]
f = interp1d(x,y, kind='cubic')