function [m,b]= linreg(X,Y)
sumxy=0;
sumy=0;
sumx=0;
sumxx=0;
n = length(X);
for ii= 1:n
    sumx = sumx + X(ii);
    sumy =sumy + Y(ii);
    sumxy = sumxy + X(ii)*Y(ii);
    sumxx = sumxx+X(ii)*X(ii);
end
m = (sumxy - sumx*sumy/n)/(sumxx - sumx*sumx/n);
b=(sumy - m*sumx)/n;