for mm in range(4):
x = t_teo[mm].reshape((-1,1)) y = t_exp[mm] regression = LinearRegression() model = regression.fit(x, y) r_sq = model.score(t_teo[mm], t_exp[m