function yhat = predar1(y,x,x0)
n=length(y);
Q = [ones(n,1) x(:,1)];
b = (Q'*Q)\(Q'*y);
yhat = x0'*b;
end