function MultipleNLEexample
clear, clc, format short g, format compact
xguess = [1. 300.]; % initial guess vector 
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
%- - - - - - - - - - - - - - - - - - - - - -
function fx = MNLEfun(x); 
CA1 = x(1); 
V = x(2); 
k = 0.075; 
v = 30; 
CA0 = 1.6; 
CA2 = 0.2 * CA0; 
fx(1,1) = k * CA1 ^ 2 - v * (CA0 - CA1) / V; 
fx(2,1) = k * CA2 ^ 2 - v * (CA1 - CA2) / V;