Partial Derivatives · applications

Applications: Partial Derivatives

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

01

Air Under a Piston: How Fast the Pressure Changes With the Temperature and With the Volume, and the Volume That Holds It Steady

methodHold the Volume Constant and Differentiate in the Temperature, Hold the Temperature Constant and Differentiate in the Volume, Then Make the Two Changes in Pressure Cancel

A cylinder closed by a sliding piston holds a fixed amount of air. When the air is at a temperature of T kelvin and fills a volume of V liters, its pressure is P(T, V) = 25TV kilopascals. At present T = 300 and V = 75. (a) Find the pressure, and the partial derivatives ∂ P∂ T and ∂ P∂ V at the present state, and say what each one means. (b) The air is then warmed slowly while the piston moves to keep the pressure the same. Use the partial derivatives to find the rate, in liters per kelvin, at which the volume must grow, and the volume after the air has warmed by 6 kelvin.

96100104727578volume V, literspressure P, kPa300 K(a) P = 25 × 300 / 75 = 100 kPa
(a) At 300 kelvin and 75 liters the pressure is 100 kilopascals. The curve is the cross-section T = 300: the pressure against the volume at that temperature.
(a) At the present state, P(300, 75) = 25 × 30075 = 750075 = 100 kilopascals, close to the pressure of the air outside.
step 1 of 5

The pressure depends on two inputs, the temperature and the volume. A partial derivative measures how fast it changes with one input while the other is held constant: ∂ P∂ T is the rise per kelvin at a fixed volume, and ∂ P∂ V the change per liter at a fixed temperature. When both inputs change a little, the change in pressure is about the sum of the two effects, so holding the pressure steady means making the two effects cancel.

  1. (a) At the present state, P(300, 75) = 25 × 30075 = 750075 = 100 kilopascals, close to the pressure of the air outside.
  2. To find ∂ P∂ T, hold V constant and differentiate with respect to T: ∂ P∂ T = 25V = 2575 = 13. With the volume held at 75 liters, the pressure rises by about 13 of a kilopascal for each kelvin of warming.
  3. To find ∂ P∂ V, hold T constant, so that P = 25T × V−1: ∂ P∂ V = −25TV2 = −75005625 = −43. With the temperature held at 300 kelvin, the pressure falls by about 43 kilopascals for each extra liter.
  4. (b) Small changes Δ T and Δ V change the pressure by about 13Δ T − 43Δ V. For the pressure to stay the same this must be 0, so 43Δ V = 13Δ T and Δ V = 14Δ T. The volume must grow at 0.25 liters per kelvin.
  5. After 6 kelvin of warming, the volume must grow by 0.25 × 6 = 1.5 liters, to 76.5 liters. Check: P(306, 76.5) = 25 × 30676.5 = 765076.5 = 100 kilopascals, so the pressure is unchanged.

answer(a) 100 kilopascals; ∂ P∂ T = 13 ≈ 0.333 kilopascals per kelvin at a fixed volume, and ∂ P∂ V = −43 ≈ −1.333 kilopascals per liter at a fixed temperature; (b) 0.25 liters per kelvin, so 76.5 liters after 6 kelvin of warming

techniqueThe Partial Derivative · Computing Partial Derivatives

Common pitfalls

  • Writing ∂ P∂ V = −25V2 and dropping the T. Holding T constant makes 25T a constant factor, 7500 here, and a constant factor stays in the derivative.
  • Answering that the volume must shrink by 0.25 liters per kelvin. Warming raises the pressure, so the volume has to grow to bring it back down; it is the negative sign of ∂ P∂ V that lets the two changes cancel.
02

A Steel Rod Cooling With Its Ends in Ice: Checking the Heat Equation, and How Fast Its Middle Cools

methodDifferentiate Once in Time and Twice Along the Rod, Holding the Other Variable Constant Each Time, and Compare the Two Results

A steel rod 100 centimeters long has both ends held in ice at 0°C. The temperature at the point x centimeters from one end, m minutes after the rod is set up, is modeled by u(x, m) = 80e−0.01msinπ x100 degrees Celsius. (a) Find ∂ u∂ m and ∂2 u∂ x2, and find the constant k for which u satisfies the heat equation ∂ u∂ m = k∂2 u∂ x2. This k is the thermal diffusivity of the steel, in square centimeters per minute. (b) Find the rate at which the middle of the rod is cooling at the start, and the time at which it is cooling at half that rate. What is its temperature then?

4080255075100x, cm from one endu, deg Cm = 0(a) rate in m: −0.8 e−0.01msin(pi x/100)
(a) At the start the rod is 80°C at the middle and 0°C at the ends. Differentiating in m with x held constant gives the rate at which each point cools.
(a) Differentiate with respect to m, holding x constant, so that sinπ x100 is a constant factor: ∂ u∂ m = 80 × (−0.01)e−0.01msinπ x100 = −0.8e−0.01msinπ x100.
step 1 of 5

The temperature depends on the position along the rod and on the time, so it has a partial derivative with respect to each. The heat equation links the rate of change in time to the second partial derivative along the rod: where the temperature profile bends downward, a point is hotter than the average of the points on either side of it and loses heat to them. Checking the equation means differentiating once in m and twice in x, holding the other variable constant each time.

  1. (a) Differentiate with respect to m, holding x constant, so that sinπ x100 is a constant factor: ∂ u∂ m = 80 × (−0.01)e−0.01msinπ x100 = −0.8e−0.01msinπ x100.
  2. Differentiate twice with respect to x, holding m constant. The first derivative is ∂ u∂ x = 80 × π100e−0.01mcosπ x100, and the second is ∂2 u∂ x2 = −80 × π210000e−0.01msinπ x100 = −π2125e−0.01msinπ x100.
  3. Both derivatives are multiples of e−0.01msinπ x100, so the heat equation holds at every point and every time when −0.8 = −π2125k, that is k = 0.8 × 125π2 = 100π2 ≈ 10.1 square centimeters per minute. The second derivative is negative all along the rod, so every point is cooling.
  4. (b) At the middle, x = 50 and sinπ2 = 1, so ∂ u∂ m = −0.8e−0.01m. At the start, m = 0, the middle is cooling at 0.8°C per minute.
  5. It cools at half that rate, 0.4°C per minute, when e−0.01m = 0.5, so −0.01m = ln 0.5 and m = 100ln 2 ≈ 69.3 minutes. The temperature of the middle is then 80 × 0.5 = 40°C. Check: at the middle ∂ u∂ m = −0.01u, and 0.01 × 40 = 0.4°C per minute.

answer(a) ∂ u∂ m = −0.8e−0.01msinπ x100 and ∂2 u∂ x2 = −π2125e−0.01msinπ x100, so the heat equation holds with k = 100π2 ≈ 10.1 square centimeters per minute; (b) 0.8°C per minute at the start, and half that after 100ln 2 ≈ 69.3 minutes, when the middle is at 40°C

techniqueHigher Partial Derivatives · Computing Partial Derivatives

Common pitfalls

  • Differentiating sinπ x100 as cosπ x100 and dropping the factor π100 that the chain rule brings. It appears once in each of the two differentiations, so the second derivative carries π210000.
  • Assuming the middle keeps cooling at 0.8°C per minute, and so is cooling at half the rate after 50 minutes. The rate is 0.01 times the temperature, so it falls as the rod cools, and the rate halves only after 69.3 minutes.
03

An Online Shop's Price and Advertising: How Each Changes the Sales, and How Each Changes the Effect of the Other

methodDifferentiate Once in Each Variable, Then Differentiate Each Result in the Other Variable and Compare the Two Mixed Derivatives

An online shop sells a phone case at p dollars and spends a hundred dollars a week on advertising. Its weekly sales are modeled by S(p, a) = 20√a(40 − p) cases, for prices from $10 to $35 and advertising from $400 to $3600 a week. At present p = 20 and a = 16. (a) Find the weekly sales, and ∂ S∂ p and ∂ S∂ a at present, and say what each one means. (b) Find both mixed partial derivatives at present: ∂2 S∂ a   ∂ p, which differentiates first in p and then in a, and ∂2 S∂ p   ∂ a, which takes the other order. Say what each one tells the shop. Then use one of them to estimate ∂ S∂ p if the advertising rises to $2000 a week, and compare the estimate with the exact value.

1600320010203040price p, dollarscases sold a weeka = 16(a) S = 20 × 4 × 20 = 1600 cases a week
(a) At a price of $20 with $1600 of advertising a week, the shop sells 1600 cases. The line is the cross-section a = 16: sales against price.
(a) At present, S(20, 16) = 20 × √16 × (40 − 20) = 20 × 4 × 20 = 1600 cases a week.
step 1 of 5

A first partial derivative is the rate of change in one variable with the other held constant. Differentiating it again in the OTHER variable gives a mixed partial derivative, which measures how one rate changes as the other variable changes: here, how advertising changes the effect of the price, and how the price changes the effect of advertising. The mixed derivative theorem says the two orders give the same result wherever the second partial derivatives are continuous.

  1. (a) At present, S(20, 16) = 20 × √16 × (40 − 20) = 20 × 4 × 20 = 1600 cases a week.
  2. Holding a constant, ∂ S∂ p = −20√a = −80: each extra dollar on the price loses about 80 sales a week. Holding p constant, ∂ S∂ a = 20(40 − p) × 12√a = 10(40 − p)√a = 2004 = 50: each extra $100 of advertising brings about 50 more sales a week.
  3. (b) Differentiate ∂ S∂ p = −20√a with respect to a: ∂2 S∂ a   ∂ p = −10√a = −104 = −2.5. Each extra $100 of advertising makes the sales lost per dollar of price about 2.5 cases larger.
  4. Differentiate ∂ S∂ a = 10(40 − p)√a with respect to p: ∂2 S∂ p   ∂ a = −10√a = −2.5. Each extra dollar on the price makes $100 of advertising bring about 2.5 fewer sales. The two mixed derivatives are equal, as the mixed derivative theorem says, since every second partial derivative of S is continuous for a > 0.
  5. Raising the advertising from a = 16 to a = 20 changes ∂ S∂ p by about −2.5 × 4 = −10, to about −90 cases per dollar. Exactly, ∂ S∂ p = −20√20 ≈ −89.4, so the estimate is within 0.6 of it.

answer(a) 1600 cases a week; ∂ S∂ p = −80 cases per dollar and ∂ S∂ a = 50 cases per $100 of advertising; (b) both mixed derivatives are −2.5; the estimate of ∂ S∂ p at $2000 a week is −90 cases per dollar, against −20√20 ≈ −89.4 exactly

techniqueThe Mixed Derivative Theorem · Higher Partial Derivatives

Common pitfalls

  • Expecting the two mixed derivatives to differ because they describe different effects for the shop. The mixed derivative theorem makes them equal wherever the second partial derivatives are continuous, which holds here for every a > 0.
  • Differentiating √a as 1√a. The derivative of a12 is 12a−12 = 12√a, so ∂ S∂ a at present is 50 cases, not 100.
04

A Kite Rising and Drifting: How Fast the Flier Must Let Out Line, and Where the Kite Goes When the Line Is Held

methodWrite the Length as a Function of the Two Distances, Multiply Each Partial Derivative by the Rate of Its Own Input, and Add

A kite is 40 meters downwind of the flier's hand, measured horizontally, and 30 meters above it. Treat the line as straight, so that its length is L(x, y) = √x2 + y2 meters when the kite is x meters downwind and y meters up. (a) At this moment the kite is rising at 0.6 meters per second and moving downwind at 1.2 meters per second. Find the rate at which the flier must let out line. (b) Later, with the kite in the same position, the flier holds the line at a fixed length while the kite still rises at 0.6 meters per second. Find the rate at which the kite moves horizontally, and whether it moves toward the flier or away.

x = 40 my = 30 mL = 50 mhandkite(a) L2= 402+ 302, so L = 50 m
(a) The kite is 40 meters downwind and 30 meters up, so the straight line is 50 meters long.
(a) At this moment L = √402 + 302 = √1600 + 900 = 50 meters. The length depends on x and y, and both depend on the time, so the chain rule gives dLds = ∂ L∂ xdxds + ∂ L∂ ydyds.
step 1 of 5

The length of the line depends on the two distances, and both distances change with time. The chain rule for partial derivatives adds the two effects: each partial derivative of the length, multiplied by the rate at which its own input changes. With s seconds as the time, dLds = ∂ L∂ xdxds + ∂ L∂ ydyds. The same equation with dLds = 0 answers part (b).

  1. (a) At this moment L = √402 + 302 = √1600 + 900 = 50 meters. The length depends on x and y, and both depend on the time, so the chain rule gives dLds = ∂ L∂ xdxds + ∂ L∂ ydyds.
  2. Holding y constant, ∂ L∂ x = x√x2 + y2 = 4050 = 0.8; holding x constant, ∂ L∂ y = y√x2 + y2 = 3050 = 0.6.
  3. So dLds = 0.8 × 1.2 + 0.6 × 0.6 = 0.96 + 0.36 = 1.32. The flier must let out line at 1.32 meters per second.
  4. (b) Now dLds = 0 and dyds = 0.6, so 0.8dxds + 0.6 × 0.6 = 0, which gives dxds = −0.360.8 = −0.45.
  5. The negative sign means that x is decreasing: the kite moves toward the flier at 0.45 meters per second horizontally, climbing along a circle of radius 50 meters about the flier's hand. Check: the velocity (−0.45, 0.6) is at right angles to the line, since 40 × (−0.45) + 30 × 0.6 = −18 + 18 = 0.

answer(a) 1.32 meters per second; (b) 0.45 meters per second horizontally, toward the flier

techniqueThe Chain Rule for Partial Derivatives · Computing Partial Derivatives

Common pitfalls

  • Adding the two speeds, 1.2 + 0.6 = 1.8 meters per second. Each rate is weighted by its partial derivative: at this position, moving downwind lengthens the line by 0.8 meters for each meter moved, and rising lengthens it by only 0.6 meters.
  • Reading dxds = −0.45 as a speed away from the flier. x is the distance downwind, so a negative rate means the distance is shrinking and the kite is moving toward the flier.
05

A Probe on a Hotplate: the Direction of Fastest Warming, and the Path of a Robot That Always Takes It

methodCollect the Partial Derivatives Into the Gradient Vector, Take Its Direction and Its Length, Then Solve the Equation That Keeps the Path Along It

A steel hotplate 12 centimeters by 8 centimeters is heated from below at its center. With the origin at the center, x centimeters along the long side and y centimeters along the short side, the temperature of the plate is T(x, y) = 100 − 0.5x2 − y2 degrees Celsius. A temperature probe rests at (4, 1). (a) Find the temperature at the probe, the gradient vector ∇ T there, the direction in which the temperature rises fastest, and that fastest rate. (b) A small robot carrying the probe always moves in the direction in which the temperature rises fastest. Find the equation of its path and the point it heads to.

979185probe(a) T = 100 − 8 − 1 = 91 deg C
(a) The isotherms are ellipses around the center of the plate; the probe at (4, 1) is on the 91°C isotherm.
(a) At the probe, T(4, 1) = 100 − 0.5 × 16 − 1 = 100 − 8 − 1 = 91°C.
step 1 of 5

The gradient vector collects the two partial derivatives, ∇ T = (∂ T∂ x, ∂ T∂ y). It points in the direction in which T increases fastest, and its length is that fastest rate. A path that always follows the gradient has, at every point, the direction of the gradient vector, and that gives a differential equation for the path.

  1. (a) At the probe, T(4, 1) = 100 − 0.5 × 16 − 1 = 100 − 8 − 1 = 91°C.
  2. The partial derivatives are ∂ T∂ x = −x and ∂ T∂ y = −2y, so ∇ T = (−x, −2y), which is (−4, −2) at the probe.
  3. The temperature rises fastest in the direction of the gradient, (−4, −2), which as a unit vector is 1√5(−2, −1). The fastest rate is the length of the gradient, √16 + 4 = √20 ≈ 4.47°C per centimeter.
  4. (b) The path always points along ∇ T = (−x, −2y), so its gradient is dydx = −2y−x = 2yx. Separating the variables, dyy = 2 dxx, so ln y = 2ln x + c and y = Cx2 for a constant C.
  5. The path passes through (4, 1), so 1 = 16C and C = 116: the robot follows y = x216 toward (0, 0), the center of the plate, where the temperature is highest, 100°C. Check: at (4, 1) the path's gradient is 2 × 416 = 0.5, and the gradient vector (−4, −2) points along a line of gradient −2−4 = 0.5 too.

answer(a) 91°C; ∇ T = (−4, −2); the temperature rises fastest in the direction 1√5(−2, −1), at √20 ≈ 4.47°C per centimeter; (b) y = x216, heading to the center (0, 0), where the plate is at 100°C

techniqueThe Gradient Vector · Directional Derivatives · Separable Differential Equations

examsH2

Common pitfalls

  • Heading straight for the center along the line y = x4. That line reaches the hottest point, but it is not the direction of fastest rise: at (4, 1) the gradient vector points along a line of gradient 0.5, not 0.25, because the temperature changes twice as fast across the short side of the plate.
  • Stopping at the unit vector 1√5(−2, −1) and giving 1 as the fastest rate. The unit vector gives only the direction; the rate is the length of the gradient itself, √20.
06

A Skier Heading for the Lift Across a Gully: Her Rate of Descent as She Sets Off, and a Gentler Traverse

methodFind the Gradient, Divide the Displacement by Its Length for a Unit Vector, Take the Dot Product, and Solve for the Angle That Gives the Rate Required

On a ski slope, the height of the snow at the point x meters east and y meters north of a marker is h(x, y) = 1800 + 0.0004x2 − 0.32y meters above sea level, for −500 ≤ x ≤ 500 and −1000 ≤ y ≤ 2000: a broad gully that runs downhill to the north. A skier stands at (300, 0). (a) She sets off toward a lift station 800 meters west and 600 meters north of her. Find the rate at which her height changes, in meters per meter, as she sets off, and compare it with the steepest rate of descent at that point. (b) She would rather traverse, descending at 0.2 meters per meter. Find the angle between such a traverse and the fall line, the direction of steepest descent, and the two directions, as unit vectors, in which she could set off.

liftskier(a) gradient (0.24, −0.32), length 0.4
(a) The curves are contours 50 meters apart, and the ground falls toward the north, which is up the page. At the skier the gradient is ∇ h = (0.24, −0.32), of length 0.4.
(a) The partial derivatives are ∂ h∂ x = 0.0008x and ∂ h∂ y = −0.32, so at (300, 0) the gradient is ∇ h = (0.24, −0.32). Its length is √0.0576 + 0.1024 = √0.16 = 0.4.
step 1 of 5

The rate of change of height in a chosen direction is the directional derivative: the dot product of the gradient with a unit vector in that direction. It equals the length of the gradient times the cosine of the angle between the gradient and the direction, so the height falls fastest along the fall line, opposite to the gradient, and does not change at all along a contour.

  1. (a) The partial derivatives are ∂ h∂ x = 0.0008x and ∂ h∂ y = −0.32, so at (300, 0) the gradient is ∇ h = (0.24, −0.32). Its length is √0.0576 + 0.1024 = √0.16 = 0.4.
  2. The lift station is √8002 + 6002 = 1000 meters away, so the unit vector toward it is u = (−0.8, 0.6). The directional derivative is ∇ h · u = 0.24 × (−0.8) + (−0.32) × 0.6 = −0.192 − 0.192 = −0.384: she sets off descending 0.384 meters for each meter she travels.
  3. The steepest descent is along the fall line, −∇ h0.4 = (−0.6, 0.8), at the length of the gradient, 0.4 meters per meter. Her direction is close to the fall line: the angle θ between them has cosθ = 0.3840.4 = 0.96, so θ ≈ 16°.
  4. (b) Along a unit vector at an angle θ to the fall line, the height falls at 0.4cosθ meters per meter. Setting 0.4cosθ = 0.2 gives cosθ = 0.5, so θ = 60°, on either side of the fall line.
  5. Turning (−0.6, 0.8) through 60° each way, with cos 60° = 0.5 and sin 60° = √32, gives (−0.3 − 0.4√3, 0.4 − 0.3√3) ≈ (−0.99, −0.12) and (−0.3 + 0.4√3, 0.4 + 0.3√3) ≈ (0.39, 0.92). Check the second: 0.24 × 0.39 − 0.32 × 0.92 = 0.0936 − 0.2944 ≈ −0.2.

answer(a) she descends 0.384 meters per meter, just under the steepest descent at that point, 0.4 meters per meter; (b) 60° from the fall line, setting off along about (−0.99, −0.12) or (0.39, 0.92)

techniqueDirectional Derivatives · The Gradient Vector

Common pitfalls

  • Using the displacement (−800, 600) in place of a unit vector and getting a rate of −384. A directional derivative needs a unit vector; dividing the displacement by the distance, 1000 meters, gives the rate per meter.
  • Taking the fall line to be the direction of the gradient, (0.6, −0.8). The gradient points the way the height rises fastest, which is uphill; the fall line is the opposite direction.
07

A Tree Measured With a Tape and a Clinometer: an Estimate From the Tangent Plane, and the Error From Each Instrument

methodFind Both Partial Derivatives at the Readings, Write the Tangent Plane With the Angle in Radians, and Add the Sizes of the Two Contributions for the Largest Error

A forester stands on level ground d meters from the foot of a tree and sights the top from ground level at an angle of elevation θ, so that the tree is H(d, θ) = dtanθ meters tall. She measures d = 30 meters and θ = 45°. (a) Find the height, and the tangent plane approximation to H near these readings, with θ in radians. A second forester stands 31 meters from the tree and reads 44°. Use the tangent plane to estimate the height from these readings, and compare the estimate with the exact value. (b) The first forester's tape is accurate to within 0.2 meters and her clinometer to within 1°. Estimate the largest error in her height, and find how accurate the clinometer must be, with the same tape, for that error to be at most 0.5 meters.

45 degd = 30 mH(a) H = 30 × tan 45 deg = 30 m
(a) From 30 meters away the top is at 45°, so the tree is 30tan 45° = 30 meters tall.
(a) At the readings, H(30, 45°) = 30tan 45° = 30 × 1 = 30 meters.
step 1 of 5

Near the readings, H is close to its tangent plane: the height at the readings plus each partial derivative times the change in its own input. The same linear approximation turns small errors in the two measurements into an error in the height. The angle must be in radians, because the derivative of tanθ is sec2θ only when θ is measured in radians.

  1. (a) At the readings, H(30, 45°) = 30tan 45° = 30 × 1 = 30 meters.
  2. Holding θ constant, ∂ H∂ d = tanθ = 1. Holding d constant, ∂ H∂ θ = dsec2θ = 30 × 2 = 60 meters per radian. So the tangent plane is H ≈ 30 + 1 × (d − 30) + 60(θ − π4).
  3. A change of −1° is −π180 ≈ −0.01745 radians, so at d = 31 and θ = 44° the plane gives H ≈ 30 + 1 − 60 × 0.01745 = 31 − 1.047 = 29.95 meters. Exactly, 31tan 44° ≈ 29.94 meters, so the estimate is out by about 0.02 meters.
  4. (b) The largest error comes when both errors push the height the same way: |Δ H| ≈ 1 × 0.2 + 60 × π180 = 0.2 + 1.047 ≈ 1.25 meters. The clinometer causes more than five times as much of it as the tape.
  5. With the tape unchanged and the clinometer accurate to within δ radians, the largest error is about 0.2 + 60δ, and 0.2 + 60δ ≤ 0.5 needs δ ≤ 0.005 radians, which is 0.005 × 180π ≈ 0.29°. Check: 0.2 + 60 × 0.005 = 0.2 + 0.3 = 0.5 meters.

answer(a) 30 meters; H ≈ 30 + (d − 30) + 60(θ − π4); the estimate for 31 meters and 44° is 29.95 meters, against 29.94 meters exactly; (b) about 1.25 meters, most of it from the clinometer, which must be accurate to within about 0.29°

techniqueTangent Planes · Computing Partial Derivatives

Common pitfalls

  • Using the angle error in degrees, 60 × 1 = 60 meters. The partial derivative dsec2θ is in meters per radian, and 1° is only 0.01745 radians.
  • Letting the two errors cancel, as in 1.047 − 0.2. The two errors are independent and either may have either sign, so the largest possible error adds their sizes.
08

A Coffee Roaster Pricing Two Blends That Customers Switch Between: the Stationary Point of the Profit, and What Kind of Point It Is

methodDifferentiate the Profit in Each Price, Solve the Two Linear Equations Together, and Classify the Point by the Sign of D and of the Second Derivative

A coffee roaster sells two blends, House at p dollars a bag and Dark at q dollars a bag. Because some customers switch between them, each week it sells 120 − 20p + 10q bags of House and 240 + 10p − 20q bags of Dark, for prices near the present ones. A bag of House costs $4 to produce and a bag of Dark costs $6. (a) Write the weekly profit P(p, q) in dollars and find the prices at which both of its partial derivatives are zero. (b) Use the second-derivative test to classify that stationary point. Find the weekly profit there, and compare it with the profit at the present prices, $8 for House and $11 for Dark.

11131581012House price p, dollarsDark price q, dollarsnow(a) P = (p − 4)(120 − 20p + 10q)+ (q − 6)(240 + 10p − 20q)
(a) The curves join prices that give the same weekly profit; the present prices, $8 and $11, lie on one of them.
(a) Each blend earns its price minus its cost on every bag sold, so P(p, q) = (p − 4)(120 − 20p + 10q) + (q − 6)(240 + 10p − 20q) dollars a week.
step 1 of 5

The profit depends on both prices, and at a stationary point both partial derivatives are zero. Solving those two equations together gives the point. The second-derivative test classifies it with D = ∂2 P∂ p2∂2 P∂ q2 − (∂2 P∂ p   ∂ q)2: if D > 0 and ∂2 P∂ p2 < 0 the point is a maximum, if D > 0 and ∂2 P∂ p2 > 0 a minimum, and if D < 0 a saddle point.

  1. (a) Each blend earns its price minus its cost on every bag sold, so P(p, q) = (p − 4)(120 − 20p + 10q) + (q − 6)(240 + 10p − 20q) dollars a week.
  2. Differentiate with respect to p, holding q constant, by the product rule: ∂ P∂ p = (120 − 20p + 10q) − 20(p − 4) + 10(q − 6) = 140 − 40p + 20q. In the same way, ∂ P∂ q = 10(p − 4) + (240 + 10p − 20q) − 20(q − 6) = 320 + 20p − 40q.
  3. Setting both to zero and dividing by 20 gives 2p − q = 7 and −p + 2q = 16. Substituting q = 2p − 7 into the second gives −p + 4p − 14 = 16, so 3p = 30, p = 10 and q = 13: House at $10 and Dark at $13.
  4. (b) The second partial derivatives are ∂2 P∂ p2 = −40, ∂2 P∂ q2 = −40 and ∂2 P∂ p   ∂ q = 20, so D = (−40)(−40) − 202 = 1600 − 400 = 1200. Since D > 0 and ∂2 P∂ p2 < 0, the stationary point is a local maximum; P is quadratic and these values hold everywhere, so it is the greatest profit.
  5. At (10, 13) the roaster sells 120 − 200 + 130 = 50 bags of House and 240 + 100 − 260 = 80 bags of Dark, for a profit of 6 × 50 + 7 × 80 = 300 + 560 = $860 a week. At the present prices it sells 120 − 160 + 110 = 70 and 240 + 80 − 220 = 100 bags, for 4 × 70 + 5 × 100 = $780, so the new prices add $80 a week.

answer(a) P = (p − 4)(120 − 20p + 10q) + (q − 6)(240 + 10p − 20q), stationary at p = 10, q = 13; (b) D = 1200 > 0 with ∂2 P∂ p2 = −40 < 0, so it is a maximum: $860 a week, $80 more than the $780 at the present prices

techniqueStationary Points of a Surface · Higher Partial Derivatives

Common pitfalls

  • Setting the price of House as if it sold alone, from House's own profit only. A higher price for House sends customers to Dark, which adds 10(q − 6) to ∂ P∂ p; leaving it out gives the wrong prices.
  • Classifying the point from ∂2 P∂ p2 < 0 and ∂2 P∂ q2 < 0 alone. A saddle point can have both negative when the mixed derivative is large, so the test needs D > 0 as well.
09

A Beam Sawn From a Round Log: the Width and Depth That Resist Bending Best, Compared With a Square Beam

methodSet the Gradient of the Strength Equal to λ Times the Gradient of the Constraint, Eliminate λ, and Substitute Back Into the Circle

A sawmill cuts a beam of rectangular cross-section, w centimeters wide and d centimeters deep, from a round log 60 centimeters across, with the corners of the rectangle on the circle, so that w2 + d2 = 3600. The beam's strength against bending under a load from above is S = kwd2, for a constant k. (a) Use a Lagrange multiplier to find the width and depth of the strongest beam. (b) Find how much stronger it is than the square beam cut from the same log, as a percentage of the square beam's strength.

204060204060width w, cmdepth d, cm(a) at the best beam the two gradients are parallel
(a) The arc is the edge of the log, w2 + d2 = 3600, and the other curves join beams of equal strength. The strongest beam is where one of them just touches the arc, and there the two gradients are parallel.
(a) ∇ S = (kd2, 2kwd) and ∇ g = (2w, 2d), so ∇ S = λ∇ g says kd2 = 2λ w and 2kwd = 2λ d.
step 1 of 5

The strength is to be as large as possible while the width and depth stay tied to the log by g(w, d) = w2 + d2 = 3600. At the best point, the level curve of the strength just touches the circle, so the two gradients are parallel: ∇ S = λ∇ g for some number λ, the Lagrange multiplier. The two components of that equation and the constraint give the width and the depth.

  1. (a) ∇ S = (kd2, 2kwd) and ∇ g = (2w, 2d), so ∇ S = λ∇ g says kd2 = 2λ w and 2kwd = 2λ d.
  2. The depth is not 0, so dividing the second equation by 2d gives λ = kw. The first then becomes kd2 = 2kw2, that is d2 = 2w2.
  3. Substituting into the constraint: w2 + 2w2 = 3600, so w2 = 1200 and w = √1200 ≈ 34.6 centimeters, and d2 = 2400, so d = √2400 ≈ 49.0 centimeters. At the ends of the arc, where w = 0 or d = 0, the strength is 0, so this point is the maximum.
  4. (b) The strongest beam has S = k × √1200 × 2400 ≈ 83138k. The square beam has w = d, so 2w2 = 3600, w = √1800 ≈ 42.43 centimeters, and S = k × √1800 × 1800 ≈ 76368k.
  5. The ratio is 8313876368 ≈ 1.089, so the strongest beam is about 8.9 percent stronger than the square one. Check: w2 + d2 = 1200 + 2400 = 3600, so the corners of the strongest beam lie on the log.

answer(a) about 34.6 centimeters wide and 49.0 centimeters deep, with w2 = 1200 and d2 = 2400; (b) about 8.9 percent stronger than the square beam

techniqueLagrange Multipliers · The Gradient Vector

Common pitfalls

  • Maximizing wd2 with nothing linking w and d. Without the constraint the strength has no greatest value; the circle of the log is what ties the width to the depth.
  • Taking the square beam to be the strongest because it has the largest cross-section of any rectangle in the circle. The strength grows with the square of the depth, so it pays to give up some width for depth.
10

A Heap of Gravel on a Rectangular Patch: Its Volume, and the Volume at One End

methodIntegrate the Depth Along the Length for the Area of a Cross-Section, Then Integrate That Area Across the Width

A heap of gravel stands on a rectangular patch of ground 12 meters long and 8 meters wide. With the origin at the center of the patch, x meters along its length and y meters across it, the gravel is h(x, y) = 1.5(1 − x236)(1 − y216) meters deep, for −6 ≤ x ≤ 6 and −4 ≤ y ≤ 4. (a) Find the volume of the heap as a double integral over the rectangle. (b) A loader first clears the part of the heap beyond the line x = 3, the last 3 meters at one end. Find the volume it clears, and compare it with a quarter of the heap.

0.51.01.5 m(a) along the length: 1 − x2/36 integrates to 8
(a) Seen from above, the curves join points where the gravel is 0.5 and 1.0 meters deep. The straight line is a slice across the heap at one value of y; integrating the depth along it gives the area of that cross-section.
(a) The volume is ∫−44∫−66 h(x, y) dx dy. Integrate first in x, holding y constant, so that 1.5(1 − y216) is a constant factor: ∫−66(1 − x236)dx = [x − x3108]−66 = (6 − 2) − (−6 + 2) = 8.
step 1 of 5

The volume under a surface over a rectangle is a double integral, worked as two single integrals. Integrating the depth along the length, with y held constant, gives the area of the cross-section of the heap at that y; integrating that area across the width adds up the slices. The depth here is a function of x times a function of y, so each integral can be worked on its own.

  1. (a) The volume is ∫−44∫−66 h(x, y) dx dy. Integrate first in x, holding y constant, so that 1.5(1 − y216) is a constant factor: ∫−66(1 − x236)dx = [x − x3108]−66 = (6 − 2) − (−6 + 2) = 8.
  2. So the cross-section of the heap at y has area 1.5 × 8 × (1 − y216) = 12(1 − y216) square meters. Integrate it across the width: ∫−4412(1 − y216)dy = 12[y − y348]−44 = 12 × 163 = 64 cubic meters.
  3. Check: over its interval, 1 − x236 has an average of 812 = 23, and 1 − y216 has an average of 163 ÷ 8 = 23. So the heap fills 23 × 23 = 49 of the box around it, 12 by 8 by 1.5 meters: 49 × 144 = 64 cubic meters.
  4. (b) The part beyond x = 3 lies over the rectangle 3 ≤ x ≤ 6, −4 ≤ y ≤ 4. The integral across the width is unchanged, 163, and ∫36(1 − x236)dx = [x − x3108]36 = (6 − 2) − (3 − 0.25) = 1.25.
  5. So the loader clears 1.5 × 1.25 × 163 = 10 cubic meters. A quarter of the heap is 16 cubic meters, and the end holds only 1064 ≈ 15.6 percent of the heap, because the gravel is shallowest near the end.

answer(a) 64 cubic meters; (b) 10 cubic meters, less than a quarter of the heap, 16 cubic meters, because the gravel is shallow near the end

techniqueThe Volume Under a Surface

Common pitfalls

  • Multiplying the area of the patch by the greatest depth, 96 × 1.5 = 144 cubic meters. That is the volume of the box around the heap; the depth falls to 0 at every edge, and the double integral adds up the actual depth over each small piece of ground.
  • Taking the end beyond x = 3 to hold a quarter of the gravel because it is a quarter of the length. The gravel there is shallow, so the end holds 10 cubic meters, not 16.
Mr. Chalk Read the guide