clc;
% Find four solutions of cirlce hyperbola intersection
% Trial 1
x0 = [3 -1.5]';
tol = 1e-8;
[root,iter] = newtonmv(x0, @newtonmv_f, @newtonmv_j, tol)
% Trial 2
x0 = [0 3]';
[root,iter] = newtonmv(x0, @newtonmv_f, @newtonmv_j, tol)
% Trial 3
x0 = [-3 1.5]';
[root,iter] = newtonmv(x0, @newtonmv_f, @newtonmv_j, tol)
% Trial 4
x0 = [0 -3]';
[root,iter] = newtonmv(x0, @newtonmv_f, @newtonmv_j, tol)