function f = fung(x) %A=[-1 -1 -2; %GLIEX % -12 -3 -15; %O % 6 2 8;%CO2 % 0 1 1];%1HB A=[1 1 1 1 1 1 1 1 1 1 1 12 12 12 12 12 12 12 12 7.5 3 3 6 6 6 6 6 6 6 6 4 2 2 0 0 0 0 0 0 0 0 0.5 1 1] %F=[-1.935913; -3.68; 4.4; 1.563028]; F=[1.97052 6.87150 4.26000 1.60392]; ro=(A*x); no=F-ro; f =norm(no)