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

    Have look the results from the same tests (table 3). For this equation, R^2 is 0.994231, which is very close to the unity, i.e. the equation as a representation of reality. However, the number represent the result of F-test is 8789, which is much smaller than the result showed in the second equation.

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

    On Autograph when you bring up the Newton-Raphson iteration object box and you type in 0.0, it says divergent in the results box, which is the opposite to convergent, when the numbers get more accurate and closer to the root.

  1. Numerical Method (Maths Investigation)

    Observe that the values obtained have increase to 1.51213, 1.51214 and so on. The numbers is either step up automatically or step down with number of steps determined by the two selected records. If you type in 2 in first row and 4 in second row under the same column,

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

    Solving 3x^5+5x�-1=0 using the "Rearrangement Method" of fixed point iteration The graph of the equation is shown below: I will rearrange f(x)=0 in the form of x=g(x). As my equation is 3x^5+5x�-1=0 , there are many possible rearrangements of this.

  1. The Gradient Fraction

    All the points at the 'x' scale have gradients of 2. I begin to comprehend the relationship between the equation and the gradient. The equation of this graph is 'y=2x'. The 2 in front of the 'x' shows the gradient in the straight line.

  2. Math Portfolio Type II - Applications of Sinusoidal Functions

    = 1.611 sin[0.017(65.157)] + 12.114 M(n) = 1.611 sin(1.107669) + 12.114 (use radian mode) M(n) = 1.441296853 + 12.114 M(n) = 13.55529685 M(n) = 13 hours and 33 minutes The number of hours and minutes of daylight in Toronto on day 144 is 15 hours and 3 minutes and the number of hours and minutes of daylight in Miami on day 144 is 13 hours and 33 minutes.

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