(a) Newton’s method for approximating a solution of an equation f(x) = 0 can be adapted to approximating a solution of a system of equations f(x, y) = 0 and g(x, y) = 0. The surfaces z = f(x, y) and z = g(x, y) intersect in a curve that intersects the xy-plane at the point (r, s), which is the solution of the system. If an initial approximation (x₁, y₁) is close to this point, then the tangent planes to the surfaces at (x₁, y₁) intersect in a straight line that intersects the xy-plane in a point (x₂, y₂), which should be closer to (r, s). Show that

where f , g, and their partial derivatives are evaluated at (x₁, y₁). If we continue this procedure, we obtain successive approximations (x_n, y_n).

(b) It was Thomas Simpson (1710 –1761) who formulated Newton’s method as we know it today and who extended it to functions of two variables as in part (a). (See the biography of Simpson in Section 7.7.) The example that he gave to illustrate the method was to solve the system of equations

x^x + y^y = 1000    x^y + y^x − 100

In other words, he found the points of intersection of the curves in the figure. Use the method of part (a) to find the coordinates of the points of intersection correct to six decimal places.

fg, - fg frgy - fy gx fag fg. fr 9, - f 9 X2 = X1 and y2 = y