reg = LinearRegression(fit_intercept=True).fit(xb, y)
# alpha und beta ausgeben
q = reg.intercept_
m = reg.coef_
y1 = reg.predict(xb)
y2 = m*x+q
print(y1)
print(y2)
sum((y1-y)**2)
sum((y2-y)**2)