• 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
  16. 16
    16
  17. 17
    17
  18. 18
    18
  19. 19
    19
  20. 20
    20
  21. 21
    21

Change of Sign Method.

Extracts from this document...

Introduction

Pure Mathematics 2                                                                                        Katie Ruck

Solution of Equations by Numerical Methods

Change of Sign Method

If we let f(x)=3x³-0.5x²-0.5x-1, in order to solve this equation and determine its roots, it is necessary for it to be written in the form 3x³-0.5x²-0.5x-1=0.

The root of the equation f(x)=0 is indicated where y=f(x) crosses the x-axis.

Roots of the equation 3x³-0.5x²-0.5x-1=0, will be found to a three decimal place accuracy.

image00.png

        Having illustrated the equation graphically using Autograph, it is evident that the equation has only root that lies between the interval [0,1].

        Before proceeding, it is necessary to check that there is a sign change in the above interval:

f(0) = 0-0-0-1= -1

f(1) = 3-0.5-0.5-1 =1

        The method that will be used for the numerical solution of the equation is the decimal search method.

        I will first take the increments in x of size 0.1 within the interval [0,1], calculating the value of the function 3x³-0.5x²-0.5x-1 for each one, until a change of sign is found.

x

0.0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

f(x)

-1

-1.052

-1.096

-1.114

-1.088

-1

-0.832

-0.566

-0.184

0.332

        The above table illustrates that there is a sign change.  It can therefore be understood that the root lies in the interval [0.8,0.9].

        Having narrowed down the interval, I will continue with the decimal search, but now using increments of 0.01 within the interval [0.8,0.9].

x

0.80

0.81

0.82

0.83

0.84

f(x)

-0.184

-0.138727

-0.092096

-0.044089

0.005312

        This indicates that the root lies in the new interval [0.83,0.84].

I will continue the process using the increments of 0.001 within the interval [0.83,0.84].

x

0.830

0.831

0.832

0.833

0.834

0.835

f(x)

-0.044089

-0.039211927

-0.034320896

-0.029415889

-0.024496888

-0.019563875

x

0.836

0.837

0.838

0.839

f(x)

-0.014616832

-0.009655741

-0.004680584

0.000308657

The root therefore lies in the interval [0.838,0.839].  I shall now use increments of 0.0001 in the interval [0.838,0.839].

x

0.8380

0.8381

0.8382

0.8383

0.8384

0.8385

f(x)

-0.004680584

-0.00418229398

-0.0036838631

-0.00318529134

-0.00268657869

-0.00218772512

x

0.8386

0.8387

0.8388

0.8389

0.8390

f(x)

-0.00168873063

-0.00118959519

-0.00069031878

-0.00019090139

0.000308657

        Therefore the root lies between 0.8389 and 0.8390.  To three decimal places, the root = 0.839.

Error Bounds

        A change of sign method such as the one used, provides bounds within which a root lies so that the maximum possible error in a result is known.

When x = 0.8385, f(0.8385) = -0.00218772512

When x = 0.8395, f(0.8395) = 0.00280856463

The error bounds of the root 0.839 are 0.839 ± 0.0005.

        However, I am able to say that I have a more accurate solution, as I know that the root lies in the interval [0.8389,0.8390].

Failure of the Change of Sign Method

        There are a number of situations that can cause problems for change of sign methods.  For example, let y=f(x)=x³-2x²-x+2.63.  

        This curve is shown graphically below.

image01.png

        With this example I shall use the decimal search method to find a change of sign and so investigate the roots of the previous equation.  Integers will be used as the x-values.

x

-2

-1

0

1

2

3

f(x)

-11.37

0.63

2.63

0.63

0.63

8.63

        From the graph it is evident that the equation x³-2x²-x+2.63=0 has more than one root, i.e. roots that lie in the intervals [-2,-1] and [1,2].  However, the above table illustrates that the method used only detects one such root that lies in the interval       [-2,-1].

        In this case, using the method of decimal search has caused an incorrect conclusion to be reached.  This is because the curve touches the x-axis between x=1 and x=2, therefore there is no change of sign and consequently all change of sign methods are doomed to failure.

Fixed Point Iteration Method

        Let y=f(x)=x3+2x2-4x-4.58.  The graph is shown below:

image12.png

        The roots of the equation can be found where f(x)=0.  From the graph, it is evident that the roots of the equation lie in the intervals [-3,-2],[-1,0] and [1,2].

        Using a sign change search verifies that the intervals within which the roots lie are indeed [-3,-2],[-1,0] and [1,2]:

x

-4

-3

-2

-1

0

1

2

3

4

f(x)

-

-

+

+

-

-

+

+

+

...read more.

Middle

x0=-1 into the iterative formula to find the value of x.

image24.png

        A line is drawn from x0to equation 1, and then ‘followed across’ to the closest point on equation 2.  Repeating the former process will cause the values where the line meets each of the two equations to come closer together until they converge at the point where the two curves intersect, and therefore at a root.  The iterative formula calculates these values as they converge.

xr+1 = ¼( x3r+2x2r-4.58)

x1 = ¼(( -1)3+2(-1)2-4.58)

=¼(-1+2-4.58)

=¼(-3.58)

x1 = -0.895

        The following table was obtained using a Microsoft Excel spreadsheet, with the formulae shown.  The sequence works using this particular iterative formula because when the gradient of y=g(x) is less that 1, y=x is used as a ‘barrier’ and makes the sequence converge thus allowing the root to be found.

image25.png

        The values that were obtained are as follows:

image26.png

Error Bounds

        Therefore the root to the equation between the interval [-1,0] is –0.9172560 correct to seven decimal places.

        To check that this is correct, substitute

x= -0.91725595 into f(x)

= -0.00000031423

and

x= -0.91725605 into f(x)

= 0.000000200258

The above calculations illustrate that there is a change of sign.  Therefore the root is    -0.9172560±0.00000005.

image27.png

Failure of the Fixed Point Iteration Method

I will now attempt to find the root that lies in the interval [1,2] of the same equation.  As before, the Fixed Point Iteration method will be applied.

image28.png

        My chosen starting value, x0, will be 2.

The following table was constructed using a Microsoft Excel spreadsheet with the formulae displayed:

image29.png

The values that were obtained are as follows:

image02.png

...read more.

Conclusion

–0.917± 0.0005.

Fixed Point Iteration Method

        The equation has already been solved to find the root that lies in the interval    [-1,0] using this method.  The following steps show how this was done.  If y=f(x)=x3+2x2-4x-4.58 then the roots of the equation can be found where f(x)=0.  It is first necessary to look for a sign change between integer values so as to choose a suitable starting value.

x

-4

-3

-2

-1

0

1

2

3

4

f(x)

-

-

+

+

-

-

+

+

+

        I have chosen to find the root that lies in the interval [-1,0] and therefore an appropriate starting value, x0, for the approximation of the root would be –1.  When using the Fixed Point Iteration method, the previous equation must be rearranged into the form x=g(x):

x3+2x2-4x-4.58=0

x3+2x2-4.58=4x

¼( x3+2x2-4.58)=x

Therefore the iterative formula for this equation is:

xr+1 = ¼( x3r+2x2r-4.58)

        For the method to succeed, the iterative formula must converge for the chosen value of x.  To ensure that this is so, the iterative formula was differentiated, and it was checked that the x0 value gave a result that was less than one:

y=0.25x3+0.5x2-1.145

dy

   dx

substituting x=-1 gives 0.75(-1)2 + (-1)

= 0.75 – 1

= -0.25

        The fact that this value is less than 1, tells us that the chosen value of x0 will lead to a root of the equation.

        The following table was obtained using a Microsoft Excel spreadsheet, with the formulae shown.  

image21.png

        The values that were obtained are as follows:

image22.png

        Therefore the root to the equation between the interval [-1,0] is –0.917 correct to three decimal places.

When x = -0.9165,  f(-0.9165) = -0.00389006712

When x = -0.9175, f(-0.9175) = 0.00125526563

The error bounds of the root –0.917 are –0.917 ± 0.0005.

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

    Here, it is a little hard to prove this in the normal sense as (x+h) is square rooted, so I will need to rationalize the numerator. V (x+h) - Vx = V (x+h) - Vx x + h -x h V (x+h)

  2. Marked by a teacher

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

    3 star(s)

    0.206148 0.1075 1.92 0.061 0.0685 LY 0.351899 0.09701 3.63 0.001 0.2083 sigma 0.0182645 RSS 0.0166796213 R^2 0.998055 F(2,50) = 1.283e+004 [0.000]** log-likelihood 138.489 DW 0.858 no. of observations 53 no. of parameters 3 mean(LC) 12.6755 var(LC) 0.161805 LC = + 0.6322*LC_1 + 0.2061 + 0.3519*LY (SE)

  1. MEI numerical Methods

    This will be significant later. However from this we can conclude that the false position method does not have a first order convergence. Method of bisection: If we apply the first order of convergence principles to this method we get the following: Unlike the secant and false position method, the method of bisection has a constant value of K, 0.5.

  2. Change of sign method --- interval bisection method

    Continue using the same method, I am going to get a graph as following: On the graph it shows that finally when x= -3.04, the value of f(x) �0, hence the root should be -3.04 with an accuracy up to two decimal places.

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

    =(B8+C8)/2 =(D8-1)*(D8+2)*(D8-3)-1 =ABS(C8-B8)/2 =A8+1 =IF(E8>0,D8,B8) =IF(E8<0,D8,C8) =(B9+C9)/2 =(D9-1)*(D9+2)*(D9-3)-1 =ABS(C9-B9)/2 =A9+1 =IF(E9>0,D9,B9) =IF(E9<0,D9,C9) =(B10+C10)/2 =(D10-1)*(D10+2)*(D10-3)-1 =ABS(C10-B10)/2 =A10+1 =IF(E10>0,D10,B10) =IF(E10<0,D10,C10) =(B11+C11)/2 =(D11-1)*(D11+2)*(D11-3)-1 =ABS(C11-B11)/2 =A11+1 =IF(E11>0,D11,B11) =IF(E11<0,D11,C11) =(B12+C12)/2 =(D12-1)*(D12+2)*(D12-3)-1 =ABS(C12-B12)/2 =A12+1 =IF(E12>0,D12,B12) =IF(E12<0,D12,C12) =(B13+C13)/2 =(D13-1)*(D13+2)*(D13-3)-1 =ABS(C13-B13)/2 =A13+1 =IF(E13>0,D13,B13) =IF(E13<0,D13,C13) =(B14+C14)/2 =(D14-1)*(D14+2)*(D14-3)-1 =ABS(C14-B14)/2 =A14+1 =IF(E14>0,D14,B14)

  2. Experimentally calculating the wavelength of an He-Ne laser by means of diffraction gratings

    Also, the distances between the central beam and the first order fringe to the left and the first order fringe to the right seem to be equal. This is the same for the second and third order fringes. However, there seems to be slightly more space between the first order

  1. maths pure

    Using the iteration formula xn+1 = (xn3-5)/3 and the starting point x0 = 2, the results for successive iterations are as follows: X0 2 X1 1 X2 -1.333333333 X3 -2.456790123 X4 -6.609579113 X5 -97.91653871 X6 -312931.4535 This particular rearrangement of the equation x3-3x-5=0 and successive iteration has failed to converge to a root.

  2. Different methods of solving equations compared. From the Excel tables of each method, we ...

    0.77929020 -0.00000271 0.77928972 -0.00000098 0.00000173 0.77928925 0.00000075 0.77928972 -0.00000098 0.77928948 -0.00000012 0.00000087 0.77928925 0.00000075 0.77928948 -0.00000012 0.77928936 0.00000032 0.00000043 0.77928936 0.00000032 0.77928948 -0.00000012 0.77928942 0.00000010 0.00000022 0.77928942 0.00000010 0.77928948 -0.00000012 0.77928945 -0.00000001 0.00000011 0.77928942 0.00000010 0.77928945 -0.00000001 0.77928944 0.00000005 0.00000005 0.77928944 0.00000005 0.77928945 -0.00000001 0.77928945 0.00000002 0.00000003 0.77928945 0.00000002 0.77928945 -0.00000001 0.77928945 0.00000001 0.00000001 Error bound: �0.000000005 (9dp)

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