You need to use the multivariate Newton-Raphson method, which involves the Jacobian.

So, let's find estimates for . and , so that is a good estimate to start with.

Then,

Plugging in gives:

Now, to apply the Newton-Raphson method, we want:

So, you need to take the inverse of the matrix.

Then, you have:

Next,

Continuing: