9 Numerical methods — revision question pack

5 specification points · notes, questions, answers and worked methods

Checked against Edexcel 9MA0 section 9. Review basis: the qualification registry sourced from the Pearson Edexcel Level 3 Advanced GCE in Mathematics (9MA0) specification; registry verification recorded 11 July 2026.

How this checking works

9.1 · Locate roots of f(x) = 0 by considering changes of sign of f(x) in an interval on which f(x) is sufficiently well behaved; understand how change of sign methods can fail.

Explanation

  • If a continuous function has values of opposite sign at the endpoints of an interval, the intermediate value theorem guarantees at least one root inside that interval.
  • Evaluate f(a)f(a) and f(b)f(b) accurately, record their signs, and conclude that a root lies in (a,b)(a,b) only after checking that ff is continuous there.
  • For example, f(1)<0f(1)<0 and f(2)>0f(2)>0 for a continuous ff, so the graph must cross the xx-axis at least once between 11 and 22.
  • A sign change across a discontinuity need not contain a root, while a repeated root can touch the axis without changing sign; a common error is to treat the sign test as an equivalence.

Worked example

For f(x)=x3+x5f(x)=x^3+x-5, show that there is exactly one root in (1,2)(1,2).

  1. 1.Calculate f(1)=1+15=3f(1)=1+1-5=-3 and f(2)=8+25=5f(2)=8+2-5=5.
  2. 2.Continuity gives at least one root in (1,2)(1,2).
  3. 3.Also f(x)=3x2+1f'(x)=3x^2+1, which is positive for every real xx, so ff is strictly increasing and cannot cross the axis more than once.

Answer: f(1)=3f(1)=-3 and f(2)=5f(2)=5, so a root lies in (1,2)(1,2).; Since f(x)=3x2+1>0f'(x)=3x^2+1>0, ff is strictly increasing and the root is unique.

Common mistakes

  • Don't conclude that an even-multiplicity root cannot exist because endpoint values have the same sign.
  • Don't claim a root from a sign change without checking continuity, or claim uniqueness without a monotonicity argument.

Exam tip

To establish exactly one root, combine a sign change on a continuous interval with a separate uniqueness argument.

Tier 1 · Easy

  1. 1.

    Let f(x)=x2x1f(x)=x^2-x-1. Use endpoint values to establish that (1,2)(1,2) brackets a zero of ff.

    (2)

    (Total for Question 1 is 2 marks)

  2. 2.

    Let p(x)=x34x+1p(x)=x^3-4x+1. Show that p(x)=0p(x)=0 has a root in the interval (0,1)(0,1).

    (2)

    (Total for Question 2 is 2 marks)

Tier 2 · Standard

  1. 1.

    The function f(x)=(x1)2f(x)=(x-1)^2 has a root in the interval [0,2][0,2]. Explain why a search based only on a change of sign between the endpoints does not detect this root.

    (3)

    (Total for Question 1 is 3 marks)

  2. 2.

    For 1x0-1\leq x\leq0, let f(x)=x+cosxf(x)=x+\cos x. Show that f(x)=0f(x)=0 has exactly one root in (1,0)(-1,0).

    (4)

    (Total for Question 2 is 4 marks)

  3. 3.

    For f(x)=x3+xaf(x)=x^3+x-a, where aa is a constant, find the range of values of aa for which the sign-change test on the interval (1,2)(1,2) guarantees a root of f(x)=0f(x)=0 in that interval.

    (4)

    (Total for Question 3 is 4 marks)

Tier 3 · Hard

  1. 1.

    Two sign tests are proposed. For h(x)=1x2h(x)=\frac{1}{x-2}, the values h(1)h(1) and h(3)h(3) have opposite signs. For k(x)=(x2)2k(x)=(x-2)^2, the values k(1.9)k(1.9) and k(2.1)k(2.1) have the same sign. Explain why neither test gives the correct conclusion about roots in the stated intervals.

    (4)

    (Total for Question 1 is 4 marks)

  2. 2.

    Let g(x)=x33x+1g(x)=x^3-3x+1. By considering changes of sign on three disjoint intervals, prove that g(x)=0g(x)=0 has exactly three real roots.

    (5)

    (Total for Question 2 is 5 marks)

  3. 3.

    Let h(x)=ex3xh(x)=e^x-3x. By considering endpoint values and the sign of h(x)h'(x), show that h(x)=0h(x)=0 has exactly two roots in the interval (0,2)(0,2).

    (5)

    (Total for Question 3 is 5 marks)

  4. 4.

    Let f(x)=x3+x1f(x)=x^3+x-1. Establish the existence and uniqueness of a root in (0.6,0.7)(0.6,0.7). By evaluating f(0.68)f(0.68) and f(0.685)f(0.685), give this root to 22 decimal places.

    (6)

    (Total for Question 4 is 6 marks)

  5. 5.

    A function gg is continuous on [3,2][-3,2]. Some of its values are shown: g(3)=2g(-3)=2, g(2)=1g(-2)=-1, g(1)=4g(-1)=4, g(0)=0g(0)=0, g(1)=3g(1)=3 and g(2)=2g(2)=-2. Determine the least number of distinct roots that g(x)=0g(x)=0 must have in [3,2][-3,2]. Identify where they occur and explain why the information does not establish the exact number of roots.

    (5)

    (Total for Question 5 is 5 marks)

9.2 · Solve equations approximately using simple iterative methods; be able to draw associated cobweb and staircase diagrams.

Explanation

  • A fixed-point iteration rewrites an equation as x=g(x)x=g(x) and generates approximations using xn+1=g(xn)x_{n+1}=g(x_n) from a chosen starting value.
  • Calculate successive values without premature rounding; on a graph of y=g(x)y=g(x) and y=xy=x, move vertically to y=g(x)y=g(x) and horizontally to y=xy=x to display each iteration.
  • Near a fixed point α\alpha, iterations usually converge when g(α)<1|g'(\alpha)|<1: a positive gradient gives a staircase pattern and a negative gradient gives a cobweb pattern.
  • Different rearrangements of the same equation can converge at different rates or diverge; a common error is to assume that obtaining x=g(x)x=g(x) automatically produces a useful recurrence.
A staircase diagram alternates between y=g(x)y=g(x) and y=xy=x to display successive fixed-point iterates.

Worked example

Use xn+1=7xnx_{n+1}=\sqrt{7-x_n} with x0=2x_0=2 to calculate x1x_1, x2x_2 and x3x_3 to 44 decimal places. State the equation satisfied by any limiting value.

  1. 1.Repeated substitution gives x1=5=2.236067x_1=\sqrt5=2.236067\ldots, x2=72.236067=2.182643x_2=\sqrt{7-2.236067\ldots}=2.182643\ldots and x3=72.182643=2.194847x_3=\sqrt{7-2.182643\ldots}=2.194847\ldots.
  2. 2.If the sequence tends to LL, then L=7LL=\sqrt{7-L}, so L2+L7=0L^2+L-7=0.

Answer: x1=2.2361x_1=2.2361, x2=2.1826x_2=2.1826 and x3=2.1948x_3=2.1948.; A limiting value satisfies x2+x7=0x^2+x-7=0.

Common mistakes

  • Don't draw horizontal moves to y=g(x)y=g(x) and vertical moves to y=xy=x, reversing the cobweb construction.
  • Don't calculate iterates and fail to state the fixed-point equation or assess convergence from the iteration graph.

Exam tip

Show each substituted iterate to the requested accuracy, then set the limiting value equal to the iteration function.

Tier 1 · Easy

  1. 1.

    The recurrence xn+1=12(xn+6xn)x_{n+1}=\frac{1}{2}\left(x_n+\frac{6}{x_n}\right) is used with x0=2x_0=2. Find x1x_1 and x2x_2.

    (2)

    (Total for Question 1 is 2 marks)

  2. 2.

    The recurrence xn+1=11xnx_{n+1}=\sqrt{11-x_n} is used with x0=3x_0=3. Find x1x_1 and x2x_2, giving each answer to 44 decimal places.

    (2)

    (Total for Question 2 is 2 marks)

Tier 2 · Standard

  1. 1.

    On a pair of axes, draw y=xy=x and y=2+14xy=2+\dfrac14x. Starting at x0=0x_0=0 on the xx-axis, draw the first three steps of the fixed-point construction. State whether the construction is a staircase or a cobweb, and find its limiting value.

    (4)

    (Total for Question 1 is 4 marks)

  2. 2.

    Starting with x0=0x_0=0, use xn+1=1+6xn+4x_{n+1}=1+\dfrac{6}{x_n+4} to find x1x_1, x2x_2 and x3x_3. Find the positive limiting value and state whether the associated construction is a staircase or a cobweb.

    (4)

    (Total for Question 2 is 4 marks)

  3. 3.

    The recurrence xn+1=ln(xn+6)x_{n+1}=\ln(x_n+6) is used with x0=1x_0=1 and is known to converge to a fixed point α\alpha. Write down an equation of the form f(x)=0f(x)=0 whose root is α\alpha. Calculate x1x_1, x2x_2 and x3x_3, giving each value to 44 decimal places. The values satisfy x0<x1<x2<x3<αx_0<x_1<x_2<x_3<\alpha. Without using derivatives, explain why the associated fixed-point path is a staircase rather than a cobweb.

    (4)

    (Total for Question 3 is 4 marks)

Tier 3 · Hard

  1. 1.

    For xn+1=30.4xnx_{n+1}=3-0.4x_n with x0=1x_0=1, find x1x_1 to x4x_4. Determine the limiting value and state whether the graphical construction is a cobweb or a staircase, giving a reason.

    (5)

    (Total for Question 1 is 5 marks)

  2. 2.

    The positive root α\alpha of x2+x4=0x^2+x-4=0 lies in (1.5,1.6)(1.5,1.6). Two proposed iterations are xn+1=4xn2x_{n+1}=4-x_n^2 and xn+1=4xnx_{n+1}=\sqrt{4-x_n}. Without carrying out a long sequence of iterations, determine which recurrence is locally convergent to α\alpha. State whether its graphical construction is a staircase or a cobweb.

    (5)

    (Total for Question 2 is 5 marks)

  3. 3.

    The recurrence xn+1=3xn+1x_{n+1}=\dfrac{3}{x_n+1} is used with x0=1x_0=1. Use exact arithmetic to find x1x_1, x2x_2, x3x_3 and x4x_4. Find the positive fixed point. Use the pattern in these four iterates to identify the type of fixed-point diagram and state whether convergence is indicated.

    (5)

    (Total for Question 3 is 5 marks)

  4. 4.

    A fixed-point iteration has the linear form xn+1=a+bxnx_{n+1}=a+bx_n. On its fixed-point diagram, the first vertical step from x0=1x_0=1 reaches (1,1.8)(1,1.8), and the next vertical step reaches (1.8,2.04)(1.8,2.04). Find aa and bb. Calculate x3x_3, find the limiting value, and classify the diagram as a convergent staircase or a convergent cobweb.

    (5)

    (Total for Question 4 is 5 marks)

  5. 5.

    The iteration xn+1=2exnx_{n+1}=2-e^{-x_n} starts from x0=1x_0=1. Use unrounded values in your working. Calculate x1x_1 to x6x_6, giving each value to 66 decimal places, and find the least value of nn for which xn+1xn<0.0005|x_{n+1}-x_n|<0.0005. Hence estimate the fixed point to 44 decimal places and state the type of fixed-point diagram.

    (6)

    (Total for Question 5 is 6 marks)

9.3 · Solve equations using the Newton-Raphson method and other recurrence relations of the form xₙ₊₁ = g(xₙ); understand how such methods can fail.

Explanation

  • Newton-Raphson replaces the curve locally by its tangent, giving xn+1=xnf(xn)f(xn)x_{n+1}=x_n-\frac{f(x_n)}{f'(x_n)} for a root of f(x)=0f(x)=0.
  • Differentiate first, substitute the current approximation into both ff and ff', and retain extra calculator digits until the requested final accuracy.
  • For instance, applying Newton-Raphson to f(x)=x2af(x)=x^2-a gives xn+1=12(xn+axn)x_{n+1}=\frac12\left(x_n+\frac{a}{x_n}\right).
  • The method can fail when f(xn)=0f'(x_n)=0, when a tangent sends the next value away from the desired root, or when the values enter a cycle; a common error is to continue without checking the iterates.
  • Set f(x)=x3+x4f(x)=x^3+x-4, so f(x)=3x2+1f'(x)=3x^2+1.
Newton–Raphson uses the tangent at xnx_n; its intercept with the xx-axis is the next iterate xn+1x_{n+1}.

Worked example

Derive a Newton-Raphson recurrence for x3+x4=0x^3+x-4=0. Starting with x0=1.3x_0=1.3, find x1x_1 and x2x_2 to 55 decimal places.

  1. 1.Set f(x)=x3+x4f(x)=x^3+x-4, so f(x)=3x2+1f'(x)=3x^2+1.
  2. 2.Substitution in the Newton-Raphson formula gives xn+1=xnxn3+xn43xn2+1x_{n+1}=x_n-\frac{x_n^3+x_n-4}{3x_n^2+1}.
  3. 3.Using x0=1.3x_0=1.3 gives x1=1.382866x_1=1.382866\ldots, and substituting this unrounded value gives x2=1.378806x_2=1.378806\ldots.

Answer: xn+1=xnxn3+xn43xn2+1x_{n+1}=x_n-\frac{x_n^3+x_n-4}{3x_n^2+1}.; x1=1.38287x_1=1.38287 and x2=1.37881x_2=1.37881.

Common mistakes

  • Don't choose a starting value where f(x0)=0f'(x_0)=0, making the first Newton step undefined.
  • Don't use the Newton formula with the derivative evaluated at the new iterate rather than the current one.

Exam tip

Write the Newton recurrence explicitly, retain guard digits between iterations and round only the reported values.

Tier 1 · Easy

  1. 1.

    Use one Newton-Raphson step on f(x)=x27f(x)=x^2-7 starting from x0=2.5x_0=2.5. Find x1x_1.

    (2)

    (Total for Question 1 is 2 marks)

  2. 2.

    Use one Newton-Raphson step on f(x)=ex3f(x)=e^x-3, starting with x0=1x_0=1. Give x1x_1 to 44 decimal places.

    (2)

    (Total for Question 2 is 2 marks)

Tier 2 · Standard

  1. 1.

    Newton-Raphson is applied to f(x)=x32x+2f(x)=x^3-2x+2 starting from x0=0x_0=0. Calculate x1x_1, x2x_2 and x3x_3, and explain why the method fails to converge to a root from this starting value.

    (4)

    (Total for Question 1 is 4 marks)

  2. 2.

    Derive a Newton-Raphson recurrence for x3+x26=0x^3+x^2-6=0. Starting with x0=1.5x_0=1.5, find x1x_1 and x2x_2, giving each answer to 44 decimal places.

    (4)

    (Total for Question 2 is 4 marks)

  3. 3.

    Newton-Raphson is used to solve lnx+x3=0\ln x+x-3=0. Derive the recurrence and, starting with x0=2x_0=2, find x1x_1 and x2x_2, giving each value to 55 decimal places. Use unrounded values in your working.

    (4)

    (Total for Question 3 is 4 marks)

Tier 3 · Hard

  1. 1.

    Newton-Raphson is applied to f(x)=x2af(x)=x^2-a, where a>0a>0, starting from x0=4x_0=4. The first iterate is x1=134x_1=\dfrac{13}{4}. Find aa. Then calculate x2x_2 exactly and x3x_3 to six decimal places, and hence give the positive root of f(x)=0f(x)=0 correct to four decimal places.

    (5)

    (Total for Question 1 is 5 marks)

  2. 2.

    Newton-Raphson is to be used to solve x33x+3=0x^3-3x+3=0. Explain why x0=1x_0=1 cannot be used. Then start with x0=2x_0=-2 and calculate x1x_1 and x2x_2, giving each value to 55 decimal places.

    (5)

    (Total for Question 2 is 5 marks)

  3. 3.

    Consider applying Newton-Raphson to f(x)=x3f(x)=\sqrt[3]{x}, whose only real root is 00. For xn0x_n\ne0, derive the recurrence in its simplest form. Starting with x0=18x_0=\dfrac18, find x1x_1, x2x_2 and x3x_3, and explain why the method does not converge to the root.

    (4)

    (Total for Question 3 is 4 marks)

  4. 4.

    Apply Newton-Raphson to cosx=x\cos x=x by first deriving a recurrence. From x0=0.7x_0=0.7, obtain the first three iterates and report each to 44 decimal places, retaining full calculator precision throughout. Verify by a change of sign between 0.739050.73905 and 0.739150.73915 that the root is 0.73910.7391 to 44 decimal places.

    (6)

    (Total for Question 4 is 6 marks)

  5. 5.

    Newton-Raphson is to be applied to 1x15=0\dfrac1x-\dfrac15=0. Derive the recurrence in polynomial form. Starting with x0=1x_0=1, find x1x_1 and x2x_2 exactly. Then show that the starting value x0=10x_0=10 makes the method fail after one step, and explain the failure.

    (5)

    (Total for Question 5 is 5 marks)

9.4 · Understand and use numerical integration of functions, including the use of the trapezium rule and estimating the approximate area under a curve and limits that it must lie between.

Explanation

  • The trapezium rule approximates an integral by replacing the curve over each strip with a straight chord and adding the resulting trapezium areas.
  • For ordinates y0,,yny_0,\ldots,y_n at equal spacing hh, use h2[y0+yn+2(y1++yn1)]\frac{h}{2}[y_0+y_n+2(y_1+\cdots+y_{n-1})] and include each interior ordinate twice.
  • A convex curve with f(x)>0f''(x)>0 lies below its chords, so the trapezium rule overestimates its integral; more, narrower strips usually improve the estimate.
  • Numerical area is not automatically the signed integral: a common error is to leave contributions below the xx-axis negative when the question asks for total area.
The trapezium rule joins consecutive ordinates with straight segments; convexity indicates whether this polygonal boundary lies above or below the curve.

Worked example

Use four equal strips and the trapezium rule to estimate 02(x2+1)dx\int_0^2(x^2+1)\,dx. State whether the estimate is an overestimate or an underestimate, with a reason.

  1. 1.Here h=0.5h=0.5 and the ordinates at x=0,0.5,1,1.5,2x=0,0.5,1,1.5,2 are 1,1.25,2,3.25,51,1.25,2,3.25,5.
  2. 2.Thus the estimate is 0.52[1+5+2(1.25+2+3.25)]=0.25(19)=4.75\frac{0.5}{2}[1+5+2(1.25+2+3.25)]=0.25(19)=4.75.
  3. 3.Since f(x)=2>0f''(x)=2>0, each chord lies above the curve, so the result is an overestimate.

Answer: 4.754.75; It is an overestimate because f(x)=2>0f''(x)=2>0, so the curve is convex.

Common mistakes

  • Don't apply weight two to the endpoints and weight one to the interior ordinates, reversing the trapezium-rule coefficients.
  • Don't use the number of ordinates as the number of strips, giving the wrong strip width and endpoint weights.

Exam tip

For the trapezium rule, calculate strip width from the interval and use half weight only on the two endpoints.

Tier 1 · Easy

  1. 1.

    The values of a function are y=1,3,9y=1,3,9 at x=0,1,2x=0,1,2 respectively. Use two strips and the trapezium rule to estimate 02ydx\int_0^2 y\,dx.

    (2)

    (Total for Question 1 is 2 marks)

  2. 2.

    At x=0,1,2x=0,1,2, the corresponding ordinates of a curve are 2,k,62,k,6. The trapezium rule with two strips gives an estimate of 1010 for the integral from 00 to 22. Find kk.

    (3)

    (Total for Question 2 is 3 marks)

Tier 2 · Standard

  1. 1.

    The values of a function at x=0,1,2,3,4x=0,1,2,3,4 are respectively 2,2.5,4,6.5,102,2.5,4,6.5,10. Use the trapezium rule to estimate 04f(x)dx\displaystyle\int_0^4 f(x)\,\mathrm dx.

    (3)

    (Total for Question 1 is 3 marks)

  2. 2.

    The table gives values of y=x+1y=\sqrt{x+1}: for x=0,0.5,1,1.5,2x=0,0.5,1,1.5,2, the corresponding values of yy are 1,1.225,1.414,1.581,1.7321,1.225,1.414,1.581,1.732. Obtain a trapezium-rule estimate for 02x+1dx\int_0^2\sqrt{x+1}\,\mathrm dx. Report the result to 33 decimal places and state, with a reason, whether it is an overestimate or an underestimate.

    (4)

    (Total for Question 2 is 4 marks)

  3. 3.

    The corresponding values of xx and yy are x=0,0.5,1.5,3x=0,0.5,1.5,3 and y=2,4,3,7y=2,4,3,7. Using a straight-line trapezium between each pair of adjacent data points, estimate 03ydx\int_0^3 y\,\mathrm dx.

    (4)

    (Total for Question 3 is 4 marks)

Tier 3 · Hard

  1. 1.

    For f(x)=1xf(x)=\frac1x on 1x21\leq x\leq2, use four equal strips to find a trapezium-rule estimate. By using left- and right-endpoint rectangles, also give lower and upper bounds for 121xdx\int_1^2\frac1x\,dx. Give decimals to 44 decimal places.

    (6)

    (Total for Question 1 is 6 marks)

  2. 2.

    A continuous curve passes through the points whose coordinates are given by x=0,1,2,3x=0,1,2,3 and y=2,1,1,4y=-2,-1,1,4 respectively. Use the trapezium rule to estimate the signed integral from 00 to 33. Then, using the straight-line trapezium approximation between each pair of adjacent points, estimate the total area between the curve and the xx-axis.

    (5)

    (Total for Question 2 is 5 marks)

  3. 3.

    At x=0,1,2,3,4x=0,1,2,3,4, the ordinates of y=f(x)y=f(x) are 8,7,9,12,148,7,9,12,14 and the ordinates of y=g(x)y=g(x) are 2,3,4,6,92,3,4,6,9. Throughout the interval, f(x)g(x)f(x)\geq g(x). Use the trapezium rule with four strips to estimate the area between the two curves. Confirm the result by separately estimating 04f(x)dx\int_0^4 f(x)\,\mathrm dx and 04g(x)dx\int_0^4 g(x)\,\mathrm dx.

    (5)

    (Total for Question 3 is 5 marks)

  4. 4.

    A solid of revolution has radius rr at distance xx along its axis. For successive unit distances x=0x=0 through x=4x=4, the radii are 1,2,3,2,11,2,3,2,1. Using V=π04r2dxV=\pi\int_0^4 r^2\,\mathrm dx, obtain a four-strip trapezium estimate of the volume. Give it exactly in terms of π\pi and show the transformed ordinates used.

    (5)

    (Total for Question 4 is 5 marks)

  5. 5.

    Five function values are equally spaced from x=0x=0 to x=2x=2 in steps of 0.50.5. The ordinates at the first, third and fifth points are 2,5,82,5,8. Using only these three ordinates, verify that the two-strip trapezium estimate is 1010. The four-strip estimate using all five ordinates is 9.59.5, and the fourth ordinate is twice the second. Find the two missing ordinates.

    (5)

    (Total for Question 5 is 5 marks)

9.5 · Use numerical methods to solve problems in context.

Explanation

  • A contextual numerical model translates a physical condition, such as a target value or accumulated quantity, into a root, recurrence or numerical integral.
  • Define the function and units before applying the method, choose an interval or starting value that is meaningful in the context, and state the requested accuracy.
  • For a threshold M(t)=CM(t)=C, solve M(t)C=0M(t)-C=0 numerically; for a varying rate q(t)q(t), approximate the accumulated change with q(t)dt\int q(t)\,dt.
  • A calculator value is not a complete contextual answer: a common error is to omit units, ignore the model's domain or round in a way that makes a safety decision invalid.

Worked example

A culture is modelled by P(t)=200e0.3tP(t)=200e^{0.3t}, where tt is in hours. Use Newton-Raphson on f(t)=200e0.3t1000f(t)=200e^{0.3t}-1000 with t0=5t_0=5 to find t1t_1 and t2t_2. Hence estimate when the culture reaches 10001000 cells.

  1. 1.Differentiate to get f(t)=60e0.3tf'(t)=60e^{0.3t}.
  2. 2.Use tn+1=tn200e0.3tn100060e0.3tnt_{n+1}=t_n-\frac{200e^{0.3t_n}-1000}{60e^{0.3t_n}}.
  3. 3.Substitution of t0=5t_0=5 gives t1=5.385502t_1=5.385502\ldots, and the next substitution gives t2=5.364857t_2=5.364857\ldots.
  4. 4.The positive time is therefore about 5.3655.365 hours.

Answer: t1=5.38550t_1=5.38550 and t2=5.36486t_2=5.36486.; The culture reaches 10001000 cells after approximately 5.3655.365 hours.

Common mistakes

  • Don't select a root that satisfies the equation but lies outside the interval allowed by the context.
  • Don't carry out a numerical iteration but give no contextual time, units or check that the root is meaningful.

Exam tip

After numerical iteration, translate the root back into the context and report a sensible accuracy with units.

Tier 1 · Easy

  1. 1.

    A greenhouse temperature is modelled by T(t)=18+4tt2T(t)=18+4t-t^2 degrees Celsius, where tt is hours after noon. Show that the temperature returns to 20C20^\circ\mathrm{C} at a time between 33 and 44 hours after noon.

    (2)

    (Total for Question 1 is 2 marks)

  2. 2.

    A sensor output is modelled by S(t)=t3+tS(t)=t^3+t volts for t0t\geq0, where tt is measured in seconds. Show that the output reaches 2020 volts at a time between 22 and 33 seconds.

    (2)

    (Total for Question 2 is 2 marks)

Tier 2 · Standard

  1. 1.

    A particle's height is modelled by h(t)=20t4.9t2h(t)=20t-4.9t^2 metres. Show that the positive time when the particle returns to the ground lies in (4.0,4.1)(4.0,4.1) seconds and, by considering the sign of h(4.05)h(4.05), show that it lies in (4.05,4.10)(4.05,4.10) seconds.

    (4)

    (Total for Question 1 is 4 marks)

  2. 2.

    A curing model predicts that the completion time tt hours satisfies t=2.4+0.35ln(t+1)t=2.4+0.35\ln(t+1). Starting with t0=3t_0=3, use tn+1=2.4+0.35ln(tn+1)t_{n+1}=2.4+0.35\ln(t_n+1) to find t1t_1, t2t_2 and t3t_3, giving each iterate to 66 decimal places. Hence give the predicted completion time to 33 decimal places.

    (4)

    (Total for Question 2 is 4 marks)

  3. 3.

    A radio tag's reliable range RR metres is modelled by R=40+4cos(R/20)R=40+4\cos(R/20). Starting with R0=45R_0=45, use Rn+1=40+4cos(Rn/20)R_{n+1}=40+4\cos(R_n/20) to find R1R_1 to R5R_5, giving each iterate to 44 decimal places. Use unrounded values in your working. Hence estimate the reliable range to the nearest 0.10.1 m.

    (4)

    (Total for Question 3 is 4 marks)

Tier 3 · Hard

  1. 1.

    Water enters a tank at rates 10,14,18,21,19,15,1110,14,18,21,19,15,11 litres per minute at times 0,5,10,15,20,25,300,5,10,15,20,25,30 minutes. The tank initially contains 350350 litres and has capacity 800800 litres. Use the trapezium rule to estimate whether the tank overflows during the 3030 minutes, and estimate any excess at the end.

    (5)

    (Total for Question 1 is 5 marks)

  2. 2.

    The volume of liquid in a vessel is modelled by V(h)=h3+4h2V(h)=h^3+4h^2 litres when the depth is hh metres, for h>0h>0. Show that the depth at which V=70V=70 litres lies between 33 m and 44 m. Starting with h0=3.2h_0=3.2, apply Newton-Raphson twice, giving each iterate to 66 decimal places, and hence give the required depth to 33 significant figures.

    (6)

    (Total for Question 2 is 6 marks)

  3. 3.

    The distance DD metres travelled by an automated trolley after tt seconds is modelled by D(t)=5t+8sin(0.5t)D(t)=5t+8\sin(0.5t) for 0t80\leq t\leq8. Show that there is a unique time in (5,6)(5,6) at which D=30D=30. Derive a Newton-Raphson recurrence and, starting with t0=5.5t_0=5.5, find t1t_1 and t2t_2 to 55 decimal places. Use unrounded values in your working. Hence give the time to 33 significant figures.

    (6)

    (Total for Question 3 is 6 marks)

  4. 4.

    A drying-process output is modelled by S(t)=t+2etS(t)=t+2e^{-t}, where tt is the time in hours. The process is complete when S=3.5S=3.5. Show that the completion time lies in (3,4)(3,4) and is unique in this interval. Starting with this interval, repeatedly test the midpoint of the current bracket and retain the half containing the sign change; carry out five such steps. State the resulting interval, use its midpoint to estimate the completion time, and give the maximum possible error in this estimate.

    (7)

    (Total for Question 4 is 7 marks)

  5. 5.

    An optical instrument has a calibration factor c>1c>1 satisfying c+lnc=4c+\ln c=4. Show that there is exactly one possible factor between 2.92.9 and 33. By testing suitable decimal values and using changes of sign, first locate cc in an interval of width 0.010.01 and then refine this to an interval of width 0.0010.001. Hence give the calibration factor to 33 significant figures.

    (6)

    (Total for Question 5 is 6 marks)

Answer key

Answers begin on a new printed page so the question pack can be completed without the solutions alongside it.

9.1 · Locate roots of f(x) = 0 by considering changes of sign of f(x) in an interval on which f(x) is sufficiently well behaved; understand how change of sign methods can fail.

Tier 1 · Easy

Mark scheme for 9.1 Tier 1 · Easy
QuestionSchemeMarks
1
  • f(1)=1f(1)=-1 and f(2)=1f(2)=1, so continuity and the change of sign give a root in (1,2)(1,2).
2
(2 marks)2
Notes
Substitute the endpoints: f(1)=111=1f(1)=1-1-1=-1 and f(2)=421=1f(2)=4-2-1=1. A polynomial is continuous, and the endpoint values have opposite signs, so the intermediate value theorem guarantees at least one root between 11 and 22.
2
  • p(0)=1p(0)=1 and p(1)=2p(1)=-2, so continuity and the change of sign give a root in (0,1)(0,1).
2
(2 marks)2
Notes
The polynomial pp is continuous. Also p(0)=1p(0)=1 and p(1)=14+1=2p(1)=1-4+1=-2. Since these values have opposite signs, the intermediate value theorem guarantees at least one root in (0,1)(0,1).

Tier 2 · Standard

Mark scheme for 9.1 Tier 2 · Standard
QuestionSchemeMarks
1
  • f(0)=f(2)=1f(0)=f(2)=1, so there is no endpoint sign change, although f(1)=0f(1)=0.
  • The graph touches the axis at the repeated (even-multiplicity) root instead of crossing it; any equivalent explanation that the sign stays positive on both sides is valid.
3
(3 marks)3
Notes
At the endpoints, f(0)=(1)2=1f(0)=(-1)^2=1 and f(2)=12=1f(2)=1^2=1, so their product is positive and the usual sign-change test gives no evidence of a root. Nevertheless, f(1)=0f(1)=0. Because the root is repeated, the non-negative graph touches the xx-axis at x=1x=1 and turns back without changing sign.
2
  • f(1)=cos11<0f(-1)=\cos1-1<0 and f(0)=1>0f(0)=1>0, so a root lies in (1,0)(-1,0); since f(x)=1sinx>0f'(x)=1-\sin x>0 on this interval, the root is unique.
4
(4 marks)4
Notes
ff is continuous on [1,0][-1,0]. Its endpoint values are f(1)=cos11<0f(-1)=\cos1-1<0 and f(0)=1>0f(0)=1>0, so there is at least one root in (1,0)(-1,0). On this interval sinx0\sin x\leq0, so f(x)=1sinx1>0f'(x)=1-\sin x\geq1>0. Thus ff is strictly increasing and can have at most one root, proving uniqueness.
3
  • 2<a<102<a<10
4
(4 marks)4
Notes
The polynomial is continuous. At the endpoints, f(1)=2af(1)=2-a and f(2)=10af(2)=10-a. The sign-change test guarantees a root when (2a)(10a)<0(2-a)(10-a)<0, which holds for 2<a<102<a<10. The endpoints are excluded because a=2a=2 or a=10a=10 places the root at an endpoint rather than inside (1,2)(1,2).

Tier 3 · Hard

Mark scheme for 9.1 Tier 3 · Hard
QuestionSchemeMarks
1
  • hh changes sign across x=2x=2 but has no root because it is discontinuous there.
  • kk has the repeated root x=2x=2 but does not change sign across it.
4
(4 marks)4
Notes
For hh, h(1)=1h(1)=-1 and h(3)=1h(3)=1, but hh is undefined at x=2x=2, so continuity fails and the sign change occurs across a vertical asymptote rather than a root. For kk, k(2)=0k(2)=0, but squaring makes k(x)0k(x)\geq0 on both sides of 22, so the graph touches the axis at a repeated root without crossing it.
2
  • There is one root in each of (2,1)(-2,-1), (0,1)(0,1) and (1,2)(1,2); these are exactly the three real roots.
5
(5 marks)5
Notes
gg is continuous. The values g(2)=1g(-2)=-1 and g(1)=3g(-1)=3 have opposite signs, as do g(0)=1g(0)=1 and g(1)=1g(1)=-1, and g(1)=1g(1)=-1 and g(2)=3g(2)=3. Hence there is a root in each of the three disjoint intervals (2,1)(-2,-1), (0,1)(0,1) and (1,2)(1,2). A non-zero cubic polynomial has at most three real roots, so these three roots are all the real roots of gg.
3
  • There is one root in (0,1)(0,1) and one root in (1,2)(1,2), and these are the only roots in (0,2)(0,2).
5
(5 marks)5
Notes
hh is continuous, with h(0)=1>0h(0)=1>0, h(1)=e3<0h(1)=e-3<0 and h(2)=e26>0h(2)=e^2-6>0. Hence sign changes give a root in each of (0,1)(0,1) and (1,2)(1,2). Also h(x)=ex3h'(x)=e^x-3, which is negative for x<ln3x<\ln3 and positive for x>ln3x>\ln3. Thus hh decreases and then increases, so it has at most one root on each side of its single turning point; the turning point itself is not a root, since h(ln3)=33ln3<0h(\ln3)=3-3\ln3<0. The two bracketed roots are therefore the only roots in (0,2)(0,2).
4
  • f(0.6)=0.184f(0.6)=-0.184 and f(0.7)=0.043f(0.7)=0.043, so a root lies in (0.6,0.7)(0.6,0.7).
  • f(x)=3x2+1>0f'(x)=3x^2+1>0, so the root is unique.
  • f(0.68)=0.005568f(0.68)=-0.005568 and f(0.685)=0.006419125f(0.685)=0.006419125, so the root is 0.680.68 to 22 decimal places.
6
(6 marks)6
Notes
The polynomial ff is continuous. Since f(0.6)=0.63+0.61=0.184f(0.6)=0.6^3+0.6-1=-0.184 and f(0.7)=0.73+0.71=0.043f(0.7)=0.7^3+0.7-1=0.043, a change of sign gives a root in (0.6,0.7)(0.6,0.7). Also f(x)=3x2+1>0f'(x)=3x^2+1>0 for every real xx, so ff is strictly increasing and the root is unique. Now f(0.68)=0.683+0.681=0.005568f(0.68)=0.68^3+0.68-1=-0.005568, while f(0.685)=0.6853+0.6851=0.006419125f(0.685)=0.685^3+0.685-1=0.006419125. Hence the root lies in (0.68,0.685)(0.68,0.685) and is 0.680.68 to 22 decimal places.
5
  • At least four distinct roots: one in each of (3,2)(-3,-2), (2,1)(-2,-1) and (1,2)(1,2), together with the root x=0x=0.
  • The data give no uniqueness or monotonicity information, so additional roots may occur between the sampled points.
5
(5 marks)5
Notes
Continuity and the sign changes between 3-3 and 2-2, between 2-2 and 1-1, and between 11 and 22 guarantee a root in each of those three disjoint open intervals. The stated value g(0)=0g(0)=0 supplies a fourth root, distinct from the three bracketed roots. Thus four is the least number forced by the data. A continuous graph could cross the axis additional times inside any sampled interval, so these values alone cannot prove that there are exactly four roots.

9.2 · Solve equations approximately using simple iterative methods; be able to draw associated cobweb and staircase diagrams.

Tier 1 · Easy

Mark scheme for 9.2 Tier 1 · Easy
QuestionSchemeMarks
1
  • x1=2.5x_1=2.5 and x2=2.45x_2=2.45
2
(2 marks)2
Notes
Substitute x0=2x_0=2 to obtain x1=12(2+6/2)=2.5x_1=\frac12(2+6/2)=2.5. Then x2=12(2.5+6/2.5)=12(4.9)=2.45x_2=\frac12(2.5+6/2.5)=\frac12(4.9)=2.45.
2
  • x1=2.8284x_1=2.8284 and x2=2.8586x_2=2.8586
2
(2 marks)2
Notes
x1=113=8=2.828427x_1=\sqrt{11-3}=\sqrt8=2.828427\ldots. Using this unrounded value, x2=112.828427=2.858596x_2=\sqrt{11-2.828427\ldots}=2.858596\ldots. Therefore x1=2.8284x_1=2.8284 and x2=2.8586x_2=2.8586 to 44 decimal places.

Tier 2 · Standard

Mark scheme for 9.2 Tier 2 · Standard
QuestionSchemeMarks
1
  • The construction moves vertically to y=2+x/4y=2+x/4 and horizontally to y=xy=x, giving x1=2x_1=2, x2=2.5x_2=2.5 and x3=2.625x_3=2.625.
  • It is a convergent staircase.
  • The limiting value is 83\dfrac83.
4
(4 marks)4
Notes
The graphical steps follow xn+1=2+xn/4x_{n+1}=2+x_n/4: from x0=0x_0=0, they give x1=2x_1=2, x2=5/2x_2=5/2 and x3=21/8x_3=21/8. The gradient g(x)=1/4g'(x)=1/4 is positive, so the path is a staircase; since 0<1/4<10<1/4<1, it converges. At the fixed point L=2+L/4L=2+L/4, so L=8/3L=8/3.
2
  • x1=52x_1=\dfrac52, x2=2513x_2=\dfrac{25}{13} and x3=15577x_3=\dfrac{155}{77}.
  • The positive limit is 22, and the construction is a convergent cobweb.
4
(4 marks)4
Notes
Successive substitution gives x1=1+6/4=2.5x_1=1+6/4=2.5, x2=1+6/6.5=25/13=1.923076x_2=1+6/6.5=25/13=1.923076\ldots and x3=1+6/(25/13+4)=155/77=2.012987x_3=1+6/(25/13+4)=155/77=2.012987\ldots. A limit LL satisfies (L1)(L+4)=6(L-1)(L+4)=6, so L2+3L10=0L^2+3L-10=0 and the positive solution is L=2L=2. Here g(x)=6/(x+4)2g'(x)=-6/(x+4)^2, so g(2)=1/6g'(2)=-1/6: the negative gradient gives a cobweb and its magnitude below 11 gives convergence.
3
  • exx6=0e^x-x-6=0
  • x1=1.9459x_1=1.9459, x2=2.0727x_2=2.0727 and x3=2.0885x_3=2.0885.
  • The path is a staircase because the iterates remain on the same side of α\alpha and increase towards it instead of alternating across it.
4
(4 marks)4
Notes
At a fixed point, α=ln(α+6)\alpha=\ln(\alpha+6), hence eα=α+6e^\alpha=\alpha+6 and exx6=0e^x-x-6=0. Using unrounded values gives x1=ln7=1.945910149x_1=\ln7=1.945910149\ldots, x2=ln(6+x1)=2.072657349x_2=\ln(6+x_1)=2.072657349\ldots and x3=ln(6+x2)=2.088482715x_3=\ln(6+x_2)=2.088482715\ldots, so the requested values are 1.94591.9459, 2.07272.0727 and 2.08852.0885. Since the iterates increase towards α\alpha while staying below it, the graphical steps approach the intersection from one side and form a staircase, not a cobweb.

Tier 3 · Hard

Mark scheme for 9.2 Tier 3 · Hard
QuestionSchemeMarks
1
  • x1=2.6x_1=2.6, x2=1.96x_2=1.96, x3=2.216x_3=2.216 and x4=2.1136x_4=2.1136.
  • The limit is 157\frac{15}{7}, and the construction is a convergent cobweb because g(x)=0.4g'(x)=-0.4.
5
(5 marks)5
Notes
Apply the recurrence successively: x1=30.4=2.6x_1=3-0.4=2.6, x2=30.4(2.6)=1.96x_2=3-0.4(2.6)=1.96, x3=30.4(1.96)=2.216x_3=3-0.4(1.96)=2.216 and x4=30.4(2.216)=2.1136x_4=3-0.4(2.216)=2.1136. A limit LL obeys L=30.4LL=3-0.4L, hence L=3/1.4=15/7L=3/1.4=15/7. Since the constant gradient is negative, the iterates alternate across the fixed point, producing a cobweb; since 0.4<1|-0.4|<1, it converges.
2
  • xn+1=4xnx_{n+1}=\sqrt{4-x_n} is locally convergent to α\alpha and gives a cobweb; xn+1=4xn2x_{n+1}=4-x_n^2 is locally divergent.
5
(5 marks)5
Notes
Both fixed-point equations rearrange to x2+x4=0x^2+x-4=0. For g1(x)=4x2g_1(x)=4-x^2, g1(α)=2α>3>1|g_1'(\alpha)|=2\alpha>3>1, so this recurrence is locally divergent. For g2(x)=4xg_2(x)=\sqrt{4-x}, g2(x)=1/(24x)g_2'(x)=-1/(2\sqrt{4-x}). At the fixed point 4α=α\sqrt{4-\alpha}=\alpha, so g2(α)=1/(2α)<1/3<1|g_2'(\alpha)|=1/(2\alpha)<1/3<1 and the recurrence is locally convergent. Its derivative is negative, so successive graphical steps alternate across the fixed point and form a cobweb.
3
  • x1=32x_1=\dfrac32, x2=65x_2=\dfrac65, x3=1511x_3=\dfrac{15}{11} and x4=3326x_4=\dfrac{33}{26}.
  • The positive fixed point is 1+132\dfrac{-1+\sqrt{13}}2; the construction is a cobweb and the iterates appear to converge.
5
(5 marks)5
Notes
Exact substitution gives x1=3/2x_1=3/2, x2=3/(5/2)=6/5x_2=3/(5/2)=6/5, x3=3/(11/5)=15/11x_3=3/(11/5)=15/11 and x4=3/(26/11)=33/26x_4=3/(26/11)=33/26. A fixed point LL satisfies L=3/(L+1)L=3/(L+1), so L2+L3=0L^2+L-3=0 and the positive solution is (1+13)/2=1.3027(-1+\sqrt{13})/2=1.3027\ldots. The iterates lie alternately above and below this value, with decreasing deviations over the four steps, so the graphical path is a converging cobweb.
4
  • a=1.5a=1.5 and b=0.3b=0.3, so xn+1=1.5+0.3xnx_{n+1}=1.5+0.3x_n.
  • x3=2.112x_3=2.112, the limit is 157\dfrac{15}{7}, and the diagram is a convergent staircase.
5
(5 marks)5
Notes
The two vertical steps show that g(1)=1.8g(1)=1.8 and g(1.8)=2.04g(1.8)=2.04. Hence a+b=1.8a+b=1.8 and a+1.8b=2.04a+1.8b=2.04. Subtracting gives 0.8b=0.240.8b=0.24, so b=0.3b=0.3 and a=1.5a=1.5. The second ordinate is x2=2.04x_2=2.04, so x3=1.5+0.3(2.04)=2.112x_3=1.5+0.3(2.04)=2.112. A limit LL satisfies L=1.5+0.3LL=1.5+0.3L, giving L=15/7L=15/7. Since 0<b<10<b<1, the steps approach the fixed point from one side as a convergent staircase.
5
  • x1=1.632121x_1=1.632121, x2=1.804485x_2=1.804485, x3=1.835441x_3=1.835441, x4=1.840457x_4=1.840457, x5=1.841255x_5=1.841255 and x6=1.841382x_6=1.841382.
  • The least value is n=5n=5, since x6x5=0.000126668<0.0005|x_6-x_5|=0.000126668\ldots<0.0005 but x5x4=0.000798258>0.0005|x_5-x_4|=0.000798258\ldots>0.0005.
  • The fixed point is approximately 1.84141.8414, and the diagram is a convergent staircase.
6
(6 marks)6
Notes
Repeated substitution, retaining the unrounded value each time, gives x1=1.632120558x_1=1.632120558\ldots, x2=1.804485465x_2=1.804485465\ldots, x3=1.835440893x_3=1.835440893\ldots, x4=1.840456855x_4=1.840456855\ldots, x5=1.841255113x_5=1.841255113\ldots and x6=1.841381782x_6=1.841381782\ldots. The first consecutive difference below 0.00050.0005 is therefore x6x5|x_6-x_5|, so n=5n=5. Continuing the iteration gives the fixed point 1.8414056601.841405660\ldots, hence 1.84141.8414 to 44 decimal places. Here g(x)=exg'(x)=e^{-x} is positive and is less than 11 for these positive iterates, so the path is a convergent staircase.

9.3 · Solve equations using the Newton-Raphson method and other recurrence relations of the form xₙ₊₁ = g(xₙ); understand how such methods can fail.

Tier 1 · Easy

Mark scheme for 9.3 Tier 1 · Easy
QuestionSchemeMarks
1
  • x1=2.65x_1=2.65
2
(2 marks)2
Notes
Here f(x)=2xf'(x)=2x. Therefore x1=2.52.5272(2.5)=2.50.755=2.65x_1=2.5-\frac{2.5^2-7}{2(2.5)}=2.5-\frac{-0.75}{5}=2.65.
2
  • x1=1.1036x_1=1.1036
2
(2 marks)2
Notes
Since f(x)=exf'(x)=e^x, Newton-Raphson gives x1=1(e3)/e=3/e=1.103638x_1=1-(e-3)/e=3/e=1.103638\ldots. Hence x1=1.1036x_1=1.1036 to 44 decimal places.

Tier 2 · Standard

Mark scheme for 9.3 Tier 2 · Standard
QuestionSchemeMarks
1
  • x1=1x_1=1, x2=0x_2=0, x3=1x_3=1
  • The iteration enters the cycle 0,1,0,1,0,1,0,1,\ldots.
4
(4 marks)4
Notes
Here f(x)=3x22f'(x)=3x^2-2. From x0=0x_0=0, x1=0f(0)/f(0)=02/(2)=1x_1=0-f(0)/f'(0)=0-2/(-2)=1. Then x2=1f(1)/f(1)=11/1=0x_2=1-f(1)/f'(1)=1-1/1=0, so x3=1x_3=1 again. The values repeat in a two-cycle and neither 00 nor 11 is a root, so this starting value does not produce convergence.
2
  • xn+1=xnxn3+xn263xn2+2xnx_{n+1}=x_n-\dfrac{x_n^3+x_n^2-6}{3x_n^2+2x_n}; x1=1.5385x_1=1.5385 and x2=1.5377x_2=1.5377.
4
(4 marks)4
Notes
Take f(x)=x3+x26f(x)=x^3+x^2-6, so f(x)=3x2+2xf'(x)=3x^2+2x. Substitution in xn+1=xnf(xn)/f(xn)x_{n+1}=x_n-f(x_n)/f'(x_n) gives the stated recurrence. With x0=1.5x_0=1.5, x1=1.538461x_1=1.538461\ldots; substituting this unrounded value gives x2=1.537656x_2=1.537656\ldots. Thus the requested values are 1.53851.5385 and 1.53771.5377.
3
  • xn+1=xnlnxn+xn31+1/xnx_{n+1}=x_n-\dfrac{\ln x_n+x_n-3}{1+1/x_n}; x1=2.20457x_1=2.20457 and x2=2.20794x_2=2.20794.
4
(4 marks)4
Notes
For f(x)=lnx+x3f(x)=\ln x+x-3, f(x)=1/x+1f'(x)=1/x+1. Substitution in the Newton-Raphson formula gives xn+1=xn(lnxn+xn3)/(1+1/xn)x_{n+1}=x_n-(\ln x_n+x_n-3)/(1+1/x_n). Starting from x0=2x_0=2 gives x1=2.204568546x_1=2.204568546\ldots. Using this unrounded value gives x2=2.207939227x_2=2.207939227\ldots, so the requested iterates are 2.204572.20457 and 2.207942.20794.

Tier 3 · Hard

Mark scheme for 9.3 Tier 3 · Hard
QuestionSchemeMarks
1
  • a=10a=10.
  • x2=329104x_2=\dfrac{329}{104} and x3=3.162278x_3=3.162278 to six decimal places.
  • The positive root is 3.16233.1623 correct to four decimal places.
5
(5 marks)5
Notes
For f(x)=x2af(x)=x^2-a, Newton-Raphson gives xn+1=12(xn+a/xn)x_{n+1}=\frac12(x_n+a/x_n). Since x0=4x_0=4 and x1=13/4x_1=13/4, 13/4=12(4+a/4)13/4=\frac12(4+a/4), so 13/2=4+a/413/2=4+a/4 and a=10a=10. Then x2=12(13/4+40/13)=329/104x_2=\frac12(13/4+40/13)=329/104. One more iteration gives x3=12(329/104+1040/329)=216401/68432=3.162277881x_3=\frac12(329/104+1040/329)=216401/68432=3.162277881\ldots. Hence the positive root is 10=3.162277660=3.1623\sqrt{10}=3.162277660\ldots=3.1623 correct to four decimal places.
2
  • x0=1x_0=1 fails because f(1)=0f'(1)=0, so the Newton step is undefined.
  • Starting from 2-2, x1=2.11111x_1=-2.11111 and x2=2.10384x_2=-2.10384.
5
(5 marks)5
Notes
Here f(x)=3x23f'(x)=3x^2-3 and xn+1=xn(xn33xn+3)/(3xn23)x_{n+1}=x_n-(x_n^3-3x_n+3)/(3x_n^2-3). At x0=1x_0=1 the denominator is zero, so no Newton iterate exists. From x0=2x_0=-2, f(2)=1f(-2)=1 and f(2)=9f'(-2)=9, giving x1=21/9=2.111111x_1=-2-1/9=-2.111111\ldots. Substitution of this unrounded value gives x2=2.103835x_2=-2.103835\ldots, hence the stated rounded values.
3
  • xn+1=2xnx_{n+1}=-2x_n; x1=14x_1=-\dfrac14, x2=12x_2=\dfrac12 and x3=1x_3=-1.
  • The iterates alternate in sign while their distance from 00 doubles, so they diverge.
4
(4 marks)4
Notes
For non-zero xx, f(x)=13x2/3f'(x)=\frac13x^{-2/3}. Newton-Raphson therefore gives xn+1=xnxn1/3/(13xn2/3)=xn3xn=2xnx_{n+1}=x_n-x_n^{1/3}/(\frac13x_n^{-2/3})=x_n-3x_n=-2x_n. From x0=1/8x_0=1/8, the next values are 1/4-1/4, 1/21/2 and 1-1. Their signs alternate and their magnitudes double, so the sequence moves away from the root 00 rather than converging to it.
4
  • xn+1=xn+cosxnxnsinxn+1x_{n+1}=x_n+\dfrac{\cos x_n-x_n}{\sin x_n+1}.
  • x1=0.7394x_1=0.7394, x2=0.7391x_2=0.7391 and x3=0.7391x_3=0.7391.
  • cos(0.73905)0.73905>0\cos(0.73905)-0.73905>0 and cos(0.73915)0.73915<0\cos(0.73915)-0.73915<0, so the root is 0.73910.7391 to 44 decimal places.
6
(6 marks)6
Notes
Take f(x)=cosxxf(x)=\cos x-x, so f(x)=sinx1f'(x)=-\sin x-1. Newton-Raphson gives xn+1=xnf(xn)/f(xn)=xn+(cosxnxn)/(sinxn+1)x_{n+1}=x_n-f(x_n)/f'(x_n)=x_n+(\cos x_n-x_n)/(\sin x_n+1). From x0=0.7x_0=0.7, unrounded calculation gives x1=0.739436497x_1=0.739436497\ldots, x2=0.739085160x_2=0.739085160\ldots and x3=0.739085133x_3=0.739085133\ldots. Also f(0.73905)=0.0000587989>0f(0.73905)=0.0000587989\ldots>0 and f(0.73915)=0.0001085633<0f(0.73915)=-0.0001085633\ldots<0. Since f(x)<0f'(x)<0 on this interval, the root is unique there and the bracket pins it as 0.73910.7391 to 44 decimal places.
5
  • xn+1=2xnxn25x_{n+1}=2x_n-\dfrac{x_n^2}{5}.
  • For x0=1x_0=1, x1=95x_1=\dfrac95 and x2=369125x_2=\dfrac{369}{125}.
  • For x0=10x_0=10, x1=0x_1=0, where the original function and the next Newton-Raphson step are undefined.
5
(5 marks)5
Notes
With f(x)=1/x1/5f(x)=1/x-1/5, f(x)=1/x2f'(x)=-1/x^2. Therefore xn+1=xn(1/xn1/5)/(1/xn2)=2xnxn2/5x_{n+1}=x_n-(1/x_n-1/5)/(-1/x_n^2)=2x_n-x_n^2/5. From x0=1x_0=1, this gives x1=21/5=9/5x_1=2-1/5=9/5 and x2=18/581/125=369/125x_2=18/5-81/125=369/125. If x0=10x_0=10, then x1=20100/5=0x_1=20-100/5=0. Since f(0)f(0) is undefined, the method cannot form another tangent step even though the equation itself has the root x=5x=5.

9.4 · Understand and use numerical integration of functions, including the use of the trapezium rule and estimating the approximate area under a curve and limits that it must lie between.

Tier 1 · Easy

Mark scheme for 9.4 Tier 1 · Easy
QuestionSchemeMarks
1
  • 88
2
(2 marks)2
Notes
The strip width is h=1h=1. The trapezium rule gives 12[1+9+2(3)]=12(16)=8\frac12[1+9+2(3)]=\frac12(16)=8.
2
  • k=6k=6
3
(3 marks)3
Notes
The strip width is 11, so the trapezium estimate is 12[2+6+2k]=4+k\frac12[2+6+2k]=4+k. Equating this to 1010 gives k=6k=6.

Tier 2 · Standard

Mark scheme for 9.4 Tier 2 · Standard
QuestionSchemeMarks
1
  • 1919
3
(3 marks)3
Notes
The strip width is h=1h=1. The trapezium rule gives 12[2+10+2(2.5+4+6.5)]=12[12+26]=19\dfrac12[2+10+2(2.5+4+6.5)]=\dfrac12[12+26]=19.
2
  • 2.7932.793; it is an underestimate because y=x+1y=\sqrt{x+1} is concave on the interval, so its chords lie below the curve.
4
(4 marks)4
Notes
Here h=0.5h=0.5. The estimate is 0.52[1+1.732+2(1.225+1.414+1.581)]=2.793\frac{0.5}{2}[1+1.732+2(1.225+1.414+1.581)]=2.793. Also y=1/[4(x+1)3/2]<0y''=-1/[4(x+1)^{3/2}]<0 for 0x20\leq x\leq2, so the graph is concave and each trapezium chord lies below the curve. The estimate is therefore an underestimate.
3
  • 12.512.5
4
(4 marks)4
Notes
The unequal widths are 0.50.5, 11 and 1.51.5, so the three trapezium areas are 0.52(2+4)=1.5\frac{0.5}{2}(2+4)=1.5, 12(4+3)=3.5\frac12(4+3)=3.5 and 1.52(3+7)=7.5\frac{1.5}{2}(3+7)=7.5. Their sum is 1.5+3.5+7.5=12.51.5+3.5+7.5=12.5.

Tier 3 · Hard

Mark scheme for 9.4 Tier 3 · Hard
QuestionSchemeMarks
1
  • Trapezium estimate =0.6970=0.6970.
  • 0.6345<121xdx<0.75950.6345<\int_1^2\frac1x\,dx<0.7595.
6
(6 marks)6
Notes
The width is h=0.25h=0.25 and the ordinates are 1,45,23,47,121,\frac45,\frac23,\frac47,\frac12. The trapezium estimate is 18[1+12+2(45+23+47)]=11711680=0.6970\frac18[1+\frac12+2(\frac45+\frac23+\frac47)]=\frac{1171}{1680}=0.6970. Since 1/x1/x is decreasing, right endpoints give the lower sum 14(45+23+47+12)=533840=0.6345\frac14(\frac45+\frac23+\frac47+\frac12)=\frac{533}{840}=0.6345, while left endpoints give the upper sum 14(1+45+23+47)=319420=0.7595\frac14(1+\frac45+\frac23+\frac47)=\frac{319}{420}=0.7595.
2
  • Estimated signed integral =1=1.
  • Estimated total area =4.5=4.5 square units.
5
(5 marks)5
Notes
With h=1h=1, the signed estimate is 12[2+4+2(1+1)]=1\frac12[-2+4+2(-1+1)]=1. The straight segment from (1,1)(1,-1) to (2,1)(2,1) crosses the axis halfway, at x=1.5x=1.5. The estimated area below the axis is 12(2+1)(1)+12(0.5)(1)=1.75\frac12(2+1)(1)+\frac12(0.5)(1)=1.75. The estimated area above it is 12(0.5)(1)+12(1+4)(1)=2.75\frac12(0.5)(1)+\frac12(1+4)(1)=2.75. Hence the total area is 1.75+2.75=4.51.75+2.75=4.5 square units.
3
  • The estimated area between the curves is 20.520.5 square units.
  • The separate estimates are 3939 and 18.518.5, whose difference is 20.520.5.
5
(5 marks)5
Notes
The ordinates of fgf-g are 6,4,5,6,56,4,5,6,5. With width 11, the trapezium estimate is 12[6+5+2(4+5+6)]=20.5\frac12[6+5+2(4+5+6)]=20.5. Separately, the estimate for 04f(x)dx\int_0^4f(x)\,\mathrm dx is 12[8+14+2(7+9+12)]=39\frac12[8+14+2(7+9+12)]=39, while that for 04g(x)dx\int_0^4g(x)\,\mathrm dx is 12[2+9+2(3+4+6)]=18.5\frac12[2+9+2(3+4+6)]=18.5. Their difference is 3918.5=20.539-18.5=20.5 square units.
4
  • The transformed ordinates are 1,4,9,4,11,4,9,4,1, giving V18πV\approx18\pi.
5
(5 marks)5
Notes
The integrand is r2r^2, not rr, so the ordinates for the trapezium rule are 12,22,32,22,121^2,2^2,3^2,2^2,1^2, namely 1,4,9,4,11,4,9,4,1. With strip width h=1h=1, 04r2dx(1/2)[1+2(4+9+4)+1]=18\int_0^4r^2\,\mathrm dx\approx(1/2)[1+2(4+9+4)+1]=18. Multiplying by π\pi gives the volume estimate 18π18\pi.
5
  • The two-strip estimate is 1010.
  • The missing ordinates are 33 at x=0.5x=0.5 and 66 at x=1.5x=1.5.
5
(5 marks)5
Notes
For two strips of width 11, the estimate is (1/2)[2+2(5)+8]=10(1/2)[2+2(5)+8]=10. Let the missing ordinates at x=0.5x=0.5 and x=1.5x=1.5 be uu and 2u2u. For four strips of width 0.50.5, the estimate is (0.5/2)[2+2(u+5+2u)+8]=5+1.5u(0.5/2)[2+2(u+5+2u)+8]=5+1.5u. Since this equals 9.59.5, 1.5u=4.51.5u=4.5 and u=3u=3. Thus the two missing ordinates are 33 and 66.

9.5 · Use numerical methods to solve problems in context.

Tier 1 · Easy

Mark scheme for 9.5 Tier 1 · Easy
QuestionSchemeMarks
1
  • T(3)20=1T(3)-20=1 and T(4)20=2T(4)-20=-2, so the return time lies between 33 and 44 hours after noon.
2
(2 marks)2
Notes
Define f(t)=T(t)20=2+4tt2f(t)=T(t)-20=-2+4t-t^2. Then f(3)=2+129=1f(3)=-2+12-9=1 and f(4)=2+1616=2f(4)=-2+16-16=-2. The model is continuous and the sign changes, so a solution occurs in (3,4)(3,4).
2
  • S(2)20=10S(2)-20=-10 and S(3)20=10S(3)-20=10, so the model reaches 2020 volts for some t(2,3)t\in(2,3) seconds.
2
(2 marks)2
Notes
Let f(t)=S(t)20=t3+t20f(t)=S(t)-20=t^3+t-20. Then f(2)=8+220=10f(2)=8+2-20=-10 and f(3)=27+320=10f(3)=27+3-20=10. The model is continuous, so the change of sign guarantees a solution between 22 and 33 seconds.

Tier 2 · Standard

Mark scheme for 9.5 Tier 2 · Standard
QuestionSchemeMarks
1
  • h(4.0)=1.6>0h(4.0)=1.6>0 and h(4.1)=0.369<0h(4.1)=-0.369<0
  • h(4.05)=0.62775>0h(4.05)=0.62775>0, so the root lies in (4.05,4.10)(4.05,4.10).
4
(4 marks)4
Notes
Calculate h(4.0)=20(4)4.9(42)=1.6h(4.0)=20(4)-4.9(4^2)=1.6 and h(4.1)=20(4.1)4.9(4.12)=0.369h(4.1)=20(4.1)-4.9(4.1^2)=-0.369. Continuity gives a root in (4.0,4.1)(4.0,4.1). The midpoint is 4.054.05, and h(4.05)=20(4.05)4.9(4.052)=0.62775>0h(4.05)=20(4.05)-4.9(4.05^2)=0.62775>0. The sign change is therefore between 4.054.05 and 4.104.10 seconds.
2
  • t1=2.885203t_1=2.885203, t2=2.875011t_2=2.875011 and t3=2.874092t_3=2.874092.
  • The predicted completion time is 2.8742.874 hours.
4
(4 marks)4
Notes
Successive substitution gives t1=2.885203t_1=2.885203\ldots, t2=2.875011t_2=2.875011\ldots and t3=2.874092t_3=2.874092\ldots. Continuing gives t4=2.874009t_4=2.874009\ldots and the values settle to 2.87400072.8740007\ldots. Therefore the model predicts a completion time of 2.8742.874 hours to 33 decimal places.
3
  • R1=37.4873R_1=37.4873, R2=38.8043R_2=38.8043, R3=38.5557R_3=38.5557, R4=38.6022R_4=38.6022 and R5=38.5935R_5=38.5935.
  • The reliable range is approximately 38.638.6 m.
4
(4 marks)4
Notes
Substituting unrounded values successively gives R1=37.4873055R_1=37.4873055\ldots, R2=38.8042885R_2=38.8042885\ldots, R3=38.5557086R_3=38.5557086\ldots, R4=38.6021810R_4=38.6021810\ldots and R5=38.5934763R_5=38.5934763\ldots. The alternating values are settling near 38.638.6, so the model gives a reliable range of approximately 38.638.6 m to the nearest 0.10.1 m.

Tier 3 · Hard

Mark scheme for 9.5 Tier 3 · Hard
QuestionSchemeMarks
1
  • Estimated inflow =487.5=487.5 litres and estimated final volume =837.5=837.5 litres.
  • The model predicts overflow, with 37.537.5 litres in excess of capacity at 3030 minutes.
5
(5 marks)5
Notes
With spacing h=5h=5, the trapezium estimate is 52[10+11+2(14+18+21+19+15)]=52(195)=487.5\frac52[10+11+2(14+18+21+19+15)]=\frac52(195)=487.5 litres. Adding the initial 350350 litres gives 837.5837.5 litres, which is 837.5800=37.5837.5-800=37.5 litres above capacity, so the tank must have overflowed by the end.
2
  • V(3)70=7V(3)-70=-7 and V(4)70=58V(4)-70=58, so the depth lies in (3,4)(3,4) m.
  • h1=3.133807h_1=3.133807, h2=3.132719h_2=3.132719, so the required depth is 3.133.13 m to 33 significant figures.
6
(6 marks)6
Notes
Set f(h)=h3+4h270f(h)=h^3+4h^2-70. Since f(3)=27+3670=7f(3)=27+36-70=-7 and f(4)=64+6470=58f(4)=64+64-70=58, continuity gives a root in (3,4)(3,4). Also f(h)=3h2+8hf'(h)=3h^2+8h, so hn+1=hn(hn3+4hn270)/(3hn2+8hn)h_{n+1}=h_n-(h_n^3+4h_n^2-70)/(3h_n^2+8h_n). From h0=3.2h_0=3.2, this gives h1=3.1338068=3.133807h_1=3.1338068\ldots=3.133807 (66 d.p.) and then h2=3.1327194=3.132719h_2=3.1327194\ldots=3.132719 (66 d.p.). Thus the required positive depth is 3.133.13 m to 33 significant figures.
3
  • D(5)<30<D(6)D(5)<30<D(6) and D(t)>0D'(t)>0, so there is exactly one required time in (5,6)(5,6).
  • tn+1=tn5tn+8sin(0.5tn)305+4cos(0.5tn)t_{n+1}=t_n-\dfrac{5t_n+8\sin(0.5t_n)-30}{5+4\cos(0.5t_n)}; t1=5.07531t_1=5.07531 and t2=5.12236t_2=5.12236.
  • The time is 5.125.12 s to 33 significant figures.
6
(6 marks)6
Notes
Let f(t)=5t+8sin(0.5t)30f(t)=5t+8\sin(0.5t)-30. Then f(5)=8sin2.55<0f(5)=8\sin2.5-5<0 and f(6)=8sin3>0f(6)=8\sin3>0, so continuity gives a root in (5,6)(5,6). Also f(t)=5+4cos(0.5t)1>0f'(t)=5+4\cos(0.5t)\geq1>0, so the root is unique. Newton-Raphson gives tn+1=tn[5tn+8sin(0.5tn)30]/[5+4cos(0.5tn)]t_{n+1}=t_n-[5t_n+8\sin(0.5t_n)-30]/[5+4\cos(0.5t_n)]. From t0=5.5t_0=5.5, unrounded substitution gives t1=5.075305513t_1=5.075305513\ldots and t2=5.122362632t_2=5.122362632\ldots. Thus the estimated time is 5.125.12 s to 33 significant figures.
4
  • For f(t)=t+2et3.5f(t)=t+2e^{-t}-3.5, f(3)<0f(3)<0 and f(4)>0f(4)>0; also f(t)=12et>0f'(t)=1-2e^{-t}>0 on (3,4)(3,4), so the root there is unique.
  • The five tested midpoints are 3.53.5, 3.253.25, 3.3753.375, 3.43753.4375 and 3.406253.40625.
  • The resulting interval is (3.40625,3.4375)(3.40625,3.4375), so the midpoint estimate is 3.4218753.421875 hours with maximum possible error 0.0156250.015625 hours.
7
(7 marks)7
Notes
Let f(t)=t+2et3.5f(t)=t+2e^{-t}-3.5. Then f(3)=0.400425<0f(3)=-0.400425\ldots<0 and f(4)=0.536631>0f(4)=0.536631\ldots>0. Since f(t)=12et>0f'(t)=1-2e^{-t}>0 when 3<t<43<t<4, there is exactly one root in the interval. The midpoint signs are f(3.5)>0f(3.5)>0, f(3.25)<0f(3.25)<0, f(3.375)<0f(3.375)<0, f(3.4375)>0f(3.4375)>0 and f(3.40625)<0f(3.40625)<0. The successive brackets therefore end at (3.40625,3.4375)(3.40625,3.4375). Its midpoint is (3.40625+3.4375)/2=3.421875(3.40625+3.4375)/2=3.421875, and half the interval width is (3.43753.40625)/2=0.015625(3.4375-3.40625)/2=0.015625, which is the maximum possible error.
5
  • For f(c)=c+lnc4f(c)=c+\ln c-4, f(2.9)<0<f(3)f(2.9)<0<f(3) and f(c)=1+1/c>0f'(c)=1+1/c>0, so there is exactly one possible factor in (2.9,3)(2.9,3).
  • f(2.92)<0<f(2.93)f(2.92)<0<f(2.93), and then f(2.926)<0<f(2.927)f(2.926)<0<f(2.927), so 2.926<c<2.9272.926<c<2.927.
  • The calibration factor is 2.932.93 to 33 significant figures.
6
(6 marks)6
Notes
Let f(c)=c+lnc4f(c)=c+\ln c-4. Then f(2.9)=0.035289f(2.9)=-0.035289\ldots and f(3)=0.098612f(3)=0.098612\ldots, so continuity gives a root in (2.9,3)(2.9,3). Since f(c)=1+1/c>0f'(c)=1+1/c>0 for c>1c>1, this root is unique. Now f(2.92)=0.008416f(2.92)=-0.008416\ldots and f(2.93)=0.005002f(2.93)=0.005002\ldots, giving the width-0.010.01 bracket (2.92,2.93)(2.92,2.93). Refining within it, f(2.926)=0.000363f(2.926)=-0.000363\ldots and f(2.927)=0.000978f(2.927)=0.000978\ldots, giving the width-0.0010.001 bracket (2.926,2.927)(2.926,2.927). Every number in this final bracket rounds to 2.932.93 to 33 significant figures.