MATH-283 (Calculus III) Curriculum

Session 13 - Optimization with Constraints

Handout Worksheet
Fall Semester 2026
Session Objectives & Overview

By the end of this session, you will be able to:

  • Formulate and locate critical points (βˆ‡f=0\nabla f = \mathbf{0}βˆ‡f=0 or undefined) for functions of two and three variables, construct the Hessian matrix and discriminant D=fxxfyyβˆ’(fxy)2D = f_{xx}f_{yy} - (f_{xy})^2D=fxx​fyyβ€‹βˆ’(fxy​)2, and rigorously classify candidates as local minima, local maxima, saddle points, or inconclusive.
  • Execute the 3-tier candidate search strategy (interior critical points, boundary edge parameterization, and corner vertices) to evaluate global extrema for continuous multivariable functions on closed and bounded planar regions under the Extreme Value Theorem.
  • Formulate and solve constrained multivariable optimization models using the method of Lagrange multipliers (βˆ‡f=Ξ»βˆ‡g\nabla f = \lambda \nabla gβˆ‡f=Ξ»βˆ‡g), resolve non-linear coordinate systems, and interpret the multiplier Ξ»\lambdaΞ» as a quantitative measure of physical system sensitivity and shadow pricing.

Micro-Lecture

1. Unconstrained Extrema & The Curvature Hessian

  • Mathematical Foundation: In multivariable modeling, stationary equilibrium states occur where the gradient vector vanishes: βˆ‡f=0\nabla f = \mathbf{0}βˆ‡f=0, creating a completely horizontal tangent plane z=f(a,b)z = f(a, b)z=f(a,b). As established in Critical Points in Multivariable Functions, stationary points represent candidates for local peaks, valleys, or directional inflection saddles. To classify these candidates, The Second Derivative Test inspects the determinant of the Hessian curvature matrix D=fxxfyyβˆ’(fxy)2D = f_{xx}f_{yy} - (f_{xy})^2D=fxx​fyyβ€‹βˆ’(fxy​)2. If D>0D > 0D>0, the surface curves uniformly in all directions (concave up for fxx>0β€…β€ŠβŸΉβ€…β€Šf_{xx} > 0 \impliesfxx​>0⟹ local min; concave down for fxx<0β€…β€ŠβŸΉβ€…β€Šf_{xx} < 0 \impliesfxx​<0⟹ local max). If D<0D < 0D<0, the surface bends upward in one cross-section and downward in a transverse cross-section, forming a saddle point.
  • Engineering Context: Aerospace engineers optimize airfoil lift-to-drag performance surfaces E(Ξ±,Ξ²)E(\alpha, \beta)E(Ξ±,Ξ²); a true local maximum (D>0,EΞ±Ξ±<0D > 0, E_{\alpha\alpha} < 0D>0,Eαα​<0) guarantees cruise stability, whereas a saddle point exposes the aircraft to violent aerodynamic pitching oscillations.

2. Compact Domains & The Multi-Tier Boundary Search

  • Mathematical Foundation: Engineering designs rarely operate across an infinite plane; physical systems are restricted to compact (closed and bounded) operating envelopes. Under the Extreme Value Theorem (explored in Extrema on Closed and Bounded Regions), a continuous multivariable function on a closed and bounded domain is guaranteed to attain an absolute maximum and minimum. Because global extrema frequently reside on perimeter boundaries or structural vertices rather than at interior stationary points, engineers must execute a systematic 3-tier search: interior critical points, 1D parameterized boundary edges, and geometric corner vertices.
  • Engineering Context: Geotechnical engineers evaluate contact soil bearing pressure P(x,y)P(x, y)P(x,y) under heavy foundation footings to guarantee that peak stresses along foundation edges or corner toes do not exceed the soil's ultimate bearing capacity, preventing structural punch-through failure.

3. Constrained Surfaces & The Method of Lagrange Multipliers

  • Mathematical Foundation: When optimization variables are strictly bound to an equation constraint g(x,y,z)=kg(x, y, z) = kg(x,y,z)=k (such as fixed material mass, trajectory pathways, or fluid flow cross-sections), interior gradients do not vanish. Instead, as proven in Lagrange Multipliers, an extreme value along the constraint path occurs when the constraint curve is tangent to a level set of the objective function. Because their tangent lines coincide, their normal gradient vectors are parallel: βˆ‡f=Ξ»βˆ‡g\nabla f = \lambda \nabla gβˆ‡f=Ξ»βˆ‡g. The scalar multiplier Ξ»\lambdaΞ» represents the marginal sensitivity (shadow price) of the system: Ξ»=dfoptdk\lambda = \frac{d f_{\text{opt}}}{d k}Ξ»=dkdfopt​​.
  • Engineering Context: Satellite mission analysts minimize orbit-raising propulsion fuel burn metrics subject to an orbital insertion velocity plane, where Ξ»\lambdaΞ» quantifies the exact propellant mass penalty per unit increase in required mission delta-V.

Skill Block 1

Worked Example A: Polynomial Gradient and Hessian Point Classification

Find and classify all critical points of the multivariable function:

f(x,y)=2x3βˆ’6x+2y2βˆ’4yf(x, y) = 2x^3 - 6x + 2y^2 - 4yf(x,y)=2x3βˆ’6x+2y2βˆ’4y

Worked Solution:
  1. Compute First-Order Partial Derivatives: State the Two-Variable Gradient Formula (defined in The Gradient Vector):

    βˆ‡f(x,y)=βŸ¨βˆ‚fβˆ‚x,βˆ‚fβˆ‚y⟩=βˆ‚fβˆ‚xi^+βˆ‚fβˆ‚yj^\nabla f(x, y) = \left\langle \frac{\partial f}{\partial x}, \frac{\partial f}{\partial y} \right\rangle = \frac{\partial f}{\partial x}\hat{i} + \frac{\partial f}{\partial y}\hat{j}βˆ‡f(x,y)=βŸ¨βˆ‚xβˆ‚f​,βˆ‚yβˆ‚fβ€‹βŸ©=βˆ‚xβˆ‚f​i^+βˆ‚yβˆ‚f​j^​
    Differentiate with respect to xxx treating yyy as a constant:
    fx(x,y)=βˆ‚βˆ‚x(2x3βˆ’6x+2y2βˆ’4y)=6x2βˆ’6f_x(x, y) = \frac{\partial}{\partial x}(2x^3 - 6x + 2y^2 - 4y) = 6x^2 - 6fx​(x,y)=βˆ‚xβˆ‚β€‹(2x3βˆ’6x+2y2βˆ’4y)=6x2βˆ’6
    Differentiate with respect to yyy treating xxx as a constant:
    fy(x,y)=βˆ‚βˆ‚y(2x3βˆ’6x+2y2βˆ’4y)=4yβˆ’4f_y(x, y) = \frac{\partial}{\partial y}(2x^3 - 6x + 2y^2 - 4y) = 4y - 4fy​(x,y)=βˆ‚yβˆ‚β€‹(2x3βˆ’6x+2y2βˆ’4y)=4yβˆ’4

  2. Locate Critical Points: State the Critical Point Condition (defined in Critical Points in Multivariable Functions):

    βˆ‡f(x,y)=0β€…β€ŠβŸΉβ€…β€Šfx(x,y)=0andfy(x,y)=0\nabla f(x, y) = \mathbf{0} \implies f_x(x, y) = 0 \quad \text{and} \quad f_y(x, y) = 0βˆ‡f(x,y)=0⟹fx​(x,y)=0andfy​(x,y)=0
    Solve the two decoupled equations:
    6x2βˆ’6=0β€…β€ŠβŸΉβ€…β€Š6(x2βˆ’1)=0β€…β€ŠβŸΉβ€…β€Šx=1orx=βˆ’16x^2 - 6 = 0 \implies 6(x^2 - 1) = 0 \implies x = 1 \quad \text{or} \quad x = -16x2βˆ’6=0⟹6(x2βˆ’1)=0⟹x=1orx=βˆ’1
    4yβˆ’4=0β€…β€ŠβŸΉβ€…β€Š4y=4β€…β€ŠβŸΉβ€…β€Šy=14y - 4 = 0 \implies 4y = 4 \implies y = 14yβˆ’4=0⟹4y=4⟹y=1
    Combining xxx and yyy solutions yields two distinct critical points:
    (1,1)and(βˆ’1,1)(1, 1) \quad \text{and} \quad (-1, 1)(1,1)and(βˆ’1,1)

  3. Compute Second-Order Partial Derivatives: State the Second-Order Partial Derivatives (defined in Second-Order Partial Derivatives):

    fxx=βˆ‚βˆ‚x(6x2βˆ’6)=12xf_{xx} = \frac{\partial}{\partial x}(6x^2 - 6) = 12xfxx​=βˆ‚xβˆ‚β€‹(6x2βˆ’6)=12x
    fyy=βˆ‚βˆ‚y(4yβˆ’4)=4f_{yy} = \frac{\partial}{\partial y}(4y - 4) = 4fyy​=βˆ‚yβˆ‚β€‹(4yβˆ’4)=4
    fxy=βˆ‚βˆ‚y(6x2βˆ’6)=0f_{xy} = \frac{\partial}{\partial y}(6x^2 - 6) = 0fxy​=βˆ‚yβˆ‚β€‹(6x2βˆ’6)=0

  4. Construct and Evaluate the Hessian Discriminant: State the Hessian Discriminant Formula (defined in Second Derivative Test for Multivariable Functions):

    D(x,y)=fxxfyyβˆ’(fxy)2=(12x)(4)βˆ’(0)2=48xD(x, y) = f_{xx}f_{yy} - (f_{xy})^2 = (12x)(4) - (0)^2 = 48xD(x,y)=fxx​fyyβ€‹βˆ’(fxy​)2=(12x)(4)βˆ’(0)2=48x

    • Evaluate at (1,1)(1, 1)(1,1):
      D(1,1)=48(1)=48>0D(1, 1) = 48(1) = 48 > 0D(1,1)=48(1)=48>0
      fxx(1,1)=12(1)=12>0f_{xx}(1, 1) = 12(1) = 12 > 0fxx​(1,1)=12(1)=12>0
      Because D>0D > 0D>0 and fxx>0f_{xx} > 0fxx​>0, the critical point (1,1)(1, 1)(1,1) is a local minimum.
      The local minimum value is:
      f(1,1)=2(1)3βˆ’6(1)+2(1)2βˆ’4(1)=2βˆ’6+2βˆ’4=βˆ’6f(1, 1) = 2(1)^3 - 6(1) + 2(1)^2 - 4(1) = 2 - 6 + 2 - 4 = -6f(1,1)=2(1)3βˆ’6(1)+2(1)2βˆ’4(1)=2βˆ’6+2βˆ’4=βˆ’6
    • Evaluate at (βˆ’1,1)(-1, 1)(βˆ’1,1):
      D(βˆ’1,1)=48(βˆ’1)=βˆ’48<0D(-1, 1) = 48(-1) = -48 < 0D(βˆ’1,1)=48(βˆ’1)=βˆ’48<0
      Because D<0D < 0D<0, the critical point (βˆ’1,1)(-1, 1)(βˆ’1,1) is a saddle point.
      The saddle point value is:
      f(βˆ’1,1)=2(βˆ’1)3βˆ’6(βˆ’1)+2(1)2βˆ’4(1)=βˆ’2+6+2βˆ’4=2f(-1, 1) = 2(-1)^3 - 6(-1) + 2(1)^2 - 4(1) = -2 + 6 + 2 - 4 = 2f(βˆ’1,1)=2(βˆ’1)3βˆ’6(βˆ’1)+2(1)2βˆ’4(1)=βˆ’2+6+2βˆ’4=2

Practice Problem 1: Concept & Reading Check

Open this Skill Block on the Memex to locate the information needed to answer the following:

  • (a) State the algebraic condition on the gradient vector βˆ‡f\nabla fβˆ‡f required for a point (a,b)(a, b)(a,b) to be classified as a critical point, and explain what this means for the surface's tangent plane.
  • (b) In Worked Example A, explain why knowing D(1,1)>0D(1, 1) > 0D(1,1)>0 alone is insufficient to conclude that (1,1)(1, 1)(1,1) is a local minimum, and state what additional derivative sign must be verified.
  • (c) In Aerospace Wing Profile Lift-to-Drag Optimization (Aerospace & Mechanical Engineering), explain why an aircraft engineer seeks an operating point with D>0D > 0D>0 and Exx<0E_{xx} < 0Exx​<0 rather than a saddle point where D<0D < 0D<0.
Solution Key:
  • (a) By the Critical Point Definition (defined in Critical Points in Multivariable Functions), all first-order partial derivatives must vanish simultaneously: βˆ‡f(a,b)=0β€…β€ŠβŸΊβ€…β€Šfx(a,b)=0\nabla f(a, b) = \mathbf{0} \iff f_x(a, b) = 0βˆ‡f(a,b)=0⟺fx​(a,b)=0 and fy(a,b)=0f_y(a, b) = 0fy​(a,b)=0. Geometrically, the tangent plane equation zβˆ’f(a,b)=fx(xβˆ’a)+fy(yβˆ’b)z - f(a, b) = f_x(x - a) + f_y(y - b)zβˆ’f(a,b)=fx​(xβˆ’a)+fy​(yβˆ’b) reduces to the constant horizontal plane z=f(a,b)z = f(a, b)z=f(a,b).
  • (b) By the Hessian Discriminant Criteria (defined in Second Derivative Test for Multivariable Functions), D>0D > 0D>0 guarantees that the principal curvatures share the same sign. To distinguish whether the surface forms an upright bowl (local minimum) or an inverted dome (local maximum), one must evaluate the sign of the pure second derivative fxxf_{xx}fxx​ (or fyyf_{yy}fyy​): fxx>0f_{xx} > 0fxx​>0 confirms upward concavity (local minimum).
  • (c) In aerospace flight mechanics, an aerodynamic efficiency peak requires negative concavity in all directions (D>0,Exx<0D > 0, E_{xx} < 0D>0,Exx​<0) so that deviations in angle of attack or camber produce graceful performance boundaries. A saddle point (D<0D < 0D<0) represents an unstable equilibrium where slight perturbations induce catastrophic loss of lift.
Show your steps and calculations here.

Practice Problem 2: Decoupled Polynomial Local Extrema

Find and classify all critical points of the function:

f(x,y)=x3βˆ’3x+y2f(x, y) = x^3 - 3x + y^2f(x,y)=x3βˆ’3x+y2

Solution Key:
  1. Compute First Partial Derivatives: State the Two-Variable Gradient Formula:

    βˆ‡f(x,y)=βŸ¨βˆ‚fβˆ‚x,βˆ‚fβˆ‚y⟩\nabla f(x, y) = \left\langle \frac{\partial f}{\partial x}, \frac{\partial f}{\partial y} \right\rangleβˆ‡f(x,y)=βŸ¨βˆ‚xβˆ‚f​,βˆ‚yβˆ‚fβ€‹βŸ©
    fx=βˆ‚βˆ‚x(x3βˆ’3x+y2)=3x2βˆ’3f_x = \frac{\partial}{\partial x}(x^3 - 3x + y^2) = 3x^2 - 3fx​=βˆ‚xβˆ‚β€‹(x3βˆ’3x+y2)=3x2βˆ’3
    fy=βˆ‚βˆ‚y(x3βˆ’3x+y2)=2yf_y = \frac{\partial}{\partial y}(x^3 - 3x + y^2) = 2yfy​=βˆ‚yβˆ‚β€‹(x3βˆ’3x+y2)=2y

  2. Locate Critical Points: Set both partial derivatives equal to zero:

    3x2βˆ’3=0β€…β€ŠβŸΉβ€…β€Š3(x2βˆ’1)=0β€…β€ŠβŸΉβ€…β€Šx=1orx=βˆ’13x^2 - 3 = 0 \implies 3(x^2 - 1) = 0 \implies x = 1 \quad \text{or} \quad x = -13x2βˆ’3=0⟹3(x2βˆ’1)=0⟹x=1orx=βˆ’1
    2y=0β€…β€ŠβŸΉβ€…β€Šy=02y = 0 \implies y = 02y=0⟹y=0
    This yields two critical points: (1,0)(1, 0)(1,0) and (βˆ’1,0)(-1, 0)(βˆ’1,0).

  3. Compute Second-Order Partials and Discriminant: State the Hessian Discriminant Formula:

    D(x,y)=fxxfyyβˆ’(fxy)2D(x, y) = f_{xx}f_{yy} - (f_{xy})^2D(x,y)=fxx​fyyβ€‹βˆ’(fxy​)2
    fxx=6x,fyy=2,fxy=0f_{xx} = 6x, \quad f_{yy} = 2, \quad f_{xy} = 0fxx​=6x,fyy​=2,fxy​=0
    D(x,y)=(6x)(2)βˆ’02=12xD(x, y) = (6x)(2) - 0^2 = 12xD(x,y)=(6x)(2)βˆ’02=12x

  4. Classify Points:

    • At (1,0)(1, 0)(1,0):
      D(1,0)=12(1)=12>0,fxx(1,0)=6(1)=6>0D(1, 0) = 12(1) = 12 > 0, \quad f_{xx}(1, 0) = 6(1) = 6 > 0D(1,0)=12(1)=12>0,fxx​(1,0)=6(1)=6>0
      Since D>0D > 0D>0 and fxx>0f_{xx} > 0fxx​>0, (1,0)(1, 0)(1,0) is a local minimum with value f(1,0)=13βˆ’3(1)+0=βˆ’2f(1, 0) = 1^3 - 3(1) + 0 = -2f(1,0)=13βˆ’3(1)+0=βˆ’2.
    • At (βˆ’1,0)(-1, 0)(βˆ’1,0):
      D(βˆ’1,0)=12(βˆ’1)=βˆ’12<0D(-1, 0) = 12(-1) = -12 < 0D(βˆ’1,0)=12(βˆ’1)=βˆ’12<0
      Since D<0D < 0D<0, (βˆ’1,0)(-1, 0)(βˆ’1,0) is a saddle point with value f(βˆ’1,0)=(βˆ’1)3βˆ’3(βˆ’1)+0=2f(-1, 0) = (-1)^3 - 3(-1) + 0 = 2f(βˆ’1,0)=(βˆ’1)3βˆ’3(βˆ’1)+0=2.
Show your steps and calculations here.

Practice Problem 3: Coupled Polynomial Extrema with Cross-Terms

Determine and classify all critical points of the function:

f(x,y)=y3βˆ’3xy+12xf(x, y) = y^3 - 3xy + 12xf(x,y)=y3βˆ’3xy+12x

Solution Key:
  1. Compute First Partial Derivatives: State the Two-Variable Gradient Formula:

    βˆ‡f(x,y)=βŸ¨βˆ‚fβˆ‚x,βˆ‚fβˆ‚y⟩\nabla f(x, y) = \left\langle \frac{\partial f}{\partial x}, \frac{\partial f}{\partial y} \right\rangleβˆ‡f(x,y)=βŸ¨βˆ‚xβˆ‚f​,βˆ‚yβˆ‚fβ€‹βŸ©
    fx=βˆ‚βˆ‚x(y3βˆ’3xy+12x)=βˆ’3y+12f_x = \frac{\partial}{\partial x}(y^3 - 3xy + 12x) = -3y + 12fx​=βˆ‚xβˆ‚β€‹(y3βˆ’3xy+12x)=βˆ’3y+12
    fy=βˆ‚βˆ‚y(y3βˆ’3xy+12x)=3y2βˆ’3xf_y = \frac{\partial}{\partial y}(y^3 - 3xy + 12x) = 3y^2 - 3xfy​=βˆ‚yβˆ‚β€‹(y3βˆ’3xy+12x)=3y2βˆ’3x

  2. Locate Critical Points: Set both equations to zero:

    βˆ’3y+12=0β€…β€ŠβŸΉβ€…β€Š3y=12β€…β€ŠβŸΉβ€…β€Šy=4-3y + 12 = 0 \implies 3y = 12 \implies y = 4βˆ’3y+12=0⟹3y=12⟹y=4
    Substitute y=4y = 4y=4 into fy=0f_y = 0fy​=0:
    3(4)2βˆ’3x=0β€…β€ŠβŸΉβ€…β€Š3(16)βˆ’3x=0β€…β€ŠβŸΉβ€…β€Š48=3xβ€…β€ŠβŸΉβ€…β€Šx=163(4)^2 - 3x = 0 \implies 3(16) - 3x = 0 \implies 48 = 3x \implies x = 163(4)2βˆ’3x=0⟹3(16)βˆ’3x=0⟹48=3x⟹x=16
    The single critical point is (16,4)(16, 4)(16,4).

  3. Compute Second-Order Partials and Discriminant: State the Hessian Discriminant Formula:

    fxx=βˆ‚βˆ‚x(βˆ’3y+12)=0f_{xx} = \frac{\partial}{\partial x}(-3y + 12) = 0fxx​=βˆ‚xβˆ‚β€‹(βˆ’3y+12)=0
    fyy=βˆ‚βˆ‚y(3y2βˆ’3x)=6yf_{yy} = \frac{\partial}{\partial y}(3y^2 - 3x) = 6yfyy​=βˆ‚yβˆ‚β€‹(3y2βˆ’3x)=6y
    fxy=βˆ‚βˆ‚y(βˆ’3y+12)=βˆ’3f_{xy} = \frac{\partial}{\partial y}(-3y + 12) = -3fxy​=βˆ‚yβˆ‚β€‹(βˆ’3y+12)=βˆ’3
    D(x,y)=fxxfyyβˆ’(fxy)2=(0)(6y)βˆ’(βˆ’3)2=0βˆ’9=βˆ’9D(x, y) = f_{xx}f_{yy} - (f_{xy})^2 = (0)(6y) - (-3)^2 = 0 - 9 = -9D(x,y)=fxx​fyyβ€‹βˆ’(fxy​)2=(0)(6y)βˆ’(βˆ’3)2=0βˆ’9=βˆ’9

  4. Classify the Point: At (16,4)(16, 4)(16,4), D=βˆ’9<0D = -9 < 0D=βˆ’9<0.
    Since D<0D < 0D<0, the critical point (16,4)(16, 4)(16,4) is a saddle point.
    Its value is f(16,4)=43βˆ’3(16)(4)+12(16)=64βˆ’192+192=64f(16, 4) = 4^3 - 3(16)(4) + 12(16) = 64 - 192 + 192 = 64f(16,4)=43βˆ’3(16)(4)+12(16)=64βˆ’192+192=64.

Show your steps and calculations here.

Practice Problem 4: Non-Linear Coordinate Elimination and Local Maxima

Find and classify all critical points of the higher-degree multivariable function:

f(x,y)=4xyβˆ’x4βˆ’y4f(x, y) = 4xy - x^4 - y^4f(x,y)=4xyβˆ’x4βˆ’y4

Solution Key:
  1. Compute First-Order Partials:

    fx=4yβˆ’4x3,fy=4xβˆ’4y3f_x = 4y - 4x^3, \quad f_y = 4x - 4y^3fx​=4yβˆ’4x3,fy​=4xβˆ’4y3

  2. Solve the Non-Linear System: Set fx=0β€…β€ŠβŸΉβ€…β€Š4y=4x3β€…β€ŠβŸΉβ€…β€Šy=x3f_x = 0 \implies 4y = 4x^3 \implies y = x^3fx​=0⟹4y=4x3⟹y=x3. Substitute y=x3y = x^3y=x3 into fy=0f_y = 0fy​=0:

    4xβˆ’4(x3)3=0β€…β€ŠβŸΉβ€…β€Š4xβˆ’4x9=0β€…β€ŠβŸΉβ€…β€Š4x(1βˆ’x8)=04x - 4(x^3)^3 = 0 \implies 4x - 4x^9 = 0 \implies 4x(1 - x^8) = 04xβˆ’4(x3)3=0⟹4xβˆ’4x9=0⟹4x(1βˆ’x8)=0
    The real solutions for xxx are x=0x = 0x=0, x=1x = 1x=1, and x=βˆ’1x = -1x=βˆ’1.

    • If x=0β€…β€ŠβŸΉβ€…β€Šy=03=0β€…β€ŠβŸΉβ€…β€Š(0,0)x = 0 \implies y = 0^3 = 0 \implies (0, 0)x=0⟹y=03=0⟹(0,0)
    • If x=1β€…β€ŠβŸΉβ€…β€Šy=13=1β€…β€ŠβŸΉβ€…β€Š(1,1)x = 1 \implies y = 1^3 = 1 \implies (1, 1)x=1⟹y=13=1⟹(1,1)
    • If x=βˆ’1β€…β€ŠβŸΉβ€…β€Šy=(βˆ’1)3=βˆ’1β€…β€ŠβŸΉβ€…β€Š(βˆ’1,βˆ’1)x = -1 \implies y = (-1)^3 = -1 \implies (-1, -1)x=βˆ’1⟹y=(βˆ’1)3=βˆ’1⟹(βˆ’1,βˆ’1)
  3. Compute Second Partials and Discriminant:

    fxx=βˆ’12x2,fyy=βˆ’12y2,fxy=4f_{xx} = -12x^2, \quad f_{yy} = -12y^2, \quad f_{xy} = 4fxx​=βˆ’12x2,fyy​=βˆ’12y2,fxy​=4
    State the Hessian Discriminant Formula:
    D(x,y)=fxxfyyβˆ’(fxy)2=(βˆ’12x2)(βˆ’12y2)βˆ’(4)2=144x2y2βˆ’16D(x, y) = f_{xx}f_{yy} - (f_{xy})^2 = (-12x^2)(-12y^2) - (4)^2 = 144x^2y^2 - 16D(x,y)=fxx​fyyβ€‹βˆ’(fxy​)2=(βˆ’12x2)(βˆ’12y2)βˆ’(4)2=144x2y2βˆ’16

  4. Classify Each Critical Point:

    • At (0,0)(0, 0)(0,0):
      D(0,0)=144(0)βˆ’16=βˆ’16<0β€…β€ŠβŸΉβ€…β€ŠSaddleΒ pointD(0, 0) = 144(0) - 16 = -16 < 0 \implies \textbf{Saddle point}D(0,0)=144(0)βˆ’16=βˆ’16<0⟹SaddleΒ point
      Value: f(0,0)=0f(0, 0) = 0f(0,0)=0.
    • At (1,1)(1, 1)(1,1):
      D(1,1)=144(1)2(1)2βˆ’16=128>0,fxx(1,1)=βˆ’12(1)2=βˆ’12<0β€…β€ŠβŸΉβ€…β€ŠLocalΒ maximumD(1, 1) = 144(1)^2(1)^2 - 16 = 128 > 0, \quad f_{xx}(1, 1) = -12(1)^2 = -12 < 0 \implies \textbf{Local maximum}D(1,1)=144(1)2(1)2βˆ’16=128>0,fxx​(1,1)=βˆ’12(1)2=βˆ’12<0⟹LocalΒ maximum
      Value: f(1,1)=4(1)(1)βˆ’14βˆ’14=2f(1, 1) = 4(1)(1) - 1^4 - 1^4 = 2f(1,1)=4(1)(1)βˆ’14βˆ’14=2.
    • At (βˆ’1,βˆ’1)(-1, -1)(βˆ’1,βˆ’1):
      D(βˆ’1,βˆ’1)=144(βˆ’1)2(βˆ’1)2βˆ’16=128>0,fxx(βˆ’1,βˆ’1)=βˆ’12(βˆ’1)2=βˆ’12<0β€…β€ŠβŸΉβ€…β€ŠLocalΒ maximumD(-1, -1) = 144(-1)^2(-1)^2 - 16 = 128 > 0, \quad f_{xx}(-1, -1) = -12(-1)^2 = -12 < 0 \implies \textbf{Local maximum}D(βˆ’1,βˆ’1)=144(βˆ’1)2(βˆ’1)2βˆ’16=128>0,fxx​(βˆ’1,βˆ’1)=βˆ’12(βˆ’1)2=βˆ’12<0⟹LocalΒ maximum
      Value: f(βˆ’1,βˆ’1)=4(βˆ’1)(βˆ’1)βˆ’(βˆ’1)4βˆ’(βˆ’1)4=2f(-1, -1) = 4(-1)(-1) - (-1)^4 - (-1)^4 = 2f(βˆ’1,βˆ’1)=4(βˆ’1)(βˆ’1)βˆ’(βˆ’1)4βˆ’(βˆ’1)4=2.

Show your steps and calculations here.

Regroup 1

  • Visualizing Curvature: Emphasize why D<0D < 0D<0 guarantees a saddle point: the quadratic expansion has eigenvalues of opposite signs, meaning the surface climbs in one direction and falls in another.
  • Algebraic Pitfall: Highlight that when solving coupled systems like 4yβˆ’4x3=04y - 4x^3 = 04yβˆ’4x3=0, students must never divide by xxx or yyy without accounting for the origin (0,0)(0, 0)(0,0) as a solution.

Skill Block 2

Worked Example A: Absolute Extrema on a Closed Rectangular Plate

Find the absolute maximum and minimum values of:

f(x,y)=x2+y2βˆ’2xβˆ’4y+5f(x, y) = x^2 + y^2 - 2x - 4y + 5f(x,y)=x2+y2βˆ’2xβˆ’4y+5
on the closed square domain D={(x,y)∈R2∣0≀x≀4, 0≀y≀4}D = \{ (x, y) \in \mathbb{R}^2 \mid 0 \le x \le 4, \, 0 \le y \le 4 \}D={(x,y)∈R2∣0≀x≀4,0≀y≀4}.

Worked Solution:
  1. Tier 1: Interior Critical Points: Compute the first partial derivatives using the Two-Variable Gradient Formula:

    fx=2xβˆ’2=0β€…β€ŠβŸΉβ€…β€Šx=1f_x = 2x - 2 = 0 \implies x = 1fx​=2xβˆ’2=0⟹x=1
    fy=2yβˆ’4=0β€…β€ŠβŸΉβ€…β€Šy=2f_y = 2y - 4 = 0 \implies y = 2fy​=2yβˆ’4=0⟹y=2
    The point (1,2)(1, 2)(1,2) lies strictly inside DDD because 0<1<40 < 1 < 40<1<4 and 0<2<40 < 2 < 40<2<4.
    Evaluate fff at (1,2)(1, 2)(1,2):
    f(1,2)=12+22βˆ’2(1)βˆ’4(2)+5=1+4βˆ’2βˆ’8+5=0f(1, 2) = 1^2 + 2^2 - 2(1) - 4(2) + 5 = 1 + 4 - 2 - 8 + 5 = 0f(1,2)=12+22βˆ’2(1)βˆ’4(2)+5=1+4βˆ’2βˆ’8+5=0

  2. Tier 2: Boundary Edge Search:

    • Bottom Edge (y=0,0≀x≀4y = 0, 0 \le x \le 4y=0,0≀x≀4):
      g1(x)=f(x,0)=x2βˆ’2x+5β€…β€ŠβŸΉβ€…β€Šg1β€²(x)=2xβˆ’2=0β€…β€ŠβŸΉβ€…β€Šx=1g_1(x) = f(x, 0) = x^2 - 2x + 5 \implies g_1'(x) = 2x - 2 = 0 \implies x = 1g1​(x)=f(x,0)=x2βˆ’2x+5⟹g1′​(x)=2xβˆ’2=0⟹x=1
      Candidate: (1,0)(1, 0)(1,0). Value: f(1,0)=12βˆ’2(1)+5=4f(1, 0) = 1^2 - 2(1) + 5 = 4f(1,0)=12βˆ’2(1)+5=4.
    • Right Edge (x=4,0≀y≀4x = 4, 0 \le y \le 4x=4,0≀y≀4):
      g2(y)=f(4,y)=16+y2βˆ’8βˆ’4y+5=y2βˆ’4y+13β€…β€ŠβŸΉβ€…β€Šg2β€²(y)=2yβˆ’4=0β€…β€ŠβŸΉβ€…β€Šy=2g_2(y) = f(4, y) = 16 + y^2 - 8 - 4y + 5 = y^2 - 4y + 13 \implies g_2'(y) = 2y - 4 = 0 \implies y = 2g2​(y)=f(4,y)=16+y2βˆ’8βˆ’4y+5=y2βˆ’4y+13⟹g2′​(y)=2yβˆ’4=0⟹y=2
      Candidate: (4,2)(4, 2)(4,2). Value: f(4,2)=22βˆ’4(2)+13=4βˆ’8+13=9f(4, 2) = 2^2 - 4(2) + 13 = 4 - 8 + 13 = 9f(4,2)=22βˆ’4(2)+13=4βˆ’8+13=9.
    • Top Edge (y=4,0≀x≀4y = 4, 0 \le x \le 4y=4,0≀x≀4):
      g3(x)=f(x,4)=x2+16βˆ’2xβˆ’16+5=x2βˆ’2x+5β€…β€ŠβŸΉβ€…β€Šg3β€²(x)=2xβˆ’2=0β€…β€ŠβŸΉβ€…β€Šx=1g_3(x) = f(x, 4) = x^2 + 16 - 2x - 16 + 5 = x^2 - 2x + 5 \implies g_3'(x) = 2x - 2 = 0 \implies x = 1g3​(x)=f(x,4)=x2+16βˆ’2xβˆ’16+5=x2βˆ’2x+5⟹g3′​(x)=2xβˆ’2=0⟹x=1
      Candidate: (1,4)(1, 4)(1,4). Value: f(1,4)=12βˆ’2(1)+5=4f(1, 4) = 1^2 - 2(1) + 5 = 4f(1,4)=12βˆ’2(1)+5=4.
    • Left Edge (x=0,0≀y≀4x = 0, 0 \le y \le 4x=0,0≀y≀4):
      g4(y)=f(0,y)=y2βˆ’4y+5β€…β€ŠβŸΉβ€…β€Šg4β€²(y)=2yβˆ’4=0β€…β€ŠβŸΉβ€…β€Šy=2g_4(y) = f(0, y) = y^2 - 4y + 5 \implies g_4'(y) = 2y - 4 = 0 \implies y = 2g4​(y)=f(0,y)=y2βˆ’4y+5⟹g4′​(y)=2yβˆ’4=0⟹y=2
      Candidate: (0,2)(0, 2)(0,2). Value: f(0,2)=22βˆ’4(2)+5=1f(0, 2) = 2^2 - 4(2) + 5 = 1f(0,2)=22βˆ’4(2)+5=1.
  3. Tier 3: Corner / Vertex Check: Evaluate fff at the four square corners:

    • At (0,0)(0, 0)(0,0): f(0,0)=5f(0, 0) = 5f(0,0)=5
    • At (4,0)(4, 0)(4,0): f(4,0)=42βˆ’2(4)+5=16βˆ’8+5=13f(4, 0) = 4^2 - 2(4) + 5 = 16 - 8 + 5 = 13f(4,0)=42βˆ’2(4)+5=16βˆ’8+5=13
    • At (4,4)(4, 4)(4,4): f(4,4)=42+42βˆ’2(4)βˆ’4(4)+5=16+16βˆ’8βˆ’16+5=13f(4, 4) = 4^2 + 4^2 - 2(4) - 4(4) + 5 = 16 + 16 - 8 - 16 + 5 = 13f(4,4)=42+42βˆ’2(4)βˆ’4(4)+5=16+16βˆ’8βˆ’16+5=13
    • At (0,4)(0, 4)(0,4): f(0,4)=42βˆ’4(4)+5=5f(0, 4) = 4^2 - 4(4) + 5 = 5f(0,4)=42βˆ’4(4)+5=5
  4. Master Comparison & Conclusion: Construct the comparison table:

    Candidate (x,y)(x, y)(x,y) Source Function Value f(x,y)f(x, y)f(x,y)
    (1,2)(1, 2)(1,2) Interior Critical Point 000 (Absolute Minimum)
    (0,2)(0, 2)(0,2) Left Boundary Edge 111
    (1,0)(1, 0)(1,0) Bottom Boundary Edge 444
    (1,4)(1, 4)(1,4) Top Boundary Edge 444
    (0,0)(0, 0)(0,0) Corner Vertex 555
    (0,4)(0, 4)(0,4) Corner Vertex 555
    (4,2)(4, 2)(4,2) Right Boundary Edge 999
    (4,0)(4, 0)(4,0) Corner Vertex 131313 (Absolute Maximum)
    (4,4)(4, 4)(4,4) Corner Vertex 131313 (Absolute Maximum)
    • Absolute Minimum: 000, occurring at (1,2)(1, 2)(1,2).
    • Absolute Maximum: 131313, occurring at (4,0)(4, 0)(4,0) and (4,4)(4, 4)(4,4).

Practice Problem 5: Concept & Reading Check

Open this Skill Block on the Memex to locate the information needed to answer the following:

  • (a) State the two geometric/topological conditions that a domain DβŠ‚R2D \subset \mathbb{R}^2DβŠ‚R2 must satisfy for the Extreme Value Theorem to guarantee the existence of absolute extrema.
  • (b) In Worked Example A, explain why calculating the Hessian discriminant D=fxxfyyβˆ’(fxy)2D = f_{xx}f_{yy} - (f_{xy})^2D=fxx​fyyβ€‹βˆ’(fxy​)2 is unnecessary when executing a global boundary-value optimization problem.
  • (c) In Triangular Foundation Slab Soil Bearing Pressure (Civil & Construction Engineering), explain why civil engineers must evaluate contact pressure at the corner vertices of the footing slab.
Solution Key:
  • (a) By the Extreme Value Theorem for Two Variables (defined in Extrema on Closed and Bounded Regions), the domain DDD must be both closed (including all boundary lines) and bounded (finite in extent). If the domain is open (missing boundary lines) or unbounded (extending infinitely), the function may grow without bound or approach an infimum on a missing edge.
  • (b) The second derivative test classifies local behavior in an open neighborhood. For global optimization on compact domains, the Extreme Value Theorem guarantees that the highest and lowest points must be among the interior critical points, boundary critical points, or corner vertices. A direct numerical evaluation of f(x,y)f(x, y)f(x,y) across this finite list automatically identifies the global extrema without computing DDD.
  • (c) In geotechnical foundation engineering, eccentric overturning loads induce linear and quadratic stress gradients. Because edges intersect sharply at vertices, mathematical boundary parameterizations terminate at corners. Physical foundation toe blowout frequently occurs at these sharp vertices.
Show your steps and calculations here.

Practice Problem 6: Absolute Extrema on a Closed Triangular Region

Find the absolute maximum and minimum values of:

f(x,y)=x2+y2f(x, y) = x^2 + y^2f(x,y)=x2+y2
on the closed triangular region bounded by the vertices (0,0)(0, 0)(0,0), (4,0)(4, 0)(4,0), and (0,3)(0, 3)(0,3).

Solution Key:
  1. Tier 1: Interior Search:

    fx=2x=0β€…β€ŠβŸΉβ€…β€Šx=0,fy=2y=0β€…β€ŠβŸΉβ€…β€Šy=0f_x = 2x = 0 \implies x = 0, \quad f_y = 2y = 0 \implies y = 0fx​=2x=0⟹x=0,fy​=2y=0⟹y=0
    The point (0,0)(0, 0)(0,0) is on the boundary (it is a vertex); thus there are no critical points strictly in the open interior.

  2. Tier 2: Boundary Edge Search:

    • Bottom Edge (y=0,0≀x≀4y = 0, 0 \le x \le 4y=0,0≀x≀4):
      f(x,0)=x2β€…β€ŠβŸΉβ€…β€Šddx[x2]=2x=0β€…β€ŠβŸΉβ€…β€Šx=0f(x, 0) = x^2 \implies \frac{d}{dx}[x^2] = 2x = 0 \implies x = 0f(x,0)=x2⟹dxd​[x2]=2x=0⟹x=0
      This reproduces vertex (0,0)(0, 0)(0,0) (f=0f = 0f=0); endpoint x=4x = 4x=4 yields (4,0)(4, 0)(4,0) (f=16f = 16f=16).
    • Left Edge (x=0,0≀y≀3x = 0, 0 \le y \le 3x=0,0≀y≀3):
      f(0,y)=y2β€…β€ŠβŸΉβ€…β€Šddy[y2]=2y=0β€…β€ŠβŸΉβ€…β€Šy=0f(0, y) = y^2 \implies \frac{d}{dy}[y^2] = 2y = 0 \implies y = 0f(0,y)=y2⟹dyd​[y2]=2y=0⟹y=0
      This reproduces vertex (0,0)(0, 0)(0,0) (f=0f = 0f=0); endpoint y=3y = 3y=3 yields (0,3)(0, 3)(0,3) (f=9f = 9f=9).
    • Hypotenuse Line (y=βˆ’34x+3,0≀x≀4y = -\frac{3}{4}x + 3, 0 \le x \le 4y=βˆ’43​x+3,0≀x≀4): Substitute y=βˆ’34x+3y = -\frac{3}{4}x + 3y=βˆ’43​x+3 into f(x,y)f(x, y)f(x,y):
      g(x)=x2+(βˆ’34x+3)2=x2+916x2βˆ’92x+9=2516x2βˆ’92x+9g(x) = x^2 + \left(-\frac{3}{4}x + 3\right)^2 = x^2 + \frac{9}{16}x^2 - \frac{9}{2}x + 9 = \frac{25}{16}x^2 - \frac{9}{2}x + 9g(x)=x2+(βˆ’43​x+3)2=x2+169​x2βˆ’29​x+9=1625​x2βˆ’29​x+9
      Differentiate and set to zero:
      gβ€²(x)=5016xβˆ’92=258xβˆ’92=0β€…β€ŠβŸΉβ€…β€Šx=9/225/8=92β‹…825=3625=1.44g'(x) = \frac{50}{16}x - \frac{9}{2} = \frac{25}{8}x - \frac{9}{2} = 0 \implies x = \frac{9/2}{25/8} = \frac{9}{2} \cdot \frac{8}{25} = \frac{36}{25} = 1.44gβ€²(x)=1650​xβˆ’29​=825​xβˆ’29​=0⟹x=25/89/2​=29​⋅258​=2536​=1.44
      Since 1.44∈[0,4]1.44 \in [0, 4]1.44∈[0,4], calculate yyy:
      y=βˆ’34(3625)+3=βˆ’2725+7525=4825=1.92y = -\frac{3}{4}\left(\frac{36}{25}\right) + 3 = -\frac{27}{25} + \frac{75}{25} = \frac{48}{25} = 1.92y=βˆ’43​(2536​)+3=βˆ’2527​+2575​=2548​=1.92
      Evaluate function at (1.44,1.92)(1.44, 1.92)(1.44,1.92):
      f(1.44,1.92)=(1.44)2+(1.92)2=2.0736+3.6864=5.76f(1.44, 1.92) = (1.44)^2 + (1.92)^2 = 2.0736 + 3.6864 = 5.76f(1.44,1.92)=(1.44)2+(1.92)2=2.0736+3.6864=5.76
  3. Tier 3: Corner Vertices:

    • At (0,0)(0, 0)(0,0): f(0,0)=0f(0, 0) = 0f(0,0)=0
    • At (4,0)(4, 0)(4,0): f(4,0)=16f(4, 0) = 16f(4,0)=16
    • At (0,3)(0, 3)(0,3): f(0,3)=9f(0, 3) = 9f(0,3)=9
  4. Comparison & Conclusion:

    • Candidate values: 000 (at (0,0)(0, 0)(0,0)), 5.765.765.76 (at (1.44,1.92)(1.44, 1.92)(1.44,1.92)), 999 (at (0,3)(0, 3)(0,3)), 161616 (at (4,0)(4, 0)(4,0)).
    • Absolute Minimum: 000 at (0,0)(0, 0)(0,0).
    • Absolute Maximum: 161616 at (4,0)(4, 0)(4,0).
Show your steps and calculations here.

Practice Problem 7: Absolute Extrema on a Closed Circular Disk

Find the absolute maximum and minimum values of:

f(x,y)=2x2+3y2βˆ’4xβˆ’5f(x, y) = 2x^2 + 3y^2 - 4x - 5f(x,y)=2x2+3y2βˆ’4xβˆ’5
on the closed circular disk D={(x,y)∈R2∣x2+y2≀4}D = \{ (x, y) \in \mathbb{R}^2 \mid x^2 + y^2 \le 4 \}D={(x,y)∈R2∣x2+y2≀4}.

Solution Key:
  1. Tier 1: Interior Search: Compute partial derivatives:

    fx=4xβˆ’4=0β€…β€ŠβŸΉβ€…β€Šx=1f_x = 4x - 4 = 0 \implies x = 1fx​=4xβˆ’4=0⟹x=1
    fy=6y=0β€…β€ŠβŸΉβ€…β€Šy=0f_y = 6y = 0 \implies y = 0fy​=6y=0⟹y=0
    Check boundary inequality: 12+02=1<41^2 + 0^2 = 1 < 412+02=1<4, so (1,0)(1, 0)(1,0) lies strictly inside DDD.
    Evaluate fff at (1,0)(1, 0)(1,0):
    f(1,0)=2(1)2+3(0)2βˆ’4(1)βˆ’5=2βˆ’4βˆ’5=βˆ’7f(1, 0) = 2(1)^2 + 3(0)^2 - 4(1) - 5 = 2 - 4 - 5 = -7f(1,0)=2(1)2+3(0)2βˆ’4(1)βˆ’5=2βˆ’4βˆ’5=βˆ’7

  2. Tier 2: Circular Boundary Search (x2+y2=4x^2 + y^2 = 4x2+y2=4 for βˆ’2≀x≀2-2 \le x \le 2βˆ’2≀x≀2): Solve for y2=4βˆ’x2y^2 = 4 - x^2y2=4βˆ’x2 and substitute directly into f(x,y)f(x, y)f(x,y):

    g(x)=2x2+3(4βˆ’x2)βˆ’4xβˆ’5=2x2+12βˆ’3x2βˆ’4xβˆ’5=βˆ’x2βˆ’4x+7g(x) = 2x^2 + 3(4 - x^2) - 4x - 5 = 2x^2 + 12 - 3x^2 - 4x - 5 = -x^2 - 4x + 7g(x)=2x2+3(4βˆ’x2)βˆ’4xβˆ’5=2x2+12βˆ’3x2βˆ’4xβˆ’5=βˆ’x2βˆ’4x+7
    Differentiate g(x)g(x)g(x) with respect to xxx and set to zero:
    gβ€²(x)=βˆ’2xβˆ’4=0β€…β€ŠβŸΉβ€…β€Š2x=βˆ’4β€…β€ŠβŸΉβ€…β€Šx=βˆ’2g'(x) = -2x - 4 = 0 \implies 2x = -4 \implies x = -2gβ€²(x)=βˆ’2xβˆ’4=0⟹2x=βˆ’4⟹x=βˆ’2
    Since x=βˆ’2x = -2x=βˆ’2 lies in the interval [βˆ’2,2][-2, 2][βˆ’2,2], find yyy:
    y2=4βˆ’(βˆ’2)2=4βˆ’4=0β€…β€ŠβŸΉβ€…β€Šy=0y^2 = 4 - (-2)^2 = 4 - 4 = 0 \implies y = 0y2=4βˆ’(βˆ’2)2=4βˆ’4=0⟹y=0
    Candidate boundary point: (βˆ’2,0)(-2, 0)(βˆ’2,0).
    Evaluate fff at (βˆ’2,0)(-2, 0)(βˆ’2,0):
    f(βˆ’2,0)=2(βˆ’2)2+3(0)2βˆ’4(βˆ’2)βˆ’5=2(4)+8βˆ’5=8+8βˆ’5=11f(-2, 0) = 2(-2)^2 + 3(0)^2 - 4(-2) - 5 = 2(4) + 8 - 5 = 8 + 8 - 5 = 11f(βˆ’2,0)=2(βˆ’2)2+3(0)2βˆ’4(βˆ’2)βˆ’5=2(4)+8βˆ’5=8+8βˆ’5=11

  3. Check Endpoints of the Interval x∈[βˆ’2,2]x \in [-2, 2]x∈[βˆ’2,2]:

    • At x=βˆ’2x = -2x=βˆ’2: y=0y = 0y=0, giving (βˆ’2,0)(-2, 0)(βˆ’2,0) already evaluated (f=11f = 11f=11).
    • At x=2x = 2x=2: y2=4βˆ’22=0β€…β€ŠβŸΉβ€…β€Šy=0y^2 = 4 - 2^2 = 0 \implies y = 0y2=4βˆ’22=0⟹y=0, giving point (2,0)(2, 0)(2,0).
      Evaluate fff at (2,0)(2, 0)(2,0):
      f(2,0)=2(2)2+3(0)2βˆ’4(2)βˆ’5=8βˆ’8βˆ’5=βˆ’5f(2, 0) = 2(2)^2 + 3(0)^2 - 4(2) - 5 = 8 - 8 - 5 = -5f(2,0)=2(2)2+3(0)2βˆ’4(2)βˆ’5=8βˆ’8βˆ’5=βˆ’5
  4. Master Comparison & Conclusion:

    Candidate Point Region Location Function Value f(x,y)f(x, y)f(x,y)
    (1,0)(1, 0)(1,0) Interior Critical Point βˆ’7-7βˆ’7 (Absolute Minimum)
    (2,0)(2, 0)(2,0) Boundary Point βˆ’5-5βˆ’5
    (βˆ’2,0)(-2, 0)(βˆ’2,0) Boundary Point 111111 (Absolute Maximum)
    • Absolute Minimum: βˆ’7-7βˆ’7 at (1,0)(1, 0)(1,0).
    • Absolute Maximum: 111111 at (βˆ’2,0)(-2, 0)(βˆ’2,0).

Show your steps and calculations here.

Regroup 2

  • Boundary Substitution Standard: Emphasize how substituting y2=4βˆ’x2y^2 = 4 - x^2y2=4βˆ’x2 reduces a 2D boundary search into single-variable Calculus 1 optimization.
  • Filtering Discarded Candidates: Remind students to always verify that interior critical points satisfy the domain inequalities before adding them to the candidate table.

Skill Block 3

Worked Example A: Constrained Distance Minimization in 3D Space

Find the point on the plane x+y+z=12x + y + z = 12x+y+z=12 that is closest to the origin, and evaluate the minimum squared distance and minimum Euclidean distance.

Worked Solution:
  1. Formulate Objective and Constraint Functions: Minimizing Euclidean distance x2+y2+z2\sqrt{x^2 + y^2 + z^2}x2+y2+z2​ is mathematically equivalent to minimizing the squared distance:

    f(x,y,z)=x2+y2+z2f(x, y, z) = x^2 + y^2 + z^2f(x,y,z)=x2+y2+z2
    subject to the linear plane constraint:
    g(x,y,z)=x+y+z=12g(x, y, z) = x + y + z = 12g(x,y,z)=x+y+z=12

  2. Compute Gradient Vectors: State the Three-Variable Gradient Formula (defined in The Gradient Vector):

    βˆ‡f=βŸ¨βˆ‚fβˆ‚x,βˆ‚fβˆ‚y,βˆ‚fβˆ‚z⟩=⟨2x,2y,2z⟩=2xi^+2yj^+2zk^\nabla f = \left\langle \frac{\partial f}{\partial x}, \frac{\partial f}{\partial y}, \frac{\partial f}{\partial z} \right\rangle = \langle 2x, 2y, 2z \rangle = 2x\hat{i} + 2y\hat{j} + 2z\hat{k}βˆ‡f=βŸ¨βˆ‚xβˆ‚f​,βˆ‚yβˆ‚f​,βˆ‚zβˆ‚fβ€‹βŸ©=⟨2x,2y,2z⟩=2xi^+2yj^​+2zk^
    βˆ‡g=βŸ¨βˆ‚gβˆ‚x,βˆ‚gβˆ‚y,βˆ‚gβˆ‚z⟩=⟨1,1,1⟩=i^+j^+k^\nabla g = \left\langle \frac{\partial g}{\partial x}, \frac{\partial g}{\partial y}, \frac{\partial g}{\partial z} \right\rangle = \langle 1, 1, 1 \rangle = \hat{i} + \hat{j} + \hat{k}βˆ‡g=βŸ¨βˆ‚xβˆ‚g​,βˆ‚yβˆ‚g​,βˆ‚zβˆ‚gβ€‹βŸ©=⟨1,1,1⟩=i^+j^​+k^

  3. Set Up the Lagrange Multiplier System: State the Lagrange Multiplier Vector System (defined in Lagrange Multipliers):

    βˆ‡f=Ξ»βˆ‡gandx+y+z=12\nabla f = \lambda \nabla g \quad \text{and} \quad x + y + z = 12βˆ‡f=Ξ»βˆ‡gandx+y+z=12
    Equate the scalar components:
    i^:2x=Ξ»(1)β€…β€ŠβŸΉβ€…β€Šx=Ξ»2\hat{i}: 2x = \lambda(1) \implies x = \frac{\lambda}{2}i^:2x=Ξ»(1)⟹x=2λ​
    j^:2y=Ξ»(1)β€…β€ŠβŸΉβ€…β€Šy=Ξ»2\hat{j}: 2y = \lambda(1) \implies y = \frac{\lambda}{2}j^​:2y=Ξ»(1)⟹y=2λ​
    k^:2z=Ξ»(1)β€…β€ŠβŸΉβ€…β€Šz=Ξ»2\hat{k}: 2z = \lambda(1) \implies z = \frac{\lambda}{2}k^:2z=Ξ»(1)⟹z=2λ​

  4. Eliminate Ξ»\lambdaΞ» and Relate Variables: Since all three variables equal Ξ»2\frac{\lambda}{2}2λ​, we obtain:

    x=y=zx = y = zx=y=z

  5. Substitute into the Constraint Equation:

    x+x+x=12β€…β€ŠβŸΉβ€…β€Š3x=12β€…β€ŠβŸΉβ€…β€Šx=4x + x + x = 12 \implies 3x = 12 \implies x = 4x+x+x=12⟹3x=12⟹x=4
    Since x=y=zx = y = zx=y=z, our candidate optimal point is (4,4,4)(4, 4, 4)(4,4,4), with multiplier:
    Ξ»=2x=2(4)=8\lambda = 2x = 2(4) = 8Ξ»=2x=2(4)=8

  6. Evaluate Minimum Distances:

    • Minimum squared distance:
      f(4,4,4)=42+42+42=16+16+16=48f(4, 4, 4) = 4^2 + 4^2 + 4^2 = 16 + 16 + 16 = 48f(4,4,4)=42+42+42=16+16+16=48
    • Minimum Euclidean distance:
      dmin⁑=f(4,4,4)=48=43β‰ˆ6.93d_{\min} = \sqrt{f(4, 4, 4)} = \sqrt{48} = 4\sqrt{3} \approx 6.93dmin​=f(4,4,4)​=48​=43β€‹β‰ˆ6.93
      The point on the plane closest to the origin is (4,4,4)(4, 4, 4)(4,4,4).

Practice Problem 8: Concept & Reading Check

Open this Skill Block on the Memex to locate the information needed to answer the following:

  • (a) Explain why the gradient vectors βˆ‡f\nabla fβˆ‡f and βˆ‡g\nabla gβˆ‡g must be parallel (βˆ‡f=Ξ»βˆ‡g\nabla f = \lambda \nabla gβˆ‡f=Ξ»βˆ‡g) at a constrained optimal point, using the geometry of level curves.
  • (b) In Worked Example A, state the value of the Lagrange multiplier Ξ»\lambdaΞ», and explain what this number predicts if the plane constant is relaxed from 121212 to 131313.
  • (c) In Orbital Insertion Satellite Delta-V Thruster Burn (Aerospace Engineering), explain how the Lagrange multiplier Ξ»\lambdaΞ» informs mission control regarding fuel mass penalties.
Solution Key:
  • (a) By the Geometric Tangency Principle (defined in Lagrange Multipliers), if a constraint curve g=kg = kg=k crosses a level curve f=cf = cf=c at an angle, the value of fff can be increased or decreased by sliding along the constraint. An extremum occurs only where the curves are tangent. Because βˆ‡f\nabla fβˆ‡f is orthogonal to the level curve of fff and βˆ‡g\nabla gβˆ‡g is orthogonal to the constraint curve, the normal vectors must be collinear, yielding βˆ‡f=Ξ»βˆ‡g\nabla f = \lambda \nabla gβˆ‡f=Ξ»βˆ‡g.
  • (b) In Worked Example A, Ξ»=8\lambda = 8Ξ»=8. By the Sensitivity Theorem Ξ»=dfoptdk\lambda = \frac{d f_{\text{opt}}}{d k}Ξ»=dkdfopt​​, relaxing the plane equation from x+y+z=12x + y + z = 12x+y+z=12 to x+y+z=13x + y + z = 13x+y+z=13 (Ξ”k=+1\Delta k = +1Ξ”k=+1) increases the minimum squared distance by Ξ”fβ‰ˆΞ»Ξ”k=8(1)=8\Delta f \approx \lambda \Delta k = 8(1) = 8Ξ”fβ‰ˆΞ»Ξ”k=8(1)=8 units (exact new value: 3Γ—(13/3)2=169/3β‰ˆ56.333 \times (13/3)^2 = 169/3 \approx 56.333Γ—(13/3)2=169/3β‰ˆ56.33, Ξ”f=8.33\Delta f = 8.33Ξ”f=8.33).
  • (c) In aerospace propulsion, propellant is mass-constrained. The multiplier Ξ»=12\lambda = 12Ξ»=12 gives the direct marginal cost in fuel metric units per unit increase in target delta-V, allowing systems engineers to evaluate aerodynamic gravity assist trade-offs.
Show your steps and calculations here.

Practice Problem 9: Constrained Distance to a 2D Line

Minimize the objective function f(x,y)=x2+y2f(x, y) = x^2 + y^2f(x,y)=x2+y2 subject to the linear constraint 4x+3y=124x + 3y = 124x+3y=12.

Solution Key:
  1. Compute Gradients:

    βˆ‡f=⟨2x,2y⟩,βˆ‡g=⟨4,3⟩\nabla f = \langle 2x, 2y \rangle, \quad \nabla g = \langle 4, 3 \rangleβˆ‡f=⟨2x,2y⟩,βˆ‡g=⟨4,3⟩

  2. Formulate Component Equations (βˆ‡f=Ξ»βˆ‡g\nabla f = \lambda \nabla gβˆ‡f=Ξ»βˆ‡g):

    i^:2x=4Ξ»β€…β€ŠβŸΉβ€…β€Šx=2Ξ»\hat{i}: 2x = 4\lambda \implies x = 2\lambdai^:2x=4λ⟹x=2Ξ»
    j^:2y=3Ξ»β€…β€ŠβŸΉβ€…β€Šy=1.5Ξ»=3Ξ»2\hat{j}: 2y = 3\lambda \implies y = 1.5\lambda = \frac{3\lambda}{2}j^​:2y=3λ⟹y=1.5Ξ»=23λ​

  3. Substitute into Constraint:

    4(2Ξ»)+3(1.5Ξ»)=12β€…β€ŠβŸΉβ€…β€Š8Ξ»+4.5Ξ»=12β€…β€ŠβŸΉβ€…β€Š12.5Ξ»=124(2\lambda) + 3(1.5\lambda) = 12 \implies 8\lambda + 4.5\lambda = 12 \implies 12.5\lambda = 124(2Ξ»)+3(1.5Ξ»)=12⟹8Ξ»+4.5Ξ»=12⟹12.5Ξ»=12
    Ξ»=1212.5=2425=0.96\lambda = \frac{12}{12.5} = \frac{24}{25} = 0.96Ξ»=12.512​=2524​=0.96

  4. Compute Optimal Coordinates and Minimum Value:

    x=2(0.96)=1.92=4825x = 2(0.96) = 1.92 = \frac{48}{25}x=2(0.96)=1.92=2548​
    y=1.5(0.96)=1.44=3625y = 1.5(0.96) = 1.44 = \frac{36}{25}y=1.5(0.96)=1.44=2536​
    Evaluate f(x,y)f(x, y)f(x,y):
    f(1.92,1.44)=(1.92)2+(1.44)2=3.6864+2.0736=5.76=14425f(1.92, 1.44) = (1.92)^2 + (1.44)^2 = 3.6864 + 2.0736 = 5.76 = \frac{144}{25}f(1.92,1.44)=(1.92)2+(1.44)2=3.6864+2.0736=5.76=25144​
    (The shortest Euclidean distance from the origin to the line is 5.76=2.4\sqrt{5.76} = 2.45.76​=2.4).

Show your steps and calculations here.

Practice Problem 10: Constrained Extrema on an Elliptical/Circular Boundary

Find the maximum and minimum values of f(x,y)=xyf(x, y) = xyf(x,y)=xy subject to the circular constraint x2+y2=8x^2 + y^2 = 8x2+y2=8.

Solution Key:
  1. Compute Gradients:

    βˆ‡f=⟨y,x⟩,βˆ‡g=⟨2x,2y⟩\nabla f = \langle y, x \rangle, \quad \nabla g = \langle 2x, 2y \rangleβˆ‡f=⟨y,x⟩,βˆ‡g=⟨2x,2y⟩

  2. Lagrange System:

    1. y=2Ξ»xy = 2\lambda xy=2Ξ»x
    2. x=2Ξ»yx = 2\lambda yx=2Ξ»y
    3. x2+y2=8x^2 + y^2 = 8x2+y2=8
  3. Eliminate Ξ»\lambdaΞ»: Multiply equation (1) by yyy and equation (2) by xxx:

    y2=2Ξ»xyandx2=2Ξ»xyβ€…β€ŠβŸΉβ€…β€Šx2=y2β€…β€ŠβŸΉβ€…β€Šy=Β±xy^2 = 2\lambda xy \quad \text{and} \quad x^2 = 2\lambda xy \implies x^2 = y^2 \implies y = \pm xy2=2Ξ»xyandx2=2Ξ»xy⟹x2=y2⟹y=Β±x

  4. Substitute into Constraint:

    x2+x2=8β€…β€ŠβŸΉβ€…β€Š2x2=8β€…β€ŠβŸΉβ€…β€Šx2=4β€…β€ŠβŸΉβ€…β€Šx=Β±2x^2 + x^2 = 8 \implies 2x^2 = 8 \implies x^2 = 4 \implies x = \pm 2x2+x2=8⟹2x2=8⟹x2=4⟹x=Β±2

    • For y=xy = xy=x, candidates are (2,2)(2, 2)(2,2) and (βˆ’2,βˆ’2)(-2, -2)(βˆ’2,βˆ’2).
    • For y=βˆ’xy = -xy=βˆ’x, candidates are (2,βˆ’2)(2, -2)(2,βˆ’2) and (βˆ’2,2)(-2, 2)(βˆ’2,2).
  5. Evaluate Objective f(x,y)=xyf(x, y) = xyf(x,y)=xy:

    • At (2,2)(2, 2)(2,2): f(2,2)=(2)(2)=4f(2, 2) = (2)(2) = 4f(2,2)=(2)(2)=4
    • At (βˆ’2,βˆ’2)(-2, -2)(βˆ’2,βˆ’2): f(βˆ’2,βˆ’2)=(βˆ’2)(βˆ’2)=4f(-2, -2) = (-2)(-2) = 4f(βˆ’2,βˆ’2)=(βˆ’2)(βˆ’2)=4
    • At (2,βˆ’2)(2, -2)(2,βˆ’2): f(2,βˆ’2)=(2)(βˆ’2)=βˆ’4f(2, -2) = (2)(-2) = -4f(2,βˆ’2)=(2)(βˆ’2)=βˆ’4
    • At (βˆ’2,2)(-2, 2)(βˆ’2,2): f(βˆ’2,2)=(βˆ’2)(2)=βˆ’4f(-2, 2) = (-2)(2) = -4f(βˆ’2,2)=(βˆ’2)(2)=βˆ’4

    Conclusion:

    • Maximum value: 444 at (2,2)(2, 2)(2,2) and (βˆ’2,βˆ’2)(-2, -2)(βˆ’2,βˆ’2).
    • Minimum value: βˆ’4-4βˆ’4 at (2,βˆ’2)(2, -2)(2,βˆ’2) and (βˆ’2,2)(-2, 2)(βˆ’2,2).
Show your steps and calculations here.

Practice Problem 11: Geometric Optimization of a Rectangular Channel

An environmental engineer is sizing a rectangular open water channel with cross-sectional area A=xy=9Β m2A = xy = 9\text{ m}^2A=xy=9Β m2, where xxx and yyy are positive dimensions. Using Lagrange multipliers, find the dimensions (x,y)(x, y)(x,y) that minimize the total perimeter P(x,y)=2x+2yP(x, y) = 2x + 2yP(x,y)=2x+2y.

Solution Key:
  1. Formulate Functions:

    • Objective: P(x,y)=2x+2yP(x, y) = 2x + 2yP(x,y)=2x+2y
    • Constraint: g(x,y)=xy=9g(x, y) = xy = 9g(x,y)=xy=9 (with x>0,y>0x > 0, y > 0x>0,y>0)
  2. Compute Gradients:

    βˆ‡P=⟨2,2⟩,βˆ‡g=⟨y,x⟩\nabla P = \langle 2, 2 \rangle, \quad \nabla g = \langle y, x \rangleβˆ‡P=⟨2,2⟩,βˆ‡g=⟨y,x⟩

  3. Set Up System:

    i^:2=Ξ»yβ€…β€ŠβŸΉβ€…β€ŠΞ»=2y\hat{i}: 2 = \lambda y \implies \lambda = \frac{2}{y}i^:2=Ξ»y⟹λ=y2​
    j^:2=Ξ»xβ€…β€ŠβŸΉβ€…β€ŠΞ»=2x\hat{j}: 2 = \lambda x \implies \lambda = \frac{2}{x}j^​:2=Ξ»x⟹λ=x2​

  4. Eliminate Ξ»\lambdaΞ» and Solve:

    2y=2xβ€…β€ŠβŸΉβ€…β€Šx=y\frac{2}{y} = \frac{2}{x} \implies x = yy2​=x2β€‹βŸΉx=y
    Substitute x=yx = yx=y into xy=9xy = 9xy=9:
    x2=9β€…β€ŠβŸΉβ€…β€Šx=3Β mx^2 = 9 \implies x = 3\text{ m}x2=9⟹x=3Β m
    Since x=yx = yx=y, y=3Β my = 3\text{ m}y=3Β m. Minimum perimeter:
    P(3,3)=2(3)+2(3)=6+6=12Β mP(3, 3) = 2(3) + 2(3) = 6 + 6 = 12\text{ m}P(3,3)=2(3)+2(3)=6+6=12Β m

Show your steps and calculations here.

Practice Problem 12: Three-Variable Box Surface Area Minimization

An open-top rectangular shipping crate must have a specified volume of V=xyz=32Β m3V = xyz = 32\text{ m}^3V=xyz=32Β m3. Using the method of Lagrange multipliers, determine the dimensions x,y,zx, y, zx,y,z that minimize the external surface area S(x,y,z)=xy+2xz+2yzS(x, y, z) = xy + 2xz + 2yzS(x,y,z)=xy+2xz+2yz.

Solution Key:
  1. Formulate Functions:

    • Objective: S(x,y,z)=xy+2xz+2yzS(x, y, z) = xy + 2xz + 2yzS(x,y,z)=xy+2xz+2yz (base plus four vertical faces)
    • Constraint: g(x,y,z)=xyz=32g(x, y, z) = xyz = 32g(x,y,z)=xyz=32, with x,y,z>0x, y, z > 0x,y,z>0
  2. Compute Gradients:

    βˆ‡S=⟨y+2z, x+2z, 2x+2y⟩\nabla S = \langle y + 2z, \, x + 2z, \, 2x + 2y \rangleβˆ‡S=⟨y+2z,x+2z,2x+2y⟩
    βˆ‡g=⟨yz, xz, xy⟩\nabla g = \langle yz, \, xz, \, xy \rangleβˆ‡g=⟨yz,xz,xy⟩

  3. Lagrange Component Equations (βˆ‡S=Ξ»βˆ‡g\nabla S = \lambda \nabla gβˆ‡S=Ξ»βˆ‡g):

    1. y+2z=Ξ»yzy + 2z = \lambda yzy+2z=Ξ»yz
    2. x+2z=Ξ»xzx + 2z = \lambda xzx+2z=Ξ»xz
    3. 2x+2y=Ξ»xy2x + 2y = \lambda xy2x+2y=Ξ»xy
  4. Eliminate Ξ»\lambdaΞ» and Relate Dimensions: Multiply equation (1) by xxx, (2) by yyy, and (3) by zzz:

    xy+2xz=Ξ»xyzxy + 2xz = \lambda xyzxy+2xz=Ξ»xyz
    xy+2yz=Ξ»xyzxy + 2yz = \lambda xyzxy+2yz=Ξ»xyz
    2xz+2yz=Ξ»xyz2xz + 2yz = \lambda xyz2xz+2yz=Ξ»xyz
    Equating the first two expressions:
    xy+2xz=xy+2yzβ€…β€ŠβŸΉβ€…β€Š2xz=2yzβ€…β€ŠβŸΉβ€…β€Šx=yxy + 2xz = xy + 2yz \implies 2xz = 2yz \implies x = yxy+2xz=xy+2yz⟹2xz=2yz⟹x=y
    Equating the first and third expressions with y=xy = xy=x:
    x2+2xz=2xz+2xzβ€…β€ŠβŸΉβ€…β€Šx2+2xz=4xzβ€…β€ŠβŸΉβ€…β€Šx2=2xzβ€…β€ŠβŸΉβ€…β€Šx=2zβ€…β€ŠβŸΉβ€…β€Šz=x2x^2 + 2xz = 2xz + 2xz \implies x^2 + 2xz = 4xz \implies x^2 = 2xz \implies x = 2z \implies z = \frac{x}{2}x2+2xz=2xz+2xz⟹x2+2xz=4xz⟹x2=2xz⟹x=2z⟹z=2x​

  5. Substitute into Volume Constraint:

    xyz=32β€…β€ŠβŸΉβ€…β€Š(x)(x)(x2)=32β€…β€ŠβŸΉβ€…β€Šx32=32β€…β€ŠβŸΉβ€…β€Šx3=64β€…β€ŠβŸΉβ€…β€Šx=4Β mxyz = 32 \implies (x)(x)\left(\frac{x}{2}\right) = 32 \implies \frac{x^3}{2} = 32 \implies x^3 = 64 \implies x = 4\text{ m}xyz=32⟹(x)(x)(2x​)=32⟹2x3​=32⟹x3=64⟹x=4Β m
    Now find yyy and zzz:
    y=x=4Β my = x = 4\text{ m}y=x=4Β m
    z=42=2Β mz = \frac{4}{2} = 2\text{ m}z=24​=2Β m

  6. Evaluate Minimum Surface Area:

    S(4,4,2)=(4)(4)+2(4)(2)+2(4)(2)=16+16+16=48Β m2S(4, 4, 2) = (4)(4) + 2(4)(2) + 2(4)(2) = 16 + 16 + 16 = 48\text{ m}^2S(4,4,2)=(4)(4)+2(4)(2)+2(4)(2)=16+16+16=48Β m2
    The optimal crate has a square base of 4Β mΓ—4Β m4\text{ m} \times 4\text{ m}4Β mΓ—4Β m and a height of 2Β m2\text{ m}2Β m, achieving a minimum surface area of 48Β m248\text{ m}^248Β m2.


Show your steps and calculations here.

Regroup 3

  • Component Equation Setup: Reiterate how multiplying by xxx, yyy, and zzz cleanly aligns terms to Ξ»xyz\lambda xyzΞ»xyz in symmetric 3D problems.
  • Physical Meaning of Constraints: Discuss how Ξ»\lambdaΞ» serves as an economic shadow price in industrial fabrication when evaluating the cost of material limits.

Synthesis Wrap-up

  • Core Takeaways:
    1. Unconstrained Extrema: Use βˆ‡f=0\nabla f = \mathbf{0}βˆ‡f=0 and classify with the Hessian discriminant D=fxxfyyβˆ’(fxy)2D = f_{xx}f_{yy} - (f_{xy})^2D=fxx​fyyβ€‹βˆ’(fxy​)2.
    2. Closed & Bounded Extrema: Use the 3-tier strategy (interior, boundary curves, vertices) and evaluate directly in a global comparison table.
    3. Constrained Extrema: Use Lagrange multipliers βˆ‡f=Ξ»βˆ‡g\nabla f = \lambda \nabla gβˆ‡f=Ξ»βˆ‡g when restricted to smooth curves or surfaces, and interpret Ξ»=βˆ‚fopt/βˆ‚k\lambda = \partial f_{\text{opt}}/\partial kΞ»=βˆ‚fopt​/βˆ‚k as sensitivity.
  • Homework Assignment (Knewton Alta):
    • Critical Points and the Second Derivative Test for Functions of Several Variables
    • Maximization and Minimization Problems for Functions of Several Variables
    • Lagrange Multipliers
  • Next Session Preview: Session 14 introduces Double Integrals over Rectangular and General Regionsβ€”accumulating volume and mass under multivariable surfaces.