function % Insert here your file name after function (Use Alphanumberic names only)
clear, clc, format short g, format compact
xguess=     % Replace this line with the xguess line(s) from Polymath report.
disp('Variable values at the initial estimate');
fguess=MNLEfun(xguess);
disp(' Variable        Value     Function Value')
for i=1:size(xguess,2);
disp([' x' int2str(i) '                 ' num2str(xguess(i)) '            ' num2str(fguess(i))]);
end
options = optimset('Diagnostics',['off'],'TolFun',[1e-9],'TolX',[1e-9]); 
xsolv=fsolve(@MNLEfun,xguess,options); 
disp('Variable values at the solution');
fsolv=MNLEfun(xsolv);
disp(' Variable        Value     Function Value')
for i=1:size(xguess,2);
disp([' x' int2str(i) '                    ' num2str(xsolv(i)) '            ' num2str(fsolv(i))])
end
%- - - - - - - - - - - - - - - - - - - - - -
% Replace this and the following line with the function copied from the Polymath report
% Do not include the xguess line(s)