function loglik = lNorm(b0, b1, X, Y, sd)

resid = Y - b0 - b1*X;
loglik = - sum(resid.^2) / (2*sd^2);