function [L] = makelike(thet,y);
m=thet(1);
s=thet(2);
T=max(size(y,1),size(y,2));
L=(T/2)*log(2*pi*s)+(ones(size(y))'*((y-m).^2))/(2*s);
%L=abs(L1);
