Using Integration · applications

Applications: Using Integration

11 question types · Pre-University · each worked step by step with a figure that follows the steps

01

A Survey Vehicle That Rises, Sinks and Rises Again: Where It Ends and How Far It Went

methodIntegrate the Velocity Once for the Displacement, Then Find the Times at Which the Velocity Is Zero, Integrate Between Them Separately and Add the Three Answers Without Their Signs for the Distance Actually Traveled

A remotely operated vehicle hangs on a tether from a survey ship. Taking upward as positive, its velocity n minutes after a maneuver begins is v = 1.5n2 − 12n + 18 meters per minute, for 0 ≤ n ≤ 8. When the maneuver begins the vehicle is 40 meters below the surface. (a) How deep is it after the eight minutes? (b) How much tether does the winch move in and out altogether during the maneuver?

−50510152002468n, minutes after the maneuver startsvelocity, m/mins = 0.5n3− 6n2+ 18n
Displacement is an antiderivative of velocity: s = 0.5n3 − 6n2 + 18n, taken as 0 at the start.
Displacement is an antiderivative of velocity. Integrating term by term, s = ∫ (1.5n2 − 12n + 18)dn = 0.5n3 − 6n2 + 18n, taking s = 0 at the start of the maneuver.
step 1 of 5

Velocity is the derivative of displacement, so displacement is an integral of velocity. That integral counts a minute of sinking against a minute of rising, which is exactly what is wanted for the depth at the end, and exactly what is not wanted for the tether. For the tether the interval has to be cut at every instant the velocity is zero, because those are the only places the vehicle can change direction, and the three pieces are then added with their signs taken off.

  1. Displacement is an antiderivative of velocity. Integrating term by term, s = ∫ (1.5n2 − 12n + 18)dn = 0.5n3 − 6n2 + 18n, taking s = 0 at the start of the maneuver.
  2. (a) Over the whole eight minutes the displacement is s(8) − s(0) = 256 − 384 + 144 = 16 meters, and it is positive, so the vehicle finishes 16 meters higher than it began. Its depth is 40 − 16 = 24 meters.
  3. The vehicle can only turn where its velocity is zero. Factorizing, v = 1.5(n − 2)(n − 6), so v = 0 at n = 2 and at n = 6, and between those two times v is negative: the vehicle is sinking.
  4. Take the three stretches separately. From 0 to 2, s(2) − s(0) = 4 − 24 + 36 = 16 meters up. From 2 to 6, s(6) − s(2) = 0 − 16 = −16, which is 16 meters down. From 6 to 8, s(8) − s(6) = 16 − 0 = 16 meters up.
  5. (b) The tether moves 16 + 16 + 16 = 48 meters in all. Check: putting the signs back gives 16 − 16 + 16 = 16 meters, the displacement of part (a), so the three pieces and the single integral agree.

answer(a) it is 24 meters below the surface, having risen 16 meters in all; (b) the winch moves 48 meters of tether

techniqueTotal Distance Traveled · Velocity and Displacement by Integration

examsO-Level

Common pitfalls

  • Integrating the velocity straight through for part (b) and answering 16 meters. That integral is the displacement: the four minutes of sinking cancel two of the rising stretches, and the winch does not care about the sign. A distance needs the interval cut at n = 2 and n = 6 and the three answers added without their signs.
  • Looking for the turning points where the velocity is least rather than where it is zero. The velocity is least at n = 4, but v = −6 there, so the vehicle is still sinking and has not turned. A change of direction happens only where v = 0.
02

A Foundry's Money In Against Its Money Out: The Profit Read as the Area Between Two Rate Curves

methodSubtract the Lower Rate From the Upper to Get the Gap, Solve for the Outputs at Which the Gap Is Zero to Fix the Limits, and Integrate the Gap Between Them for the Money the Area Stands For

A foundry can cast up to 900 fittings a day. At an output of x hundred fittings a day, money comes in at a rate of r = 30 + 6x − x2 thousand dollars for each extra hundred fittings, and goes out at a rate of c = 37 − 2x thousand dollars for each extra hundred. (a) Find the two outputs at which the two rates are equal. (b) Find the output that makes the foundry the most money, and how much more it makes there than at an output of 100 fittings a day.

01020304001479x, hundred fittings a daythousand dollars per hundredmoney inmoney outgap = money in − money out = −x2+ 8x − 7
The gap between the rates is r − c = (30 + 6x − x2) − (37 − 2x) = −x2 + 8x − 7.
The gap between the two rates is r − c = (30 + 6x − x2) − (37 − 2x) = −x2 + 8x − 7, in thousand dollars per hundred fittings.
step 1 of 5

Each curve is a rate, so the money each one stands for over a stretch of output is the area under it, and the profit gained over that stretch is the area between them. The gap between the curves is itself a function, and where that gap is zero the two rates are equal: those outputs are the limits of the region, and they are also the outputs where the profit stops rising and starts falling.

  1. The gap between the two rates is r − c = (30 + 6x − x2) − (37 − 2x) = −x2 + 8x − 7, in thousand dollars per hundred fittings.
  2. (a) The rates are equal where the gap is zero: x2 − 8x + 7 = 0, so (x − 1)(x − 7) = 0 and x = 1 or x = 7. The two outputs are 100 fittings a day and 700 fittings a day.
  3. Between those outputs the gap is positive, since at x = 4 it is −16 + 32 − 7 = 9. So every extra hundred fittings up to x = 7 adds money, and past x = 7 each one takes money away: the best output is 700 a day.
  4. (b) The money gained between the two outputs is the area between the curves, ∫17(−x2 + 8x − 7)dx = [−x33 + 4x2 − 7x]17.
  5. At x = 7 that is −3433 + 196 − 49 = 983, and at x = 1 it is −13 + 4 − 7 = −103. The difference is 1083 = 36 thousand dollars a day. Check: the gap is a parabola with roots 1 and 7 reaching 9 at the middle, and the mean height of such a parabola is two thirds of its greatest height, so the area is 6 × 6 = 36.

answer(a) the rates are equal at x = 1 and x = 7, that is at 100 and at 700 fittings a day; (b) 700 fittings a day, which makes 36 thousand dollars a day more than 100 fittings a day

techniqueArea Between Two Curves

examsH2

Common pitfalls

  • Integrating the gap from x = 0 instead of from x = 1. Below an output of 100 fittings the costs are the greater, so that stretch of the region lies below the axis and takes 103 thousand dollars off the total, leaving 983 instead of 36. The limits of a region between two curves are the outputs where they cross.
  • Multiplying the widest gap by the width of the region, 9 × 6 = 54. That is the area of a rectangle, and the gap is a curve which is 9 only at x = 4 and 0 at both ends. The integral is what adds the varying gap up.
03

A Vase Turned From Beech on a Lathe: Its Capacity by Discs and the Wood in Its Wall by Washers

methodTake a Slice at Right Angles to the Axis of Revolution: for the Cavity It Is a Disc of Area Pi Times the Inside Radius Squared, and for the Wall It Is a Washer Whose Area Is Pi Times the Difference of the Two Squared Radii, Then Integrate Each Up the Height

A vase is turned from beech on a lathe, so that it is a solid of revolution about a vertical axis. Measured upward from the inside of its base, at a height of y centimeters the inside radius is √y + 4 centimeters and the outside radius is √y + 9 centimeters, for 0 ≤ y ≤ 12. (a) Find the capacity of the vase. (b) Find the volume of wood in the wall between those two heights.

a discdisc area = pi × (y + 4)
A slice of the cavity at height y is a disc of radius √y + 4, so its area is π(y + 4).
A slice of the cavity at height y is a disc of radius √y + 4, so its area is π(√y + 4)2 = π(y + 4) square centimeters.
step 1 of 5

A solid of revolution is built from slices taken at right angles to its axis, and each slice is a circle whose radius the question gives. Where the solid is filled the slice is a disc; where it has been hollowed out the slice is a washer, a disc with a smaller disc removed, and its area is the difference of the two areas rather than the area of the difference of the radii. Integrating the slice area up the height adds the slices back into a volume.

  1. A slice of the cavity at height y is a disc of radius √y + 4, so its area is π(√y + 4)2 = π(y + 4) square centimeters.
  2. (a) The capacity is π∫012(y + 4)dy = π[y22 + 4y]012 = π(72 + 48) = 120π, which is 377 cubic centimeters to the nearest cubic centimeter.
  3. A slice of the wall at the same height is a washer: a disc of radius √y + 9 with a disc of radius √y + 4 taken out. Its area is π[(y + 9) − (y + 4)] = 5π square centimeters, the same at every height.
  4. (b) The wood in the wall is therefore π∫0125 dy = 5π × 12 = 60π, which is 188 cubic centimeters to the nearest cubic centimeter.
  5. Check: the block the vase was turned from is π∫012(y + 9)dy = π(72 + 108) = 180π, and the cavity and the wall add back to it, since 120π + 60π = 180π.

answer(a) the capacity is 120π, which is 377 cubic centimeters to the nearest cubic centimeter; (b) the wall holds 60π, which is 188 cubic centimeters of wood

techniqueVolumes of Revolution · Volumes by Washers

examsH2

Common pitfalls

  • Writing the washer's area as π(√y + 9 − √y + 4)2, the area of a circle whose radius is the wall's thickness. At the base that gives π(3 − 2)2 = π, against the true 5π: five times too small. A washer's area is the difference of the two areas, not the area built on the difference of the radii.
  • Integrating the radius rather than its square. ∫012√y + 4 dy has the units of an area, not a volume, and leaves out the π that makes a radius into a circle. Every slice of a solid of revolution is a circle, so its area carries π r2.
04

A Marquee Pitched Over a Round Dance Floor: The Air Under the Canvas From Its Ribs

methodWrite the Area of the Cross-Section as a Function of the Distance Along the Floor, Using the Shape the Question Gives the Ribs, and Integrate That Area Across the Whole Floor Instead of Multiplying One Cross-Section by a Length

A marquee is pitched over a circular dance floor of radius 3 meters. Its canvas is carried on parallel ribs set across the floor; each rib is an isosceles triangle standing on the floor, with both sides rising at 45°. Measure x meters from the center of the floor, at right angles to the ribs. (a) Find the volume of air under the canvas. (b) Ribs of the same design, running parallel to one side, are pitched over a square floor 6 meters by 6 meters instead. How much air does that hold, and how much more is it?

the sides rise at 45 deg, so the height is half the base
At a distance x from the center the floor is 2√9 − x2 meters wide, and the sides rise at 45°.
At a distance x from the center, the round floor is b = 2√9 − x2 meters wide, because the floor's edge is the circle x2 + y2 = 9.
step 1 of 5

The solid is not a prism, because its cross-section changes from rib to rib, so no single cross-section can simply be multiplied by a length. What is fixed is the shape of each rib, and the floor fixes how wide each rib has to be. Writing the area of a rib in terms of its position turns the solid into a stack of slabs, and integrating that area across the floor adds the slabs up.

  1. At a distance x from the center, the round floor is b = 2√9 − x2 meters wide, because the floor's edge is the circle x2 + y2 = 9.
  2. The rib there is an isosceles triangle on that base with both sides rising at 45°, so its height is half its base, √9 − x2. Its area is A = 12 × 2√9 − x2 × √9 − x2 = 9 − x2 square meters. At the center that is 9 square meters, a rib 6 meters wide and 3 meters high.
  3. (a) The air under the canvas is ∫−33(9 − x2)dx = [9x − x33]−33 = (27 − 9) − (−27 + 9) = 36 cubic meters.
  4. (b) Over a square floor 6 meters wide every rib has the same base of 6 meters, so A = 9 square meters throughout and the marquee is a prism: 9 × 6 = 54 cubic meters.
  5. That is 54 − 36 = 18 cubic meters more, a half as much again. Check: A = 9 − x2 is a parabola with roots at ± 3 reaching 9 in the middle, and the mean value of such a parabola is two thirds of its greatest value, so the mean cross-section is 6 square meters and the volume is 6 × 6 = 36.

answer(a) 36 cubic meters of air; (b) the square floor holds 54 cubic meters, which is 18 cubic meters more

techniqueVolumes with Known Cross Sections

Common pitfalls

  • Multiplying the largest rib by the width of the floor, 9 × 6 = 54 cubic meters. That is the answer for the square floor of part (b), where every rib really is the largest one. Over a round floor the ribs near the edges are both narrower and lower, and the integral is what accounts for them.
  • Taking the rib's height to be equal to its base rather than half of it. Sides rising at 45° meet above the middle of the base at half the base's width, so a 6 meter rib stands 3 meters high, not 6; using 6 doubles every slice and gives 72 cubic meters.
05

A Roller Coaster's First Drop: The Length of the Graded Track Against the Straight Line Under It

methodTake Pythagoras' Theorem on a Short Piece of the Curve, So That the Element of Length Is the Square Root of One Plus the Gradient Squared, and Integrate That Along the Run Rather Than Integrating the Gradient Itself

A roller coaster's first drop is graded so that, x meters from the crest, it falls with gradient √x96: level as it leaves the crest, and steepening all the way down. The graded track runs 54 meters from the crest. (a) Find the length of the graded surface. (b) Find how far the track drops over that run, and how much longer the graded surface is than a straight ramp between the same two points.

−27−18−900183654x, meters from the crestmeters fallena short piecethe gradient squared is x/96
A short piece of the track has sides dx and dydxdx, so its length is √1 + x96 dx.
A short piece of the surface has horizontal side dx and vertical side dydxdx, so by Pythagoras' theorem its length is √1 + (dydx)2 dx. Here (dydx)2 = x96, so the length is L = ∫054√1 + x96 dx.
step 1 of 5

A short piece of any curve is near enough a straight line, with a horizontal side and a vertical side, so Pythagoras' theorem gives its length as the horizontal step times the square root of one plus the gradient squared. Adding those pieces up is an integral along the run. The gradient on its own integrates to the drop, which is a different question, and the two are easy to confuse.

  1. A short piece of the surface has horizontal side dx and vertical side dydxdx, so by Pythagoras' theorem its length is √1 + (dydx)2 dx. Here (dydx)2 = x96, so the length is L = ∫054√1 + x96 dx.
  2. Integrate: ∫√1 + x96 dx = 96 × 23(1 + x96)32 = 64(1 + x96)32.
  3. (a) At x = 54, 1 + 5496 = 2516, whose three-halves power is 12564, so the bracket is 64 × 12564 = 125; at x = 0 it is 64. The graded surface is 125 − 64 = 61 meters long.
  4. (b) The drop is the integral of the gradient itself: ∫054√x96 dx = 23 × 5432√96 = 23 × 54 × √5496 = 23 × 54 × 34 = 27 meters.
  5. A straight ramp from the lip to the landing point is √542 + 272 = √3645 = 27√5 = 60.37 meters, so the graded surface is 61 − 60.37 = 0.63 meters longer. Check: the gradient nowhere exceeds 34, so the surface cannot be longer than 54 × √1 + 916 = 54 × 1.25 = 67.5 meters, nor shorter than the straight ramp, and 61 lies between.

answer(a) the graded surface is 61 meters long; (b) it drops 27 meters, and a straight ramp between the same two points is 27√5 = 60.37 meters, so the graded surface is only 0.63 meters longer

techniqueArc Length

Common pitfalls

  • Answering part (a) with the straight distance √542 + 272 = 60.37 meters. That is the chord from the lip to the landing point and runs under the hill the whole way; the graded surface bends away from it and is 61 meters.
  • Integrating the gradient instead of the element of length, ∫054√x96 dx = 27. That number is the drop, not the length: the 1 under the root is the horizontal side of the little triangle, and dropping it throws the horizontal run away.
06

A Grain Store's Hooped Roof: The Mean Height Across the Floor Against the Average of Five Measurements

methodTake the Mean Value of a Function as Its Integral Divided by the Width It Covers, and Set That Beside the Plain Average of Equally Spaced Readings, Which Gives the Two End Readings the Same Weight as the Ones in the Middle

A grain store has a hooped roof. Its height above the floor, x meters from the center line, is h = 8 − x28 meters, and the store is 16 meters wide, so the roof meets the floor at each side. (a) Find the mean height of the roof across the width of the store. (b) A worker measures the height every 4 meters across and averages the five readings. What does that average come to, and what do the same five readings give by the trapezium rule?

02468−8−4048x, meters from the center lineheight, mmean value = the area divided by the width
The mean value of a function is its integral divided by the width: h = 116∫−88(8 − x28)dx.
The mean value of a function over an interval is its integral divided by the width: h = 116∫−88(8 − x28)dx.
step 1 of 5

The mean value of a function over an interval is its integral divided by the width of that interval: it is the height of the rectangle with the same area. A handful of readings averaged is a different thing, because each reading is treated as though it stood for an equal share of the width, and the two readings at the ends stand for only half a share each. The trapezium rule is what puts that right without any more measuring.

  1. The mean value of a function over an interval is its integral divided by the width: h = 116∫−88(8 − x28)dx.
  2. Integrating, ∫−88(8 − x28)dx = [8x − x324]−88 = (64 − 643) − (−64 + 643) = 128 − 1283 = 2563, which is the area of the cross-section in square meters.
  3. (a) Dividing by the width, h = 2563 × 116 = 163 = 5.33 meters to two decimal places. That is two thirds of the 8 meter crown, which is what the mean of a parabola over the whole of its arch must be.
  4. (b) The five readings, at x = −8, −4, 0, 4 and 8, are 0, 6, 8, 6 and 0 meters, and their average is 205 = 4 meters. That is well short, because the two readings at the walls stand for only half a strip of floor each and yet carry a full share of the average.
  5. The trapezium rule gives each end reading half weight, which is the share it really stands for: the area is 4(02 + 6 + 8 + 6 + 02) = 4 × 20 = 80 square meters, so the mean height is 8016 = 5 meters. Check: 80 is short of the true 2563 = 85.33 only because each chord cuts inside a roof that curves.

answer(a) the mean height is 5.33 meters, that is 16/3 meters; (b) the average of the five readings is 4 meters, while the trapezium rule on the same readings gives an area of 80 square meters and a mean height of 5 meters

techniqueThe Mean Value of a Function · Numerical Integration

Common pitfalls

  • Averaging the five readings and calling it the mean height. That gives 4 meters, over a meter short, because it weights the two zero readings at the walls as heavily as the 8 meter crown. The mean value of a function shares the width out evenly along the floor, not evenly among the readings.
  • Taking the mean height as halfway between the lowest and the highest, 0 + 82 = 4 meters. Halving the two extremes is right for a straight line and wrong for anything that bends: the roof spends more of its width near the crown than a straight slope would.
07

Two Ion Thrusters Left Running Forever: The Total Impulse of Each

methodIntegrate to a Finite Upper Limit First, Read the Answer as a Function of That Limit, and Then Let the Limit Grow

A probe's ion thruster gives a thrust of F = 240e−0.002n millinewtons, where n is the number of seconds since it was lit, and the mission plans to leave it burning without end. The total impulse is the integral of the thrust over all time. (a) Find the total impulse, and the time by which nine tenths of it has been delivered. (b) A second design would give a thrust of F = 2401 + 0.002n millinewtons, which also falls away to nothing. Find its total impulse to a million seconds and to a hundred million seconds, and say whether it has a total impulse at all.

0100200011511500n, secondsthrust, mN060120180011511500T, secondstotal, thousand mN sto a time T the total is 120000(1 − e−0.002T)
Integrate to a finite time first: ∫0T240e−0.002ndn = 120000(1 − e−0.002T).
Integrate the first thrust to a finite time T: ∫0T240e−0.002ndn = [−120000e−0.002n]0T = 120000(1 − e−0.002T) millinewton-seconds.
step 1 of 6

An integral to infinity is not a sum of infinitely many things but a limit: integrate to a finite time, keep that time as a letter, and watch what the answer does as the time grows. Both thrusts here fall away to nothing, so neither is obviously hopeless, and the thing that decides the matter is how fast each one falls. An exponential falls fast enough that the tail adds almost nothing; a reciprocal does not.

  1. Integrate the first thrust to a finite time T: ∫0T240e−0.002ndn = [−120000e−0.002n]0T = 120000(1 − e−0.002T) millinewton-seconds.
  2. (a) As T grows, e−0.002T → 0, so the total settles at 120000 millinewton-seconds, that is 120 newton-seconds. At T of a million seconds and again at T of a hundred million the figure reads 120000 to every place a calculator shows: the tail beyond a million seconds is smaller than any of them.
  3. Nine tenths of that total has arrived when 1 − e−0.002n = 0.9, that is when e−0.002n = 0.1. Taking logarithms, 0.002n = ln 10, so n = 500ln 10 = 1151 seconds, about nineteen minutes.
  4. (b) Integrate the second thrust to a finite time in the same way: ∫0T2401 + 0.002ndn = [120000ln(1 + 0.002n)]0T = 120000ln(1 + 0.002T).
  5. At T = 106 that is 120000ln 2001 = 912000 millinewton-seconds, and at T = 108 it is 120000ln 200001 = 1465000. The answer has not settled; it has gone up by half as much again. A logarithm has no limit, so the second design has no total impulse: whatever total is asked for, it is reached eventually.
  6. Check: both thrusts fall to nothing, so falling away is not what decides it. At T = 1016 the second design would have delivered about 3680000 millinewton-seconds and would still be climbing.

answer(a) the total impulse settles at 120000 millinewton-seconds, that is 120 newton-seconds, and nine tenths of it has arrived by 1151 seconds; (b) the second design gives about 912000 millinewton-seconds by a million seconds and about 1465000 by a hundred million, so it has no total impulse at all

techniqueImproper Integrals

Common pitfalls

  • Writing e−∞ = 0 straight into the antiderivative without the limit. The working comes out right here, but the same short cut applied to the second design gives ln∞, which is not a number, and a reader who has never written the limit down has nothing to look at when that happens. The limit is what tells the two cases apart.
  • Arguing that the second thrust falls to zero, so its total must settle. Falling to zero is necessary and nowhere near enough: 2401 + 0.002n falls to zero and yet its total passes every figure you name. What matters is how fast the thrust falls, and a reciprocal falls far too slowly.
08

A River Gauged From Soundings: The Flow by Simpson's Rule and What the Trapezium Rule Leaves Out

methodWith Readings Instead of a Formula, Fit a Parabola Through Each Pair of Strips by Simpson's Rule, Weighting the Odd Soundings Four and the Even Interior Ones Two, and Compare It With the Trapezium Rule, Whose Straight Chords Cut Inside a Bed That Curves

A river is 24 meters wide at a gauging station. A survey boat sounds the depth every 3 meters from one bank, and reads 0, 1.75, 3.00, 3.75, 4.00, 3.75, 3.00, 1.75 and 0 meters. The water is moving at a mean speed of 0.5 meters per second. (a) Use Simpson's rule on the nine soundings to find the area of the cross-section, and hence the flow. (b) Work the same soundings by the trapezium rule, and say which answer is the smaller and why.

−4−3−2−1006121824meters from the bankdepth, meight strips, each 3 m wide
Eight strips of 3 meters span the width, and eight is even, so Simpson's rule applies to the nine soundings.
Eight strips of 3 meters span the 24 meter width, and eight is even, so Simpson's rule applies. It weights the two bank soundings 1, the odd-numbered soundings 4 and the interior even-numbered ones 2, and multiplies by a third of the strip width.
step 1 of 5

There is no formula for a river bed, only the soundings, so the area has to be estimated from them. The trapezium rule joins neighboring soundings by a straight chord; Simpson's rule takes the soundings two strips at a time and lays a parabola through the three, which follows a curving bed far more closely. The two rules use exactly the same nine numbers, so any difference between them is a difference of method, not of measurement.

  1. Eight strips of 3 meters span the 24 meter width, and eight is even, so Simpson's rule applies. It weights the two bank soundings 1, the odd-numbered soundings 4 and the interior even-numbered ones 2, and multiplies by a third of the strip width.
  2. (a) The area is 33[0 + 4(1.75 + 3.75 + 3.75 + 1.75) + 2(3.00 + 4.00 + 3.00) + 0] = 4 × 11 + 2 × 10 = 44 + 20 = 64 square meters.
  3. The flow is the area of the cross-section times the mean speed of the water: 64 × 0.5 = 32 cubic meters per second.
  4. (b) The trapezium rule gives every interior sounding weight 1 and halves the two at the banks: 3(02 + 1.75 + 3.00 + 3.75 + 4.00 + 3.75 + 3.00 + 1.75 + 02) = 3 × 21 = 63 square meters, a flow of 31.5 cubic meters per second.
  5. The trapezium rule is the smaller by 1 square meter. Between soundings the bed dips away below the straight chord drawn across it, so every trapezium leaves out a sliver, and the eight slivers come to that square meter. Check: Simpson's mean depth is 6424 = 2.67 meters, two thirds of the deepest sounding of 4.00 meters, which is what a bed shaped like a parabola gives, and Simpson's rule is exact on a parabola.

answer(a) the cross-section is 64 square meters and the flow is 32 cubic meters per second; (b) the trapezium rule gives 63 square meters and a flow of 31.5 cubic meters per second, 1 square meter the smaller, because every straight chord cuts inside a bed that dips away beneath it

techniqueSimpson's Rule · Numerical Integration

Common pitfalls

  • Giving the two bank soundings a weight of 4 or 2 along with the rest. The first and last readings are counted once each in both rules; here they are zero so nothing changes, but on a river with a steep bank the same slip adds several square meters. The pattern is 1, 4, 2, 4, 2, 4, 2, 4, 1.
  • Using Simpson's rule on an odd number of strips. The rule fits one parabola to each pair of strips, so the strips have to pair up; with seven strips one is left over and the rule cannot be applied to it. Either take one more sounding or use the trapezium rule on the odd strip.
09

A Water Butt Emptying Through Its Tap: A Sentence About a Rate Turned Into an Equation and Solved

methodTurn the Words Into a Differential Equation With a Constant of Proportionality and a Sign That Matches the Falling Depth, Separate the Variables and Integrate Both Sides, Then Use the Two Readings, One to Fix the Constant of Integration and One to Fix the Constant of Proportionality

A water butt with straight sides empties through a tap at its foot. Torricelli's law says that the depth of water falls at a rate proportional to the square root of the depth. The depth is 64 centimeters when the tap is opened, and 49 centimeters five minutes later. (a) Find the depth after 20 minutes. (b) Find when the butt is empty.

0164964052040n, minutesdepth h, cm0478052040n, minutesthe root of hdh/dn = −k√h, the minus because the depth falls
The depth falls at a rate proportional to √h, so dhdn = −k√h, the minus sign because it falls.
Let h be the depth in centimeters and n the number of minutes since the tap was opened. The depth falls at a rate proportional to √h, so dhdn = −k√h with k > 0, the minus sign because the depth is falling.
step 1 of 5

A sentence about a rate becomes a differential equation as soon as the rate is written as a derivative, the words after proportional to are written out, and a sign is chosen to match what is happening. This equation separates, because every h can be moved to one side and every n to the other. Two readings are then needed and both are given: one fixes the constant of integration, the other fixes the constant of proportionality.

  1. Let h be the depth in centimeters and n the number of minutes since the tap was opened. The depth falls at a rate proportional to √h, so dhdn = −k√h with k > 0, the minus sign because the depth is falling.
  2. Separate the variables and integrate: ∫ h−12dh = −k∫ dn, so 2√h = −kn + C, which is tidier written as √h = A − k2n.
  3. At n = 0 the depth is 64, so √64 = 8 = A. At n = 5 the depth is 49, so 7 = 8 − 5k2, giving k2 = 0.2. Hence √h = 8 − 0.2n, that is h = (8 − 0.2n)2.
  4. (a) At n = 20, √h = 8 − 4 = 4, so the depth is h = 42 = 16 centimeters.
  5. (b) The butt is empty when h = 0, that is when 8 − 0.2n = 0, so n = 40 minutes. Check: the rule gives (8 − 1)2 = 49 at n = 5, as it must, and the emptying slows down as it should, losing 15 centimeters in the first five minutes and only 1 centimeter in the last five.

answer(a) the depth is 16 centimeters; (b) the butt is empty 40 minutes after the tap was opened

techniqueSeparable Differential Equations · Forming Differential Equations

examsH2

Common pitfalls

  • Reading the sentence as ordinary decay, dhdn = −kh, and solving it as h = 64e−kn. The rate is proportional to the square root of the depth, not to the depth, and the difference is not a detail: an exponential never reaches zero, so that model says the butt never empties, while the real one is dry after 40 minutes.
  • Reading √h = 8 − 0.2n at n = 20 as h = 4. The 4 is the square root of the depth, so the depth is 42 = 16 centimeters. Squaring is the last step and it is easy to leave out when the number that comes out looks like an answer.
10

An Aeration Tank Whose Bacteria Saturate: A Concentration Stepped Out Where No Formula Can Be Read Off

methodWhen Separating the Variables Leaves an Equation That Cannot Be Rearranged for the Unknown, Draw the Slope Field and Step Along It by Euler's Method, Then Halve the Step and Compare the Two Answers to See Which Way the Method Errs and by How Much

In an aeration tank the bacteria take up ammonia at a rate that levels off while ammonia is plentiful and falls away as it runs short. With c milligrams of ammonia per liter after n hours, dcdn = −6c4 + c, and the tank starts at c = 8. (a) Use Euler's method with four steps of half an hour to estimate the concentration two hours later. (b) Do it again with eight steps of a quarter of an hour, and say which way Euler's method errs here and why.

0246800.511.52n, hoursammonia, mg per literseparating gives 4 ln c + c = K − 6n
Separating gives 4ln c + c = K − 6n, a true relation that cannot be rearranged to give c.
Separating gives ∫4 + ccdc = −∫ 6 dn, that is 4ln c + c = K − 6n. It is a true relation between c and n, but no rearrangement of it gives c on its own, so a value has to be stepped out instead of read off.
step 1 of 6

Separating the variables here gives a relation between c and n that is perfectly correct and perfectly useless for reading a value off, because it cannot be rearranged to give c. The slope field still shows what every solution does, and Euler's method walks along it: from the point you are at, follow the slope the equation gives for one step, and start again from wherever that lands you. Halving the step and comparing is how the size of the error is judged.

  1. Separating gives ∫4 + ccdc = −∫ 6 dn, that is 4ln c + c = K − 6n. It is a true relation between c and n, but no rearrangement of it gives c on its own, so a value has to be stepped out instead of read off.
  2. Euler's method steps along the slope field: cnext = c + h × f, where f = −6c4 + c is the slope the equation gives at the point reached so far.
  3. (a) With h = 0.5: at c = 8, f = −4812 = −4, so c → 8 − 2 = 6; at c = 6, f = −3610 = −3.6, so c → 6 − 1.8 = 4.2; at c = 4.2, f = −3.073, so c → 2.663; at c = 2.663, f = −2.398, so c → 1.464. After two hours the estimate is 1.46 milligrams per liter.
  4. (b) Eight steps of h = 0.25, worked the same way, give 1.667, so the estimate is 1.67 milligrams per liter.
  5. A step small enough for the answer to settle gives 1.852. Euler's method reads low both times, because it holds the slope of the start of each step for the whole of that step, and the slope eases as the ammonia falls: each step takes off more than it should. Halving the step roughly halves the error, from 0.39 to 0.19, which is what a first-order method does.
  6. Check: at n = 0 the relation gives K = 4ln 8 + 8 = 16.318, so at n = 2 it asks for 4ln c + c = 16.318 − 12 = 4.318. Putting c = 1.852 in gives 4ln 1.852 + 1.852 = 4.318, so the settled value is right and it is Euler's two estimates that are short.

answer(a) about 1.46 milligrams per liter; (b) about 1.67 milligrams per liter, and Euler's method reads low both times, since a step small enough to settle gives 1.852

techniqueEuler’s Method · Slope Fields

Common pitfalls

  • Using the slope at the start for every step. The slope is −4 only while c is 8; by the last step it is nearer −2.4. Holding −4 throughout gives 8 − 4 × 2 = 0 after two hours, which says the tank is clean when it is not.
  • Quoting Euler's answer as though it were exact. It is an estimate from a first-order method: halving the step only halves the error, so 1.46 and 1.67 are both short of 1.852, and neither is worth more than two significant figures at these step sizes.
11

A Car Driven Onto a Stretch of Undulating Road: The Body's Motion From the First Moment, and the Swing It Settles Into

methodSolve the Equation With Zero on the Right for the Complementary Function, Try a Cosine and a Sine at the Road's Frequency for the Particular Integral, Add the Two Before Fitting the Starting Conditions, and Let the Exponential in the Complementary Function Fall Away for the Long Run

One corner of a car carries 400 kilograms of its body on a spring of stiffness 16000 newtons per meter and a damper that resists with 1600 newtons for each meter per second of speed. The car has been driving on level road, with the body resting on the spring, when it reaches a stretch whose slabs have settled into long, regular undulations. From then on the road pushes up on that corner with an extra force of 300sin 5t newtons, where t is the number of seconds since the car reached the stretch. Let x be the height of the body in millimeters above where it rested. (a) Find x in terms of t. (b) Find the amplitude of the motion once it has settled down, and how many times as large that is as the height at which a steady upward push of 300 newtons would hold the body.

−300300123time, secondsheight x, mmat rest when the push beginsdivide by 400, in mm: 4, 40 and 750
Newton's second law, in millimeters and divided by 400: d2xdt2 + 4dxdt + 40x = 750sin 5t, starting from rest.
Let X be the height in meters. Newton's second law, with the spring and the damper both pulling back, is 400d2Xdt2 = 300sin 5t − 1600dXdt − 16000X. Putting x = 1000X millimeters and dividing every term by 400 gives d2xdt2 + 4dxdt + 40x = 750sin 5t, with x = 0 and dxdt = 0 at t = 0, since the body was resting.
step 1 of 5

The body on its spring and damper obeys a second-order linear equation, and the road's push is its right side. Such an equation is solved in two parts that are added: the complementary function makes the left side zero and carries the two constants, and the particular integral is one solution built to give the right side. The starting conditions can fix the constants only once the two parts are added, and in the long run the damper wears the complementary function away, leaving the particular integral as the motion the road keeps up.

  1. Let X be the height in meters. Newton's second law, with the spring and the damper both pulling back, is 400d2Xdt2 = 300sin 5t − 1600dXdt − 16000X. Putting x = 1000X millimeters and dividing every term by 400 gives d2xdt2 + 4dxdt + 40x = 750sin 5t, with x = 0 and dxdt = 0 at t = 0, since the body was resting.
  2. The complementary function solves the equation with 0 on the right. Its auxiliary equation is r2 + 4r + 40 = 0, so r = −4 ± √16 − 1602 = −2 ± 6i, and the complementary function is e−2t(Acos 6t + Bsin 6t): a swing of 6 radians a second inside a factor that shrinks by e−2 every second.
  3. For the particular integral, copy the shape of the right side with a cosine as well as a sine, because the damper turns one into the other: try x = pcos 5t + qsin 5t. Then d2xdt2 = −25x and dxdt = −5psin 5t + 5qcos 5t, so the left side is (15p + 20q)cos 5t + (15q − 20p)sin 5t. Matching it to 750sin 5t gives 15p + 20q = 0 and 15q − 20p = 750, so q = −34p and −1254p = 750, that is p = −24 and q = 18.
  4. (a) The general solution is the two parts added: x = e−2t(Acos 6t + Bsin 6t) − 24cos 5t + 18sin 5t. At t = 0, x = A − 24 = 0, so A = 24; and dxdt = −2A + 6B + 90 = 0, so 6B = 48 − 90 and B = −7. Hence x = e−2t(24cos 6t − 7sin 6t) − 24cos 5t + 18sin 5t millimeters.
  5. (b) By t = 3 the factor e−2t is below 0.003, so the complementary function has died away and the body moves as −24cos 5t + 18sin 5t, whose amplitude is √242 + 182 = √900 = 30 millimeters. A steady push of 300 newtons puts 750 on the right in place of 750sin 5t; its particular integral is the constant x = 75040 = 18.75 millimeters. So the undulations swing the body 3018.75 = 1.6 times as far as the same force held steady. Check: putting the particular integral back gives 15(−24) + 20(18) = 0 and 15(18) − 20(−24) = 750, as it must.

answer(a) x = e−2t(24cos 6t − 7sin 6t) − 24cos 5t + 18sin 5t millimeters; (b) the amplitude settles at 30 millimeters, 1.6 times the 18.75 millimeters at which a steady push of 300 newtons would hold the body

techniqueComplementary Function and Particular Integral · The Auxiliary Equation

Common pitfalls

  • Fitting the starting conditions to the complementary function before the particular integral is added. Applied to e−2t(Acos 6t + Bsin 6t) alone, x = 0 and dxdt = 0 make both constants 0, and the answer is the particular integral by itself, which starts at −24 millimeters rather than at rest. The conditions describe the whole motion, so they go on the whole solution.
  • Trying x = qsin 5t for the particular integral because the right side is a sine. The damping term differentiates the sine into a cosine, so the cosine terms must cancel as well, and with no cosine in the trial that asks for 20q = 0: there is then no particular integral at all. The trial needs both pcos 5t and qsin 5t.
Mr. Chalk Read the guide