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

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

    0.688896179 -6.65815E-06 0.688892365 -3.03145E-07 3.8147E-06 19 0.68888855 6.05186E-06 0.688892365 -3.03145E-07 0.688890457 2.87436E-06 1.90735E-06 20 0.688890457 2.87436E-06 0.688892365 -3.03145E-07 0.688891411 1.28561E-06 9.53674E-07 Newton Raphson's method n xn f(xn) f'(xn) f(xn)/f'(xn) xn+1 1 1 -0.5 -1.5 0.333333333 0.666666667 2 0.666666667 0.037037037 -1.66667 -0.022222222 0.688888889 3 0.688888889 5.48697E-06 -1.66593 -3.29364E-06 0.688892183 4 0.688892183

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

    0.0038113 0.001953125 10 0.634765625 0.00381134 0.6367188 -0.0192111 0.635742188 -0.007702 0.000976563 11 0.634765625 0.00381134 0.6357422 -0.0077017 0.635253906 -0.001946 0.000488281 12 0.634765625 0.00381134 0.6352539 -0.0019456 0.635009766 0.0009327 0.000244141 13 0.635009766 0.00093274 0.6352539 -0.0019456 0.635131836 -0.000506 0.00012207 14 0.635009766 0.00093274 0.6351318 -0.0005065 0.635070801 0.0002131 6.10352E-05 15 0.635070801 0.00021312 0.6351318 -0.0005065 0.635101318 -0.000147 3.05176E-05

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

    The use of polynomials: The mid-point rule is a second order polynomial which overestimates convex curves. The function I am integrating also follows a convex curve. As a result I expect the consecutive values of the mid-point to decrease as this polynomial overestimates.

  2. C3 Coursework: Numerical Methods

    This means that x=2 would be a suitable starting point. The table below shows the result to the Newton Raphson formula when using 2 as the starting point. I shall now use the Change of Sign method to establish the error bounds.

  1. Solving Equations Using Numerical Methods

    I will use the equation y=1.676x� + 3.018x� -1.851x + 0.2427. Here is the equation. I will try to find the middle root of the equation using the decimal search method. I will start with the integer immediately below the root (0)

  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