%EXS 3.8
clear, clc, close all;

ab0=[0.6;1];
delta=inf;
k1=0.12;
k2=0.22;
while norm(delta./ab0)>0.5e-3
   F=[ab0(1)+k1*ab0(1)*ab0(2)+k2*ab0(1)*ab0(2)^2-0.6
       ab0(2)+k1*ab0(1)*ab0(2)+2*k2*ab0(1)*ab0(2)^2-1];
   J=[1+k1*ab0(2)+k2*ab0(2)^2 k1*ab0(1)+2*k2*ab0(1)*ab0(2);
       k1*ab0(2)+2*k2*ab0(2)^2 1+k1*ab0(1)+4*ab0(1)*ab0(2)];
   delta=-J\F;
   ab0=ab0+delta;
end
ab0
F=[ab0(1)+k1*ab0(1)*ab0(2)+k2*ab0(1)*ab0(2)^2-0.6
    ab0(2)+k1*ab0(1)*ab0(2)+2*k2*ab0(1)*ab0(2)^2-1]
