Newton's Method generalized to 3 dimensions

Click For Summary
Newton's Method can be generalized to find the roots of systems of equations in two variables, f(x,y) = 0 and g(x,y) = 0, by using the intersection of their corresponding surfaces. The tangent planes at an initial approximation (x1,y1) lead to updated coordinates (x2,y2) through specific formulas involving partial derivatives. The discussion highlights the equivalence of choosing either function for z in the equations, as both yield the same results when evaluated at the intersection point. The confusion arises from the aesthetic pairing of functions with their respective derivatives, which ultimately does not affect the outcome. The thread concludes that the method remains valid regardless of the choice made between f and g.
kostoglotov
Messages
231
Reaction score
6

Homework Statement



Prelim: my question is about a very specific part of a question whereby the student is asked to derive the final formula for the general solution in two vars, but I will post the entire question for clarify.

Newton's Method for approximating the roots of an equation f(x)=0 can be adapted to approximating solutions 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 (x1,y1) is close to the point (r,s), then the tangent planes to the surface at (x1,y1) intersect in a straight line that should intersect the xy-plane closer to (r,s).

Show that

x_2 = x_1 - \frac{fg_y-f_yg}{f_xg_y-f_yg_x}

and

y_2 = y_1 - \frac{f_xg-fg_x}{f_xg_y-f_y-g_x}

Homework Equations



\vec{v} = \vec{n_f} \times \vec{n_g}: direction vector of line of intersection of the two planes, basically the cross product of the normal vectors to the two planes

x_n = x_{n-1} - \frac{f(x_{n-1})}{f'(x_{n-1})}: single var form of Newton's Method to be adapted.

The Attempt at a Solution



So I consider this line of intersection of two planes from two perspectives. One looking at it "down" the y-axis to consider just the relationship between z and x, and then by looking "down" the x-axis to consider just the relationship between z and y.

I can see that \frac{\partial z}{\partial x} = \frac{f_xg_y-f_yg_x}{g_y-f_y} and \frac{\partial z}{\partial y} = \frac{f_xg_y-f_yg_x}{f_x-g_x} by taking the z component of \vec{v} (the direction vector of the line formed from the cross product) and dividing it by the x and y components of \vec{v} respectively.

I get to here with it:

x_2 = x_1 - z_1 \frac{g_y-f_y}{f_xg_y-f_yg_x}

Now, I can see that since all the points we are working with here are points of intersection, that the x,y and z's are all common for all the surfaces, planes and lines concerned. So z1 = f(x1,y1) = g(x1,y1)...so I can accept that we could pick either f or g to be z in the equation...but the text specifically states

x_2 = x_1 - \frac{fg_y-f_yg}{f_xg_y-f_yg_x}

The text shows the z = f when multiplied by the gy and z = g when multiplied by the fy.

Why is this?

A follow up question to this one is to find the points of intersection of x^x+y^y=1000 and x^y + y^x = 100 so I can see how it might be important to choose z = f for multiplication with gy and z = g for multiplication with fy, but I can't spy the reasoning behind it.
 
Physics news on Phys.org
I don't quite get the logic of your solution - but going back to the "relevant equations" part, if you calculate the vector ## \vec v ## then you are searching for the intersection of the line ## \{(x,y,z)+t\vec v,t\in\mathbb{R}\} ## with the plane z=0. Would that work ?
 
wabbit said:
I don't quite get the logic of your solution - but going back to the "relevant equations" part, if you calculate the vector ## \vec v ## then you are searching for the intersection of the line ## \{(x,y,z)+t\vec v,t\in\mathbb{R}\} ## with the plane z=0. Would that work ?

Actually, I realized that I was just confusing myself. Provided f(x,y) = 0 and g(x,y) = 0, z = f = g when evaluated at the same x,y, so it doesn't matter which partial deriv you pair the z,f, or g with, they'll always be the same value anyway, so I think the text was just pairing f with gy and g with fy for aesthetic purposes, because regardless of which you choose that second problem I mentioned x^x+y^y=1000 and x^y + y^x = 100 will have the same solutions, which look just fine when you graph them.
 
Well there is only one Newton formula you can arrive at, no aestethic choice involved - so I am not sure what you mean here.

Also my question in the previous post was rhetorical. It will work.
 
  • Like
Likes kostoglotov
Question: A clock's minute hand has length 4 and its hour hand has length 3. What is the distance between the tips at the moment when it is increasing most rapidly?(Putnam Exam Question) Answer: Making assumption that both the hands moves at constant angular velocities, the answer is ## \sqrt{7} .## But don't you think this assumption is somewhat doubtful and wrong?

Similar threads

Replies
2
Views
2K
Replies
4
Views
1K
  • · Replies 8 ·
Replies
8
Views
2K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 2 ·
Replies
2
Views
1K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
6
Views
2K
  • · Replies 12 ·
Replies
12
Views
2K
Replies
6
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K