Skip to main content
Study NotesCELE · Engineering Mathematics

CELE Engineering Mathematics — Differential CalculusStudy Notes

Study notes covering Differential Calculus for CELE reviewers. Use them for your first full pass through the chapter; the shorter formats suit revision later.

Super Tutor's review for GELE (as “Mathematics”) uses these same chapters, so they live here in the CELE section.

Differential Calculus - Study Notes

Differential calculus is the mathematical foundation for analyzing rates of change, optimizing structural designs, and solving motion problems essential to civil engineering practice. This study module covers limits, continuity, derivatives, extrema, and related rates—topics that recur annually on the MSTE (Mathematics, Science and Technology) section of the PRC Licensure Examination. Mastery of these concepts enables engineers to solve real-world problems such as minimizing material cost in beam design, finding critical buckling loads, and analyzing traffic flow dynamics. This material is pitched at the professional review level, emphasizing board-style worked examples in SI units and practical applications aligned with NSCP 2015, AISC 360, and ACI 318 standards.

1. Limits and Continuity — Foundation of the Derivative

A limit describes the value that a function approaches as the independent variable moves toward a specified point. Formally, lim⁡x→af(x)=L\lim_{x\to a}f(x) = L means that for any small positive number ϵ\epsilon, there exists a small positive number δ\delta such that whenever 0<∣x−a∣<δ0 < |x - a| < \delta, we have ∣f(x)−L∣<ϵ|f(x) - L| < \epsilon.

Continuity at a point: A function $f(x)$ is continuous at x=ax = a if and only if three conditions hold:

  1. $f(a)$ is defined
  2. lim⁡x→af(x)\lim_{x\to a}f(x) exists
  3. lim⁡x→af(x)=f(a)\lim_{x\to a}f(x) = f(a)

Indeterminate forms arise when direct substitution yields 00\frac{0}{0}, ∞∞\frac{\infty}{\infty}, 0⋅∞0 \cdot \infty, ∞−∞\infty - \infty, 000^0, 1∞1^{\infty}, or ∞0\infty^0. For these cases, L'Hôpital's Rule provides a systematic solution: if lim⁡x→af(x)g(x)\lim_{x\to a}\frac{f(x)}{g(x)} yields an indeterminate form and both $f'(x)$ and $g'(x)$ exist near aa, then

lim⁡x→af(x)g(x)=lim⁡x→af′(x)g′(x)\lim_{x\to a}\frac{f(x)}{g(x)} = \lim_{x\to a}\frac{f'(x)}{g'(x)}

provisioning the latter limit exists. You may apply L'Hôpital's Rule repeatedly if necessary.

Limit laws (for limits that exist):

  • Sum rule: lim⁡[f(x)+g(x)]=lim⁡f(x)+lim⁡g(x)\lim[f(x) + g(x)] = \lim f(x) + \lim g(x)
  • Product rule: lim⁡[f(x)⋅g(x)]=[lim⁡f(x)]⋅[lim⁡g(x)]\lim[f(x) \cdot g(x)] = [\lim f(x)] \cdot [\lim g(x)]
  • Quotient rule: lim⁡f(x)g(x)=lim⁡f(x)lim⁡g(x)\lim\frac{f(x)}{g(x)} = \frac{\lim f(x)}{\lim g(x)} (provided lim⁡g(x)≠0\lim g(x) \neq 0)
  • Power rule: lim⁡[f(x)]n=[lim⁡f(x)]n\lim[f(x)]^n = [\lim f(x)]^n

Important standard limits:

lim⁡x→0sin⁡xx=1,lim⁡x→01−cos⁡xx2=12,lim⁡x→∞(1+1x)x=e\lim_{x\to 0}\frac{\sin x}{x} = 1, \quad \lim_{x\to 0}\frac{1 - \cos x}{x^2} = \frac{1}{2}, \quad \lim_{x\to \infty}\left(1 + \frac{1}{x}\right)^x = e

These limits appear frequently in engineering applications—for example, when analyzing small-angle approximations in structural dynamics or wave propagation.

Examples

Evaluate lim⁡x→2x2−4x−2\lim_{x\to 2}\frac{x^2 - 4}{x - 2}.

Solution

Direct substitution gives 00\frac{0}{0}, an indeterminate form. Factor the numerator:

lim⁡x→2x2−4x−2=lim⁡x→2(x−2)(x+2)x−2=lim⁡x→2(x+2)=4\lim_{x\to 2}\frac{x^2 - 4}{x - 2} = \lim_{x\to 2}\frac{(x-2)(x+2)}{x-2} = \lim_{x\to 2}(x+2) = 4

Alternatively, apply L'Hôpital's Rule:

lim⁡x→2x2−4x−2=lim⁡x→22x1=4\lim_{x\to 2}\frac{x^2 - 4}{x - 2} = \lim_{x\to 2}\frac{2x}{1} = 4

Use L'Hôpital's Rule to evaluate lim⁡x→0ex−1−xx2\lim_{x\to 0}\frac{e^x - 1 - x}{x^2}.

Solution

Direct substitution: e0−1−002=00\frac{e^0 - 1 - 0}{0^2} = \frac{0}{0}, indeterminate. Apply L'Hôpital's:

lim⁡x→0ex−1−xx2=lim⁡x→0ex−12x\lim_{x\to 0}\frac{e^x - 1 - x}{x^2} = \lim_{x\to 0}\frac{e^x - 1}{2x}

This is still 00\frac{0}{0}, so apply L'Hôpital's again:

lim⁡x→0ex−12x=lim⁡x→0ex2=12\lim_{x\to 0}\frac{e^x - 1}{2x} = \lim_{x\to 0}\frac{e^x}{2} = \frac{1}{2}

Note: Each application of L'Hôpital's Rule requires differentiation of the numerator and denominator separately.

Determine where f(x)=x2−1x−1f(x) = \frac{x^2 - 1}{x - 1} is continuous.

Solution

The function has a removable discontinuity at x=1x = 1 because the denominator vanishes there. For x≠1x \neq 1:

f(x)=(x−1)(x+1)x−1=x+1f(x) = \frac{(x-1)(x+1)}{x-1} = x + 1

Thus lim⁡x→1f(x)=2\lim_{x\to 1}f(x) = 2, but $f(1)$ is undefined. If we redefine f(1)=2f(1) = 2, the function becomes continuous everywhere. For all other points, the function is continuous by the limit laws.

Key points

  • A limit is the value a function approaches, not necessarily the value it attains at a point
  • Continuity requires the limit to exist and equal the function value; check all three conditions
  • Indeterminate forms must be resolved algebraically or via L'Hôpital's Rule before evaluating
  • L'Hôpital's Rule applies only to indeterminate forms and requires differentiation of numerator and denominator separately
  • Standard limits like sin⁡x/x→1\sin x / x \to 1 as x→0x \to 0 are essential for engineering calculations

2. The Derivative — Definition and Geometric Interpretation

The derivative of $f(x)$ at point xx measures the instantaneous rate of change of ff with respect to xx. Geometrically, it is the slope of the tangent line to the curve y=f(x)y = f(x) at that point.

Definition (limit form):

f′(x)=lim⁡h→0f(x+h)−f(x)hf'(x) = \lim_{h\to 0}\frac{f(x+h) - f(x)}{h}

Alternatively, at a specific point x=ax = a:

f′(a)=lim⁡x→af(x)−f(a)x−af'(a) = \lim_{x\to a}\frac{f(x) - f(a)}{x - a}

Notation: The derivative of y=f(x)y = f(x) is denoted as $f'(x)$, dydx\frac{dy}{dx}, dfdx\frac{df}{dx}, or $y'$.

Equation of the tangent line at point (x0,f(x0))(x_0, f(x_0)):

y−f(x0)=f′(x0)(x−x0)y - f(x_0) = f'(x_0)(x - x_0)

Differentiability: A function is differentiable at x=ax = a if $f'(a)$ exists. Differentiability implies continuity, but continuity does not guarantee differentiability (e.g., y=∣x∣y = |x| is continuous at x=0x = 0 but not differentiable there due to a sharp corner).

Physical interpretation in civil engineering: If yy represents displacement and xx represents time, then dydx\frac{dy}{dx} is velocity. If stress σ\sigma varies with strain ϵ\epsilon, then dσdϵ\frac{d\sigma}{d\epsilon} is the elastic modulus. Derivatives quantify how quickly structural responses change under loading—critical for analyzing beam deflections, column buckling, and dynamic behavior.

Examples

Using the limit definition, find the derivative of f(x)=x2f(x) = x^2 at x=3x = 3.

Solution

Apply the definition:

f′(3)=lim⁡h→0f(3+h)−f(3)h=lim⁡h→0(3+h)2−9hf'(3) = \lim_{h\to 0}\frac{f(3+h) - f(3)}{h} = \lim_{h\to 0}\frac{(3+h)^2 - 9}{h}

Expand (3+h)2=9+6h+h2(3+h)^2 = 9 + 6h + h^2:

f′(3)=lim⁡h→09+6h+h2−9h=lim⁡h→06h+h2h=lim⁡h→0(6+h)=6f'(3) = \lim_{h\to 0}\frac{9 + 6h + h^2 - 9}{h} = \lim_{h\to 0}\frac{6h + h^2}{h} = \lim_{h\to 0}(6 + h) = 6

Therefore, f′(3)=6f'(3) = 6.

Find the equation of the tangent line to y=x2−2xy = x^2 - 2x at the point where x=1x = 1.

Solution

First, find the point on the curve: y(1)=1−2=−1y(1) = 1 - 2 = -1, so the point is $(1, -1)$.
Next, find the derivative: y′=2x−2y' = 2x - 2, so y′(1)=2(1)−2=0y'(1) = 2(1) - 2 = 0.
The tangent line has slope m=0m = 0 and passes through $(1, -1)$:

y−(−1)=0(x−1)⇒y=−1y - (-1) = 0(x - 1) \Rightarrow y = -1

The tangent line is horizontal, which makes sense because x=1x = 1 is the vertex of the parabola y=(x−1)2−1y = (x-1)^2 - 1.

Analyze the differentiability of f(x)=∣x∣f(x) = |x| at x=0x = 0.

Solution

Check the left and right derivatives:

  • Left derivative: lim⁡h→0−f(0+h)−f(0)h=lim⁡h→0−∣h∣−0h=lim⁡h→0−−hh=−1\lim_{h\to 0^-}\frac{f(0+h) - f(0)}{h} = \lim_{h\to 0^-}\frac{|h| - 0}{h} = \lim_{h\to 0^-}\frac{-h}{h} = -1
  • Right derivative: lim⁡h→0+f(0+h)−f(0)h=lim⁡h→0+∣h∣−0h=lim⁡h→0+hh=1\lim_{h\to 0^+}\frac{f(0+h) - f(0)}{h} = \lim_{h\to 0^+}\frac{|h| - 0}{h} = \lim_{h\to 0^+}\frac{h}{h} = 1

Since the left and right derivatives differ (−1≠1-1 \neq 1), the derivative does not exist at x=0x = 0. The function f(x)=∣x∣f(x) = |x| is continuous but not differentiable at x=0x = 0 due to the sharp corner (cusp). This illustrates that a function can be continuous everywhere but differentiable only where it is smooth.

Key points

  • The derivative is the limit of the difference quotient; it measures instantaneous rate of change
  • Geometrically, f′(x0)f'(x_0) is the slope of the tangent line at (x0,f(x0))(x_0, f(x_0))
  • Differentiability implies continuity, but not vice versa; check for sharp corners or discontinuities
  • The tangent line equation is y−y0=f′(x0)(x−x0)y - y_0 = f'(x_0)(x - x_0)
  • In engineering: derivatives represent rates of change (velocity, strain-stress relationships, reaction rates)

3. Derivative Rules — Power, Product, Quotient, and Chain Rules

Rather than use the limit definition repeatedly, we apply systematic rules to find derivatives efficiently.

Power Rule: For any real constant nn,

ddxxn=nxn−1\frac{d}{dx}x^n = nx^{n-1}

This applies even when nn is negative or fractional (e.g., ddxx−2=−2x−3\frac{d}{dx}x^{-2} = -2x^{-3}, ddxx=ddxx1/2=12x−1/2\frac{d}{dx}\sqrt{x} = \frac{d}{dx}x^{1/2} = \frac{1}{2}x^{-1/2}).

Constant Multiple Rule:

ddx[cf(x)]=cdfdx\frac{d}{dx}[cf(x)] = c\frac{df}{dx}

Sum and Difference Rule:

ddx[f(x)±g(x)]=dfdx±dgdx\frac{d}{dx}[f(x) \pm g(x)] = \frac{df}{dx} \pm \frac{dg}{dx}

Product Rule: For the product of two functions $u(x)$ and $v(x)$,

(uv)′=u′v+uv′(uv)' = u'v + uv'

Mnemonic: "first times derivative of second, plus second times derivative of first."

Quotient Rule: For the quotient uv\frac{u}{v} where v≠0v \neq 0,

(uv)′=u′v−uv′v2\left(\frac{u}{v}\right)' = \frac{u'v - uv'}{v^2}

Mnemonic: "low dd-high minus high dd-low, over low low" (denominator-derivative high minus high-derivative low, divided by low squared).

Chain Rule: If y=f(u)y = f(u) and u=g(x)u = g(x), then

dydx=dydu⋅dudx=f′(g(x))⋅g′(x)\frac{dy}{dx} = \frac{dy}{du} \cdot \frac{du}{dx} = f'(g(x)) \cdot g'(x)

Alternatively, if y=f(g(x))y = f(g(x)), then y′=f′(g(x))⋅g′(x)y' = f'(g(x)) \cdot g'(x). The chain rule is crucial for composite functions and appears constantly in engineering calculations.

Derivatives of trigonometric functions:

ddxsin⁡x=cos⁡x,ddxcos⁡x=−sin⁡x,ddxtan⁡x=sec⁡2x\frac{d}{dx}\sin x = \cos x, \quad \frac{d}{dx}\cos x = -\sin x, \quad \frac{d}{dx}\tan x = \sec^2 x
ddxcot⁡x=−csc⁡2x,ddxsec⁡x=sec⁡xtan⁡x,ddxcsc⁡x=−csc⁡xcot⁡x\frac{d}{dx}\cot x = -\csc^2 x, \quad \frac{d}{dx}\sec x = \sec x \tan x, \quad \frac{d}{dx}\csc x = -\csc x \cot x

Derivatives of exponential and logarithmic functions:

ddxex=ex,ddxax=axln⁡a(a>0,a≠1)\frac{d}{dx}e^x = e^x, \quad \frac{d}{dx}a^x = a^x \ln a \quad (a > 0, a \neq 1)
ddxln⁡x=1x,ddxlog⁡ax=1xln⁡a\frac{d}{dx}\ln x = \frac{1}{x}, \quad \frac{d}{dx}\log_a x = \frac{1}{x \ln a}

Inverse trigonometric function derivatives:

ddxarcsin⁡x=11−x2,ddxarccos⁡x=−11−x2,ddxarctan⁡x=11+x2\frac{d}{dx}\arcsin x = \frac{1}{\sqrt{1-x^2}}, \quad \frac{d}{dx}\arccos x = -\frac{1}{\sqrt{1-x^2}}, \quad \frac{d}{dx}\arctan x = \frac{1}{1+x^2}

These rules form the toolkit for all derivative calculations. Mastering them enables rapid solution of board-exam problems.

Examples

Differentiate y=(3x2+2x−1)(x3−4)y = (3x^2 + 2x - 1)(x^3 - 4).

Solution

Use the product rule with u=3x2+2x−1u = 3x^2 + 2x - 1 and v=x3−4v = x^3 - 4.
First, find the derivatives:

  • u′=6x+2u' = 6x + 2
  • v′=3x2v' = 3x^2

Apply the product rule:

y′=u′v+uv′=(6x+2)(x3−4)+(3x2+2x−1)(3x2)y' = u'v + uv' = (6x + 2)(x^3 - 4) + (3x^2 + 2x - 1)(3x^2)

Expand the first term: (6x+2)(x3−4)=6x4−24x+2x3−8(6x + 2)(x^3 - 4) = 6x^4 - 24x + 2x^3 - 8
Expand the second term: (3x2+2x−1)(3x2)=9x4+6x3−3x2(3x^2 + 2x - 1)(3x^2) = 9x^4 + 6x^3 - 3x^2
Combine:

y′=6x4−24x+2x3−8+9x4+6x3−3x2=15x4+8x3−3x2−24x−8y' = 6x^4 - 24x + 2x^3 - 8 + 9x^4 + 6x^3 - 3x^2 = 15x^4 + 8x^3 - 3x^2 - 24x - 8

Differentiate y=2x2+3x−1y = \frac{2x^2 + 3}{x - 1}.

Solution

Use the quotient rule with u=2x2+3u = 2x^2 + 3 and v=x−1v = x - 1.
Find the derivatives:

  • u′=4xu' = 4x
  • v′=1v' = 1

Apply the quotient rule:

y′=u′v−uv′v2=4x(x−1)−(2x2+3)(1)(x−1)2=4x2−4x−2x2−3(x−1)2=2x2−4x−3(x−1)2y' = \frac{u'v - uv'}{v^2} = \frac{4x(x-1) - (2x^2 + 3)(1)}{(x-1)^2} = \frac{4x^2 - 4x - 2x^2 - 3}{(x-1)^2} = \frac{2x^2 - 4x - 3}{(x-1)^2}

Differentiate y=(5x3−2)7y = (5x^3 - 2)^7.

Solution

Use the chain rule. Let u=5x3−2u = 5x^3 - 2, so y=u7y = u^7.
Apply the chain rule:

dydx=dydu⋅dudx=7u6⋅15x2=7(5x3−2)6⋅15x2=105x2(5x3−2)6\frac{dy}{dx} = \frac{dy}{du} \cdot \frac{du}{dx} = 7u^6 \cdot 15x^2 = 7(5x^3 - 2)^6 \cdot 15x^2 = 105x^2(5x^3 - 2)^6

Differentiate y=esin⁡xy = e^{\sin x}.

Solution

Use the chain rule. Let u=sin⁡xu = \sin x, so y=euy = e^u.
Apply the chain rule:

dydx=dydu⋅dudx=eu⋅cos⁡x=esin⁡xcos⁡x\frac{dy}{dx} = \frac{dy}{du} \cdot \frac{du}{dx} = e^u \cdot \cos x = e^{\sin x} \cos x

Differentiate y=ln⁡(x2+3x)y = \ln(x^2 + 3x).

Solution

Use the chain rule. Let u=x2+3xu = x^2 + 3x, so y=ln⁡uy = \ln u.
Apply the chain rule:

dydx=1u⋅dudx=1x2+3x⋅(2x+3)=2x+3x2+3x\frac{dy}{dx} = \frac{1}{u} \cdot \frac{du}{dx} = \frac{1}{x^2 + 3x} \cdot (2x + 3) = \frac{2x + 3}{x^2 + 3x}

Key points

  • Power rule: ddxxn=nxn−1\frac{d}{dx}x^n = nx^{n-1} applies to all real exponents, including fractions and negatives
  • Product rule and quotient rule have specific, easy-to-confuse forms; practice to avoid sign errors
  • Chain rule is essential for nested functions; apply from outside to inside
  • Memorize standard derivatives of trig, exponential, and log functions—they recur constantly
  • Combine rules systematically: often a single problem requires power, product, and chain rules together

4. Higher-Order Derivatives and Concavity

The second derivative $f''(x)$ (or d2ydx2\frac{d^2y}{dx^2}) is the derivative of the first derivative. It measures how quickly the rate of change is changing—i.e., the curvature or concavity of the graph.

Concavity interpretation:

  • Concave up (convex): $f''(x) > 0$. The graph curves upward; tangent lines lie below the curve. Geometrically, the function is "holding water."
  • Concave down: $f''(x) < 0$. The graph curves downward; tangent lines lie above the curve.
  • Inflection point: A point where the concavity changes, where f′′(x)=0f''(x) = 0 (and the sign of $f''$ changes).

Second Derivative Test for Extrema: At a critical point x=cx = c where f′(c)=0f'(c) = 0:

  • If $f''(c) > 0$, then ff has a local minimum at cc.
  • If $f''(c) < 0$, then ff has a local maximum at cc.
  • If f′′(c)=0f''(c) = 0, the test is inconclusive; use the first derivative test instead.

Third and higher derivatives: Successive differentiation yields $f'''(x)$, f(4)(x)f^{(4)}(x), etc. In engineering applications, the third derivative (often called "jerk" in motion analysis) represents the rate of change of acceleration and is important in vehicle dynamics and seismic analysis. Higher derivatives are rare in basic optimization but essential in Taylor series approximations.

Application in civil engineering: When analyzing beam deflection under load, the deflection function $y(x)$ has a physical meaning:

  • $y'(x)$ = slope of the beam (angle of rotation)
  • $y''(x)$ = curvature (proportional to bending moment MM; the beam equation is M=EId2ydx2M = EI\frac{d^2y}{dx^2} where EE is elastic modulus and II is second moment of area)
  • $y'''(x)$ = shear force (related to internal shear VV; V=EId3ydx3V = EI\frac{d^3y}{dx^3})
  • y(4)(x)y^{(4)}(x) = distributed load qq (the load equation is q=EId4ydx4q = EI\frac{d^4y}{dx^4})

This direct relationship between derivatives and physical quantities makes differential calculus indispensable in structural mechanics.

Examples

For f(x)=x3−3x2+2f(x) = x^3 - 3x^2 + 2, find the second derivative and identify concavity regions.

Solution

Find the first derivative: f′(x)=3x2−6xf'(x) = 3x^2 - 6x.
Find the second derivative: f′′(x)=6x−6=6(x−1)f''(x) = 6x - 6 = 6(x - 1).

Set f′′(x)=0f''(x) = 0: x=1x = 1 is the inflection point.

  • For x<1x < 1: $f''(x) < 0$ → concave down
  • For x>1x > 1: $f''(x) > 0$ → concave up

At x=1x = 1, the concavity changes from down to up, confirming the inflection point.

A cantilever beam of length L=4L = 4 m has deflection y(x)=−w024EIx2(x2−4Lx+6L2)y(x) = -\frac{w_0}{24EI}x^2(x^2 - 4Lx + 6L^2), where w0w_0 is the distributed load. Using beam theory, find the bending moment at x=2x = 2 m. (Given: E=200E = 200 GPa, I=8×10−5I = 8 \times 10^{-5} m⁴, w0=10w_0 = 10 kN/m.)

Solution

The bending moment is M(x)=EId2ydx2M(x) = EI\frac{d^2y}{dx^2}. Rather than expand the full deflection formula, recognize that the deflection equation encodes the load distribution. Differentiate twice:

First derivative: dydx=−w024EIddx[x2(x2−4Lx+6L2)]\frac{dy}{dx} = -\frac{w_0}{24EI}\frac{d}{dx}[x^2(x^2 - 4Lx + 6L^2)]

Using the product rule: dydx=−w024EI[2x(x2−4Lx+6L2)+x2(2x−4L)]\frac{dy}{dx} = -\frac{w_0}{24EI}[2x(x^2 - 4Lx + 6L^2) + x^2(2x - 4L)]

Simplifying: dydx=−w024EI[2x3−8Lx2+12L2x+2x3−4Lx2]=−w024EI[4x3−12Lx2+12L2x]\frac{dy}{dx} = -\frac{w_0}{24EI}[2x^3 - 8Lx^2 + 12L^2x + 2x^3 - 4Lx^2] = -\frac{w_0}{24EI}[4x^3 - 12Lx^2 + 12L^2x]

Second derivative: d2ydx2=−w024EI[12x2−24Lx+12L2]=−w02EI[x2−2Lx+L2]=−w02EI(x−L)2\frac{d^2y}{dx^2} = -\frac{w_0}{24EI}[12x^2 - 24Lx + 12L^2] = -\frac{w_0}{2EI}[x^2 - 2Lx + L^2] = -\frac{w_0}{2EI}(x - L)^2

Therefore, M(x)=EI⋅(−w02EI(x−L)2)=−w02(x−L)2M(x) = EI \cdot \left(-\frac{w_0}{2EI}(x - L)^2\right) = -\frac{w_0}{2}(x - L)^2

At x=2x = 2 m with L=4L = 4 m:

M(2)=−10×1032(2−4)2=−100002×4=−20,000 N⋅m=−20 kN⋅mM(2) = -\frac{10 \times 10^3}{2}(2 - 4)^2 = -\frac{10000}{2} \times 4 = -20,000\text{ N}\cdot\text{m} = -20\text{ kN}\cdot\text{m}

The negative sign indicates the beam is in hogging (negative moment) at x=2x = 2 m from the cantilever support.

Key points

  • Second derivative $f''(x)$ measures concavity; $f'' > 0$ is concave up, $f'' < 0$ is concave down
  • Inflection points occur where f′′(x)=0f''(x) = 0 and the sign of $f''$ changes
  • Second derivative test: at critical point, $f'' > 0$ → minimum, $f'' < 0$ → maximum
  • In beam theory, M=EId2ydx2M = EI\frac{d^2y}{dx^2} relates bending moment to deflection curvature
  • Higher derivatives encode physical information: slope, curvature, shear, and load in structural systems

5. Maxima, Minima, and Optimization — Critical Points and Extrema

Finding extrema (maximum and minimum values of a function) is among the most important applications of calculus in engineering. Optimal structural designs minimize cost or weight while meeting safety constraints—a fundamental principle in NSCP 2015 and modern engineering practice.

Critical points: A critical point occurs at x=cx = c if:

  1. f′(c)=0f'(c) = 0 (stationary point), or
  2. $f'(c)$ does not exist (e.g., sharp corner, cusp).

Every local extremum in the interior of a domain must occur at a critical point.

First Derivative Test (always reliable):

  1. Find all critical points by solving f′(x)=0f'(x) = 0 (or locating where $f'$ is undefined).
  2. Examine the sign of $f'(x)$ on intervals around each critical point:

    • If $f'$ changes from positive to negative as xx increases through cc, then cc is a local maximum.
    • If $f'$ changes from negative to positive as xx increases through cc, then cc is a local minimum.
    • If $f'$ does not change sign, then cc is neither a maximum nor a minimum (saddle point or inflection).

Second Derivative Test (convenient, sometimes inconclusive):
At a critical point cc where f′(c)=0f'(c) = 0:

  • If $f''(c) > 0$, then cc is a local minimum.
  • If $f''(c) < 0$, then cc is a local maximum.
  • If f′′(c)=0f''(c) = 0, the test is inconclusive; use the first derivative test.

Global extrema on a closed interval $[a,b]$:
Evaluate ff at:

  • All critical points in $(a, b)$
  • The endpoints $f(a)$ and $f(b)$

The largest value is the global maximum, and the smallest is the global minimum. This is essential for constrained optimization in engineering design.

Optimization strategy:

  1. Identify the quantity to optimize and express it as a function of a single variable using constraint equations.
  2. Find the domain (physical bounds on the variable).
  3. Differentiate and set equal to zero to find critical points.
  4. Check endpoints if the domain is closed.
  5. Verify using the second derivative test or by examining sign changes in the first derivative.
  6. Interpret the result physically—does it make sense given the context?

Example context in civil engineering: Minimizing the cost of concrete in a rectangular foundation of fixed volume, or maximizing the load capacity of a beam of given cross-section—these are prototypical optimization problems solved via calculus.

Examples

Find the absolute maximum and minimum of f(x)=x3−3x2−9x+5f(x) = x^3 - 3x^2 - 9x + 5 on the closed interval $[-2, 4]$.

Solution

Step 1: Find critical points.

f′(x)=3x2−6x−9=3(x2−2x−3)=3(x−3)(x+1)f'(x) = 3x^2 - 6x - 9 = 3(x^2 - 2x - 3) = 3(x - 3)(x + 1)

Setting f′(x)=0f'(x) = 0: x=3x = 3 or x=−1x = -1. Both lie in $[-2, 4]$.

Step 2: Evaluate ff at critical points and endpoints.

f(−2)=(−2)3−3(−2)2−9(−2)+5=−8−12+18+5=3f(-2) = (-2)^3 - 3(-2)^2 - 9(-2) + 5 = -8 - 12 + 18 + 5 = 3
f(−1)=(−1)3−3(−1)2−9(−1)+5=−1−3+9+5=10f(-1) = (-1)^3 - 3(-1)^2 - 9(-1) + 5 = -1 - 3 + 9 + 5 = 10
f(3)=(3)3−3(3)2−9(3)+5=27−27−27+5=−22f(3) = (3)^3 - 3(3)^2 - 9(3) + 5 = 27 - 27 - 27 + 5 = -22
f(4)=(4)3−3(4)2−9(4)+5=64−48−36+5=−15f(4) = (4)^3 - 3(4)^2 - 9(4) + 5 = 64 - 48 - 36 + 5 = -15

Step 3: Compare values.

  • Absolute maximum: f(−1)=10f(-1) = 10
  • Absolute minimum: f(3)=−22f(3) = -22

A rectangular garden is to be enclosed by 80 m of fencing, with one side bounded by an existing wall (no fence needed). Maximize the enclosed area.

Solution

Step 1: Set up variables. Let xx = width (perpendicular to wall), yy = length (parallel to wall).

Step 2: Constraint (fence equation): The three free sides require 2x+y=802x + y = 80, so y=80−2xy = 80 - 2x.

Step 3: Express area in one variable: A(x)=xy=x(80−2x)=80x−2x2A(x) = xy = x(80 - 2x) = 80x - 2x^2.

Step 4: Domain: Since x>0x > 0 and y>0y > 0, we need 80−2x>080 - 2x > 0, so 0<x<400 < x < 40.

Step 5: Find critical points.

A′(x)=80−4x=0⇒x=20A'(x) = 80 - 4x = 0 \Rightarrow x = 20

Step 6: Verify using second derivative test.

A′′(x)=−4<0⇒local maximum at x=20A''(x) = -4 < 0 \Rightarrow \text{local maximum at } x = 20

Step 7: Find corresponding dimensions and area.

x=20 m,y=80−2(20)=40 mx = 20\text{ m}, \quad y = 80 - 2(20) = 40\text{ m}
Amax⁡=20×40=800 m2A_{\max} = 20 \times 40 = 800\text{ m}^2

Interpretation: To maximize the garden area with 80 m of fencing and one side against a wall, make the width perpendicular to the wall 20 m and the length parallel to the wall 40 m.

Design a cylindrical storage tank (open top) with volume V=1000V = 1000 m³. Minimize the surface area of the metal used (base plus lateral surface). Assume radius rr and height hh.

Solution

Step 1: Volume constraint: V=πr2h=1000V = \pi r^2 h = 1000, so h=1000πr2h = \frac{1000}{\pi r^2}.

Step 2: Surface area (base + lateral, no top): S=πr2+2πrhS = \pi r^2 + 2\pi r h (base area + lateral area).

Step 3: Substitute hh:

S(r)=πr2+2πr⋅1000πr2=πr2+2000rS(r) = \pi r^2 + 2\pi r \cdot \frac{1000}{\pi r^2} = \pi r^2 + \frac{2000}{r}

Step 4: Find critical points.

S′(r)=2πr−2000r2=0S'(r) = 2\pi r - \frac{2000}{r^2} = 0
2πr=2000r2⇒2πr3=2000⇒r3=1000π⇒r=1000π3≈6.83 m2\pi r = \frac{2000}{r^2} \Rightarrow 2\pi r^3 = 2000 \Rightarrow r^3 = \frac{1000}{\pi} \Rightarrow r = \sqrt[3]{\frac{1000}{\pi}} \approx 6.83\text{ m}

Step 5: Verify with second derivative test.

S′′(r)=2π+4000r3>0(always positive)⇒minimumS''(r) = 2\pi + \frac{4000}{r^3} > 0 \quad \text{(always positive)} \Rightarrow \text{minimum}

Step 6: Find height.

h=1000πr2=1000π(6.83)2≈6.83 mh = \frac{1000}{\pi r^2} = \frac{1000}{\pi (6.83)^2} \approx 6.83\text{ m}

Note that the optimal design has h=rh = r—the height equals the radius. This is a classic result: a cylinder of fixed volume has minimum lateral + base surface area when height equals diameter.

Minimum surface area: S=π(6.83)2+20006.83≈146.4+292.8=439.2 m2S = \pi(6.83)^2 + \frac{2000}{6.83} \approx 146.4 + 292.8 = 439.2\text{ m}^2.

Key points

  • Critical points occur where f′(x)=0f'(x) = 0 or $f'(x)$ does not exist; all interior extrema occur at critical points
  • First derivative test (sign change) is always valid; second derivative test is convenient but sometimes inconclusive
  • Always check endpoints and domain boundaries—the global maximum/minimum may not occur at a critical point
  • Optimization strategy: express the quantity in one variable via constraint, differentiate, solve f′=0f' = 0, verify
  • Physical interpretation: verify that the critical point makes physical sense in the engineering context

6. Related Rates — Time-Dependent Optimization

Related rates problems involve quantities that change with respect to time and are connected by constraint equations. The goal is to find the rate of change of one quantity given the rates of others.

Strategy:

  1. Draw and label a diagram with all quantities and given rates.
  2. Identify the variables and constants. Variables usually change with time; constants do not.
  3. Write the constraint equation relating all variables (geometric formula, volume, etc.).
  4. Differentiate both sides with respect to time tt, treating all variables as functions of tt (use chain rule systematically).
  5. Substitute the known values (rates, positions) to solve for the unknown rate.
  6. Verify units and sign—positive rates mean increasing, negative means decreasing.

Key principle: Differentiate before substituting numerical values. Substitution after differentiation prevents errors from missing rate terms.

Common constraint equations in related rates:

  • Pythagorean theorem: x2+y2=c2x^2 + y^2 = c^2 (distances)
  • Circle area: A=πr2A = \pi r^2; circumference: C=2πrC = 2\pi r
  • Sphere volume: V=43πr3V = \frac{4}{3}\pi r^3
  • Cylinder volume: V=πr2hV = \pi r^2 h
  • Cone volume: V=13πr2hV = \frac{1}{3}\pi r^2 h
  • Similar triangles: x1y1=x2y2\frac{x_1}{y_1} = \frac{x_2}{y_2} (proportional dimensions)

Application in civil engineering: In traffic flow analysis, if the number of vehicles $N(t)$ passing a point increases at a known rate dNdt\frac{dN}{dt}, and the total distance traveled DD is related to NN by D=L⋅ND = L \cdot N (where LL is the route length), then the rate of distance increase is dDdt=LdNdt\frac{dD}{dt} = L \frac{dN}{dt}. Monitoring such rates helps engineers optimize traffic signal timing and estimate congestion dynamics.

Examples

A spherical balloon is inflated at a constant rate of dVdt=50\frac{dV}{dt} = 50 cm³/s. At what rate is the radius increasing when r=10r = 10 cm?

Solution

Step 1: Identify variables and rates.

  • $V(t)$ = volume (cm³)
  • $r(t)$ = radius (cm)
  • Given: dVdt=50\frac{dV}{dt} = 50 cm³/s
  • Find: drdt\frac{dr}{dt} when r=10r = 10 cm

Step 2: Write the volume constraint.

V=43πr3V = \frac{4}{3}\pi r^3

Step 3: Differentiate with respect to tt.

dVdt=43π⋅3r2drdt=4πr2drdt\frac{dV}{dt} = \frac{4}{3}\pi \cdot 3r^2 \frac{dr}{dt} = 4\pi r^2 \frac{dr}{dt}

Step 4: Solve for drdt\frac{dr}{dt}.

drdt=14πr2dVdt\frac{dr}{dt} = \frac{1}{4\pi r^2}\frac{dV}{dt}

Step 5: Substitute r=10r = 10 cm and dVdt=50\frac{dV}{dt} = 50 cm³/s.

drdt=14π(10)2×50=50400π=18π≈0.0398 cm/s\frac{dr}{dt} = \frac{1}{4\pi (10)^2} \times 50 = \frac{50}{400\pi} = \frac{1}{8\pi} \approx 0.0398\text{ cm/s}

Interpretation: When the radius is 10 cm, it is increasing at approximately 0.04 cm/s (or about 0.4 mm/s). As the balloon inflates, the radius grows more slowly as rr increases (because the surface area 4πr24\pi r^2 increases with rr).

A 5 m ladder leans against a vertical wall. The base slides away from the wall at $1$ m/s. When the base is 3 m from the wall, how fast is the top of the ladder sliding down the wall?

Solution

Step 1: Draw a diagram. Let xx = distance from base to wall, yy = height of ladder on wall. The ladder has fixed length ℓ=5\ell = 5 m.

Step 2: Identify rates.

  • Given: dxdt=1\frac{dx}{dt} = 1 m/s (base moving away)
  • Find: dydt\frac{dy}{dt} when x=3x = 3 m (rate at top moving down)

Step 3: Write the Pythagorean constraint.

x2+y2=25x^2 + y^2 = 25

Step 4: Differentiate with respect to tt.

2xdxdt+2ydydt=02x\frac{dx}{dt} + 2y\frac{dy}{dt} = 0
xdxdt+ydydt=0x\frac{dx}{dt} + y\frac{dy}{dt} = 0

Step 5: Solve for dydt\frac{dy}{dt}.

dydt=−xydxdt\frac{dy}{dt} = -\frac{x}{y}\frac{dx}{dt}

Step 6: Find yy when x=3x = 3 m.

32+y2=25⇒y2=16⇒y=4 m3^2 + y^2 = 25 \Rightarrow y^2 = 16 \Rightarrow y = 4\text{ m}

Step 7: Substitute.

dydt=−34×1=−0.75 m/s\frac{dy}{dt} = -\frac{3}{4} \times 1 = -0.75\text{ m/s}

Interpretation: The top of the ladder is sliding down at 0.75 m/s (negative sign indicates downward motion). Note that as the base gets farther from the wall, xx increases and yy decreases, so dydt\frac{dy}{dt} is negative. The rate increases in magnitude as xx approaches 5 m (the top approaches the ground faster near the end).

Water flows into a conical tank (vertex down) at $2$ m³/min. The tank has a fixed height of $4$ m and base radius of $2$ m. How fast is the water level rising when the water depth is $1$ m?

Solution

Step 1: Identify variables. Let $h(t)$ = depth of water, $r(t)$ = radius of water surface, $V(t)$ = volume of water.

  • Given: dVdt=2\frac{dV}{dt} = 2 m³/min
  • Find: dhdt\frac{dh}{dt} when h=1h = 1 m

Step 2: Set up the cone. The tank has fixed height 4 m and base radius 2 m. By similar triangles, the water radius at depth hh is:

rh=24=12⇒r=h2\frac{r}{h} = \frac{2}{4} = \frac{1}{2} \Rightarrow r = \frac{h}{2}

Step 3: Write the volume equation (cone).

V=13πr2h=13π(h2)2h=πh312V = \frac{1}{3}\pi r^2 h = \frac{1}{3}\pi \left(\frac{h}{2}\right)^2 h = \frac{\pi h^3}{12}

Step 4: Differentiate with respect to tt.

dVdt=π12⋅3h2dhdt=πh24dhdt\frac{dV}{dt} = \frac{\pi}{12} \cdot 3h^2 \frac{dh}{dt} = \frac{\pi h^2}{4}\frac{dh}{dt}

Step 5: Solve for dhdt\frac{dh}{dt}.

dhdt=4πh2dVdt\frac{dh}{dt} = \frac{4}{\pi h^2}\frac{dV}{dt}

Step 6: Substitute h=1h = 1 m and dVdt=2\frac{dV}{dt} = 2 m³/min.

dhdt=4π(1)2×2=8π≈2.546 m/min\frac{dh}{dt} = \frac{4}{\pi (1)^2} \times 2 = \frac{8}{\pi} \approx 2.546\text{ m/min}

Interpretation: When the water is 1 m deep, the water level rises at approximately 2.55 m/min. As the tank fills, hh increases, and h2h^2 increases, so dhdt\frac{dh}{dt} decreases—the water level rises more slowly as the tank fills (the cone widens).

Key points

  • Draw a clear diagram and label all quantities, including rates of change
  • Write the constraint equation before differentiating
  • Differentiate with respect to time tt, applying the chain rule to all variables
  • Substitute numerical values after differentiating—never before
  • Check that units match and the sign makes physical sense (positive = increasing, negative = decreasing)

7. Tangent Lines and Linear Approximation

The tangent line to a curve y=f(x)y = f(x) at a point (x0,f(x0))(x_0, f(x_0)) is a straight line that touches the curve at that point and has the same slope as the curve there.

Equation of the tangent line:

y−f(x0)=f′(x0)(x−x0)y - f(x_0) = f'(x_0)(x - x_0)

Alternatively, in slope-intercept form:

y=f(x0)+f′(x0)(x−x0)y = f(x_0) + f'(x_0)(x - x_0)

The tangent line is the best linear approximation to ff near x0x_0. For values of xx close to x0x_0:

f(x)≈f(x0)+f′(x0)(x−x0)f(x) \approx f(x_0) + f'(x_0)(x - x_0)

This is called linear approximation or differential approximation. The error E=f(x)−[f(x0)+f′(x0)(x−x0)]E = f(x) - [f(x_0) + f'(x_0)(x - x_0)] is small when xx is near x0x_0 and becomes negligible for very small intervals (e.g., small measurement uncertainties).

Differential notation: The differential dydy represents the change in the tangent line over an interval dxdx:

dy=f′(x) dxdy = f'(x)\,dx

This is useful for estimating small changes in yy due to small changes in xx. For example, if a structural parameter changes slightly, engineers can estimate the resulting change in response without full recalculation.

Normal line: The line perpendicular to the tangent at (x0,f(x0))(x_0, f(x_0)) has slope −1f′(x0)-\frac{1}{f'(x_0)} (negative reciprocal) and equation:

y−f(x0)=−1f′(x0)(x−x0)y - f(x_0) = -\frac{1}{f'(x_0)}(x - x_0)

Application in engineering: In structural analysis, if a load parameter varies slightly from its design value, the engineer can use linear approximation to quickly estimate the change in deflection or stress without re-solving the full problem. This enables sensitivity analysis and uncertainty quantification.

Examples

Find the equation of the tangent line to y=x2−3xy = x^2 - 3x at the point $(2, -2)$.

Solution

Step 1: Verify the point is on the curve.

y(2)=(2)2−3(2)=4−6=−2✓y(2) = (2)^2 - 3(2) = 4 - 6 = -2 \quad \checkmark

Step 2: Find the derivative.

y′=2x−3y' = 2x - 3

Step 3: Evaluate the slope at x=2x = 2.

y′(2)=2(2)−3=1y'(2) = 2(2) - 3 = 1

Step 4: Write the tangent line.

y−(−2)=1(x−2)y - (-2) = 1(x - 2)
y+2=x−2y + 2 = x - 2
y=x−4y = x - 4

Use linear approximation to estimate 25.1\sqrt{25.1} based on f(x)=xf(x) = \sqrt{x} at x0=25x_0 = 25.

Solution

Step 1: Identify the function and point.

f(x)=x,x0=25,f(25)=5f(x) = \sqrt{x}, \quad x_0 = 25, \quad f(25) = 5

Step 2: Find the derivative.

f′(x)=12x,f′(25)=1225=110f'(x) = \frac{1}{2\sqrt{x}}, \quad f'(25) = \frac{1}{2\sqrt{25}} = \frac{1}{10}

Step 3: Use linear approximation with x=25.1x = 25.1, so Δx=0.1\Delta x = 0.1.

f(x)≈f(x0)+f′(x0)Δx=5+110(0.1)=5+0.01=5.01f(x) \approx f(x_0) + f'(x_0)\Delta x = 5 + \frac{1}{10}(0.1) = 5 + 0.01 = 5.01

Step 4: Check with a calculator.

25.1≈5.00999≈5.01\sqrt{25.1} \approx 5.00999 \approx 5.01

The linear approximation is accurate to two decimal places. The small error arises because the second derivative (curvature) slightly bends the curve.

A cylindrical pipe has radius r=10r = 10 cm. Due to corrosion, the radius decreases by dr=0.1dr = 0.1 cm. Estimate the decrease in cross-sectional area using a differential.

Solution

Step 1: Set up the area function.

A(r)=πr2A(r) = \pi r^2

Step 2: Find the differential.

dA=dAdr dr=2πr drdA = \frac{dA}{dr}\,dr = 2\pi r \, dr

Step 3: Substitute r=10r = 10 cm and dr=−0.1dr = -0.1 cm (decrease).

dA=2π(10)(−0.1)=−2π≈−6.28 cm2dA = 2\pi (10)(-0.1) = -2\pi \approx -6.28\text{ cm}^2

Step 4: Interpret.
The cross-sectional area decreases by approximately 6.28 cm² when the radius decreases by 0.1 cm.

Verification (exact): A(10)=100πA(10) = 100\pi, A(9.9)=(9.9)2π=98.01πA(9.9) = (9.9)^2\pi = 98.01\pi. The exact decrease is 100π−98.01π=1.99π≈6.25100\pi - 98.01\pi = 1.99\pi \approx 6.25 cm². The differential approximation gives 6.28 cm², a very close estimate (error < 0.5%).

Key points

  • The tangent line is the best linear approximation near a point; use it to estimate function values close to x0x_0
  • Tangent line equation: y=f(x0)+f′(x0)(x−x0)y = f(x_0) + f'(x_0)(x - x_0)
  • Differential dy=f′(x) dxdy = f'(x)\,dx estimates the change in yy for a small change dxdx
  • Linear approximation error is small for small ∣x−x0∣|x - x_0| and shrinks as x→x0x \to x_0
  • Normal line: perpendicular to tangent; slope is −1f′(x0)-\frac{1}{f'(x_0)}

8. Curve Sketching and Radius of Curvature

Complete curve sketching synthesizes derivatives, critical points, concavity, and asymptotes to create an accurate graph. This skill is valuable for visualizing structural responses, force distributions, and optimization problems.

Curve-sketching procedure:

  1. Domain: Identify where the function is defined.
  2. Intercepts: Find xx-intercepts (set y=0y = 0) and yy-intercepts (set x=0x = 0).
  3. Asymptotes:

    • Vertical: Denominator = 0 (check left and right limits)
    • Horizontal: lim⁡x→±∞f(x)\lim_{x \to \pm\infty} f(x)
    • Oblique (slant): Polynomial long division when degree of numerator exceeds denominator by 1
  4. First derivative: Find critical points (f′(x)=0f'(x) = 0), determine intervals of increase/decrease.
  5. Second derivative: Find inflection points (f′′(x)=0f''(x) = 0), determine concavity.
  6. Sketch: Plot intercepts, critical points, inflection points, and asymptotes; connect with smooth curves respecting concavity and monotonicity.

Radius of curvature: The curvature κ\kappa measures how sharply a curve bends. The radius of curvature RR is the reciprocal:

κ=∣f′′(x)∣[1+(f′(x))2]3/2,R=1κ=[1+(f′(x))2]3/2∣f′′(x)∣\kappa = \frac{|f''(x)|}{[1 + (f'(x))^2]^{3/2}}, \quad R = \frac{1}{\kappa} = \frac{[1 + (f'(x))^2]^{3/2}}{|f''(x)|}

Geometrically, RR is the radius of the osculating circle (best-fitting circle) at a point on the curve.

Physical meaning:

  • Large RR → gentle curve (small curvature)
  • Small RR → sharp curve (high curvature)
  • At inflection points where f′′(x)=0f''(x) = 0, the radius becomes infinite (curvature is zero).

Application in civil engineering: Road and railway alignment design use curvature formulas. For a highway curve, the minimum radius of curvature is constrained by speed and friction to prevent skidding: Rmin⁡=v2g(e+f)R_{\min} = \frac{v^2}{g(e + f)}, where vv is speed, ee is banking angle, ff is friction coefficient, and gg is gravity. The radius of curvature of the designed centerline must exceed this minimum everywhere.

Examples

Sketch the curve y=xx2−1y = \frac{x}{x^2 - 1}.

Solution

Step 1: Domain. The denominator x2−1=0x^2 - 1 = 0 at x=±1x = \pm 1, so domain is (−∞,−1)∪(−1,1)∪(1,∞)(-\infty, -1) \cup (-1, 1) \cup (1, \infty).

Step 2: Intercepts.

  • yy-intercept: x=0⇒y=0x = 0 \Rightarrow y = 0. Point $(0, 0)$.
  • xx-intercept: y=0⇒x=0y = 0 \Rightarrow x = 0. Same point.

Step 3: Asymptotes.

  • Vertical: x=−1x = -1 and x=1x = 1 (denominator zero)

    • As x→−1−x \to -1^-: y→−∞y \to -\infty (negative numerator, small positive denom)
    • As x→−1+x \to -1^+: y→+∞y \to +\infty (negative numerator, small negative denom)
    • As x→1−x \to 1^-: y→+∞y \to +\infty (positive numerator, small negative denom)
    • As x→1+x \to 1^+: y→−∞y \to -\infty (positive numerator, small positive denom)
  • Horizontal: lim⁡x→±∞xx2−1=lim⁡x→±∞1/x1−1/x2=0\lim_{x \to \pm\infty}\frac{x}{x^2-1} = \lim_{x \to \pm\infty}\frac{1/x}{1 - 1/x^2} = 0. So y=0y = 0 is a horizontal asymptote.

Step 4: First derivative.

y′=(x2−1)−x(2x)(x2−1)2=−x2−1(x2−1)2=−(x2+1)(x2−1)2y' = \frac{(x^2-1) - x(2x)}{(x^2-1)^2} = \frac{-x^2 - 1}{(x^2-1)^2} = \frac{-(x^2+1)}{(x^2-1)^2}

The numerator −(x2+1)<0-(x^2+1) < 0 always, and the denominator (x2−1)2>0(x^2-1)^2 > 0 always. Thus $y' < 0$ everywhere the function is defined: the function is decreasing on each connected interval.

Step 5: Second derivative (for concavity). Differentiate y′=−(x2+1)(x2−1)2y' = \frac{-(x^2+1)}{(x^2-1)^2}. Using the quotient rule:

y′′=−2x(x2−1)2−(−(x2+1))⋅2(x2−1)⋅2x(x2−1)4y'' = \frac{-2x(x^2-1)^2 - (-(x^2+1)) \cdot 2(x^2-1) \cdot 2x}{(x^2-1)^4}

Simplify (factor (x2−1)(x^2-1) from numerator):

y′′=−2x(x2−1)+4x(x2+1)(x2−1)3=−2x3+2x+4x3+4x(x2−1)3=2x3+6x(x2−1)3=2x(x2+3)(x2−1)3y'' = \frac{-2x(x^2-1) + 4x(x^2+1)}{(x^2-1)^3} = \frac{-2x^3 + 2x + 4x^3 + 4x}{(x^2-1)^3} = \frac{2x^3 + 6x}{(x^2-1)^3} = \frac{2x(x^2+3)}{(x^2-1)^3}

Sign analysis:

  • x2+3>0x^2 + 3 > 0 always
  • Sign of $y''$ depends on xx and (x2−1)3(x^2-1)^3:

    • $x < -1$: x<0x < 0, (x2−1)3<0(x^2-1)^3 < 0 → $y'' > 0$ (concave up)
    • $-1 < x < 0$: x<0x < 0, (x2−1)3>0(x^2-1)^3 > 0 → $y'' < 0$ (concave down)
    • 0<x<10 < x < 1: x>0x > 0, (x2−1)3<0(x^2-1)^3 < 0 → $y'' < 0$ (concave down)
    • x>1x > 1: x>0x > 0, (x2−1)3>0(x^2-1)^3 > 0 → $y'' > 0$ (concave up)

Inflection point at x=0x = 0 (where the concavity changes).

Step 6: Sketch.

  • Plot the point $(0, 0)$.
  • Draw vertical asymptotes at x=−1x = -1 and x=1x = 1.
  • Draw horizontal asymptote at y=0y = 0.
  • On (−∞,−1)(-\infty, -1): decreasing, concave up, approaches y=0y = 0 as x→−∞x \to -\infty, approaches −∞-\infty as x→−1−x \to -1^-.
  • On $(-1, 1)$: passes through $(0, 0)$, decreasing, concave down for x<0x < 0, concave down for 0<x<10 < x < 1. Approaches +∞+\infty as x→−1+x \to -1^+ and +∞+\infty as x→1−x \to 1^-.
  • On (1,∞)(1, \infty): decreasing, concave up, approaches −∞-\infty as x→1+x \to 1^+, approaches $0$ from below as x→∞x \to \infty.

Find the radius of curvature of y=x2y = x^2 at x=0x = 0 and at x=1x = 1.

Solution

Step 1: Find derivatives.

y′=2x,y′′=2y' = 2x, \quad y'' = 2

Step 2: Apply the radius of curvature formula.

R=[1+(y′)2]3/2∣y′′∣=[1+(2x)2]3/22=[1+4x2]3/22R = \frac{[1 + (y')^2]^{3/2}}{|y''|} = \frac{[1 + (2x)^2]^{3/2}}{2} = \frac{[1 + 4x^2]^{3/2}}{2}

At x=0x = 0:

R(0)=[1+0]3/22=12=0.5R(0) = \frac{[1 + 0]^{3/2}}{2} = \frac{1}{2} = 0.5

The osculating circle at the origin has radius 0.5 (the parabola is sharply curved at the vertex).

At x=1x = 1:

R(1)=[1+4]3/22=(5)3/22=552≈11.182≈5.59R(1) = \frac{[1 + 4]^{3/2}}{2} = \frac{(5)^{3/2}}{2} = \frac{5\sqrt{5}}{2} \approx \frac{11.18}{2} \approx 5.59

At $(1, 1)$, the curve is gentler; the osculating circle has radius ≈ 5.59. This makes intuitive sense: as we move away from the vertex of the parabola, it becomes less sharply curved.

Key points

  • Follow a systematic curve-sketching procedure: domain, intercepts, asymptotes, derivatives, concavity, then sketch
  • Critical points (where f′(x)=0f'(x) = 0) and inflection points (where f′′(x)=0f''(x) = 0) are key features
  • Radius of curvature R=[1+(y′)2]3/2∣y′′∣R = \frac{[1+(y')^2]^{3/2}}{|y''|} quantifies the tightness of a curve
  • High curvature (small RR) indicates sharp bending; low curvature (large RR) indicates gentle bending
  • In road design, minimum radius of curvature is determined by vehicle dynamics and friction

9. Partial Derivatives and Multivariable Optimization

Many engineering problems involve functions of multiple independent variables. For example, the deflection of a beam depends on load, span length, material properties, and geometry—all varying independently.

Partial derivative notation: The partial derivative of $f(x, y)$ with respect to xx, denoted ∂f∂x\frac{\partial f}{\partial x} or fxf_x, is the derivative of ff with respect to xx alone, treating all other variables as constants.

∂f∂x=lim⁡h→0f(x+h,y)−f(x,y)h\frac{\partial f}{\partial x} = \lim_{h \to 0}\frac{f(x+h, y) - f(x, y)}{h}

Similarly, ∂f∂y=fy\frac{\partial f}{\partial y} = f_y is the derivative with respect to yy, treating xx as a constant.

Second-order partial derivatives:

  • ∂2f∂x2=fxx\frac{\partial^2 f}{\partial x^2} = f_{xx}: differentiate twice with respect to xx
  • ∂2f∂y2=fyy\frac{\partial^2 f}{\partial y^2} = f_{yy}: differentiate twice with respect to yy
  • ∂2f∂x∂y=fxy\frac{\partial^2 f}{\partial x \partial y} = f_{xy}: differentiate first with respect to xx, then yy
  • ∂2f∂y∂x=fyx\frac{\partial^2 f}{\partial y \partial x} = f_{yx}: differentiate first with respect to yy, then xx

Schwarz's theorem: If ff and its partial derivatives are continuous, then fxy=fyxf_{xy} = f_{yx} (mixed partials are equal).

Critical points for extrema: A point (x0,y0)(x_0, y_0) is a critical point if fx(x0,y0)=0f_x(x_0, y_0) = 0 and fy(x0,y0)=0f_y(x_0, y_0) = 0 simultaneously.

Second Derivative Test (2D Hessian): At a critical point (x0,y0)(x_0, y_0), define

D=fxx(x0,y0)⋅fyy(x0,y0)−[fxy(x0,y0)]2D = f_{xx}(x_0, y_0) \cdot f_{yy}(x_0, y_0) - [f_{xy}(x_0, y_0)]^2

Then:

  • If D>0D > 0 and fxx>0f_{xx} > 0: local minimum
  • If D>0D > 0 and fxx<0f_{xx} < 0: local maximum
  • If D<0D < 0: saddle point (neither maximum nor minimum)
  • If D=0D = 0: test is inconclusive

Application in civil engineering: In foundation design optimization, the settlement SS depends on soil cohesion cc and angle of internal friction ϕ\phi. Minimizing cost C(c,ϕ)C(c, \phi) subject to settlement limits requires finding critical points via partial derivatives. This approach scales to problems with dozens of design variables.

Examples

Find the critical points and classify them for f(x,y)=x2+y2−4x+6y+5f(x, y) = x^2 + y^2 - 4x + 6y + 5.

Solution

Step 1: Find first partial derivatives.

fx=2x−4,fy=2y+6f_x = 2x - 4, \quad f_y = 2y + 6

Step 2: Solve fx=0f_x = 0 and fy=0f_y = 0.

2x−4=0⇒x=22x - 4 = 0 \Rightarrow x = 2
2y+6=0⇒y=−32y + 6 = 0 \Rightarrow y = -3

The critical point is $(2, -3)$.

Step 3: Find second partial derivatives.

fxx=2,fyy=2,fxy=0f_{xx} = 2, \quad f_{yy} = 2, \quad f_{xy} = 0

Step 4: Apply the second derivative test.

D=fxxfyy−(fxy)2=(2)(2)−02=4>0D = f_{xx} f_{yy} - (f_{xy})^2 = (2)(2) - 0^2 = 4 > 0
fxx=2>0f_{xx} = 2 > 0

Since D>0D > 0 and fxx>0f_{xx} > 0, the point $(2, -3)$ is a local minimum.

Step 5: Find the minimum value.

f(2,−3)=4+9−8−18+5=−8f(2, -3) = 4 + 9 - 8 - 18 + 5 = -8

Interpretation: The paraboloid f(x,y)=(x−2)2+(y+3)2−8f(x, y) = (x-2)^2 + (y+3)^2 - 8 has its vertex (global minimum) at $(2, -3, -8)$.

A rectangular box is to be designed with volume V=100V = 100 m³. The bottom and top cost $5/m², the four sides cost $3/m². Find the dimensions that minimize the total cost.

Solution

Step 1: Set up variables. Let xx = length, yy = width, hh = height.

Step 2: Volume constraint: V=xyh=100V = xyh = 100, so h=100xyh = \frac{100}{xy}.

Step 3: Cost function. Bottom and top: 2xy2xy at $5/m² = $10xy$. Four sides: 2(xh+yh)=2h(x+y)2(xh + yh) = 2h(x+y) at $3/m² = $6h(x+y)$.

C=10xy+6h(x+y)=10xy+6⋅100xy⋅(x+y)=10xy+600(x+y)xyC = 10xy + 6h(x+y) = 10xy + 6\cdot\frac{100}{xy}\cdot(x+y) = 10xy + \frac{600(x+y)}{xy}

Simplify:

C(x,y)=10xy+600y+600xC(x,y) = 10xy + \frac{600}{y} + \frac{600}{x}

Step 4: Find critical points.

∂C∂x=10y−600x2=0⇒10y=600x2⇒y=60x2\frac{\partial C}{\partial x} = 10y - \frac{600}{x^2} = 0 \Rightarrow 10y = \frac{600}{x^2} \Rightarrow y = \frac{60}{x^2}
∂C∂y=10x−600y2=0⇒10x=600y2⇒x=60y2\frac{\partial C}{\partial y} = 10x - \frac{600}{y^2} = 0 \Rightarrow 10x = \frac{600}{y^2} \Rightarrow x = \frac{60}{y^2}

Step 5: Solve simultaneously. Substitute the first into the second:

x=60(60/x2)2=60⋅x43600=x460x = \frac{60}{(60/x^2)^2} = \frac{60 \cdot x^4}{3600} = \frac{x^4}{60}
60x=x4⇒x3=60⇒x=603≈3.91 m60x = x^4 \Rightarrow x^3 = 60 \Rightarrow x = \sqrt[3]{60} \approx 3.91\text{ m}

From y=60x2y = \frac{60}{x^2}:

y=60(603)2=60602/3=601−2/3=601/3≈3.91 my = \frac{60}{(\sqrt[3]{60})^2} = \frac{60}{60^{2/3}} = 60^{1-2/3} = 60^{1/3} \approx 3.91\text{ m}

From h=100xyh = \frac{100}{xy}:

h=100(603)2=100602/3=10036003≈6.57 mh = \frac{100}{(\sqrt[3]{60})^2} = \frac{100}{60^{2/3}} = \frac{100}{\sqrt[3]{3600}} \approx 6.57\text{ m}

Step 6: Verify with second derivative test (or note that cost is convex for positive dimensions). The critical point is a minimum by physical reasoning (cost increases as dimensions approach 0 or ∞).

Optimal dimensions: x≈y≈3.91x \approx y \approx 3.91 m, h≈6.57h \approx 6.57 m. Note that the optimal base is square; this is a common feature in many optimization problems due to symmetry.

Key points

  • Partial derivatives treat one variable as changing; hold other variables constant during differentiation
  • Find critical points by solving ∂f∂x=0\frac{\partial f}{\partial x} = 0 and ∂f∂y=0\frac{\partial f}{\partial y} = 0 simultaneously
  • Second Derivative Test (Hessian): D=fxxfyy−fxy2D = f_{xx} f_{yy} - f_{xy}^2 classifies critical points (local min, max, or saddle)
  • Schwarz's theorem: mixed partial derivatives are equal if the function is sufficiently smooth
  • Multivariable optimization: reduce constraints and solve systems of partial derivative equations

Diagrams

Differential Calculus Concept Map

This mind map organizes the entire differential calculus curriculum into six major branches: foundational concepts (limits and the derivative), computational rules, practical applications (optimization and related rates), curve sketching, advanced derivatives, and multivariable extension. The hierarchical structure helps reviewees see how topics interconnect and prioritize study time.

Optimization Problem-Solving Workflow

This flowchart provides a step-by-step decision tree for solving optimization problems. It emphasizes the key decision points (finding critical points, classifying extrema, comparing candidates) and the importance of verifying that the mathematical solution makes physical sense. This workflow is applicable to virtually all board-exam optimization questions.

Related Rates Problem-Solving Strategy

This flowchart emphasizes the critical principle: differentiate before substituting numerical values. Common errors arise when students substitute too early, missing rate terms. The diagram shows the iterative verification steps (checking units and physical interpretation) that distinguish confident solutions from rushed ones.

Concavity and Inflection Point Analysis State Diagram

This state diagram illustrates the logic flow when analyzing concavity and identifying inflection points. It emphasizes that an inflection point requires both f′′(x)=0f''(x) = 0 and a sign change in $f''$. The diagram helps reviewees avoid the common error of assuming every zero of $f''$ is an inflection point.

In summary

Differential calculus is the cornerstone of quantitative analysis in civil engineering. This study module covered:

  1. Limits and Continuity — The foundation of calculus, enabling the definition of derivatives and reliable function behavior analysis. L'Hôpital's Rule resolves indeterminate forms efficiently.
  2. The Derivative — The instantaneous rate of change, geometrically interpreted as the slope of a tangent line. Differentiability implies continuity but is stronger (requires smoothness).
  3. Derivative Rules — Power, product, quotient, and chain rules provide systematic computation without the limit definition. Memorizing standard derivatives of trig, exponential, and logarithmic functions is essential for board exams.
  4. Higher-Order Derivatives — The second derivative measures concavity; in structural mechanics, it directly relates to bending moment (M=EIy′′M = EI y''). Inflection points and the second derivative test classify extrema efficiently.
  5. Optimization and Extrema — Critical to engineering design: express the objective in one variable via constraints, set f′(x)=0f'(x) = 0, verify with the second derivative test, and check domain endpoints. This workflow solves problems ranging from minimum-cost materials to maximum load capacity.
  6. Related Rates — Time-dependent quantities linked by constraint equations. Differentiate before substituting; this avoids missing rate terms and is the most common source of student error.
  7. Tangent Lines and Linear Approximation — The tangent line y=f(x0)+f′(x0)(x−x0)y = f(x_0) + f'(x_0)(x - x_0) enables rapid estimation of function values near a known point. The differential dy=f′(x) dxdy = f'(x)\,dx quantifies small changes and is crucial for sensitivity analysis.
  8. Curve Sketching and Curvature — Systematic graphing (intercepts, asymptotes, critical points, concavity) synthesizes all derivative information. The radius of curvature R=[1+(y′)2]3/2∣y′′∣R = \frac{[1+(y')^2]^{3/2}}{|y''|} is essential for road and railway alignment design.
  9. Partial Derivatives and Multivariable Optimization — For functions of multiple variables, partial derivatives generalize single-variable calculus. The Hessian matrix (second derivative test in 2D) classifies critical points. This framework extends to high-dimensional optimization in structural and geotechnical design.

Exam Strategy:

  • Master the standard derivatives and derivative rules; they appear in nearly every problem.
  • For optimization: reduce to one variable using constraints, solve f′(x)=0f'(x) = 0, verify with $f''$, and check boundaries.
  • For related rates: sketch a diagram, write the constraint, differentiate with respect to tt, then substitute.
  • For indeterminate limits: try factoring or L'Hôpital's Rule (differentiate numerator and denominator separately).
  • Always verify units, interpret the sign (increasing vs. decreasing), and check physical reasonableness.

Key Formulas to Memorize:

  • Power rule: ddxxn=nxn−1\frac{d}{dx}x^n = nx^{n-1}
  • Product rule: (uv)′=u′v+uv′(uv)' = u'v + uv'
  • Quotient rule: (u/v)′=u′v−uv′v2(u/v)' = \frac{u'v - uv'}{v^2}
  • Chain rule: ddxf(g(x))=f′(g(x))⋅g′(x)\frac{d}{dx}f(g(x)) = f'(g(x)) \cdot g'(x)
  • Second derivative test: $f''(c) > 0$ → minimum, $f''(c) < 0$ → maximum
  • Radius of curvature: R=[1+(y′)2]3/2∣y′′∣R = \frac{[1+(y')^2]^{3/2}}{|y''|}
  • Tangent line: y=f(x0)+f′(x0)(x−x0)y = f(x_0) + f'(x_0)(x - x_0)

These concepts appear repeatedly on the MSTE section of the PRC Licensure Examination. Proficiency in differential calculus—particularly optimization and related rates—distinguishes strong candidates and prepares engineers for real-world design challenges.

Differential Calculus is chapter 5 of 10 in the Engineering Mathematics part of Super Tutor's CELE review. The Civil Engineer Licensure Examination is given by Professional Regulation Commission (PRC) — Board of Civil Engineering; its dates, format and passing rules are on the CELE page. See the CELE guide.

Study Differential Calculus in the Super Tutor app

Super Tutor's CELE review has every chapter in 11 study formats, practice quizzes at three difficulty levels and an AI tutor, plus full CELE mock exams. Start free — 5-day trial, no card.

Reviewing for another exam? Start free for: