Newton-Raphson Method used in a real engineering/physical/mathematical situation

I've been using the Newton-Raphson Method in my Numerical Methods course for a while now, blindly solving non-linear equations and systems of equations . This makes me somehow lose motivation, as I can't manage to find a real problem in which this tool is applied .

Can anyone suggest such a problem in an engineering/physical or mathematical environment ? I'm interested in seeing a meaningful equation related to a possible problem and the motivation behind finding its roots.

• The method solves the vast majority of all variable equations. The problem is not that they don't exist, you just are looking right past them. (equations make up 99% of physics and 1% of engineering (just kidding!)) – Zach466920 Apr 24 '15 at 21:49
• I wasn't implying that such problems don't exist, it's pretty obvious this is not the case since they're teaching this method centuries later. I'm just saying that I haven't encountered such a situation on my own and I would like to see one . Feel free to make a suggestion :) – Victor Apr 24 '15 at 21:53
• One powerful way to get nonlinear algebraic equations is from taking Fourier transforms of nonlinear differential equations. For instance, in my undergrad work I dealt with an equation like $(u_{xx} + a F(u) + b u)_{xx}=0$. In my setting $F$ was a nonlinear function. In Fourier space this is an algebraic equation $K^4 \hat{u} + K^2 a \hat{F}(\hat{u}) + K^2 b \hat{u}=0$. Here $K$, at least when dealing with a finite Fourier expansion (for a numerical approximation) is a diagonal matrix. – Ian Apr 24 '15 at 21:53
• You can use them with Kirchhoff's Laws for electronics and for Triangulation for determining location. – Zach466920 Apr 24 '15 at 22:16

Just to give you a simple one, consider Colebrook equation $$\frac 1{\sqrt f}=-2\log_{10}\Big(\frac{\epsilon}{3.71D_h}+\frac{2.51}{Re \sqrt f}\Big)$$ in which $f$ is the friction factor, $\epsilon$ the roughness height, $D_h$ the hydraulic diameter and $Re$ the Reynolds number. This friction factor determines the pressure drop in a pipeline and I suppose that you perceive the impact of that when we speak about gas or oil pipelines.
Even if we know an explicit solution of this equation in terms of Lambert function, most of the time, this equation is solved using Newton method. But, in order to save iterations and ensure stability, we use to find proper transforms and "good" starting values. For example, changing variables $f=\frac 1{x^2}$ leads to $$\Phi(x)=x+2\log_{10}\Big(\frac{\epsilon}{3.71D_h}+\frac{2.51}{Re} x\Big)=0$$ Dozens of papers have been published proposing unexpensive ways for the generation of the initial estimate; these make the convergence of Newton method very fast to high accuracy.
More difficult is the case of Underwood equations $$\Phi(x)=\sum_{i=1}^n \frac {a_i}{b_i-x}-c=0$$ for which all the roots must be found ($n$ can be several hundreds). This equation is also solved zillions of time in chemical process simulations. Hundreds of papers addressed how to estimate the starting points, how to transform the equation to ensure proper convexity but by the end, Newton--Raphson does the job.