# Numerical Methods/Equation Solving

## Contents

## Solution of Algebraic and Transcendental Equations[edit]

An equation of the type is either algebraic or transcendental.

E.g, these equations are algebraic.

and these are transcendental

While roots can be found directly for algebraic equations of fourth order or lower, and for a few special transcendental equations, in practice we need to solve equations of higher order and also arbitrary transcendental equations.

As analytic solutions are often either too cumbersome or simply do not exist, we need to find an approximate method of solution. This is where numerical analysis comes into the picture.

### Some Useful Observations[edit]

- The total number of roots an algebraic equation can have is the same as its degree.
- An algebraic equation can have at most as many positive roots as the number of changes of sign in .
- An algebraic equation can have at most as many negative roots as the number of changes of sign in .
- In an algebraic equation with real coefficients, complex roots occur in conjugate pairs
- If with roots then the following hold good:
- If is continuous in the interval and then a root must exist in the interval

### Initial Approximation[edit]

The last point about the interval is one of the most useful properties numerical methods use to find the roots. All of them have in common the requirement that we need to make an initial guess for the root. Practically, this is easy to do graphically. Simply plot the equation and make a rough estimate of the solution. Analytically, we can usually choose any point in an interval where a change of sign takes place. However, this is subject to certain conditions that vary from method to method.

### Convergence[edit]

A numerical method to solve equations will be a long process. We would like to know, if the method will lead to a solution (close to the exact solution) or will lead us away from the solution. If the method, leads to the solution, then we say that the method is convergent. Otherwise, the method is said to be divergent.i.e, in case of linear and non linear interpolation convergence means tends to 0.

### Rate of Convergence[edit]

Various methods converge to the root at different rates. That is, some methods are slow to converge and it takes a long time to arrive at the root, while other methods can lead us to the root faster. This is in general a compromise between ease of calculation and time.

For a computer program however, it is generally better to look at methods which converge quickly. The rate of convergence could be linear or of some higher order. The higher the order, the faster the method converges.

If is the magnitude of the error in the th iteration, ignoring sign, then the order is if is approximately constant.

It is also important to note that the chosen method will converge only if .

### Bisection Method[edit]

This is one of the simplest methods and is strongly based on the property of intervals. To find a root using this method, the first thing to do is to find an interval such that . Bisect this interval to get a point . Choose one of or so that the sign of is opposite to the ordinate at that point. Use this as the new interval and proceed until you get the root within desired accuracy.

##### Example[edit]

**Solve correct up to 2 decimal places.**

#### Error Analysis[edit]

The maximum error after the th iteration using this process will be given as

As the interval at each iteration is halved, we have . Thus this method converges linearly.

If we are interested in the number of iterations the *Bisection Method* needs to converge to a root within a certain tolerance than we can use the formula for the maximum error.

##### Example[edit]

**How many iterations do you need to get the root if you start with a = 1 and b = 2 and the tolerance is 10^{−4}?**

The error needs to be smaller than 10^{−4}. Use the formula for the maximum error:

Solve for *i* using log rules

Hence 14 iterations will ensure an approximation to the root accurate to . Note: the error analysis only gives a bound approximation to the error; the actual error may be much smaller.

### False Position Method[edit]

The false position method (sometimes called the regula falsi method) is essentially same as the bisection method -- except that instead of bisecting the interval, we find where the chord joining the two points meets the X axis. The roots are calculated using the equation of the chord, i.e. putting in

- .

The rate of convergence is still linear but faster than that of the bisection method. Both these methods will fail if *f* has a double root.

#### Example[edit]

Consider *f*(*x*)=*x*^{2}-1. We already know the roots of this equation, so we can easily check how fast the regula falsi method converges.

For our initial guess, we'll use the interval [0,2].

Since *f* is concave upwards and increasing, a quick sketch of the geometry shows that the chord will always intersect the *x*-axis to the left of the solution. This can be confirmed by a little algebra.

We'll call our *n*^{th} iteration of the interval [*a*_{n}, 2]

The chord intersects the *x*-axis when

Rearranging and simplifying gives

Since this is always less than the root, it is also *a*_{n+1}

The difference between *a*_{n} and the root is *e*_{n}=*a*_{n}-1, but

This is always smaller than *e*_{n} when *a*_{n} is positive. When *a*_{n} approaches 1, each extra iteration reduces the error by two-thirds, rather than one-half as the bisection method would.

The order of convergence of this method is 2/3 and is linear.

In this case, the lower end of the interval tends to the root, and the minimum error tends to zero, but the upper limit and maximum error remain fixed. This is not uncommon.

### Fixed Point Iteration (or Staircase method)[edit]

If we can write *f*(*x*)=0 in the form *x*=*g*(*x*), then the point *x* would be a __fixed point__ of the function *g* (that is, the input of *g* is also the output). Then an obvious sequence to consider is

If we look at this on a graph we can see how this could converge to the intersection.

Any function can be written in this form if we define *g*(*x*)=*f*(*x*)+*x*, though in some cases other rearrangements may prove more useful.

#### Error analysis[edit]

We define the error at the *n*^{th} step to be

Then we have

So, when |*g*′(*x*)|<l, this sequence converges to a root, and the error will be approximately proportional to *n*.

Because the relationship between *e*_{n+1} and *e*_{n} is *linear*, we say that this method *converges linearly*, if it converges at all.

When *g*(*x*)=*f*(*x*)+*x* this means that if

- .

converges to a root, *x*, of *f* then

Note that this convergence will only happen for a certain range of *x*. If the first estimate is outside that range then no solution will be found.

Also note that although this is a necessary condition for convergence, it does not guarantee convergence. In the error analysis we neglected higher powers of *e*_{n}, but we can only do this if *e*_{n} is small. If our initial error is large, the higher powers may prevent convergence, even when the condition is satisfied.

If |*g*′(*x*)|<1 is true at the root, the iteration sequence will converge in some interval around the root, which may be smaller than the interval where |*g*′(*x*)|<1. If |*g*′(*x*)| isn't smaller than one at the root, the iteration will not converge to that root.

#### Example[edit]

Lets consider , which we can see has a single root at *x*=1. There are several ways *f*(*x*)=0 can be written in the desired form, *x*=*g*(*x*).

The simplest is

1):

In this case, , and the convergence condition is

Since this is never true, this doesn't converge to the root.

2)An alternate rearrangement is

This converges when

Since this range does not include the root, this method won't converge either.

3)Another obvious rearrangement is

In this case the convergence condition becomes

Again, this region excludes the root.

4)Another possibility is obtained by dividing by *x*^{2}+1

In this case the convergence condition becomes

Consideration of this inequality shows it is satisfied if *x*>1, so if we start with such an *x*, this will converge to the root.

Clearly, finding a method of this type which converges is not always straightforwards.

### Newton-Raphson[edit]

In numerical analysis, Newton's method (also known as the Newton–Raphson method or the Newton–Fourier method) is an efficient algorithm for finding approximations to the zeros (or roots) of a real-valued function. As such, it is an example of a root-finding algorithm.

Any zero-finding method (Bisection Method, False Position Method, Newton-Raphson, etc.) can also be used to find a minimum or maximum of such a function, by finding a zero in the function's first derivative, see Newton's method as an optimization algorithm.

#### Description of the method[edit]

The idea of the Newton-Raphson method is as follows: one starts with an initial guess which is reasonably close to the true root, then the function is approximated by its tangent line (which can be computed using the tools of calculus), and one computes the x-intercept of this tangent line (which is easily done with elementary algebra). This x-intercept will typically be a better approximation to the function's root than the original guess, and the method can be iterated. Suppose f : [a, b] → R is a differentiable function defined on the interval [a, b] with values in the real numbers R. The formula for converging on the root can be easily derived. Suppose we have some current approximation xn. Then we can derive the formula for a better approximation, xn+1 by referring to the diagram on the right. We know from the definition of the derivative at a given point that it is the slope of a tangent at that point.

We can get better convergence if we know about the function's derivatives. Consider the tangent to the function:

Near any point, the tangent at that point is approximately the same as *f*('x*) itself, so we can use the tangent to approximate the function.*

The tangent through the point (*x*_{n}, *f*(*x*_{n})) is

The next approximation, *x*_{n+1}, is where the tangent line intersects the axis, so where *y*=0. Rearranging, we find

#### Error analysis[edit]

Again, we define the root to be *x*, and the error at the *n*^{th} step to be

Then the error at the next step is

- (1)

where we've written *f* as a Taylor series round its root, *x*. Rearranging this, and using *f*(*x*)=0, we get

where we've neglected cubic and higher powers of the error, since they will be much smaller than the squared term, when the error itself is small.

Notice that the error is squared at each step. This means that the number of correct decimal places *doubles* with each step, much faster than linear convergence.

This sequence will converge if

If *f*′ isn't zero at the root, then there will always be a range round the root where this method converges.

If *f*′ is zero at the root, then on looking again at (1) we see that we get

and the convergence becomes merely linear.

Overall, this method works well, provided *f* does not have a minimum near its root, but it can only be used if the derivative is known.

#### Example[edit]

Let's consider *f*(*x*)=*x*^{2}-*a*. Here, we know the roots exactly, so we can see better just how well the method converges.

We have

This method is easily implemented, even with just pen and paper, and has been used to rapidly estimate square roots since long before Newton.

The *n*^{th} error is *e*_{n}=*x*_{n}-√*a*, so we have

If *a*=0, this simplifies to *e*_{n}/2, as expected.

If *a*>0, *e*_{n+1} will be positive, provided *e*_{n} is greater than -√*a*, i.e provided *x*_{n} is positive. Thus, starting from any positive number, all the errors, except perhaps the first will be positive.

The method converges when,

so, assuming *e*_{n} is positive, it converges when

which is always true.

This method converges to the square root, *starting from any positive number*, and it does so quadratically.

### Higher order methods[edit]

There are methods that converge even faster than Newton-Raphson. e.g

which converge cubicly, tripling the number of correct digits at each iteration, which is 50% faster than Newton-Raphson.

However, if iterating each step takes 50% longer, due to the more complex formula, there is no net gain in speed. For this reason, methods such as this are seldom used.