% Session p. 113/114

format long e;

u=[0.014142; 1.028102; 0.636408; -0.596767; [65;113]*pi/180];
while true,
   f=p5func(u); [u(1)-u(2)+u(3)-u(4), max(abs(f))]
   if abs(f)<5e-15, break, end;
   u=u-jac('p5func',u)\p5func(u);
end
u