• Join over 1.2 million students every month
  • Accelerate your learning by 29%
  • Unlimited access from just £6.99 per month
Page
  1. 1
    1
  2. 2
    2
  3. 3
    3
  4. 4
    4
  5. 5
    5
  6. 6
    6
  7. 7
    7
  8. 8
    8
  9. 9
    9
  10. 10
    10
  11. 11
    11
  12. 12
    12
  13. 13
    13
  14. 14
    14

Numerical solutions of equations

Extracts from this document...

Introduction

Kosturi Ash 12/3

A2 Mathematics Coursework C3

Year 12

Numerical solutions of equations

Solving 0 = x5+x-5 using the “Change Of Sign” Method

The method I will use to solve 0 = x5+x-5 is the Change of Sign Method involving the Decimal Search method. I have drawn this graph using the Autograph Software, and the print screen of this is below:

image00.png

image01.png

From my graph above, I can see that the root of this equation is between x =1 and x = 1.5. The table of x values and f(x) values is shown below. I can work out the f(x) values by substituting the x-values into the equation.

x

1

1.1

1.2

1.3

1.4

1.5

f(x)

-3

-2.28949

-1.31168

0.01293

1.77824

4.09375

From my table of values above, it is clear that the change of sign from negative to positive occurs between x = 1.2 and x = 1.3. So, I can narrow these values down further to find another change of sign.

x

f(x)

1.21

-1.19626

1.22

-1.07729

1.23

-0.95469

1.24

-0.82837

1.25

-0.69824

1.26

-0.56420

1.27

-0.42616

1.28

-0.28403

1.29

-0.13769

1.30

0.01293

I can see that the change of sign is between x = 1.29 and x = 1.30.

x

f(x)

1.291

-0.12283

1.292

-0.10792

1.293

-0.09296

1.294

-0.07797

1.295

-0.06293

1.296

-0.04784

1.297

-0.03271

1.298

-0.01754

1.299

-0.00233

1.300

0.01293

The change of sign is in the interval [1.299, 1.300]

x

f(x)

1.2991

-0.000805

1.2992

0.000720

1.2993

0.002244

The root of this equation lies in the interval [1.2991, 1.2992]. This means that I can take the root to be x = 1.29915(5 decimal places).

...read more.

Middle

3+6x2+3

The iterative formula is: xn+1 = xn- f(xn)

                                       f’(xn)

So the iterative formula for this equation is: xn+1 = xn- xn4+2xn3+3xn-4

                                                              4xn3+6xn2+3

On my diagram of this function, I have taken my x1 value to be about –1.6

x1 = -1.6

x2 = 3.682591

x3 = 2.658976

x4 = 1.899865

x5 = 1.353797

x6 = 1.003365

x7 = 0.847116

x8 = 0.819172

x9 = 0.818390

x10 = 0.818389

x11 = 0.818389

From that, I can quite clearly see that there is a convergence beginning to show towards the positive root. However, my aim was to find the negative root, and I could not find it. Therefore, the Newton-Raphson method has failed to find that particular root, even though the starting value was close to it.

Another example where the Newton-Raphson would not work is the function (x+21/2)1/2 ln (x+21/2) = 0. The graph and calculations of this is below:

image07.png

The function is: (x+21/2)1/2 ln (x+21/2) = 0

f’(x) = ((x+21/2)1/2X 1                 ) + (ln (x+21/2) X ½ (x+21/2) -1/2)

x + 2 1/2

(x+21/2)1/2   + ln (x+21/2) (x+21/2) -1/2

(x+21/2)1                2

        1         + ln(x+21/2)

                  (x+21/2)1/2               2(x+21/2)1/2

=        2+ ln((x+21/2)

                2(x+21/2)1/2

The iterative formula for the Newton-Raphson method is:

        xn+1 = xn = f’(xn)

                      f’(xn)

        x n+1 = xn  - 2(x+21/2)ln(x+21/2)

 2 + ln(x+21/2)

I choose my starting value of x (x1) to beabout –1.3.

x1 = -1.3

x2 = -4.22079

x3 = not defined

x4 =

Another example where the Newton-Raphson method would fail is with the function x3-5x+0.1=0.

...read more.

Conclusion

Decimal Search method.

Both the Newton-Raphson method and the Rearrangementmethod were fixed point estimates and involved an iterative process. Therefore, these methods were very similar. However, these two methods differ because there is a specific formula for the Newton-Raphson method. Although the Newton-Raphson method was more complicated than the Rearrangement method, it gave a much more rapid rate of convergence. In fact, the Newton-Raphson method gave the most rapid rate of convergence to 6 decimal places, whereas, for the Decimal Search method, it took me very many calculations to converge to 5 decimal places. So, the Decimal Search method gave the slowest rate of convergence. These two methods were not as time-consuming as the Decimal Search method because it did not involve having to find a change of sign in many particular intervals; there was already a convergence towards the particular root I wanted. However, it was really easy to make mistakes on the calculator due the order of the terms in the iterative formulae.

I used the Autograph software to draw these graphs. This was so much quicker than drawing the graphs by hand, and Excel does not have the facilities to draw such advanced graphs. The Autograph software was very easy to use and it also helped me to find out which rearrangement of a function were suitable (when carrying out the Rearrangement method) and in which interval the root was in (when using the Decimal Search method).

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

    C3 Coursework - different methods of solving equations.

    5 star(s)

    Fixed point Iteration This method is probably the hardest method to use out of the three. It involves rearranging a y = f(x) formula into x = g(x), by using any one of the x values in the f(x) formula.

  2. Marked by a teacher

    The Gradient Function

    5 star(s)

    Therefore this proves that all positive values of n in the equation nx n-1 apply to this formula. After I try to prove negative powers are correct, I shall thoroughly investigate, in depth, the overall gradient function for the curves y =ax n-1, where "a" is not constant.

  1. Marked by a teacher

    Estimate a consumption function for the UK economy explaining the economic theory and statistical ...

    3 star(s)

    The same conclusion can be gained from the other tests. Such as the p-test, the higher the value of t-test the lower the value of p-test. Thus, the consumption function which incorporating inflation and error correction is successful. The predictions from this equation are shown in Figure 7(a).

  2. C3 Coursework: Numerical Methods

    It was not as tedious as when using the Newton Raphson typing in this formula.

  1. MEI numerical Methods

    Assuming it does we will be able to approximate the root for another value of K. It will only be a approximation as the value of our roots are approximations. To prove that there is a common factor, C we do the following: (assume k=(k+1)

  2. In my coursework I will be using three equations to investigate their solutions using ...

    6.292479904 6.292514666 6.292514666 6.292537349 6.292537349 6.29255215 The error bound is x= 6.29257994 ± 0.000000005 The magnitude of g'(x)

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

    + f(?????h?]?h??????[f(?+2h) + f(?????h?]?h??? + [f(?+3h) + f(?????h?]?h?????? [f(??????[f(??h) + f(?????h????f(?????h?] + f(?????h?] Thus the general approximation to for n trapezia is: Tn= (h/2) [f0 + 2(f1 + f2 + f3 +...+ fn-1) + fn] where f0 is the first strip and is the last fn strip. Notice that f1, f2, f3 ect.

  2. I am going to solve equations by using three different numerical methods in this ...

    f'(Xr) =the value after differentiated f (Xr) By using the Excel, we can easily find the x value in the many terms satisfy my required degree of accuracy which is answer to 7 decimal places. Root is -1.694532 to 6d.p. Error bounds is -1.69453±0.000005. Root bounds is -1.694535<x<-1.694525. Check X Y -1.694535 -0.00011128 (negative) -1.694525 0.000147178 (positive)

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