• Join over 1.2 million students every month
  • Accelerate your learning by 29%
  • Unlimited access from just £6.99 per month

Numerical Method Course Work

Extracts from this document...


Numerical Method Course Work

  1. Problem Specification

Finding all the roots of equation (x+1) sin (1/x)=0 for 0.05<x<1.The reason for why I use numerical methods is I can’t solve this equation by any analytical method. And it lies beyond my mathematics knowledge to solve by an algebraic method.

  1. Strategy

I am not sure which method I will use for finding all roots of the equation. So I want to see which method is most efficient method for finding a root, efficient method means can find the roots most quickly. The methods are shown below,

 Bisection Method with the formula: If a root lies between X = X0 and X = X1, then X2 = 0.5 ( X0 + X1 ) will give an approximation to the roots;

 Secant Method     with the formula: If X = X0 and X1 are approximation to a root of X = X0. A better approximation to the root will usually be given by X = X1 – (X1 - X0 )*f (X1 )/[f(X1 )-f(X0 )];

 Linear Iteration    If x = a and x = b are approximations to a root of f (x)

...read more.


), X0is given. This will usually converge to a root of f (x) = 0 near to x = x0. The method is called the Newton-Raphson method and has second order convergence. (Appendix 5)

In my course work, the equation is  (x+1) sin (1/x) = 0.

So differentiate this function by using chain rule, product rule, and trig. Differentiation.

Y =F (x) = (x+1) sin (1/x))

Y= u * v

        U’ = 1

        V’ = -(1/x^2)* cos. (1/x)   (Chain rule and Trig. Differentiation)

        F ‘(x) = u v’ + v u’  (Product Rule)

                  = Sin (1/x) –  (x+1)/x^2* cos. (1/x)

Then I put the function f (x) = (x+1) sin (1/x) into the graphic calculator, using the table to store the information with the range 0.05 as starting point, 1 as the end point, 0.01 as the pitch. I use the calculator to see how many roots are there roughly. The calculator lists a table  and draws a graph for me. (Appendix 3 and 4)

From the table and graph, we can see that there may be 6 roots in the interval (0.05, 1). There are between 0.05 and 0.06, 0.06 and 0.07, etc. Then I use the same method to check whether or not there is exactly only 1 root in the interval. Instead, I use 0.05 as starting value, 0.

...read more.


0 is correct to 1 decimal place, X1 will probably be correct to 2 decimal places, X2 to 4 decimal places, and so on.

Because Ms Excel spreadsheet just can go to the accuracy as far as 9 decimal places, so my solution just can achieve 8 decimal places, because of rounding in the computer program

  1. Interpretation

The roots of the equation (x+1) sin (1/x) =0 are

* 0.053051648,

* 0.063661977,

* 0.106103295,

* 0.159154943,

* 0.318309886,

* 0.079577472.  (8 decimal places).

Because the Successive errors of all the roots become smaller and turn to 0 at the end, and the En/ En-1^2 ratio becomes constant, I can say my solution is valid.

So the limitation of Newton-Raphson method is I need to choose the roots very carefully otherwise I can’t find the roots I want to find. Particularly to find roots at 0.0795774472, it is very important to choose the correct value.

 And the disadvantage of Secant method is that, like the Newton-Raphson method, can lead to divergence in some cases. (Appendix 7)

And for the fix-iteration, this method seem not so good to the equation (x+1)sin (1/X) =0

The bisection method is good, reliable method for solving an equation f (x) = 0, it always go to the convergence but it reduces the rate of going to convergence, because it halves the interval width.

...read more.

This student written piece of work is one of many that can be found in our AS and A Level Core & Pure Mathematics section.

Found what you're looking for?

  • Start learning 29% faster today
  • 150,000+ documents available
  • Just £6.99 a month

Not the one? Search for your essay title...
  • Join over 1.2 million students every month
  • Accelerate your learning by 29%
  • Unlimited access from just £6.99 per month

See related essaysSee related essays

Related AS and A Level Core & Pure Mathematics essays

  1. Marked by a teacher

    The Gradient Function

    5 star(s)

    the overall gradient function, but I must find far more evidence than one solitary a value to give it any meaning. If this is the case, however, I think that when the co-efficient is 3, the whole gradient will be 3(nxn-1).

  2. The open box problem

    then I will know that this is right. x x 24 24-2x x 24-2x 24 Once again I will use the same equation but change the 15 to a 24 this time so the equation will be V=x(24-2x)^2. Again I will construct tables and graphs to find out the maximum volume and the value of x.

  1. MEI numerical Methods

    If we therefore do the function of the upper and lower bounds, where the function is, y=x+tanx-1, if the root is a correct approximation one function will be negative and the other positive. f(0.4797310075)= 0.0000000005 f(0.4797310065)= -0.0000000018 Therefore this interval obeys the change in sign law, the fact that

  2. Solving Equations. Three numerical methods are discussed in this investigation. There are advantages and ...

    Other 2 root are also found: Target root 2: n xn f(xn) f'(x) f'(xn) f(xn)/f'(xn) xn+1 1 -1 -0.5 3x2+6x+0.5 -2.5 0.2 -1.2 2 -1.2 -0.008 3x2+6x+0.5 -2.38 0.003361345 -1.203361345 3 -1.203361 -6.81716E-06 3x2+6x+0.5 -2.37593 2.86926E-06 -1.203364214 4 -1.203364 -5.02265E-12 3x2+6x+0.5 -2.37593 2.11397E-12 -1.203364 5 -1.203364 0 3x2+6x+0.5 -2.37593 0

  1. Numerical solution of equations, Interval bisection---change of sign methods, Fixed point iteration ---the Newton-Raphson ...

    r Xr f(Xr) f'(Xr) f(Xr)/f'(Xr) 1 2 29 47 0.617021277 2 1.382978723 7.671007388 23.2775011 0.329546 3 1.053432723 1.678766068 13.4149463 0.125141468 4 0.928291256 0.205236732 10.1818519 0.020157112 5 0.908134144 0.004995226 9.68744176 0.000515639 6 0.907618504 3.23625E-06 9.67489018 3.345E-07 7 0.90761817 1.36247E-12 9.67488204 1.40825E-13 Spreadsheet 2.3 From spreadsheet 2.3 above, rounding up the solution to 5 decimal places, I get the solution of 0.90762.

  2. The method I am going to use to solve x&amp;amp;#8722;3x-1=0 is the Change ...

    see that I can see that there is really no convergence to that particular root I am looking for(the one in Rearrangement B and my observation from the graph is correct.

  1. Numerical integration can be described as set of algorithms for calculating the numerical value ...

    The formula can be derived analytically and has the advantage of involving the same amount of arithmetic as the trapezium rule with the addition of providing more accurate results. Below is a diagram of incorporating Simpsons rule: The interval a to b is divided into 2n strips, each of width h.

  2. Solving Equations Using Numerical Methods

    Here is the table of results. I can now see that the root lies between -1.526 and -1.525. As I only intend to find the root to 3 decimal places I don't need to continue the method. To be certain that the root is between -1.526 and -1.525, I will calculate f(x)

  • Over 160,000 pieces
    of student written work
  • Annotated by
    experienced teachers
  • Ideas and feedback to
    improve your own work