• Join over 1.2 million students every month
  • Accelerate your learning by 29%
  • Unlimited access from just £6.99 per month

numerical solutions-Comparison of the three methods and Newton Raphson

Extracts from this document...

Introduction

Numerical Solutions of equations

  1. Newton Raphson method

Equation to be solved is x³−5x−1=0  The function f(x)= x³−5x−1is shown below

image00.png

There are 3 roots. I will first find the root in the interval [2, 3]

I will do the first few lines of calculation manually.

The formula to use is:  xn+1= xn –f(xn)/ f’(xn)

Therefore I must first differentiate x³−5x−1 which is 3x2-5

Using x1=3

X2:   3 – [(33-5 x 3-1)/(3 x 32 -5)] = 2.5

X3:   2.5 – [(2.53-5 x 2.5-1)/(3 x 2.52 -5)] = 2.3455

I will now work out all 3 roots using autograph until 5 significant figures are guaranteedimage01.pngimage06.pngimage07.png

The 3 boxes above show how I obtained the 3 roots of the equation x³−5x−1=0    

Which are -0.20164, -2.1284 and 2.3301

Below is the function f(x)= x³−5x−1 showing where I applied the

...read more.

Middle

image10.png

image11.png

Error bounds for interval [2, 3]

2.3301 is the root in this interval to 5 significant figures.

Therefore the error is 2.3301 ± 0.00005 I will now perform the change of sign test to confirm it is within these limits.

Lower limit is 2.33005 then f (2.33005) = -0.000098648

Upper limit is 2.33015 then f (2.33015) = 0.0010302

There is a change of sign which confirms root in interval is 2.3301 ± 0.00005

When does this method fail?

The Newton Raphson method does not always work, I will show this Using the equation (5x+4)1/7 =0

This function is shown below with autographs attempt to find the root using tangents to the curve starting at x= -1

The function f(x) = (5x+4)1/7

image12.png

As you can see from the graph and the Values of x in the table, it was unable to find the root in the interval [-1, 0].image13.png

This is because of the gradient of the curve, where my starting value is very close to a turning point. Then rather than converging towards the root, it is diverging further and further away on the x axis, therefore overflows and is unable to find the root.

Comparison of methods

...read more.

Conclusion

Comparing Ease of use of Hardware and software

It’s very obvious that the hardware and software speeds up the process of working out the solutions dramatically, as manually is very time consuming and easier to make mistakes.

Autograph is very easy to use and therefore Newton Raphson and rearrangement methods are very easy to perform, as you simply have to add the equation and perform the functions required which only require very few steps and choosing the starting point.  However change of sign requires more knowledge of Excel, having to input the formula yourself and choosing appropriate values for x to converge to the root. However once you know how to use excel you are able to drag down and it copies the formula for each value of x without having to retype it. For the manual calculations I was able to use a calculator, although this was time consuming, setting the starting point to “ans” allowed me to quickly get all the values of x by simply pressing the equals button. Autograph also allowed me to zoom in on roots, changing axes, etc, this helped make the process of finding roots and showing failure much easier.

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

    I shall try and investigate this. 3x� x y second value x second value y gradient 4 48 4.1 50.43 24.6 4 48 4.01 48.2403 24.06 4 48 4.001 48.024003 24.006 3 27 3.1 28.83 18.6 3 27 3.01 27.1803 18.06 3 27 3.001 27.018003 18.006 2 12 2.1 13.23

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

    "T256" 0.93285210652957800 0.00000062026201424 0.25000000000000 "T512" 0.93285272679159200 0.00000015506550356 0.25000000000000 "T1024" 0.93285288185709500 0.00000003876637589 0.25000000000000 "T2048" 0.93285292062347100 0.00000000969159397 0.25000000000000 In this coursework I made used the extrapolation method. Table 1.1 is an extract of the spreadsheet I used in this coursework. The values of T up to T128 are worked out using Tn= (h/2)

  1. MEI numerical Methods

    of excel; I should be able to produce a formula on excel for the secant method thus making this method suitable for this equation. False position method: The false position method is similar to the secant method in the sense that it uses a straight line to approximate the function of the area of interest.

  2. Numerical Method (Maths Investigation)

    It allows me to design what I would like to be seen in my Coursework as it has the function of WHAT YOU SEE IS WHAT YOU GET (WYSIWYG). I can copy tables from Microsoft Excel and graphs from Graphmatica and organise them neatly in my Coursework.

  1. Numerical solutions of equations

    My negative root is x = -1.220744 (6 decimal places) Now I will find the positive root of this function. The starting value (x1) is approximately 1.5 when f(x)=0. x1 = 1.5 x2 = 1.116379 x3 = 0.862068 x4 = 0.745761 x5 = 0.725049 x6 = 0.724492 x7 = 0.724491959 x8 = 0.724491959 I can see some convergence from x6.

  2. Solving Equations Using Numerical Methods

    and this will be the lower bound. I will find the values of f(x) for this and consecutive increments of 0.1 until a change of sign is found. Here is my table of values. When looking at this table of values, no change of sign occurs.

  1. newton raphson

    0,680760 -0,000056 0.67 -0,06939 0,687 0,040909 0,6807 -0,00045 0,680770 0,000010 0.68 -0,00501 0,688 0,047526 0,6808 0,000205 0,680780 0,000075 0.69 0,060803 0,689 0,054157 0,6809 0,000857 0,680790 0,000140 0.7 0,12807 0,69 0,060803 0,681 0,00151 0,680800 0,000205 I can express this interval using the Error bounds.

  2. C3 COURSEWORK - comparing methods of solving functions

    Also, if the too roots are too close together, we might miss the other two after finding the first one. Consequently, we would be unlikely to seach further in the intervals. Despite the fact that x=g(x) method also involved some calculations, it?s not the quickest way to find all the

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