• 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

OCR MEI C3 Coursework - Numerical Methods

Extracts from this document...

Introduction

C3 Coursework                Raoul Harris

C3 Coursework

Change of sign (decimal search)

Finding a root of an equation with graphical illustration

f(x)=x3-7x2+2x+1

This is graph of y=x3-7x2+2x+1

image00.png

The graph shows that roots of f(x)=0 exist in the intervals [-1,0]; [0,1]; [6,7]

We shall test for a root in the interval [6,7]

image01.png

x

f(x)

6.0

-23

6.1

-20.289

6.2

-17.352

6.3

-14.183

6.4

-10.776

6.5

-7.125

6.6

-3.224

6.7

0.933

Change of sign indicates root exists in interval [6.6,6.7]

This means that x=6.65±0.05image12.png

x=7 (0d.p.)

x

f(x)

6.60

-3.224

6.61

-2.81992

6.62

-2.41327

6.63

-2.00405

6.64

-1.59226

6.65

-1.17787

6.66

-0.7609

6.67

-0.34134

6.68

0.080832


Change of sign indicates root in interval [6.67,6.68]

This means that x=6.675±0.005

x=6.7 (1d.p.)

image13.png

x

f(x)

6.670

-0.34134

6.671

-0.29924

6.672

-0.25711

6.673

-0.21496

6.674

-0.17278

6.675

-0.13058

6.676

-0.08835

6.677

-0.04609

6.678

-0.00381

6.679

0.038498

Change of sign indicates root in interval [6.678,6.679]

This means that x=6.6785±0.0005

x=6.68 (2d.p.)image14.png

X

f(x)

6.6780

-0.00381

6.6781

0.000419

Change of sign indicates root in interval [6.6780,6.6781]

This means that x=6.67805±0.00005

x = 6.678 (3d.p.)

Failure of the decimal search

Two of the roots of f(x)=0 where f(x)=5x3-20x2+2x-0.05 could not be found with this method. This is because they are so close together that there is no sign change between f(0) and f(0.1). The graph on the next page shows the graph of y=f(x)

...read more.

Middle

Finding the remaining root

The other root can be found using x1=-2.

xn

f(xn)

f'(xn)

x1

-2

-2.12500

30.27344

x2

-1.92981

-0.87732

10.43071

x3

-1.84570

-0.29110

4.65311

x4

-1.78314

-0.05792

2.98408

x5

-1.76373

-0.00339

2.64527

x6

-1.76245

-0.00001

2.62496

x7

-1.76244

0.00000

2.62488

x8

-1.76244

0.00000

2.62488

x=2.62488 (5d.p.)

Failure of the Newton-Raphson method

The method fails for the same function if x1=0

xn

f(xn)

f'(xn)

x1

0

0.60317

-0.14172

x2

4.25600

1.05148

-0.02495

x3

46.40573

1.00046

-0.00002

x4

Value could not be found

1.00000

0.00000

The value diverges, despite the starting value being close to the root. This is because the tangent crosses the asymptote, as shown below:

image04.png


Fixed point iteration with x=g(x)

Finding a root

f(x)=x5-4x+3

The graph of y=f(x):

image05.png

We shall rearrange the equation into the form x=g(x) to find roots of f(x)=0.

0= x5-4x+3

x5=4x-3

x=(4x-3)1/5

xn+1=(4xn-3)(1/5)

We shall take a starting value of x1=-1.

x1

-1

x2

-1.47577

x3

-1.54849

x4

-1.55848

x5

-1.55983

x6

-1.56001

x7

-1.56004

x8

-1.56004

x=-1.56004 (5d.p.)

Graphical illustration

The graph below shows the lines y=x and y=g(x). The lines intersect at the roots of f(x)=0.

image06.pngimage07.png

image08.png

The importance of the magnitude of g’(x)

For a root of f(x)

...read more.

Conclusion

The Newton-Raphson method is more difficult to automate because f'(x) must be found. Some computer programs, such as Autograph, can carry out the iterations for you, which, if they are available to you, can make it easier to use than a decimal search. Care must be taken when choosing a starting value and asymptotes can cause the method to fail. In some cases you may not be able to differentiate f(x).

Fixed point iteration with x=g(x) can be carried using programs like Autograph. It must still be rearranged manually though, and a large proportion of rearrangements fail. Because of this it is the most difficult method to use, especially if you do not have software to automate the iterations.

Overall the change of sign method is the easiest as it is purely numerical calculation. Newton-Raphson is the second simplest because it fails far less often than x=g(x).

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

    hx(x+h) hx(x+h) x(x+0)* x2 *h tends to 0. Voila, this gradient function is in accordance with nx n-1. By getting this far, I have nearly convinced myself that all values work, whether they are integers or not, or positive or not.

  2. Methods of Advanced Mathematics (C3) Coursework.

    I find the solutions numerically it only shows one change of sign being the first one and the second is over looked. The routes are visible on the graph but the table below shows only the one change of sign in red.

  1. Borders Coursework

    (4/3) 102/3 (32/3) 36 (108/3) 851/3 (256/3) 1662/3 (500/3) Difference between sequence and 11/3 n3 / 4/3 n3 -1/3 (-1/3) -32/3 (-11/3) -11 (-33/3) -221/3 (-67/3) -372/3 (-113/3) New sequence: -1/3 -11/3 -33/3 -67/3 -113/3 -10/3 -22/3 -34/3 -46/3 -12/3 -12/3 -12/3 General Quadratic equation: bn2 + cn + d 2b= -12/3 (-4)

  2. C3 Coursework: Numerical Methods

    This gives a second estimate of 1.63. It is possible to iterate again to find a more accurate estimate of the root at point C. The tangent to the point x=1.63 crosses the x axis at the point x=1.472. This gives a third estimate of 1.474. The tangent to the point 1.474 crosses the x axis at the point 1.442.

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

    3.61711E-13 -1.66593 -2.17123E-13 0.688892183 5 0.688892183 0 -1.66593 0 0.688892183 6 0.688892183 0 -1.66593 0 0.688892183 , Rearrange into form n xn g(xn) g'(xn) 1 0.500000 0.812500 -0.625 2 0.812500 0.608032 -0.63477 3 0.608032 0.742693 -0.66151 4 0.742693 0.653240 -0.658 5 0.653240 0.712654 -0.6664 6 0.712654 0.673094 -0.66349 7 0.673094

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

    Comparing those three methods Speed of convergence After applying those three methods, I then know that the Bisectional method takes 18 steps, shown on spreadsheet 4.2, which is the longest procedure, to get the answer. On the other hand, Newton-Raphson method proceeds more quickly than that of bisection, which takes only 6 steps, shown on spreadsheet 4.3.

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

    #NUM! 7 #NUM! #NUM! #NUM! 8 #NUM! The root found is not correct, as the gradient is greater than 1; therefore, it is impossible to use this method to find the root. Comparison of methods: I have finished investigating three methods. In the whole process I found I can solve the non-trivial equation by these three methods properly.

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

    + f2n-1) + 2(f2 + f4 + f6 + .... + f2n-2)] However, seeing as I am carrying out the mid-point rule then the trapezium rule, I will be using a different version of Simpson's rule (2Mn + Tn)/3. This is simply due to the fact that it involves less arithmetic than using than using the formula stated in the page above.

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