Using the "Newton's Method"

Write a MATLAB script to solve for the following nonlinear system of equations:

x² + y² + z² = 3

x² + y² - z = 1

x + y + z =3

using the initial guess (x,y,z) = (1,0,1), tolerance tol = 1e-7, and maximum number of iterations maxiter = 20.