f= @(x)[((93-x(1))^2)+((63-x(2))^2)-(55.1)^2;((6-x(1))^2)+((16-x(2))^2)-(46.2)^2];
J= @(x)[-2*(93-x(1)),-2*(63-x(2));-2*(6-x(1)),-2*(16-x(2))];
h = J(x)\f(x);
TOL = 1e-10;
x = [0;0];
while norm(h)>TOL

x = x-h
end
h