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

    4 0.25 4.1 0.243902439 0.059488 4 0.25 4.01 0.249376559 0.062189 4 0.25 4.001 0.249937516 0.062469 2 0.5 2.1 0.476190476 0.226757 2 0.5 2.01 0.497512438 0.247519 2 0.5 2.001 0.499750125 0.24975 1 1 1.1 0.909090909 0.826446 1 1 1.01 0.99009901 0.980296 1 1 1.001 0.999000999 0.998003 x x-1 Gradient 1 1

  1. Marked by a teacher

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

    3 star(s)

    Thus the residual in the Figure 7(b) starts show a pattern. There are many reasons, which used to explain the fail prediction of model. The most common two are house price rise and credit liberalization, and reduced uncertainty. As the credit liberalization, which made household easier to borrow money to buy house.

  2. C3 Coursework: Numerical Methods

    This is also because the formula was more complex. Like the previous methods, I used Microsoft Excel when using the x=g(x) method. This method, like the others, was simple to use with Microsoft Excel. It was simple to type in the formula I used as well as the replicate the same formula.

  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. Numerical integration can be described as set of algorithms for calculating the numerical value ...

    Notice that f1, f2, f3 ect. to fn-1 do double duty as right-hand and left-hand body. Note that is the dependent on the number of trapezia. Thus T16 means that area under the curve is split in 16 rectangles whereas T32 means its split up in 32.

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

    x y 4 4.687548148 4.687548148 5.195135269 5.195135269 5.554033303 5.554033303 5.800790441 5.800790441 5.967384127 5.967384127 6.078530709 6.078530709 6.152113534 6.152113534 6.200582955 6.200582955 6.232405254 6.232405254 6.253253236 6.253253236 6.266892437 6.266892437 6.275807357 6.275807357 6.281630902 6.281630902 6.285433572 6.285433572 6.287916018 6.287916018 6.289536331 6.289536331 6.290593809 6.290593809 6.291283911 6.291283911 6.291734245 x y 6.29255215 6.292561807 6.292561807 6.292568109 6.292568109 6.29257222 6.29257222 6.292574903

  2. Methods of Advanced Mathematics (C3) Coursework.

    -0.04277 -2.493 -0.02912 -2.492 -0.01548 -2.491 -0.00186 -2.49 0.011751 Newton-Raphson xn f(xn) f'(xn) xn+1 -2 5 7 -2.714285714 -2.71428571 -3.425655977 17.10204082 -2.513978861 -2.51397886 -0.318677638 13.96026914 -2.49115139 -2.49115139 -0.003918158 13.61750575 -2.490863661 -2.49086366 -6.18691E-07 13.61320533 -2.490863615 The Newton raphson method shows the route to be -2.490863615 ?

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