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

Numerical Methods used to solve those equations which cannot be solved analytically.

Extracts from this document...

Introduction

Mathematics Coursework

Numerical Methods can be used to solve those equations which cannot be solved analytically.

One such method is the ‘Change of Sign Method’, focusing on Decimal Search.

Decimal search is based on the principal that f(x) changes sign as a curve passes through the x axis at a root.

image00.png

Above the x axis, f(x) is positive and below the x axis f(x) is negative.

Using this principal, values of f(x) are initially found and consecutive values are found in increments of 0.1 until a change in sign is found. If further accuracy is required, consecutives values could then be found in increments of 0.01 and 0.001, until a change of accuracy is found, therefore giving one root of the equation.

Evidence

Equation: y=0.6938x3 – 0.9157x2 – 1.421x + 1.671

Graph:

image01.png

Table of results

First Interval

x

f(x)

1

0.0281

1.1

-0.07665

1.2

-0.1539

1.3

-0.1996

1.4

-0.2094

1.5

-0.1792

1.6

-0.105

1.7

0.01757

1.8

0.1926

1.9

0.4242

The table shows that the sign changes within the interval 1.6 to 1.7. This tell us that the roots lie somewhere within this interval;

Therefore the root is 1.65 ± 0.05

Second Interval

X

f(x)

1.6

-0.105

1.61

-0.09497

1.62

-0.08447

1.63

-0.07348

1.64

-0.06199

1.65

-0.05001

1.66

-0.03752

1.67

-0.02452

1.68

-0.01101

1.69

0.003019

The table of values shows that the root lies in the interval [1.68, 1.69]. Therefore the root of the equation is 1.685 ± 0.05

Third Interval

x

f(x)

1.68

-0.01101

1.681

-0.009628

1.682

-0.008243

1.683

-0.006854

1.684

-0.005459

1.685

-0.004059

1.686

-0.002654

1.687

-0.001243

1.688

0.0001725

1.689

0.001593

1.69

0.003019

...read more.

Middle

Using this method, a first approximation, x1 is made. The intercept of the tangent at (x1,y1) with the x axis is found as the next approximation, x2. This process is repeated until the desired accuracy is achieved. This method is also known as the ‘Tangent Sliding’ method.

This method can be demonstrated manually, using the equation: y= 20x3 + 15x2 – 15x – 12

F (x) = 20x3 + 15x2 – 15x – 12

F’(x) = 60x2 + 30x – 15

First estimate x1 will be 1

First iteration

F (1) = 20(13) + 15(12) -15(1) – 12 = 8

F’(1) = 60(12) + 30(1) – 15 = 75

X2 = 1 – 8/75

X2 = 0.8933

Second Iteration

F (0.8933) = 20(0.89333) + 15(0.89332) – 15(0.8933) – 12 = 0.82906074

F’(0.8933) = 60(0.89332) + 30(0.8933) – 15 = 59.68266667

X3 = 0.8933 – 0.829/59.68266667

X3 = 0.879442165

X3 = 0.8794

Third iteration

F (0.8794) = 20(0.87943) + 15(0.87942) -15(0.8794) – 12 = 0.01318348

F’(0.8794) = 60(0.87942) + 30(0.8794) – 15 =  57.78837624

X4 = 0.8794 – 0.013/57.79

X4 = 0.879214031

X4 = 0.8792

This can also be demonstrated graphically using the computer programme autograph:

y= 20x3 + 15x2 – 15x – 12

image09.png

First iteration

image02.png

 X = 0.8933

Second iteration

image02.png

X = 0.8794

Third Iteration

image10.png

X = 0.8792

The results have converged to 0.8792

To confirm that this is a root by looking for a change of sign

F (x) = 20x3 + 15x2 – 15x – 12

F (0.8791) = -6.5827 x 10-3 = negative

F (0.8792) = positive

Therefore the root is 0.87925 ±0.00005

Where the method cannot work

If the initial estimate is too far from the required root, the method can fail.

This can be shown using the equation y=20x3 + 15x2 – 8x – 19

image11.png

Take the initial estimate as x1 = -1

First Iteration

image02.png

 x = - 0.2727

Second Iteration

image02.png

x = -1.647

Third Iteration

image02.png

x=-1.13

Evidently, this method fails in this instance.

...read more.

Conclusion

Whilst the Newton Raphson method allows quick convergence, and is easy to perform on autograph, the method could be criticised, as it requires calculus, unlike decimal search. It can also be subject to failure, for example, if the initial estimate is made too far from the required root. Furthermore, to find an error bound for increased accuracy, decimal search needs to be carried out, again adding to the amount of work required. This disadvantage can also be applied to the graphical convergence method. In this method, there is also a high failure rate, as it can diverge, oscillate or converge to the wrong root, as shown in the calculations. Also, you may not be able to rearrange f(x) =0, in which case this method could not be used.

Therefore, in conclusion, whilst each method can be used to find roots, a d=sketch of the graph should be drawn, before deciding upon which method will work the best, and the level of accuracy that is desired.

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

    4.001 768.768288 768.5761 3 243 3.1 277.0563 357.492 3 243 3.01 246.256236 327.2508 3 243 3.001 243.324162 324.3241 2 48 2.1 58.3443 111.132 2 48 2.01 48.96722403 97.44721 2 48 2.001 48.09607202 96.14407 1 3 1.1 4.3923 15.972 1 3 1.01 3.12181203 12.36361 1 3 1.001 3.012018012 12.03604 x x4

  2. MEI numerical Methods

    I will use one root a as 0 and the other root b as ?/4. Bearing this in mind we can substitute the values into the formula: The formula is the same as the secant method however remember the false position method roots approximation must support the sign change argument.

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

    it had an error bound of, in other words it had a solution bound of with a maximum error. This method can also be shown graphically by computer: This method can be drawn in software 'autograph'. By using this program the equation can be enter and there is a function

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

    16 0.635070801 0.00021312 0.6351013 -0.0001467 0.63508606 3.322E-05 1.52588E-05 17 0.63508606 3.3219E-05 0.6351013 -0.0001467 0.635093689 -5.67E-05 7.62939E-06 18 0.63508606 3.3219E-05 0.6350937 -5.673E-05 0.635089874 -1.18E-05 3.8147E-06 Spreadsheet 4.2 The bisection method turns out to achieve the root 0.63509 with error bounds of Again, by applying these into the spreadsheet 4.3, the Newton-Raphson method gives us the following result: r Xr f(Xr)

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

    Check X Y -0.34500825 -0.345008213 -0.34500815 -0.34500824 X g'(x) -0.34500825 -0.271173069 -0.34500815 -0.271173032 Below shows the formulae for using in the Excel: x x=(2x�+3x�-2)/5 n g'(x) 0 =(2*A4^3+3*A4^2-2)/5 0 =(6*A4^2+6*A4)/5 =B4 =(2*A4^3+3*A4^2-2)/6 1 =(6*A4^2+6*A4)/6 =B5 =(2*A4^3+3*A4^2-2)/7 2 =(6*A4^2+6*A4)/7 =B6 =(2*A4^3+3*A4^2-2)/8 3 =(6*A4^2+6*A4)/8 =B7 =(2*A4^3+3*A4^2-2)/9 4 =(6*A4^2+6*A4)/9 =B8 =(2*A4^3+3*A4^2-2)/10 5 =(6*A4^2+6*A4)/10 =B9 =(2*A4^3+3*A4^2-2)/11 6 =(6*A4^2+6*A4)/11 =B10 =(2*A4^3+3*A4^2-2)/12 7 =(6*A4^2+6*A4)/12 =B11

  2. Best shape for gutter and further alegbra - using Excel to solve some mathematical ...

    Total interest paid in this case = �201.89 To clear the debt off in half the time the payments will have to increase to approximately �130. This is an increase of approximately 85% of the original monthly payment. To clear off the debt in double the time the payments will be lower, decreasing to around �40.

  1. The method I am going to use to solve x−3x-1=0 is the Change ...

    see that I can see that there is really no convergence to that particular root I am looking for(the one in Rearrangement B and my observation from the graph is correct.

  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.

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