function t = normalMLE(t0,Y)

s0 =[t0(1),log(t0(2))]'; 

options = optimset('LargeScale','off','DerivativeCheck','off',...
    'GradObj','off','TolX',1e-6,'Display','off',...
    'Diagnostics','off','MaxIter',1000000);

s=fminunc('normalLLF',s0,options,Y);

t = [s(1),exp(s(2))]';