x = np.array([0.1, 0.2, 0.4, 0.6, 0.8, 1.0, 2.0, 4.0, 6.0, 8.0, 10.0, 20.0,
    40.0, 60.0, 80.0])
y = np.array([0.5050533250540701, 1.1207373784533172, 2.198184471902, 
    3.174620900339869, 4.290548247126004, 6.281622667807696, 
    11.073788414382639, 23.24847977054601, 32.12046230136718, 
    44.036117671229206, 54.009003143831116, 102.7077685684846, 
    185.72880217806673, 256.1218314554581, 301.97120103079675])
x = x[:, (np.newaxis)]
a, _, _, _ = np.linalg.lstsq(x, y)
plt.plot(x, y, 'bo')
plt.plot(x, a * x, 'r-')
plt.show()
Как заставить нулевой перехват в линейной регрессии?