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);