• 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

    The Gradient Function

    5 star(s)

    -1 2 0.5 -0.25 3 0.333333 -0.111111111 4 0.25 -0.0625 I predict that the gradient function for when n=-1 = -x -2, in correspondence to the formula nx n-1. I can tell whether these data do in fact agree with this gradient function, but I shall binomially expand this to verify this is the case: 1 .

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

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

    This can be shown graphically on the Graph 3.7. Not knowing the magnitude of the gradient, the value of x taken is near to the approximation---2. The slope of the graph at x = 4 is g' (4) which is 3.692307692..., which is larger than 1. Spreadsheet 3.6 The magnitude of the rearrangement is quite good as you have

  2. MEI numerical Methods

    -4 -0.31838296 -5 -0.24381388 -6 -0.19689847 The relationship between the negative integer K and the root ? is: As k-> -? then ? -> 0 Essentially this means that there is a positive correlation, as the value of K increases the value of ?

  1. In this piece of coursework, I will use three different methods (Change of Sign, ...

    and so this method will not work. The gradient may be found by entering the value of '2' into the equation found when differentiating the equation used in Rearrangement 1: Rearrangement 1 equation: y = 4 - x3 3 Differentiated: dy/dx = -x2 Putting x = 2 into the equation:

  2. Numerical Method (Maths Investigation)

    It also gives my variety of fonts and sizes that I may prefer to. Neat Alignment can help my Coursework to look neat and tidy. I can even get an eye of my coursework before I get it print out by PREVIEW it, making sure that I am satisfies with it.

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

    -1.735649917 -1.735649917 -1.09912868 27.62106745 -1.695856789 -1.695856789 -0.03431212 25.90287655 -1.694532143 -1.69453214 -3.745E-05 25.8463345 -1.69453069 -1.694530694 -4.48E-11 25.84627268 -1.694530694 -1.694530694 3.55271E-15 25.84627268 -1.694530694 -1.694530694 -7.1054E-15 25.84627268 -1.694530694 -1.694530694 3.55271E-15 25.84627268 -1.694530694 In this spreadsheet, Xr = estimated value Xr. f'(Xr) =the value after differentiated f (Xr)

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

    from this fix point and eventually it will converge to the roots. This method requires us to start with a point which is close to the target root and we call it. Then we will find the gradient - at Cartesian co-ordinate and we can find the better estimation -

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