function [A,B]=lsline(X,Y)
% Linear regression
n=length(X);
xmean=mean(X);
ymean=mean(Y);
sumx=sum(X);
sumy=sum(Y);
sumxy=sum(X.*Y);
sumx2=sum(X.*X);
A=(n*sumxy-sumx*sumy)/(n*sumx2-sumx^2);
B=ymean-A*xmean;