function loglike=ll(phi,y,x1,x2)
n=length(y);
loglike=n*log(2*pi)/2 + n*log(phi(3))/2 +((y-phi(1)*(x1+phi(2)*x2))'*(y-phi(1)*(x1+phi(2)*x2))/(2*phi(3)));