function [p1,p2]=lsm(x,y)
%least squares method for linear interpolating

X=[x ones(length(x),1)];
p=inv(X'*X)*X'*y;
p1=p(1);
p2=p(2);