Interpolation

import numpy as np
from scipy.interpolate import interp1d

x1 = np.array([1, 2, 3, 4, 5, 45])
y1 = np.array([1, 8, 27, 64, 125, 91125])

# Create an interpolation function
f = interp1d(x1, y1, kind='linear', fill_value='extrapolate')

# Define x_new
x_new = np.linspace(1, 45, 10)

# Interpolate y values
y_interp = f(x_new)

# Print the result
print(y_interp)