python import numpy as np import statsmodels.api as sm x = np.array([1, 2, 3, 4, 5]) y = np.array([2, 4, 6, 8, 10]) X = sm.add_constant(x) model = sm.OLS(y, X) result = model.fit() print(result.summary())