g,t=var('g t')
x=vector([0,1,2,3,4,5,6])
y=vector([3,5,7,8,10,10,22])
n=len(x)
summe_x=sum(x)
summe_y=sum(y)
summe_xy=x.dot_product(y)
summe_xquad=x.dot_product(x)
summe_yquad=y.dot_product(y)
b=(n*summe_xy-summe_x*summe_y)/(n*summe_xquad-(summe_x)^2)
a=(summe_y-summe_x*b)/n
g=a+b*t
p = plot(g, 0, 10, axes_labels=['x','y'],color='red')
d = list_plot(zip(x,y))
r=(n*summe_xy-summe_x*summe_y)/(sqrt(n*summe_xquad-summe_x^2)*sqrt(n*summe_yquad-summe_y^2))