Finding common roots of two multivariate equations - Computational Science Stack Exchange most recent 30 from scicomp.stackexchange.com 2019-09-20T07:53:37Z https://scicomp.stackexchange.com/feeds/question/24483 https://creativecommons.org/licenses/by-sa/4.0/rdf https://scicomp.stackexchange.com/q/24483 2 Finding common roots of two multivariate equations user2272296 https://scicomp.stackexchange.com/users/7837 2016-07-19T16:28:03Z 2016-07-19T21:13:47Z <p>I am trying to find points of a surface whose normal is parallel with the 3 dimensional axes (i.e. local maxima, minima and saddle points in x, y and z dimensions).</p> <p>The function of the surface is defined as $f(u,v)$ which is continuous and gives a 3D vector indicating the position of the surface for a given set of parameters u and v.</p> <p>To find the points I require I understand that I need to find the derivatives of the surface with respect to u and v. The relevant points will occur when:</p> <p>$\frac{\delta f_i(u,v)}{\delta u}=0$ &amp; $\frac{\delta f_i(u,v)}{\delta v}=0$</p> <p>where $i$ is the respective dimension (x,y or z) and $0\leq u\leq1$ &amp; $0\leq v\leq1$.</p> <p>What kind of strategy should I take to solving this problem? I was thinking it may be possible to use a numerical method with the partial derivatives to iteratively find solutions. But then how would I remove the solution as you would in univariate root finding?</p> <p>Any insights would be really useful.</p> https://scicomp.stackexchange.com/questions/24483/-/24484#24484 1 Answer by Futurologist for Finding common roots of two multivariate equations Futurologist https://scicomp.stackexchange.com/users/20964 2016-07-19T21:13:47Z 2016-07-19T21:13:47Z <p>The derivatives $\partial_u f(u,v)$ and $\partial_v f(u,v)$ are two vectors tangent to yor surface that in general span the tangent space to the surface at the point $f(u,v)$. Then the normal vector to the surface at that point is $n(u,v)=\partial_u f(u, v)\times\partial_v f(u,v)$ and you want it parallel to say the coordinate direction $e_z=(0,0,1)^T$. Then yo want to solve the equation $n(u,v)\times e_z=(\partial_u f(u, v)\times\partial_v f(u,v))\times e_z=0.$ After applying a simplification formula for a double cross-product like this one, you get equations like the ones you have written. To solve the system of equations numerically maybe you can simply use Newton's method or something of that kind. </p>