• 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

Decimal search.

Extracts from this document...

Introduction

Decimal search

The equation f (x) = 0, where f (x) = x3-x2+x-2, has only one real root, but there is no simple analytical method of finding it. Therefore, a spreadsheet had been used to solve the equations numerically using decimal search.

x

f(x)

-3

-41

-2

-16

-1

-5

0

-2

1

-1

2

4

3

19

image16.pngimage00.png

The table and graph above illustrated the first approximations to the roots of the equation  x3-x2+x-2=0.

As the curve crosses the x-axis, f (x) changes sign, so provided that f (x) is a continuous function, once you have located an interval in which f(x) changes sign, you know that that interval must contain a root.

In the table, you first take increments in x of size 1 within the interval 1<x<2, working out the value of the function x5-5x+3 for each one. You do this until you find a change of sign of f(x) between the value of x = 1 and 2 , so you should know that there is a root lying in the interval 1<x<2 on the graph.

Having narrow down the interval, you can now continue 1.3<x<1.4, you could now homing in on such root using decimal search.

x

f(x)

1

-1

1.1

-0.779

1.2

-0.512

1.3

-0.193

1.4

0.184

1.5

0.625

1.6

1.136

1.7

1.723

1.8

2.392

1.9

3.149

2

4

  • Tabulate f (x) for 1<x<2 with increments in x of 0.1, a change of sign reveals that 1.3<α<1.4.image07.png

  • Tabulate f (x) for 1.3<x<1.4 with increments in x of 0.01, a change of sign gives 1.35<α<1.36image17.png

x

f(x)

1.3

-0.193

1.31

-0.158009

1.32

-0.122432

1.33

-0.086263

1.34

-0.049496

1.35

-0.012125

1.36

0.025856

1.37

0.064453

1.38

0.103672

1.39

0.143519

1.4

0.184

image15.pngimage21.png

image00.png

  • Tabulate f (x) for 1.35<x<1.36 with increments in x of 0.001, a change of sign gives 1.353<α<1.354

image22.png

x

f(x)

1.35

-0.012125

1.351

-0.008354449

1.352

-0.004577792image15.png

1.353

-0.000795023image00.png

1.354

0.002993864

1.355

0.006788875

1.356

0.010590016

1.357

0.014397293

1.358

0.018210712

1.359

0.022030279

1.36

0.025856

...read more.

Middle

0.000341

1.3534

0.000719797

1.3535

0.001098655

1.3536

0.001477575

1.3537

0.001856555

1.3538

0.002235597

1.3539

0.0026147

1.354

0.002993864


Having established that 1.3532 <α<1.3533, since f(1.3532) <0 and f(1.3533) >0, the maximum error is 1.3533 –1.3532 = 0.0001, which is more than sufficient to give the root correct to 3

decimal places as α = 1.353

Example for problems with change of sign method

If the curve touches the x-axis without crossing it, there will be no change of sign, so change of sign methods are doomed to failure. For example,               has a repeated rational root, which can be expressed as a recurring decimal i.e      = 0.3333……image01.png

Sketch graph of f (x) =(3x-1)2(2x+4)  is drawn below, you would notice that there is no change of sign on the graph as well as the table. It seems that there is repeated root between the interval 1<x<2 , it touches the x-axis when you looking at the graph, however, the table shows that all the values of f(x) within the interval 1<x<2 are all positive, there is no change of sign. You are unable to do any further by decimal search for a root without a change of sign.

image24.png

x

f(x)

1

24

1.1

32.798

1.2

43.264

1.3

55.506

1.4

69.632

1.5

85.75

1.6

103.968

1.7

124.394

1.8

147.136

1.9

172.302

2

200


Newton-Raphson Iteration

This method is based on the iteration.

image02.png

xn+1 = xn -              , n = 0,1,2,3,…….with initial approximation x0.

The Newton – Raphson iterative formula is based on evaluating the gradient of the tangent to the curve y = f(x) at image25.png

x = x0

The gradient of the tangent at (x0, f (x0) is f ′ (x)

image03.png

f ′ (x0) =

Rearranging this,

x1 = x0 –

To solve the equation ex-7x-3, f (x) = ex-7x-3⇒ f ′ (x) = ex-7. This will give rise to the Newton Raphson iteration formula, that is,

xn+1 = xn

image26.png

image27.png

You will then be able to find the upper root α, let the initial approximation x0 = 4

image04.png

image05.png

image06.png

image08.png

                                                                                                                                  etc

The sequence converge rapidly towards the upper root to give α=3.2478 to 5 s.f.

In this case, Newton-Raphson Iteration gives an extremely rapid rate of convergence. This is the case for examples, even when the first approximation is not particularly good. For manual caluations it is almost always the most efficient method.

You may try to solve the equation using spreadsheet, it gives a

Iteration: xn+1 = x0 - (ex-7x-3)/(ex-7)

Convergence towards upper root

n

xn

0

4

1

3.504221277

2

3.286134403

3

3.248830171

4

3.247850646

5

3.247849987

6

3.247849987

7

3.247849987

8

3.247849987

9

3.247849987

10

3.247849987

...read more.

Conclusion

When the formula for a rearrangement of an equation g(x) is entered correctly with an input of initial approximation(i.e. x0 in this case). Then the approximation would converge towards the root required by simply copy and paste the formula until there is no any changing result.

Sometimes, you have to try which rearrangement is suitable for a particular roots. Plot the line y =x and the curve g (x), then find the gradient near the root, if the gradient is greater than 1 (i.e. the gradient of the line y = x ), illustrated that this is a unsuitable rearrangement of finding the root, and vice versa.

It requires the same techniques as you used in Fixed Point Iteration.

With the first approximation, it then converges gradually to the root required.

Sometimes, there is a problem of plotting a graph by using spreadsheet. For example,      +5  ,

it is impossible to plot the point where x = 0, it ends up with a continuous function which is not the expected discontinuous function.

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

    x + h - x h = h (5x4 +10x�h + 10x�h + 5xh� + h4) = 5x4 +10x�h + 10x�h + 5xh� + h4 h Every term apart from 5x4 contains an h term, leaving every term to tend to 0.

  2. Marked by a teacher

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

    3 star(s)

    graph are relate to another, as they display a recurrent cyclical pattern. Even the Figure 2(a) shows that the predict value fits the actual value well. Different way to measure permanent income will lead to different equation. If the permanent income has been treated as the lagged income, the following

  1. The open box problem

    I will now draw table sand graphs to find out the volume and the value for x. X 0.5 1 1.5 2 2.5 V 18 24 21 12 0 We can see here that the maximum volume is between 0.5 and 1.5 so I will draw another table focusing in on the data.

  2. Investigate the solution of equations, comparing the following methods, Systematic search for change of ...

    Divergence in this case is when the tangent never reaches the x-axis. Rearranging f(x)=0 in the form x=g(x) I will be using the equation y=x�-5x+1. This is equal to 0. After looking at the graph I can see that it cuts through between -3 and -2, 0 and 1, 2 and 3.

  1. The method I am going to use to solve x&amp;amp;#8722;3x-1=0 is the Change ...

    �)] = -1.380277569 x8 = x7 - [(x7 ^4+x7 �-1=0)/( 4x7 �+3x7 �)] = -1.380277569 I can see some convergence from x6. There has been no change in the x-values between x7 and x8 for this number of decimal place.

  2. Math Portfolio Type II - Applications of Sinusoidal Functions

    Explain the method you used to determine these dates. The predicted start date and end date of the growing season in Toronto would be Day 143 and Day 200. The method that is used to determine these dates is through using the TI-84 Plus graphing calculator. The sinusoidal function T(n)

  1. Solving Equations Numerically

    We can carry on using this method to get to many decimal places, but it is very slow work and can lead to severe depression and boredom. You can see on the next page that I have produced tables to find a root to 6 decimal places.

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

    = 0. x -2 -1 0 1 2 3 f (x) -29 -1 5 1 -1 11 Table1.2 From the table1.2 above, it can be seen that the roots lie in the interval of [-1, 0], [1, 2] and [2, 3] Here, I can take a=1 and b=2, whereas the mid-value is taken as (a+b)

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