• 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...

Introduction

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.

Middle

), 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.

Conclusion

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)

    3(x+h)� - x� = 3(x� +3x�h + 3h�x+h�) -3x� = 3h(3x� +3hx + h�) = x + h - x h h 3(3x� +3hx + h�) = 9x� + 9hx +3h� = 9x�* - "h" limits to 0. axn value gradient function x3 3x2 2x3 6x2 3x3 9x2 This pattern

  2. MEI numerical Methods

    So therefore the higher the rate of convergence the quicker the method reaches the final approximation. First order of convergence: The formula to work out the first order of convergence is: Absolute error in x1 / Absolute error in x2 In order to find absolute error, we usually do (x(approximation)

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

    f(Xr)/ f'(Xr) 1 -1 20 -10 -2 2 1 -4 -10 0.4 3 0.6 0.416 -11.92 -0.034899329 4 0.634899329 0.00223484 -11.79070853 -0.000189542 5 0.635088871 6.84356E-08 -11.78998638 -5.80456E-09 6 0.635088877 0 -11.78998635 0 7 0.635088877 0 -11.78998635 0 Spreadsheet 4.3 Newton-Raphson has found the root 0.63509 as well.

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

    4x2 �+3x2 �)] = -1.454113738 x4 = x3 - [(x3 ^4+x �-1=0)/( 4x3 �+3x3 �)] = -1.387577585 x5 = x4 - [(x4 ^4+x4 �-1=0)/( 4x4 �+3x4 �)] = -1.380357406 x6 = x5 - [(x5 ^4+x5 �-1=0)/( 4x5 �+3x5 �)] = -1.380277579 x7 = x6 - [(x6 ^4+x6 �-1=0)/( 4x6 �+3x6

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

    This is because the mid-point overestimates whilst the trapezium rule underestimates. Therefore, by using both methods we can confine the exact value between the two different polynomials. This means by obtaining values of Tn and Mn where n is a large number, we can find an accurate approximation to the

  2. C3 Coursework: Numerical Methods

    The root at the point B appears to be between 0 and 1. As a result of this x=0 would be a suitable starting point. The table below shows the result to the Newton Raphson Formula when using 0 as the starting point.

  1. Solving Equations Using Numerical Methods

    Although I have been successful in finding this root for this equation, sometimes the decimal search method fails to find the root. I will now show an example of the decimal search method failing and find out why it fails.

  2. Numerical Differentiation

    Or using the error formula: Error = X - x ?= 2.3921495 the error in this approximation is 5.0x107 When change of sign does not work If I translate my curve that I have used in the last part vertically 4 units, I get the new equation: y=x�-5.5x�+2x+17 As you

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