import numpy as np
def interp(power, time, masses):
"""
This module creates the polynomial interpolation function for the mass/time data
"""
logtime = np.log(time)
logmasses = np.log(masses)
coefs = np.polyfit(logtime,logmasses, power)
polynomial = np.poly1d(coefs) # creates a polynomial function
return polynomial