FP1-6 Further numerical methods (Further Pure 1) — revision question pack

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

Checked against Edexcel 9FM0 section FP1-6. Review basis: the qualification registry sourced from the Pearson Edexcel Level 3 Advanced GCE in Further Mathematics (9FM0) specification; registry verification recorded 17 July 2026.

How this checking works

FP1-6.1 · Numerical solution of first order and second order differential equations.

Explanation

  • Three finite-difference approximations are quoted in the specification and must be known. The forward difference (dydx)nyn+1ynh\left(\dfrac{\mathrm{d}y}{\mathrm{d}x}\right)_n\approx\dfrac{y_{n+1}-y_n}{h} steps forward from a single known value.
  • The central difference (dydx)nyn+1yn12h\left(\dfrac{\mathrm{d}y}{\mathrm{d}x}\right)_n\approx\dfrac{y_{n+1}-y_{n-1}}{2h} is more accurate but needs two consecutive starting values.
  • The second-derivative formula is (d2ydx2)nyn+12yn+yn1h2\left(\dfrac{\mathrm{d}^2y}{\mathrm{d}x^2}\right)_n\approx\dfrac{y_{n+1}-2y_n+y_{n-1}}{h^2}.
  • In every case, rearrange to make yn+1y_{n+1} the subject and then step forward, keeping more decimal places in the working than the question asks for in the answer.
  • The forward difference has error of order hh, whereas both central formulae have error of order h2h^2, which is why halving the step size reduces their error by about a factor of four.

Worked example

Given dydx=x+y\dfrac{\mathrm{d}y}{\mathrm{d}x}=x+y and y=1y=1 at x=0x=0, use the approximation (dydx)nyn+1ynh\left(\dfrac{\mathrm{d}y}{\mathrm{d}x}\right)_n\approx\dfrac{y_{n+1}-y_n}{h} with h=0.1h=0.1 to estimate yy at x=0.2x=0.2.

  1. 1.Rearranged, yn+1=yn+h(xn+yn)y_{n+1}=y_n+h\left(x_n+y_n\right).
  2. 2.y1=1+0.1(0+1)=1.1y_1=1+0.1(0+1)=1.1.
  3. 3.y2=1.1+0.1(0.1+1.1)=1.1+0.12=1.22y_2=1.1+0.1(0.1+1.1)=1.1+0.12=1.22.

Answer: y(0.2)1.22y(0.2)\approx1.22.

Common mistakes

  • Don't fall into the trap of using xn+1x_{n+1} instead of xnx_n inside the forward-difference step.
  • Don't fall into the trap of starting a central-difference calculation from a single value, when it needs two.
  • Don't fall into the trap of rounding the intermediate values to the accuracy demanded in the final answer, which loses a digit at every step.

Exam tip

Write the rearranged recurrence on its own line before any arithmetic, then tabulate nn, xnx_n and yny_n; the method mark is for the recurrence.

Tier 1 · Easy

  1. 1.

    Given dydx=x+y\dfrac{\mathrm{d}y}{\mathrm{d}x}=x+y with y=1y=1 at x=0x=0, use the approximation (dydx)nyn+1ynh\left(\dfrac{\mathrm{d}y}{\mathrm{d}x}\right)_n\approx\dfrac{y_{n+1}-y_n}{h} with h=0.1h=0.1 to estimate yy at x=0.1x=0.1 and at x=0.2x=0.2.

    (4)

    (Total for Question 1 is 4 marks)

  2. 2.

    Given dydx=x+y\dfrac{\mathrm{d}y}{\mathrm{d}x}=\sqrt{x+y} with y=2y=2 at x=1x=1, use the forward-difference approximation with h=0.1h=0.1 to estimate y(1.2)y(1.2), giving your answer to 44 decimal places.

    (5)

    (Total for Question 2 is 5 marks)

Tier 2 · Standard

  1. 1.

    For dydx=x2y\dfrac{\mathrm{d}y}{\mathrm{d}x}=x^2-y it is given that y(1)=0.5y(1)=0.5 and y(1.1)=0.55y(1.1)=0.55. Use the approximation (dydx)nyn+1yn12h\left(\dfrac{\mathrm{d}y}{\mathrm{d}x}\right)_n\approx\dfrac{y_{n+1}-y_{n-1}}{2h} with h=0.1h=0.1 to estimate y(1.2)y(1.2).

    (5)

    (Total for Question 1 is 5 marks)

  2. 2.

    For d2ydx2=xy\dfrac{\mathrm{d}^2y}{\mathrm{d}x^2}=x-y it is given that y(0)=1y(0)=1 and y(0.1)=1.05y(0.1)=1.05. Use the approximation (d2ydx2)nyn+12yn+yn1h2\left(\dfrac{\mathrm{d}^2y}{\mathrm{d}x^2}\right)_n\approx\dfrac{y_{n+1}-2y_n+y_{n-1}}{h^2} with h=0.1h=0.1 to estimate y(0.2)y(0.2).

    (5)

    (Total for Question 2 is 5 marks)

  3. 3.

    Use the forward-difference approximation with h=0.1h=0.1 to estimate y(0.2)y(0.2) for dydx=y\dfrac{\mathrm{d}y}{\mathrm{d}x}=y with y(0)=1y(0)=1, and compare your answer with the exact value.

    (5)

    (Total for Question 3 is 5 marks)

Tier 3 · Hard

  1. 1.

    Given dydx=x+y2\dfrac{\mathrm{d}y}{\mathrm{d}x}=x+y^2 with y=0.5y=0.5 at x=0x=0, use the forward-difference approximation with h=0.05h=0.05 to estimate y(0.1)y(0.1), giving your answer to 44 decimal places.

    (6)

    (Total for Question 1 is 6 marks)

  2. 2.

    For d2ydx2+dydx+y=0\dfrac{\mathrm{d}^2y}{\mathrm{d}x^2}+\dfrac{\mathrm{d}y}{\mathrm{d}x}+y=0 it is given that y(0)=1y(0)=1 and y(0.1)=0.995y(0.1)=0.995. Using h=0.1h=0.1 with the central-difference approximations for both derivatives, estimate y(0.2)y(0.2).

    (7)

    (Total for Question 2 is 7 marks)

  3. 3.

    Use Taylor's series to derive the approximation (d2ydx2)nyn+12yn+yn1h2\left(\dfrac{\mathrm{d}^2y}{\mathrm{d}x^2}\right)_n\approx\dfrac{y_{n+1}-2y_n+y_{n-1}}{h^2} and show that its error is of order h2h^2.

    (6)

    (Total for Question 3 is 6 marks)

  4. 4.

    Given dydx=11+x2\dfrac{\mathrm{d}y}{\mathrm{d}x}=\dfrac{1}{1+x^2} with y(0)=0y(0)=0, use the forward-difference approximation with h=0.25h=0.25 to estimate y(1)y(1), and compare your answer with the exact value.

    (7)

    (Total for Question 4 is 7 marks)

  5. 5.

    Explain why the central-difference approximation (dydx)nyn+1yn12h\left(\dfrac{\mathrm{d}y}{\mathrm{d}x}\right)_n\approx\dfrac{y_{n+1}-y_{n-1}}{2h} requires two starting values, and for dydx=y\dfrac{\mathrm{d}y}{\mathrm{d}x}=y with y(0)=1y(0)=1 and y(0.1)=1.1051709y(0.1)=1.1051709 use it with h=0.1h=0.1 to estimate y(0.2)y(0.2), comparing the error with that of the forward difference.

    (7)

    (Total for Question 5 is 7 marks)

FP1-6.2 · Simpson's rule.

Explanation

  • Simpson's rule estimates abydx\displaystyle\int_a^b y\,\mathrm{d}x by fitting a quadratic through each consecutive group of three ordinates.
  • With nn strips of width h=banh=\dfrac{b-a}{n} and ordinates y0,y1,,yny_0,y_1,\dots,y_n, the rule is abydxh3[(y0+yn)+4(y1+y3+)+2(y2+y4+)]\displaystyle\int_a^b y\,\mathrm{d}x\approx\frac{h}{3}\left[\left(y_0+y_n\right)+4\left(y_1+y_3+\cdots\right)+2\left(y_2+y_4+\cdots\right)\right].
  • The number of strips must be even, because the ordinates are consumed three at a time with the outer ones shared, so the strips pair up.
  • The end ordinates take coefficient 11, the odd-numbered ones take 44 and the even-numbered interior ones take 22; a quick check is that the coefficients must total 3n3n.
  • The rule is exact for every polynomial of degree at most three, and its error is of order h4h^4, so halving the strip width reduces the error by roughly a factor of sixteen.

Worked example

Use Simpson's rule with four strips to estimate 12lnxdx\displaystyle\int_1^2\ln x\,\mathrm{d}x, giving your answer to 55 decimal places.

  1. 1.h=0.25h=0.25 and the ordinates are 00, 0.22314360.2231436, 0.40546510.4054651, 0.55961580.5596158 and 0.69314720.6931472.
  2. 2.The bracket is (0+0.6931472)+4(0.2231436+0.5596158)+2(0.4054651)(0+0.6931472)+4(0.2231436+0.5596158)+2(0.4054651).
  3. 3.=0.6931472+3.1310376+0.8109302=4.635115=0.6931472+3.1310376+0.8109302=4.635115.
  4. 4.Multiplying by h3=0.0833333\dfrac{h}{3}=0.0833333 gives 0.38625960.3862596.

Answer: 12lnxdx0.38626\displaystyle\int_1^2\ln x\,\mathrm{d}x\approx0.38626.

Common mistakes

  • Don't fall into the trap of using an odd number of strips, for which the rule is not defined.
  • Don't fall into the trap of attaching the coefficient 44 to the even-numbered ordinates instead of the odd-numbered ones.
  • Don't fall into the trap of using h2\dfrac{h}{2} from the trapezium rule instead of h3\dfrac{h}{3}.

Exam tip

Tabulate nn, xnx_n, yny_n and the coefficient in four columns and check that the coefficients sum to 3n3n before multiplying.

Tier 1 · Easy

  1. 1.

    Use Simpson's rule with four strips to estimate 011+x3dx\displaystyle\int_0^1\sqrt{1+x^3}\,\mathrm{d}x, giving your answer to 44 decimal places.

    (4)

    (Total for Question 1 is 4 marks)

  2. 2.

    Use Simpson's rule with four strips to estimate 0π/2sinxdx\displaystyle\int_0^{\pi/2}\sin x\,\mathrm{d}x, giving your answer to 66 decimal places, and state the error.

    (5)

    (Total for Question 2 is 5 marks)

Tier 2 · Standard

  1. 1.

    Use Simpson's rule with four strips to estimate 01ex2dx\displaystyle\int_0^1\mathrm{e}^{-x^2}\,\mathrm{d}x, giving your answer to 66 decimal places.

    (5)

    (Total for Question 1 is 5 marks)

  2. 2.

    Use Simpson's rule with four strips to estimate 13lnxxdx\displaystyle\int_1^3\frac{\ln x}{x}\,\mathrm{d}x, giving your answer to 44 decimal places, and compare it with the exact value 12(ln3)2\tfrac12\left(\ln3\right)^2.

    (6)

    (Total for Question 2 is 6 marks)

  3. 3.

    Explain why Simpson's rule requires an even number of strips, and state the coefficient pattern used for six strips.

    (4)

    (Total for Question 3 is 4 marks)

Tier 3 · Hard

  1. 1.

    Show that Simpson's rule is exact for 02hx3dx\displaystyle\int_0^{2h}x^3\,\mathrm{d}x.

    (5)

    (Total for Question 1 is 5 marks)

  2. 2.

    Use Simpson's rule with six strips to estimate 0π/3secxdx\displaystyle\int_0^{\pi/3}\sec x\,\mathrm{d}x, giving your answer to 44 decimal places, and compare it with the exact value ln(2+3)\ln\left(2+\sqrt3\right).

    (7)

    (Total for Question 2 is 7 marks)

  3. 3.

    Use Simpson's rule with four strips and then with eight strips to estimate 12dxx\displaystyle\int_1^2\frac{\mathrm{d}x}{x}, and use the two errors to confirm that the error behaves like h4h^4.

    (8)

    (Total for Question 3 is 8 marks)

  4. 4.

    The speed vm s1v\,\text{m s}^{-1} of a particle at times t=0t=0, 0.50.5, 11, 1.51.5 and 22 seconds is 00, 2.42.4, 4.14.1, 5.05.0 and 5.35.3 respectively. Use Simpson's rule to estimate the distance travelled in the first 22 seconds.

    (5)

    (Total for Question 4 is 5 marks)

  5. 5.

    A student applies Simpson's rule to 01f(x)dx\displaystyle\int_0^1 f(x)\,\mathrm{d}x with five strips. Explain why the calculation is invalid, and state the smallest number of strips greater than five that could be used together with the coefficient pattern that results.

    (4)

    (Total for Question 5 is 4 marks)

Answer key

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

FP1-6.1 · Numerical solution of first order and second order differential equations.

Tier 1 · Easy

Mark scheme for FP1-6.1 Tier 1 · Easy
QuestionSchemeMarks
1
  • yn+1=yn+h(xn+yn)y_{n+1}=y_n+h\left(x_n+y_n\right)
  • y1=1+0.1(1)=1.1y_1=1+0.1(1)=1.1
  • y2=1.1+0.1(1.2)=1.22y_2=1.1+0.1(1.2)=1.22
4
(4 marks)4
Notes
Rearrange the forward difference and step twice. Independent check against the exact solution y=2exx1y=2\mathrm{e}^{x}-x-1: y(0.1)=1.11034y(0.1)=1.11034 and y(0.2)=1.24281y(0.2)=1.24281, so the estimates are low, as expected from a first-order method.
2
  • yn+1=yn+hxn+yny_{n+1}=y_n+h\sqrt{x_n+y_n}
  • y1=2+0.13=2.1732051y_1=2+0.1\sqrt{3}=2.1732051
  • y2=2.1732051+0.13.2732051=2.3541251y_2=2.1732051+0.1\sqrt{3.2732051}=2.3541251
  • y(1.2)2.3541y(1.2)\approx2.3541
5
(5 marks)5
Notes
Apply the recurrence twice, keeping seven decimal places in the working. Independent check by halving the step to h=0.05h=0.05 and taking four steps: the estimate becomes 2.35802.3580, closer to the fourth-order value 2.36202.3620, and both step-wise estimates fall short of it as expected for a first-order method on an increasing solution.

Tier 2 · Standard

Mark scheme for FP1-6.1 Tier 2 · Standard
QuestionSchemeMarks
1
  • At x=1.1x=1.1: y(1.2)y(1.0)0.21.120.55\dfrac{y(1.2)-y(1.0)}{0.2}\approx1.1^2-0.55
  • =1.210.55=0.66=1.21-0.55=0.66
  • y(1.2)0.5+0.2(0.66)y(1.2)\approx0.5+0.2(0.66)
  • =0.632=0.632
5
(5 marks)5
Notes
Apply the central difference at the middle point x=1.1x=1.1, where both neighbouring ordinates are involved, and rearrange for the unknown one. Independent check by the forward difference from x=1.1x=1.1: y(1.2)0.55+0.1(0.66)=0.616y(1.2)\approx0.55+0.1(0.66)=0.616, and the central-difference value lies on the more accurate side of it.
2
  • y2=2y1y0+h2(x1y1)y_2=2y_1-y_0+h^2\left(x_1-y_1\right)
  • =2(1.05)1+0.01(0.11.05)=2(1.05)-1+0.01(0.1-1.05)
  • =1.10.0095=1.1-0.0095
  • y(0.2)1.0905y(0.2)\approx1.0905
5
(5 marks)5
Notes
Rearrange the second-derivative formula to make yn+1y_{n+1} the subject and substitute at n=1n=1. Independent check: the second derivative at x=0.1x=0.1 is 0.11.05=0.950.1-1.05=-0.95, so the sequence should be slightly concave, and 1,1.05,1.09051,1.05,1.0905 has second difference 0.0095=h2×(0.95)-0.0095=h^2\times(-0.95) exactly as the formula demands.
3
  • yn+1=yn(1+h)y_{n+1}=y_n(1+h)
  • y1=1.1y_1=1.1 and y2=1.21y_2=1.21
  • The exact solution is y=exy=\mathrm{e}^{x}, so y(0.2)=1.2214y(0.2)=1.2214
  • The estimate is low by about 0.01140.0114
5
(5 marks)5
Notes
The recurrence multiplies by 1+h1+h at each step, giving (1.1)2=1.21(1.1)^2=1.21. Independent check by the series: e0.2=1+0.2+0.02+0.00133+=1.22140\mathrm{e}^{0.2}=1+0.2+0.02+0.00133+\cdots=1.22140, so the forward difference has lost the quadratic term, an error of order hh per step as expected.

Tier 3 · Hard

Mark scheme for FP1-6.1 Tier 3 · Hard
QuestionSchemeMarks
1
  • yn+1=yn+h(xn+yn2)y_{n+1}=y_n+h\left(x_n+y_n^{2}\right)
  • y1=0.5+0.05(0+0.25)=0.5125y_1=0.5+0.05(0+0.25)=0.5125
  • y2=0.5125+0.05(0.05+0.51252)y_2=0.5125+0.05\left(0.05+0.5125^{2}\right)
  • =0.5125+0.05(0.3126563)=0.5281328=0.5125+0.05(0.3126563)=0.5281328
  • y(0.1)0.5281y(0.1)\approx0.5281
6
(6 marks)6
Notes
Two forward-difference steps with the squared term evaluated at the current value of yy. Independent check by the Taylor series method: y(0)=0.25y'(0)=0.25, y=1+2yyy''=1+2yy' gives y(0)=1.25y''(0)=1.25, so y(0.1)0.5+0.025+0.00625=0.53125y(0.1)\approx0.5+0.025+0.00625=0.53125, and the step-wise estimate lies just below it, as a first-order method should.
2
  • At x=0.1x=0.1: y22y1+y0h2+y2y02h+y1=0\dfrac{y_2-2y_1+y_0}{h^2}+\dfrac{y_2-y_0}{2h}+y_1=0
  • 100(y21.99+1)+5(y21)+0.995=0100\left(y_2-1.99+1\right)+5\left(y_2-1\right)+0.995=0
  • 100y299+5y25+0.995=0100y_2-99+5y_2-5+0.995=0
  • 105y2=103.005105y_2=103.005
  • y(0.2)0.9810y(0.2)\approx0.9810
7
(7 marks)7
Notes
Substitute both finite-difference approximations at the interior point x=0.1x=0.1, which leaves a single linear equation in the unknown y2y_2. Independent check against the exact solution ex/2(cos3x2+13sin3x2)\mathrm{e}^{-x/2}\left(\cos\tfrac{\sqrt3x}{2}+\tfrac{1}{\sqrt3}\sin\tfrac{\sqrt3x}{2}\right), which gives y(0.1)=0.99517y(0.1)=0.99517 and y(0.2)=0.98133y(0.2)=0.98133, within 0.00040.0004 of the estimate.
3
  • y(x+h)=y+hy+h22y+h36y+h424y(4)+y(x+h)=y+hy'+\dfrac{h^2}{2}y''+\dfrac{h^3}{6}y'''+\dfrac{h^4}{24}y^{(4)}+\cdots
  • y(xh)=yhy+h22yh36y+h424y(4)y(x-h)=y-hy'+\dfrac{h^2}{2}y''-\dfrac{h^3}{6}y'''+\dfrac{h^4}{24}y^{(4)}-\cdots
  • Adding, the odd-order terms cancel: y(x+h)+y(xh)=2y+h2y+h412y(4)+y(x+h)+y(x-h)=2y+h^2y''+\dfrac{h^4}{12}y^{(4)}+\cdots
  • y(x+h)2y(x)+y(xh)h2=y+h212y(4)+\dfrac{y(x+h)-2y(x)+y(x-h)}{h^2}=y''+\dfrac{h^2}{12}y^{(4)}+\cdots
  • The leading error term is h212y(4)\dfrac{h^2}{12}y^{(4)}, so the error is of order h2h^2
6
(6 marks)6
Notes
Expanding both neighbouring values and adding cancels every odd derivative, leaving h2yh^2y'' as the first surviving term. Independent check with y=x4y=x^4 at x=1x=1, h=0.1h=0.1: the formula gives 1.46412+0.65610.01=12.02\dfrac{1.4641-2+0.6561}{0.01}=12.02, against the exact y=12y''=12, an error of 0.02=h212×240.02=\tfrac{h^2}{12}\times24.
4
  • yn+1=yn+0.251+xn2y_{n+1}=y_n+\dfrac{0.25}{1+x_n^{2}}
  • y1=0.25y_1=0.25, y2=0.4852941y_2=0.4852941, y3=0.6852941y_3=0.6852941
  • y4=0.8452941y_4=0.8452941
  • The exact solution is y=arctanxy=\arctan x, so y(1)=π4=0.7853982y(1)=\dfrac{\pi}{4}=0.7853982
  • The estimate is high by about 0.060.06
7
(7 marks)7
Notes
Four forward-difference steps using the left-hand ordinate at each stage. Because the integrand is decreasing, using its value at the left of each interval overestimates the area, so the estimate exceeds arctan1\arctan1. Independent check: the same calculation is the left-hand rectangle rule for 01dx1+x2\displaystyle\int_0^1\frac{\mathrm{d}x}{1+x^2}, and the trapezium rule with the same ordinates gives 0.78280.7828, much closer to π4\tfrac{\pi}{4}.
5
  • The formula links three consecutive ordinates, so two must be known before the third can be found
  • y2=y0+2hy1=1+0.2(1.1051709)y_2=y_0+2hy_1=1+0.2(1.1051709)
  • =1.2210342=1.2210342
  • The exact value is e0.2=1.2214028\mathrm{e}^{0.2}=1.2214028, an error of 0.000370.00037
  • The forward difference gave 1.211.21, an error of 0.01140.0114, about 3030 times larger
7
(7 marks)7
Notes
Rearranging the central difference gives yn+1=yn1+2hf(xn,yn)y_{n+1}=y_{n-1}+2hf\left(x_n,y_n\right), which involves three ordinates at once and so cannot start from a single value. Independent check of the order of accuracy: the central difference has error of order h2h^2 and the forward difference error of order hh, so at h=0.1h=0.1 the ratio of errors should be roughly h1=10h^{-1}=10 or more, consistent with the observed factor of about 3030.

FP1-6.2 · Simpson's rule.

Tier 1 · Easy

Mark scheme for FP1-6.2 Tier 1 · Easy
QuestionSchemeMarks
1
  • h=0.25h=0.25 and the ordinates are 11, 1.00778221.0077822, 1.06066021.0606602, 1.19242401.1924240 and 1.41421361.4142136
  • The bracket is 2.4142136+4(2.2002062)+2(1.0606602)=13.33635882.4142136+4(2.2002062)+2(1.0606602)=13.3363588
  • The estimate is 0.253×13.3363588=1.1114\dfrac{0.25}{3}\times13.3363588=1.1114
4
(4 marks)4
Notes
Tabulate the five ordinates, apply the 1,4,2,4,11,4,2,4,1 pattern and multiply by h3\tfrac{h}{3}. Independent check with eight strips: the estimate becomes 1.11141.1114 to four decimal places as well, confirming convergence.
2
  • h=π8h=\dfrac{\pi}{8} and the ordinates are 00, 0.38268340.3826834, 0.70710680.7071068, 0.92387950.9238795 and 11
  • The bracket is 1+4(1.3065629)+2(0.7071068)=7.64046521+4(1.3065629)+2(0.7071068)=7.6404652
  • The estimate is π24×7.6404652=1.000135\dfrac{\pi}{24}\times7.6404652=1.000135
  • The exact value is 11, so the error is 0.0001350.000135
5
(5 marks)5
Notes
Apply the rule with h=π8h=\tfrac{\pi}{8} and compare with the exact integral [cosx]0π/2=1\left[-\cos x\right]_0^{\pi/2}=1. Independent check of the fourth-order error: with eight strips the estimate is 1.00000831.0000083, about sixteen times closer, exactly as an h4h^4 error predicts.

Tier 2 · Standard

Mark scheme for FP1-6.2 Tier 2 · Standard
QuestionSchemeMarks
1
  • h=0.25h=0.25 and the ordinates are 11, 0.93941310.9394131, 0.77880080.7788008, 0.56978280.5697828 and 0.36787940.3678794
  • The bracket is 1.3678794+4(1.5091959)+2(0.7788008)=8.96226461.3678794+4(1.5091959)+2(0.7788008)=8.9622646
  • The estimate is 0.253×8.9622646=0.746855\dfrac{0.25}{3}\times8.9622646=0.746855
5
(5 marks)5
Notes
The integrand has no elementary antiderivative, which is why a numerical method is required; apply the standard coefficient pattern. Independent check: the true value of this integral is 0.74682410.7468241, so Simpson's rule with only four strips is already accurate to five significant figures.
2
  • h=0.5h=0.5 and the ordinates are 00, 0.27031010.2703101, 0.34657360.3465736, 0.36651630.3665163 and 0.36620410.3662041
  • The bracket is 0.3662041+4(0.6368264)+2(0.3465736)=3.60665690.3662041+4(0.6368264)+2(0.3465736)=3.6066569
  • The estimate is 0.53×3.6066569=0.6011\dfrac{0.5}{3}\times3.6066569=0.6011
  • The exact value is 12(ln3)2=0.6035\tfrac12\left(\ln3\right)^2=0.6035, so the estimate is low by about 0.00240.0024
6
(6 marks)6
Notes
Apply the rule with h=0.5h=0.5; the relatively large strip width and the turning point of the integrand at x=ex=\mathrm{e} make the error visible at the third decimal place. Independent check by substitution: u=lnxu=\ln x gives 0ln3udu=12(ln3)2=0.6034745\displaystyle\int_0^{\ln3}u\,\mathrm{d}u=\tfrac12\left(\ln3\right)^2=0.6034745.
3
  • The rule fits a quadratic through each consecutive group of three ordinates
  • Each quadratic spans two strips, so the strips must pair up exactly
  • An odd number of strips would leave one strip without a partner
  • For six strips the pattern is 1,4,2,4,2,4,11,4,2,4,2,4,1
4
(4 marks)4
Notes
Each parabolic arc covers two adjacent strips and shares its end ordinates with the neighbouring arcs, which is why an even number is needed. Independent check on the coefficients: they sum to 1+4+2+4+2+4+1=18=3n1+4+2+4+2+4+1=18=3n with n=6n=6, the required total.

Tier 3 · Hard

Mark scheme for FP1-6.2 Tier 3 · Hard
QuestionSchemeMarks
1
  • The ordinates are y0=0y_0=0, y1=h3y_1=h^3 and y2=8h3y_2=8h^3
  • Simpson's rule gives h3[0+4h3+8h3]=h3(12h3)=4h4\dfrac{h}{3}\left[0+4h^3+8h^3\right]=\dfrac{h}{3}\left(12h^3\right)=4h^4
  • The exact integral is [x44]02h=4h4\left[\dfrac{x^4}{4}\right]_0^{2h}=4h^4
  • The two agree, so the rule is exact for this cubic
5
(5 marks)5
Notes
Evaluate the three ordinates of x3x^3 at 00, hh and 2h2h, apply the rule and compare with the exact integral. Independent check of the general claim: the rule is exact for 11, xx and x2x^2 because it integrates the interpolating quadratic exactly, and the error term for a general function is proportional to the fourth derivative, which vanishes for every cubic.
2
  • h=π18h=\dfrac{\pi}{18} and the ordinates are 11, 1.01542661.0154266, 1.06417781.0641778, 1.15470051.1547005, 1.30540731.3054073, 1.55572381.5557238 and 22
  • The bracket is 3+4(3.7258510)+2(2.3695851)=22.64257403+4(3.7258510)+2(2.3695851)=22.6425740
  • The estimate is π54×22.6425740=1.3173\dfrac{\pi}{54}\times22.6425740=1.3173
  • The exact value is ln(2+3)=1.3170\ln\left(2+\sqrt3\right)=1.3170, so the estimate is high by about 0.00030.0003
7
(7 marks)7
Notes
Apply the rule with seven ordinates and the 1,4,2,4,2,4,11,4,2,4,2,4,1 pattern. Independent check of the exact value: secxdx=lnsecx+tanx\displaystyle\int\sec x\,\mathrm{d}x=\ln\left|\sec x+\tan x\right|, and at x=π3x=\tfrac{\pi}{3} this is ln(2+3)=1.3169579\ln\left(2+\sqrt3\right)=1.3169579, while at x=0x=0 it is ln1=0\ln1=0.
3
  • With four strips the estimate is 0.69325400.6932540
  • With eight strips the estimate is 0.69315450.6931545
  • The exact value is ln2=0.6931472\ln2=0.6931472
  • The two errors are 1.068×1041.068\times10^{-4} and 7.35×1067.35\times10^{-6}
  • Their ratio is about 14.514.5, close to the factor of 1616 that halving hh predicts for an h4h^4 error
8
(8 marks)8
Notes
Apply the rule twice with h=0.25h=0.25 and h=0.125h=0.125 and compare each estimate with ln2\ln2. Independent check of the theory: an error proportional to h4h^4 should shrink by 24=162^4=16 when the strip width is halved, and the small shortfall from 1616 is the usual effect of the higher-order terms at these step sizes.
4
  • h=0.5h=0.5 with five ordinates
  • The bracket is (0+5.3)+4(2.4+5.0)+2(4.1)(0+5.3)+4(2.4+5.0)+2(4.1)
  • =5.3+29.6+8.2=43.1=5.3+29.6+8.2=43.1
  • The distance is 0.53×43.1=7.18m\dfrac{0.5}{3}\times43.1=7.18\,\text{m}
5
(5 marks)5
Notes
Distance is the integral of speed with respect to time, so apply the rule to the tabulated ordinates with h=0.5h=0.5. Independent check by the trapezium rule on the same data: 0.52[0+5.3+2(2.4+4.1+5.0)]=7.08m\tfrac{0.5}{2}\left[0+5.3+2(2.4+4.1+5.0)\right]=7.08\,\text{m}, slightly lower as expected for a concave speed-time graph.
5
  • Simpson's rule pairs the strips into parabolic arcs, so an even number is required
  • With five strips one strip is left over and no quadratic can be fitted to it
  • The smallest valid number greater than five is six
  • The pattern is 1,4,2,4,2,4,11,4,2,4,2,4,1, whose coefficients sum to 18=3×618=3\times6
4
(4 marks)4
Notes
The rule integrates a quadratic through each consecutive triple of ordinates, and each such quadratic covers exactly two strips, so an odd count cannot be partitioned. Independent check: with nn strips there are n+1n+1 ordinates, of which n2\tfrac{n}{2} carry coefficient 44 and n21\tfrac{n}{2}-1 carry coefficient 22; both counts are whole numbers only when nn is even.