# Newton & Square Roots

The web is an interesting place, it humbles you down and shows how little you actually know. While hitchhiking through the hyperlinks, I had one such realization recently.

What if I ask you to calculate the square root of a number? Depending upon your background you are likely to come up with a variety of solutions. Let me talk about that specific fraction of you who will be interested in churning out code for doing so. *[No, sqrt()]*

*Binary Search* seems to be a valid and intuitive approach, after all, we can easily discard half of the search space in each iteration. But I want you to hold that thought and glance over this small snippet below.

```
double sq(double num)
{
double xn = num/2;
while(abs((xn*xn)-num)>eps) //eps is precision
{
xn = (xn + num/xn)/2;
}
return xn;
}
```

Magical isn't it? When I first looked at this, I was really startled by the elegance of this solution. So what is it and why does it work?

## Sweet Mathematics

You just witnessed the glamour of ** Newton Raphson** method.

In simple terms, it is a way to quickly find a good approximation for the roots of a real-valued function

**f(x) = 0**.

This is our cue, square root for a number can be minimally represented by solving for the roots of **f(x) = x^2- N**, where *N* is the number we have to find the root of.

That is all cool! But how does it actually work? Let's look at the simple math behind this:-

Consider a general function:

**y = f(x)**We start by choosing an arbitrary value**x1**as our starting point[1]. For our specific case, we can choose**N/2**as our starting point and can safely discard values greater than that from the candidature of**√N**.Now let's draw a tangent (Differential) to the curve at that point and mark its intersection with the x-axis. Let's call this

**x2**.

Repeat

*step 2*for point**x2**to obtain**x3**. You must have noticed that with each iteration we are nearing the roots of the quadratic equation. Slow and steady!*Depending upon the precision, you can decide upon the number of iterations required.*

## Formal Expression

Let's try to workout an elegant formula for **f(x) = x^2 - N** which will get us the compact code we saw in the beginning. Please note that we are aiming to get a relation between **N** and **(N+1)** iteration.

- For any point
**xn**, by using the equation of line we can get:**f'(x[n]) = (f(x[n+1]) - f(x[n]))/(x[n+1]- x[n])**

[*f'(x)*being the slope at point*x*]

- Putting
**0**for y coordinate at point**x[n+1]**where tangent crosses x-axis we get:**f'(x[n]) = (- f(x[n]))/(x[n+1]- x[n])**

- Rearranging the terms & putting
**f'(x) = 2x[n]**we have our final expression:**x[n+1] = (x[n]+ N/x[n])/2**

Visualising Newton Raphson method for N = 10000

## Notes

[1] For the sake of simplicity, I have omitted the details about choosing the initial value which plays a crucial role in more complex functions.

*Special thanks to AMSI for facilitating illustrations*

If you liked this post, you can find more by:

Following me on

Connecting on

Visiting my

website: https://abhimanyutimes.com

Thanks for reading!