function [ llv ] = lln(mu)
% computes the normal loglikelihood value
% 

global data sigmahat;

llv=0;
n=size(data,1);

for i=1:n;
   
    llv=llv-.5*log(2*pi)-log(sigmahat)-((data(i)-mu)^2)/(2*sigmahat^2);
   
end

  llv=-llv;

end