import scipy.interpolate as intp # SciPy has more interpolation functions than NumPy. N = 2**n # Fix nr. of knots. x = # The locations of the knots. y = # The y-values at the knots. S = ...