x= [stuff]
y=[stuff]
n=length(x)
xy=x.*y
xx=x.*x
sx=sum(x)
sy=sum(y)
sxy=sum(xy)
sxx=sum(xx)
m=(n*sxy - (sx)*(sy)) / (nsxx - (sx)^2)
b=(sy - m*(sx)) / n
newy=m*x+b