Calculus I

One continuous manuscript. Worked solutions and resource URLs are included when printing; no interactive control is required to understand a chapter.

Study plan →

Mathematics Project · James Howard

Calculus I

Limits, derivatives, and accumulated change

First-pass textbook manuscript · November 2026

Original exposition and exercises. Independent mathematical and classroom review pending.

How to use this book

Read with pencil and paper. Reproduce the worked examples, then attempt the exercises before consulting the solutions. A proof should identify the assumptions that make each step legal. A modeled quantity should have units and a domain. The external readings and videos offer additional explanations; they do not replace the chapter.

The two online laboratories illustrate secant limits and a quadratic epsilon–delta bound. Their necessary mathematics is developed in the text, so this print edition stands independently. The manuscript assumes algebra, elementary functions, logarithms, and trigonometry in radians.

Chapter 1

Functions, Measurements, and Change

Domains, units, average rates, and the algebra that an instantaneous rate will require.

Prerequisites: Algebraic manipulation, graphs of elementary functions, and the coordinate plane.

  • Determine domains before manipulating a formula.
  • Calculate and interpret a difference quotient with units.
  • Distinguish average speed, average velocity, and an instantaneous rate.

A position sensor records where a cart is at successive times. It can tell us how far the cart moved between two readings. A speedometer appears to answer a more demanding question: how fast the cart is moving at one instant. Dividing distance by elapsed time seems to require two instants. Calculus will make that tension precise rather than disguise it with a sufficiently small number.

For a deliberately simple model, let the cart’s position be s(t)=t2s(t)=t^2 metres, where tt is measured in seconds and t≥0t\geq0. The coefficient of t2t^2 therefore has units metres per second squared. The formula is a model, not a claim about every cart. Between t=2t=2 and t=5t=5, the displacement is 25−4=2125-4=21 metres and the elapsed time is 33 seconds. The average velocity is 77 metres per second. It does not follow that the cart moved at that velocity throughout the interval.

Functions include their domains

A real-valued function assigns exactly one real output to each allowed input. Its domain is part of its specification. The formulas x2x^2 on R\mathbb R and t2t^2 on [0,∞)[0,\infty) use the same algebra but represent different functions when domains matter.

Three restrictions recur throughout this book. A denominator cannot be zero. The argument of a real square root must be nonnegative. The argument of a real logarithm must be positive. For example,

f(x)=x+2x−1has domain[−2,1)∪(1,∞).f(x)=\frac{\sqrt{x+2}}{x-1} \quad\hbox{has domain}\quad [-2,1)\cup(1,\infty).

Writing x+2\sqrt{x+2} requires x≥−2x\geq-2; division requires x≠1x\ne1. Both restrictions apply. A graphing program that silently omits the forbidden point has not removed the restriction.

Composition also carries restrictions. If g(x)=log⁡xg(x)=\log x and f(u)=1/(u−2)f(u)=1/(u-2), then (f∘g)(x)=1/(log⁡x−2)(f\circ g)(x)=1/(\log x-2) requires x>0x>0 and x≠e2x\ne e^2. We must check the input to each operation, in the order in which it is performed.

The difference quotient

For distinct inputs aa and a+ha+h, the average rate of change is

f(a+h)−f(a)h,h≠0.\frac{f(a+h)-f(a)}h,\qquad h\ne0.

The numerator is the output change and the denominator the input change. On a graph, this is the slope of the line through (a,f(a))(a,f(a)) and (a+h,f(a+h))(a+h,f(a+h)), called a secant line. If the output is position and the input is time, its units are distance per time. If output is cost and input is production, the units are currency per item. A slope without its variable definitions can be mathematically correct and practically useless.

For s(t)=t2s(t)=t^2, expand before simplifying:

s(a+h)−s(a)h=a2+2ah+h2−a2h=2a+h(h≠0).\frac{s(a+h)-s(a)}h =\frac{a^2+2ah+h^2-a^2}h=2a+h\qquad(h\ne0).

At a=2a=2, intervals of length 11, 0.10.1, and 0.010.01 give average velocities 55, 4.14.1, and 4.014.01. Backward intervals use negative hh: h=−0.1h=-0.1 compares times 1.91.9 and 22, giving 3.93.9. Both directions suggest 44. The original quotient remains undefined at h=0h=0, however. The simplified expression describes its behavior where the quotient is defined; it does not authorize division by zero.

Direction matters

Average velocity is displacement divided by elapsed time. Average speed is total distance travelled divided by elapsed time. A person who walks three metres east and three metres west in four seconds has zero average velocity and average speed 6/4=1.56/4=1.5 metres per second. The endpoint positions alone determine displacement; they do not determine total distance.

For a graph that is a straight line f(x)=mx+bf(x)=mx+b, every secant slope is mm. For a curved graph the interval matters. Calculus asks whether there is a well-defined local rate as the interval contracts. It does not begin by assuming that every curve has one.

Algebraic preparation

The identity (u+v)2=u2+2uv+v2(u+v)^2=u^2+2uv+v^2 explains the cancellation above. Factoring and conjugates provide two more tools. Since x2−a2=(x−a)(x+a)x^2-a^2=(x-a)(x+a), a quotient can simplify away from x=ax=a. Since (x−a)(x+a)=x−a(\sqrt{x}-\sqrt a)(\sqrt{x}+\sqrt a)=x-a, multiplying by a conjugate can expose a hidden factor. In either case, preserve the original domain while manipulating the expression.

A useful check is to calculate the quotient at one ordinary permitted input before and after simplifying. Agreement at one input is not proof of an identity, but disagreement immediately exposes an error.

Exercises

  1. Determine the real domain of 9−x2/(x−2)\sqrt{9-x^2}/(x-2).
  2. For p(t)=3t2−2t+1p(t)=3t^2-2t+1 metres, calculate average velocity from t=1t=1 to t=3t=3 seconds.
  3. Simplify [p(a+h)−p(a)]/h[p(a+h)-p(a)]/h and state its restriction.
  4. A tank contains V(t)=100−4tV(t)=100-4t litres for 0≤t≤250\leq t\leq25 minutes. Interpret its average rate over any nonzero interval in this domain.
  5. Explain why cancelling x−1x-1 in (x2−1)/(x−1)(x^2-1)/(x-1) does not define the original function at 11.
  6. A car returns to its starting position after an hour. What can be concluded about its average velocity and average speed?
Worked solutions
  1. The square root requires −3≤x≤3-3\leq x\leq3 and the denominator excludes 22, giving [−3,2)∪(2,3][-3,2)\cup(2,3].
  2. p(1)=2p(1)=2 and p(3)=22p(3)=22, so the average is (22−2)/(3−1)=10(22-2)/(3-1)=10 metres per second.
  3. Expansion gives 6a+3h−26a+3h-2, valid for h≠0h\ne0 and permitted times. The terms 3a2−2a+13a^2-2a+1 cancel, leaving h(6a+3h−2)h(6a+3h-2).
  4. The rate is −4-4 litres per minute. The negative sign indicates loss of volume. The model stops when the tank empties; extrapolating to negative volume is outside its domain.
  5. The two expressions agree for x≠1x\ne1. The expression x+1x+1 has a value at 11 that the original quotient lacks. Extending a function is a separate definition.
  6. Average velocity is zero. Average speed is nonnegative and positive if any nonzero distance was travelled. Its value requires the total distance.

Open this chapter with navigation and any associated laboratory →

Chapter 2

Limits and One-Sided Behavior

Separate a function's value from its nearby behavior, and recognize when a limit fails.

Prerequisites: Functions, domains, interval notation, and difference quotients.

  • Distinguish a limit from evaluation at a point.
  • Use one-sided limits to test existence of a two-sided limit.
  • Distinguish finite limits, unbounded behavior, and oscillation.

The secant calculation gave a quantity that was undefined at the input we wanted. This is not an unusual defect to patch. It is the reason we need a new concept. We want to describe what outputs do when inputs approach a point through other allowed inputs.

Consider

f(x)=x2−4x−2,x≠2.f(x)=\frac{x^2-4}{x-2},\qquad x\ne2.

Factoring shows that f(x)=x+2f(x)=x+2 at every input in its domain. As xx approaches 22, the outputs approach 44. We write lim⁡x→2f(x)=4\lim_{x\to2}f(x)=4. We have made a statement about values near 22, although f(2)f(2) is not defined.

A punctured neighborhood

The informal meaning of a finite limit LL is that outputs can be made arbitrarily close to LL by keeping inputs sufficiently close to aa, while excluding aa itself. The exclusion matters. A punctured neighborhood has the form 0<∣x−a∣<δ0<|x-a|<\delta. For the two-sided limits used here, our functions are defined throughout some such neighborhood. More general domain-relative limits are possible, but we will state explicitly when an endpoint requires a one-sided approach.

Now define a new function gg that agrees with ff for x≠2x\ne2 but has g(2)=100g(2)=100. Its limit at 22 is still 44. Altering one value does not change nearby behavior. On the other hand, altering all values to the right of 22 can change the limit even if the value at 22 stays fixed.

The assertion lim⁡x→af(x)=L\lim_{x\to a}f(x)=L also does not require f(x)f(x) to equal LL near aa. For f(x)=xf(x)=x, the limit at 00 is 00, but no nonzero input has output 00. Nor does it forbid equality: a constant function can equal its limit everywhere.

Left and right must agree

The notation x→a−x\to a^- restricts the approach to x<ax<a; x→a+x\to a^+ restricts it to x>ax>a. A finite two-sided limit exists exactly when the two finite one-sided limits exist and agree, under the neighborhood assumption above.

For a concrete jump, let q(x)=x+1q(x)=x+1 when x<0x<0 and q(x)=x+3q(x)=x+3 when x≥0x\geq0. Then

lim⁡x→0−q(x)=1,lim⁡x→0+q(x)=3.\lim_{x\to0^-}q(x)=1,\qquad \lim_{x\to0^+}q(x)=3.

No two-sided finite limit exists at 00. Choosing q(0)=1q(0)=1, 22, or 33 cannot repair this disagreement. There is no single output target that describes both sides.

At a domain endpoint, a one-sided question is often the appropriate one. For x\sqrt{x} on [0,∞)[0,\infty), the right-hand limit at 00 is 00. There is no real-valued graph to its left. We should state the right-hand limit rather than pretend we have checked both sides.

A table is evidence, not a guarantee

For f(x)=(x2−4)/(x−2)f(x)=(x^2-4)/(x-2), inputs 1.91.9, 1.991.99, 2.012.01, and 2.12.1 give 3.93.9, 3.993.99, 4.014.01, and 4.14.1. The factorization establishes what every permitted nearby input does; the table only illustrates a few inputs. A finite collection of samples cannot rule out a narrow spike or rapid oscillation between sampled points.

For example, sin⁡(1/x)\sin(1/x) has no limit as x→0x\to0. Along xn=1/(π/2+2πn)x_n=1/(\pi/2+2\pi n) its values are 11. Along yn=1/(3π/2+2πn)y_n=1/(3\pi/2+2\pi n) they are −1-1. Both input sequences tend to zero. A common finite limit would force both output sequences toward the same value, which is impossible.

Contrast this with xsin⁡(1/x)x\sin(1/x). Its oscillations are trapped in an envelope whose height tends to zero. That function does have limit zero. Oscillation by itself does not settle the question; its amplitude matters.

Infinity describes behavior

Writing lim⁡x→01/x2=+∞\lim_{x\to0}1/x^2=+\infty means the values eventually exceed any prescribed positive bound as nonzero inputs approach zero. It does not assign a real value called infinity. For 1/x1/x, the right-hand behavior is +∞+\infty and the left-hand behavior is −∞-\infty, so even their unbounded directions disagree.

Limits at infinity concern large inputs instead: lim⁡x→∞1/x=0\lim_{x\to\infty}1/x=0. The graph approaches the horizontal axis while never reaching it. A vertical asymptote concerns behavior near a finite input; a horizontal asymptote concerns behavior as inputs become unbounded. They answer different questions.

Exercises

  1. Let f(x)=x+5f(x)=x+5 for x≠1x\ne1 and f(1)=−8f(1)=-8. Find the limit at 11 and the function value there.
  2. Find both one-sided limits of ∣x∣/x|x|/x at zero. Does a two-sided limit exist?
  3. Explain why 1/x21/x^2 has no finite real limit at zero.
  4. Give two functions that have the same limit at 33 but different values at 33.
  5. Explain why a graph drawn from 1,000 sample points cannot establish a limit for every function.
  6. Find two sequences showing that cos⁡(1/x)\cos(1/x) has no limit at zero.
Worked solutions
  1. The limit is 66 because the nearby formula is x+5x+5; f(1)=−8f(1)=-8.
  2. The left limit is −1-1 and the right limit is 11, so no two-sided limit exists.
  3. Values exceed every real bound near zero. No finite candidate can remain close to those arbitrarily large outputs.
  4. Use f(x)=xf(x)=x everywhere and g(x)=xg(x)=x except g(3)=7g(3)=7. Both limits are 33.
  5. Sampling leaves unexamined inputs. A function can agree at all sampled points but behave differently between them or closer to the target.
  6. For positive integers nn, use xn=1/(2πn)x_n=1/(2\pi n) and yn=1/((2n+1)π)y_n=1/((2n+1)\pi). Outputs are 11 and −1-1 while both sequences approach zero.

Open this chapter with navigation and any associated laboratory →

Chapter 3

Computing Limits with Reasons

Apply limit laws, factoring, conjugates, squeezing, and dominant terms without discarding their assumptions.

Prerequisites: Finite and one-sided limits; algebraic factoring; trigonometric functions in radians.

  • Apply quotient and composition rules only when their hypotheses hold.
  • Resolve removable algebraic indeterminacies.
  • Use bounding arguments and distinguish growth at infinity.

Substitution is often the last line of a limit calculation, but it is not the definition of a limit. It works for polynomials because their limits obey algebraic rules; it works for a rational expression at points where the denominator has a nonzero limit. Recognizing those conditions lets us distinguish an ordinary calculation from one that needs more work.

The algebra of finite limits

Suppose lim⁡x→af(x)=L\lim_{x\to a}f(x)=L and lim⁡x→ag(x)=M\lim_{x\to a}g(x)=M, with L,ML,M finite. Then limits commute with sums, differences, constant multiples, and products. Thus lim⁡(f+g)=L+M\lim(f+g)=L+M and lim⁡(fg)=LM\lim(fg)=LM. For quotients we additionally require M≠0M\ne0, giving lim⁡(f/g)=L/M\lim(f/g)=L/M. Positive integer powers follow by repeated products.

For composition, if f(x)→Lf(x)\to L and FF is continuous at LL, then F(f(x))→F(L)F(f(x))\to F(L), provided the composition is defined nearby. Continuity will receive a careful definition shortly. Root functions are continuous on their real domains, with one-sided care at endpoints. Hence

lim⁡x→23x2−x+1x+4=116.\lim_{x\to2}\frac{3x^2-x+1}{x+4}=\frac{11}{6}.

The denominator approaches 66, not zero. That single check justifies the final division.

Zero over zero is a prompt to investigate

For (x2−9)/(x−3)(x^2-9)/(x-3) at 33, substitution gives the meaningless quotient 0/00/0. It does not give an answer of zero, infinity, or “undefined limit.” For x≠3x\ne3,

x2−9x−3=x+3,\frac{x^2-9}{x-3}=x+3,

so the limit is 66. The two expressions agree throughout a punctured neighborhood. Limits depend on exactly that nearby behavior.

Compare x/xx/x, x2/xx^2/x, and x/x2x/x^2 at zero. All display 0/00/0 on substitution, yet their behaviors are respectively a constant 11, a limit 00, and incompatible unbounded one-sided limits. The notation indeterminate form records insufficient information, not a result.

Conjugates handle some square-root cases. For x≠4x\ne4 near 44,

x−2x−4=1x+2,\frac{\sqrt{x}-2}{x-4} =\frac1{\sqrt{x}+2},

because multiplying numerator and denominator by x+2\sqrt{x}+2 exposes x−4x-4. The limit is 1/41/4. Multiplying by the conjugate is legal nearby because x+2\sqrt{x}+2 is nonzero.

Squeezing every nearby value

The squeeze theorem says that if u(x)≤f(x)≤v(x)u(x)\leq f(x)\leq v(x) for all sufficiently near allowed inputs, and both bounding functions approach LL, then f(x)f(x) also approaches LL. The inequalities must hold throughout the relevant neighborhood, not only at selected table entries.

Since ∣sin⁡(1/x)∣≤1|\sin(1/x)|\leq1, we have

−∣x∣≤xsin⁡(1/x)≤∣x∣.-|x|\leq x\sin(1/x)\leq |x|.

Both outer expressions tend to zero. This supplies a conclusion even though the middle function oscillates infinitely often near zero.

An important trigonometric limit is lim⁡x→0sin⁡x/x=1\lim_{x\to0}\sin x/x=1, with angles measured in radians. On 0<x<π/20<x<\pi/2, the unit-circle area comparison gives sin⁡x<x<tan⁡x\sin x<x<\tan x. Dividing yields cos⁡x<sin⁡x/x<1\cos x<\sin x/x<1. As x→0+x\to0^+, both bounds tend to 11; symmetry gives the left side. This argument uses elementary geometric continuity of cosine, not a derivative formula for sine. Also,

1−cos⁡xx=sin⁡xxsin⁡x1+cos⁡x⟶0.\frac{1-\cos x}{x}=\frac{\sin x}{x}\frac{\sin x}{1+\cos x}\longrightarrow0.

These results will justify trigonometric derivatives. Degrees would introduce a conversion factor; the radian convention is mathematical infrastructure, not a cosmetic preference.

Unbounded inputs and dominant terms

To study a rational expression as x→∞x\to\infty, divide numerator and denominator by a suitable highest power. For example,

5x2+2x−12x2+7=5+2/x−1/x22+7/x2⟶52.\frac{5x^2+2x-1}{2x^2+7} =\frac{5+2/x-1/x^2}{2+7/x^2}\longrightarrow\frac52.

If the denominator degree is larger, the ratio tends to zero. If the numerator degree is larger, polynomial growth remains; its sign and direction must still be checked. A memorized degree rule should summarize this calculation rather than replace it.

For radicals, x2=∣x∣\sqrt{x^2}=|x|, not always xx. Therefore x2+1/x\sqrt{x^2+1}/x tends to 11 as x→+∞x\to+\infty and −1-1 as x→−∞x\to-\infty. Forgetting the absolute value changes the second answer.

Exercises

  1. Evaluate lim⁡x→2(x2+x−6)/(x−2)\lim_{x\to2}(x^2+x-6)/(x-2).
  2. Evaluate lim⁡x→0(1+x−1)/x\lim_{x\to0}(\sqrt{1+x}-1)/x.
  3. Find lim⁡x→0sin⁡(5x)/x\lim_{x\to0}\sin(5x)/x and justify the factor 55.
  4. Use squeezing to evaluate lim⁡x→0x2cos⁡(1/x)\lim_{x\to0}x^2\cos(1/x).
  5. Find the limits of (3x3−x)/(x3+2)(3x^3-x)/(x^3+2) at positive and negative infinity.
  6. Explain why the quotient law cannot establish lim⁡x→0sin⁡x/x\lim_{x\to0}\sin x/x directly.
  7. Find lim⁡x→−∞4x2+1/x\lim_{x\to-\infty}\sqrt{4x^2+1}/x.
  8. Give an example showing that lim⁡(fg)\lim(fg) can exist even when neither factor has a finite limit.
Worked solutions
  1. Factor the numerator as (x+3)(x−2)(x+3)(x-2); the limit is 55.
  2. The conjugate gives 1/(1+x+1)1/(\sqrt{1+x}+1) near zero, with limit 1/21/2.
  3. Write 5[sin⁡(5x)/(5x)]5[\sin(5x)/(5x)]. The bracket tends to 11 in radians, giving 55.
  4. The absolute value is bounded by x2x^2, so the limit is zero.
  5. Dividing by x3x^3 gives 33 in either direction. The nonzero denominator limit justifies the quotient step.
  6. The denominator limit is zero, violating that law’s nonzero-denominator hypothesis.
  7. Write ∣x∣4+1/x2/x|x|\sqrt{4+1/x^2}/x. For negative xx, ∣x∣/x=−1|x|/x=-1, giving −2-2.
  8. Let f(x)=1/xf(x)=1/x and g(x)=xg(x)=x for x≠0x\ne0. Their product has limit 11 at zero, although the first factor has no finite limit there. To make neither factor converge, take f(x)=2+sin⁡(1/x)f(x)=2+\sin(1/x) and g(x)=1/f(x)g(x)=1/f(x); their product is identically 11 while both factors oscillate.

Open this chapter with navigation and any associated laboratory →

Chapter 4

The Epsilon–Delta Definition

Turn arbitrarily close into a quantified claim, construct delta, and prove that a proposed limit fails.

Prerequisites: Absolute-value inequalities, punctured neighborhoods, and elementary limit calculations.

  • State the quantifiers in the definition of a finite limit in their correct order.
  • Construct epsilon–delta proofs for linear and quadratic functions.
  • Negate the definition and establish uniqueness of a limit.

A numerical table can make a limit plausible. A proof must handle every input sufficiently near the target, including inputs that no table contains. We therefore need a statement that says exactly what “sufficiently near” must accomplish.

Let ff be defined throughout a punctured neighborhood of aa. We say lim⁡x→af(x)=L\lim_{x\to a}f(x)=L if

for every ε>0 there exists δ>0 such that 0<∣x−a∣<δ  ⟹  ∣f(x)−L∣<ε.\text{for every }\varepsilon>0\text{ there exists }\delta>0 \text{ such that }0<|x-a|<\delta \implies |f(x)-L|<\varepsilon.

Here ε\varepsilon is the requested output tolerance and δ\delta is the promised input tolerance. We receive an arbitrary positive ε\varepsilon, then choose a suitable δ\delta, then guarantee the implication for every permitted xx. A choice of δ\delta may depend on ε\varepsilon, aa, and the fixed function. It must not depend on an input xx selected afterward.

A linear proof from start to finish

Prove lim⁡x→2(3x−1)=5\lim_{x\to2}(3x-1)=5. The output error is

∣(3x−1)−5∣=3∣x−2∣.|(3x-1)-5|=3|x-2|.

For a given ε>0\varepsilon>0, choose δ=ε/3\delta=\varepsilon/3. If 0<∣x−2∣<δ0<|x-2|<\delta, then ∣(3x−1)−5∣=3∣x−2∣<3δ=ε|(3x-1)-5|=3|x-2|<3\delta=\varepsilon. Since this works for every positive tolerance, the proof is complete.

Notice the distinction between scratch work and the proof. We discovered δ\delta by working backward from the desired inequality. We presented the proof forward, starting with an arbitrary tolerance and the explicit choice. A successful proof does not require the largest possible δ\delta. Any positive choice that ensures the implication works.

A quadratic needs a local bound

To prove lim⁡x→2x2=4\lim_{x\to2}x^2=4, factor the error:

∣x2−4∣=∣x−2∣∣x+2∣.|x^2-4|=|x-2||x+2|.

The second factor depends on xx. We first arrange ∣x−2∣<1|x-2|<1, which gives 1<x<31<x<3 and hence ∣x+2∣<5|x+2|<5. Now choose

δ=min⁡{1,ε/5}.\delta=\min\{1,\varepsilon/5\}.

If 0<∣x−2∣<δ0<|x-2|<\delta, then ∣x2−4∣<5∣x−2∣<5δ≤ε|x^2-4|<5|x-2|<5\delta\leq\varepsilon. The final strict inequality comes from the strict bound on ∣x−2∣|x-2|. The minimum enforces two requirements simultaneously: remain in a neighborhood where the auxiliary factor is controlled, and make the final error small.

For a general fixed real aa, ∣x−a∣<1|x-a|<1 implies ∣x+a∣<2∣a∣+1|x+a|<2|a|+1. Thus δ=min⁡{1,ε/(2∣a∣+1)}\delta=\min\{1,\varepsilon/(2|a|+1)\} proves lim⁡x→ax2=a2\lim_{x\to a}x^2=a^2.

The largest symmetric tolerance for a square

The laboratory uses an exact bound to make the geometry visible. Put x=a+hx=a+h. For ∣h∣<δ|h|<\delta,

∣x2−a2∣=∣2ah+h2∣≤2∣a∣∣h∣+h2<2∣a∣δ+δ2.|x^2-a^2|=|2ah+h^2| \leq2|a||h|+h^2 <2|a|\delta+\delta^2.

The supremum of the output error on this open neighborhood is 2∣a∣δ+δ22|a|\delta+\delta^2, approached from the outward direction. A sufficient and, for the symmetric neighborhood, maximal choice is

δ∗=a2+ε−∣a∣=εa2+ε+∣a∣.\delta_* = \sqrt{a^2+\varepsilon}-|a| =\frac{\varepsilon}{\sqrt{a^2+\varepsilon}+|a|}.

The second form avoids subtracting nearly equal numbers in computation. At a=0a=0, this becomes ε\sqrt\varepsilon. At a=2a=2 and ε=0.1\varepsilon=0.1, it is approximately 0.02484570.0248457. Setting δ=δ∗\delta=\delta_* is valid: the boundary supremum equals ε\varepsilon, but boundary inputs are excluded and every interior error is strictly smaller. This is a useful example of why strict and nonstrict inequalities cannot be interchanged casually.

How to show a candidate fails

Negating the definition gives: there exists some ε0>0\varepsilon_0>0 such that for every δ>0\delta>0, there is a permitted xx with 0<∣x−a∣<δ0<|x-a|<\delta and ∣f(x)−L∣≥ε0|f(x)-L|\geq\varepsilon_0.

For f(x)=∣x∣/xf(x)=|x|/x, the candidate L=1L=1 fails at zero. Choose ε0=1\varepsilon_0=1. Whatever δ>0\delta>0 is offered, choose x=−δ/2x=-\delta/2. Then f(x)=−1f(x)=-1 and ∣f(x)−1∣=2≥1|f(x)-1|=2\geq1. To rule out every candidate, the one-sided disagreement established previously is more economical: any two-sided limit would equal both −1-1 and 11.

A finite limit is unique

Suppose two distinct finite limits LL and MM existed. Let ε=∣L−M∣/3\varepsilon=|L-M|/3. There would be neighborhoods in which outputs are within ε\varepsilon of each limit. Choose a noncentral input inside both neighborhoods. The triangle inequality would give

∣L−M∣≤∣L−f(x)∣+∣f(x)−M∣<2∣L−M∣/3,|L-M|\leq|L-f(x)|+|f(x)-M|<2|L-M|/3,

a contradiction. The existence of nearby domain points matters; our punctured-neighborhood assumption supplies them. This proof establishes uniqueness without relying on a picture.

Exercises

  1. Prove lim⁡x→−1(5x+2)=−3\lim_{x\to-1}(5x+2)=-3 with an explicit δ\delta.
  2. Prove lim⁡x→0x2=0\lim_{x\to0}x^2=0 using the largest symmetric tolerance.
  3. Prove lim⁡x→3x2=9\lim_{x\to3}x^2=9 using a local bound of radius one.
  4. Explain why choosing δ=ε/∣x+2∣\delta=\varepsilon/|x+2| is not an acceptable construction in the quadratic proof.
  5. For a=1a=1, ε=0.25\varepsilon=0.25, does δ=0.2\delta=0.2 guarantee the square-function tolerance? Find the supremum error.
  6. Show that the candidate limit L=0L=0 fails for the constant function f(x)=2f(x)=2 at zero.
Worked solutions
  1. Choose δ=ε/5\delta=\varepsilon/5. The output error is 5∣x+1∣<ε5|x+1|<\varepsilon.
  2. Choose δ=ε\delta=\sqrt\varepsilon. Then 0<∣x∣<δ0<|x|<\delta implies x2<εx^2<\varepsilon.
  3. If ∣x−3∣<1|x-3|<1, then ∣x+3∣<7|x+3|<7. Choose min⁡{1,ε/7}\min\{1,\varepsilon/7\} and factor ∣x2−9∣|x^2-9|.
  4. This choice depends on the later input xx. We must provide a single radius that works for every input inside it; bounding ∣x+2∣|x+2| by a fixed constant accomplishes that.
  5. The supremum is 2(1)(0.2)+0.22=0.442(1)(0.2)+0.2^2=0.44, greater than 0.250.25, so the guarantee fails. For example, x=1.19x=1.19 is allowed and produces error 0.41610.4161.
  6. Choose ε0=1\varepsilon_0=1. Every nonzero input within every proposed positive radius has error 2≥12\geq1.

Open this chapter with navigation and any associated laboratory →

Chapter 5

Continuity and Existence

Diagnose discontinuities, apply the intermediate and extreme value theorems, and distinguish existence from uniqueness.

Prerequisites: Finite limits, one-sided limits, and the epsilon–delta definition.

  • Check the three conditions for continuity at a point.
  • Apply IVT and EVT with all interval hypotheses stated.
  • Use bisection to turn an existence argument into a controlled approximation.

A thermostat reading can pass from below a target to above it. If temperature varies continuously, it must equal the target somewhere between the readings. This conclusion concerns an unobserved instant. It follows from a mathematical assumption about the entire interval, not from the two readings alone. A sensor that jumps between quantized readings need not display every intermediate number.

Matching nearby behavior to the actual value

A function is continuous at aa when f(a)f(a) exists, the finite limit lim⁡x→af(x)\lim_{x\to a}f(x) exists, and that limit equals f(a)f(a). For an interior point these are two-sided conditions. At an included endpoint of an interval, use the appropriate one-sided limit. Continuity on [a,b][a,b] means continuity at all interior points, right continuity at aa, and left continuity at bb.

For f(x)=(x2−1)/(x−1)f(x)=(x^2-1)/(x-1), the limit at 11 is 22 but the function is undefined there. Defining f(1)=2f(1)=2 creates a continuous extension. Defining f(1)=7f(1)=7 leaves a removable discontinuity: changing one value can repair it. A jump with different left and right limits cannot be repaired that way. Neither can the unbounded behavior of 1/x1/x at zero.

Polynomials are continuous everywhere on the real line. Rational functions are continuous wherever their denominators are nonzero. Compositions of continuous functions are continuous where the composition is defined. These are consequences of limit laws. They let us verify continuity without repeating a tolerance proof at every point.

Joining two formulas

Let f(x)=kx+1f(x)=kx+1 for x<2x<2 and f(x)=x2−1f(x)=x^2-1 for x≥2x\geq2. Each piece is continuous on its own interval. At the join, the left limit is 2k+12k+1, the right limit is 33, and f(2)=3f(2)=3. Continuity requires 2k+1=32k+1=3, hence k=1k=1. Checking only that both formulas are polynomials misses the join, which is the only point that needs extra work.

The intermediate value theorem

If ff is continuous on [a,b][a,b] and NN lies between f(a)f(a) and f(b)f(b), then some c∈[a,b]c\in[a,b] satisfies f(c)=Nf(c)=N. If NN lies strictly between the endpoint values, such a cc lies in (a,b)(a,b). The theorem guarantees existence. It does not guarantee one solution, a formula for it, or that the function is monotone.

For p(x)=x3+x−1p(x)=x^3+x-1, continuity holds on [0,1][0,1], and p(0)=−1p(0)=-1, p(1)=1p(1)=1. Therefore a zero lies in (0,1)(0,1). It is unique because for v>uv>u,

p(v)−p(u)=(v−u)(v2+uv+u2+1)>0.p(v)-p(u)=(v-u)(v^2+uv+u^2+1)>0.

The second factor is positive since v2+uv+u2=(v+u/2)2+3u2/4≥0v^2+uv+u^2=(v+u/2)^2+3u^2/4\geq0. We supplied a separate monotonicity argument; IVT alone would not establish uniqueness.

The continuity hypothesis is essential. A function equal to −1-1 for x<0x<0 and 11 for x≥0x\geq0 changes sign across [−1,1][-1,1] without ever equalling zero. A sign change is not a substitute for continuity.

Bisection with an error bound

If a continuous function has opposite signs at the endpoints, evaluate its midpoint. Keep the half whose endpoints still bracket a sign change, unless the midpoint is already a root. Repeating gives nested intervals. After nn halvings, the bracket length is (b−a)/2n(b-a)/2^n, and its midpoint is within (b−a)/2n+1(b-a)/2^{n+1} of at least one root remaining in the bracket.

For pp above, p(0.5)=−0.375p(0.5)=-0.375, while p(0.75)=0.171875p(0.75)=0.171875. A root therefore lies in (0.5,0.75)(0.5,0.75). The midpoint 0.6250.625 has a guaranteed error at most 0.1250.125 at this stage. To guarantee midpoint error no greater than 0.0010.001 from the initial unit interval, nine halvings suffice because 1/210<0.0011/2^{10}<0.001. A small function residual alone is a different notion of error; a very flat function can have a small residual far from a root.

The extreme value theorem

If ff is continuous on a closed bounded interval [a,b][a,b], it attains an absolute maximum and an absolute minimum on that interval. The conclusion is stronger than having a finite upper or lower bound: actual inputs achieve the extreme outputs.

Both interval conditions matter. On (0,1)(0,1), the continuous function f(x)=xf(x)=x attains neither its supremum 11 nor its infimum 00. On R\mathbb R, the continuous function x2x^2 has no absolute maximum. On [0,1][0,1], a discontinuous function can be unbounded: set f(0)=0f(0)=0 and f(x)=1/xf(x)=1/x for x>0x>0.

IVT and EVT are foundational real-number theorems. Their full proofs depend on completeness of the real numbers. Here we use their stated conclusions and identify the hypotheses each application needs. Later, EVT will ensure optimization problems have candidates worth seeking.

Exercises

  1. Find kk making f(x)=kx+2f(x)=kx+2 for x<1x<1, f(x)=x2+4f(x)=x^2+4 for x≥1x\geq1 continuous.
  2. Prove that x3−2x^3-2 has a zero in (1,2)(1,2). Is uniqueness supplied by IVT?
  3. How many bisection halvings reduce an initial bracket of length 44 to midpoint error at most 0.010.01?
  4. Give a continuous function on an open bounded interval with no maximum.
  5. Does EVT imply that a maximum occurs at an interior point?
  6. For f(x)=1/xf(x)=1/x on [−1,1]∖{0}[-1,1]\setminus\{0\}, explain why opposite endpoint signs do not establish a zero.
Worked solutions
  1. The right value and limit at 11 are 55; k+2=5k+2=5 gives k=3k=3.
  2. The polynomial is continuous and has endpoint values −1-1 and 66. IVT establishes a root; uniqueness follows separately from strict increase of the cube function.
  3. Require 4/2n+1≤0.014/2^{n+1}\leq0.01. Eight halvings give 4/512=0.00781254/512=0.0078125; seven give 0.0156250.015625, so eight suffice.
  4. f(x)=xf(x)=x on (0,1)(0,1) has supremum 11 that is never attained.
  5. No. For f(x)=xf(x)=x on [0,1][0,1], the maximum is at 11.
  6. The function is not continuous on the full closed interval; it is undefined at zero. IVT does not apply.

Open this chapter with navigation and any associated laboratory →

Chapter 6

The Derivative at a Point

Derive instantaneous rates and tangent lines, and prove why a corner can be continuous without being differentiable.

Prerequisites: Difference quotients, two-sided limits, and continuity.

  • Calculate a derivative from its defining limit.
  • Construct a tangent line and interpret derivative units.
  • Prove differentiability implies continuity and explain why the converse fails.

For the cart with position s(t)=t2s(t)=t^2, the average velocity from aa to a+ha+h was 2a+h2a+h. Both positive and negative increments approach the same value 2a2a. We can now define the instantaneous velocity without assigning a meaning to division by zero: take the limit of the nonzero-interval average velocities.

For a function defined on an open interval containing aa, the derivative at aa is the finite limit

f′(a)=lim⁡h→0f(a+h)−f(a)h,f'(a)=\lim_{h\to0}\frac{f(a+h)-f(a)}h,

if it exists. An equivalent form uses x=a+hx=a+h and lets x→ax\to a. Both formulations compare the output at the target with outputs at nearby inputs. We use a two-sided limit for interior differentiability. Endpoint derivatives, when useful, will be named one-sided derivatives.

A quadratic from the definition

Let f(x)=2x2−3x+1f(x)=2x^2-3x+1. Expanding gives

f(a+h)−f(a)h=4ah+2h2−3hh=4a+2h−3,h≠0.\frac{f(a+h)-f(a)}h =\frac{4ah+2h^2-3h}h=4a+2h-3,\qquad h\ne0.

Therefore f′(a)=4a−3f'(a)=4a-3. At a=2a=2, the function value is 33 and derivative is 55. The tangent line has equation

y=f(2)+f′(2)(x−2)=3+5(x−2).y=f(2)+f'(2)(x-2)=3+5(x-2).

A tangent is specified by its local limiting slope. It need not stay on one side of the curve or touch it at only one point. The tangent y=0y=0 to x3x^3 at zero crosses the curve there.

Derivatives carry units

If ff measures temperature in degrees Celsius and xx measures hours, f′(a)f'(a) has units degrees Celsius per hour. A derivative of −2-2 means temperature is locally decreasing at that rate; it does not mean the temperature is negative. If position is measured in metres and time in seconds, velocity has units metres per second and acceleration, the derivative of velocity, has units metres per second squared.

For s(t)=t3−3ts(t)=t^3-3t on a specified time interval, v(t)=3t2−3v(t)=3t^2-3. Positive velocity indicates motion in the positive coordinate direction. Speed is ∣v(t)∣|v(t)|, so negative velocity is not negative speed. Later we will integrate speed to compute total distance.

Continuous does not mean differentiable

For f(x)=∣x∣f(x)=|x| at zero, the difference quotient is ∣h∣/h|h|/h. It equals −1-1 for negative hh and 11 for positive hh. Thus the two-sided derivative does not exist even though ∣x∣→0=f(0)|x|\to0=f(0). The graph has a corner.

The function x1/3x^{1/3} at zero has quotient h−2/3h^{-2/3}, which grows without bound. This is often described geometrically as a vertical tangent. It is not a finite derivative under our definition. A discontinuity also prevents differentiability, but a smooth-looking numerical plot is not proof that discontinuities or corners are absent.

Why differentiability forces continuity

Suppose f′(a)f'(a) exists and is finite. For h≠0h\ne0,

f(a+h)−f(a)=hf(a+h)−f(a)h.f(a+h)-f(a)=h\frac{f(a+h)-f(a)}h.

As h→0h\to0, the first factor tends to zero and the second tends to f′(a)f'(a). The product tends to zero. Hence f(a+h)→f(a)f(a+h)\to f(a), establishing continuity at aa. Finiteness of the derivative matters because the ordinary product limit law needs finite limits.

This proof gives a useful diagnostic: if continuity fails, a finite derivative cannot exist. Its converse remains false, as the absolute-value example shows. A theorem and its converse are separate claims.

The local linear model

Rearrange the defining quotient into

f(a+h)=f(a)+f′(a)h+r(h),r(h)h⟶0.f(a+h)=f(a)+f'(a)h+r(h),\qquad \frac{r(h)}h\longrightarrow0.

The derivative supplies a linear change, and the residual is small relative to the increment, not merely small in absolute size. For f(x)=x2f(x)=x^2, r(h)=h2r(h)=h^2, so r(h)/h=h→0r(h)/h=h\to0. This relative-error statement distinguishes a tangent approximation from any arbitrary line passing through the point.

The secant laboratory displays exact formulas for the square, cube, and absolute-value functions. A finite list of small increments cannot prove a derivative. The associated algebra explains the limit; the table and graph help us inspect its consequences. For polynomials the program evaluates algebraically simplified quotients to avoid catastrophic subtraction at small increments.

Exercises

  1. Use the definition to differentiate f(x)=3x2+2xf(x)=3x^2+2x at a general aa.
  2. Find the tangent to x2x^2 at x=−1x=-1.
  3. Use first principles to find the derivative of 1/x1/x at a≠0a\ne0.
  4. Show that the tangent to x3x^3 at zero crosses its graph.
  5. A temperature derivative is zero at noon. Must the temperature be zero? Must it remain constant afterward?
  6. Explain the different reasons that ∣x∣|x| and x1/3x^{1/3} lack a finite derivative at zero.
Worked solutions
  1. Expansion gives [f(a+h)−f(a)]/h=6a+3h+2[f(a+h)-f(a)]/h=6a+3h+2 for nonzero hh, hence f′(a)=6a+2f'(a)=6a+2.
  2. The point is (−1,1)(-1,1) and slope −2-2, so y=1−2(x+1)=−2x−1y=1-2(x+1)=-2x-1.
  3. The quotient simplifies to −1/[a(a+h)]-1/[a(a+h)], provided h≠0h\ne0 and a+h≠0a+h\ne0. Its limit is −1/a2-1/a^2.
  4. The derivative at zero is lim⁡h2=0\lim h^2=0, so the tangent is the horizontal axis. Cubes have opposite signs on opposite sides of zero.
  5. Neither. The derivative describes an instantaneous rate, not the level or a future interval. For example, T(t)=20+t2T(t)=20+t^2 has derivative zero at t=0t=0 without being constant.
  6. Absolute value has two distinct finite one-sided slopes; the cube root has a quotient unbounded in the positive direction. Neither gives the required finite two-sided limit.

Open this chapter with navigation and any associated laboratory →

Chapter 7

The Derivative as a Function

Derive power and reciprocal examples, track derivative domains, and distinguish a function from its rate function.

Prerequisites: Derivatives at a point; binomial expansion and conjugates.

  • Construct a derivative function from pointwise limits.
  • Derive the positive-integer power rule.
  • Explain how the derivative's domain can be smaller than the original domain.

A derivative at one point is a number. A derivative function assigns such numbers to all inputs where the required finite limit exists. For f(x)=x2f(x)=x^2, the calculation f′(a)=2af'(a)=2a holds for every real aa, so the derivative function is f′(x)=2xf'(x)=2x. The variable has changed its role from a fixed base point in one proof to the argument of a new function.

The notation dy/dxdy/dx also denotes the derivative when y=f(x)y=f(x). At this stage it is notation for a limit-defined rate, not permission to cancel arbitrary symbols. Later, the chain rule and substitution will explain why certain manipulations that resemble fraction arithmetic work.

Cubes and integer powers

For f(x)=x3f(x)=x^3, the difference quotient is

(x+h)3−x3h=3x2+3xh+h2,h≠0.\frac{(x+h)^3-x^3}h=3x^2+3xh+h^2,\qquad h\ne0.

Letting h→0h\to0 gives f′(x)=3x2f'(x)=3x^2. Both positive and negative nonzero increments yield the same expression’s limit.

For a positive integer nn, the binomial theorem gives

(x+h)n=xn+nxn−1h+∑k=2n(nk)xn−khk.(x+h)^n=x^n+nx^{n-1}h+\sum_{k=2}^n {n\choose k}x^{n-k}h^k.

Subtract xnx^n and divide by nonzero hh. Every term of the sum contains a positive power of hh and tends to zero. The result is the power rule d(xn)/dx=nxn−1d(x^n)/dx=nx^{n-1} for positive integers. The case n=1n=1 is simply a constant quotient 11. A constant function has derivative zero because its output difference is always zero. Extending the power rule to all real exponents needs additional results; the binomial argument alone does not supply that extension.

Reciprocal and square root

The reciprocal calculation gives

(1x)′=−1x2,x≠0.\left(\frac1x\right)'=-\frac1{x^2},\qquad x\ne0.

The derivative cannot include zero because the original function has no value there. For f(x)=xf(x)=\sqrt{x} at x>0x>0, conjugates give

x+h−xh=1x+h+x⟶12x.\frac{\sqrt{x+h}-\sqrt{x}}h =\frac1{\sqrt{x+h}+\sqrt{x}} \longrightarrow\frac1{2\sqrt{x}}.

The function’s domain is [0,∞)[0,\infty), while its ordinary two-sided derivative has domain (0,∞)(0,\infty). Even the right-hand difference quotient at zero grows without bound, so it does not supply a finite endpoint derivative. A formula and its domain must travel together.

For ∣x∣|x|, the derivative is −1-1 on (−∞,0)(-\infty,0) and 11 on (0,∞)(0,\infty), and is undefined at zero. The derivative function records a change of slope and the missing derivative at the corner. Drawing a line between the two branches would add information the function does not contain.

Higher derivatives and interpretation

If f′f' is differentiable, its derivative is f′′f'', the second derivative. For f(x)=x3f(x)=x^3, f′=3x2f'=3x^2, f′′=6xf''=6x, f′′′=6f'''=6, and every derivative of order four or greater is zero. The nnth derivative is written f(n)f^{(n)}, while fnf^n normally means a power; these symbols are not interchangeable.

For position s(t)s(t), the first derivative is velocity and the second acceleration. Positive acceleration means velocity is increasing, but speed need not increase. If a vehicle is moving in the negative direction with v=−10v=-10 and a=2a=2, its velocity moves toward zero and its speed initially decreases. The signs of velocity and acceleration together determine this behavior.

Reading a derivative graph

Where f′>0f'>0, tangent slopes are positive; where f′<0f'<0, they are negative. A theorem in a later chapter will justify conclusions about increase on an interval. A derivative equal to zero at one point does not itself identify a maximum or minimum: x3x^3 has derivative zero at zero and continues increasing through it.

Two different functions can have the same derivative. For instance x2x^2 and x2+7x^2+7 both have derivative 2x2x. Differentiation discards a vertical offset. The mean value theorem will show that, on an interval, functions with the same derivative can differ only by a constant. This observation will become the foundation of antiderivatives.

Exercises

  1. Differentiate x5x^5 directly by expansion, identifying the term that survives the limit.
  2. Give the domains of x−2\sqrt{x-2} and its derivative.
  3. Find the first three derivatives of x4x^4.
  4. Describe the derivative of ∣x−3∣|x-3| at every real input where it exists.
  5. Give two functions with derivative 4x34x^3.
  6. Explain why negative velocity with negative acceleration produces increasing speed locally when both signs persist.
Worked solutions
  1. The quotient is 5x4+10x3h+10x2h2+5xh3+h45x^4+10x^3h+10x^2h^2+5xh^3+h^4. Only 5x45x^4 survives.
  2. The function is defined on [2,∞)[2,\infty); the finite ordinary derivative 1/(2x−2)1/(2\sqrt{x-2}) is defined on (2,∞)(2,\infty).
  3. 4x34x^3, 12x212x^2, and 24x24x.
  4. It equals −1-1 for x<3x<3, 11 for x>3x>3, and is undefined at 33 because the one-sided slopes disagree.
  5. x4x^4 and x4+10x^4+10 are examples.
  6. Velocity becomes more negative; its magnitude therefore increases. Speed is −v-v while v<0v<0, so its derivative is −a>0-a>0.

Open this chapter with navigation and any associated laboratory →

Chapter 8

Sum, Product, and Quotient Rules

Derive the rules that make differentiation practical and learn to recognize the operation at the outside of an expression.

Prerequisites: Derivative functions and finite limit laws.

  • Derive the product rule by adding and subtracting an intermediate term.
  • Apply sum, product, and quotient rules with their domain restrictions.
  • Verify a result by an independent algebraic representation.

First principles establish what a derivative means. They are less convenient for differentiating a complicated expression repeatedly. We can prove rules once, then use them whenever their hypotheses hold. Every rule in this chapter assumes the component functions are differentiable at the point in question. The quotient rule also requires a nonzero denominator.

Linearity comes from the limit laws

For constants c,dc,d, expand the difference quotient of cf+dgcf+dg and split it into two quotients. Taking limits gives

(cf+dg)′=cf′+dg′.(cf+dg)'=cf'+dg'.

Thus (4x3−2x2+7)′=12x2−4x(4x^3-2x^2+7)'=12x^2-4x. The constant contributes zero. Differentiating the square of a function is not covered by linearity; squaring is multiplication, not addition.

Why both factors contribute

For F=fgF=fg, write the numerator as

f(x+h)g(x+h)−f(x)g(x)=f(x+h)[g(x+h)−g(x)]+g(x)[f(x+h)−f(x)].f(x+h)g(x+h)-f(x)g(x) =f(x+h)[g(x+h)-g(x)]+g(x)[f(x+h)-f(x)].

After division by hh, take limits. Differentiability implies f(x+h)→f(x)f(x+h)\to f(x), so

(fg)′=fg′+gf′.(fg)'=fg'+gf'.

The extra term that a naive “differentiate both factors and multiply” approach loses records the change in one factor while the other remains present. For x⋅xx\cdot x, the rule gives x+x=2xx+x=2x, whereas multiplying the two derivatives gives 11, which fails even this basic check.

For F(x)=(x2+1)(x3−2)F(x)=(x^2+1)(x^3-2),

F′=2x(x3−2)+(x2+1)3x2=5x4+3x2−4x.F'=2x(x^3-2)+(x^2+1)3x^2=5x^4+3x^2-4x.

Expanding first gives F=x5+x3−2x2−2F=x^5+x^3-2x^2-2 and the same derivative. Independent representations are useful error checks. They are often more informative than repeating the same calculation with slightly different notation.

The quotient rule

For a differentiable nonzero gg, the difference quotient for 1/g1/g simplifies to

−g(x+h)−g(x)h g(x+h)g(x).-\frac{g(x+h)-g(x)}{h\,g(x+h)g(x)}.

Its limit is −g′/g2-g'/g^2. Continuity and g(x)≠0g(x)\ne0 ensure g(x+h)g(x+h) stays nonzero for sufficiently small increments. Applying the product rule to f(1/g)f(1/g) gives

(fg)′=f′g−fg′g2.\left(\frac fg\right)'=\frac{f'g-fg'}{g^2}.

The order in the numerator matters. For Q(x)=(x2+1)/(x−1)Q(x)=(x^2+1)/(x-1),

Q′(x)=2x(x−1)−(x2+1)(x−1)2=x2−2x−1(x−1)2,x≠1.Q'(x)=\frac{2x(x-1)-(x^2+1)}{(x-1)^2} =\frac{x^2-2x-1}{(x-1)^2},\quad x\ne1.

Polynomial division gives Q=x+1+2/(x−1)Q=x+1+2/(x-1), whose derivative is 1−2/(x−1)21-2/(x-1)^2, algebraically identical to the result. The forbidden input 11 remains forbidden in both versions.

Negative integer powers

For a positive integer nn, apply the reciprocal rule to xnx^n:

ddxx−n=−nxn−1x2n=−nx−n−1,x≠0.\frac d{dx}x^{-n}=-\frac{nx^{n-1}}{x^{2n}}=-nx^{-n-1},\qquad x\ne0.

The power rule therefore holds for all integer exponents on their relevant domains. At exponent zero the function is constant wherever it is defined. This extension came from a new argument, not from silently expanding the earlier binomial proof beyond its assumptions.

Choosing a method

First identify the outermost operation. An expression can contain sums inside a product inside a quotient. Start with the outside, retain its component functions intact, then differentiate those components. Simplify when simplification reduces error; leaving a correctly factored derivative is often preferable to an unnecessarily long polynomial.

Numerical difference quotients at ordinary inputs can detect a sign error. Very small increments may instead amplify floating-point cancellation. A numerical agreement supports a check at that input; the proved rules establish the symbolic result across the domain.

Exercises

  1. Differentiate (2x−1)(x2+3)(2x-1)(x^2+3) by the product rule and by expansion.
  2. Differentiate (x+1)/(x+2)(x+1)/(x+2) and state its domain.
  3. Differentiate x−3+5x2−8x^{-3}+5x^2-8.
  4. Explain why (f/g)′=f′/g′(f/g)'=f'/g' cannot be a general rule.
  5. Derive the rule for a product of three differentiable functions using the two-factor rule twice.
  6. Find the tangent slope of (x2+1)/(x−1)(x^2+1)/(x-1) at x=0x=0.
Worked solutions
  1. 2(x2+3)+(2x−1)2x=6x2−2x+62(x^2+3)+(2x-1)2x=6x^2-2x+6. Expanding first gives 2x3−x2+6x−32x^3-x^2+6x-3 and the same derivative.
  2. [(x+2)−(x+1)]/(x+2)2=1/(x+2)2[(x+2)-(x+1)]/(x+2)^2=1/(x+2)^2, for x≠−2x\ne-2.
  3. −3x−4+10x-3x^{-4}+10x, for x≠0x\ne0.
  4. Take f=g=xf=g=x away from zero. The quotient is constant with derivative zero, while f′/g′=1f'/g'=1.
  5. (fgh)′=f′gh+fg′h+fgh′(fgh)'=f'gh+fg'h+fgh'.
  6. Substitute into the verified derivative: (−1)/1=−1(-1)/1=-1.

Open this chapter with navigation and any associated laboratory →

Chapter 9

Composition and Trigonometric Derivatives

Justify the chain rule without dividing by a possibly zero increment, then derive sine and cosine derivatives.

Prerequisites: Product and quotient rules; the fundamental trigonometric limits in radians.

  • Differentiate nested functions using the chain rule.
  • Derive sine and cosine derivatives from angle-addition formulas.
  • Track the units and radian convention in oscillatory models.

If a quantity depends on temperature and temperature depends on time, its time rate should reflect both changes. The chain rule makes that statement exact. If gg is differentiable at xx and ff at g(x)g(x), then

ddxf(g(x))=f′(g(x))g′(x).\frac d{dx}f(g(x))=f'(g(x))g'(x).

The outer derivative is evaluated at the actual inner value. Writing f′(x)g′(x)f'(x)g'(x) usually differentiates a different expression.

A justification that handles zero inner changes

Differentiability gives f(u+k)−f(u)=f′(u)k+kη(k)f(u+k)-f(u)=f'(u)k+k\eta(k), where η(k)→0\eta(k)\to0 as k→0k\to0 and we may set η(0)=0\eta(0)=0. Put u=g(x)u=g(x) and k=g(x+h)−g(x)k=g(x+h)-g(x). Divide by hh, not by kk:

f(g(x+h))−f(g(x))h=[f′(g(x))+η(k)]g(x+h)−g(x)h.\frac{f(g(x+h))-f(g(x))}h =[f'(g(x))+\eta(k)]\frac{g(x+h)-g(x)}h.

Continuity of gg gives k→0k\to0, and the second quotient tends to g′(x)g'(x). The product has the claimed limit. This avoids the common gap in a proof that divides by g(x+h)−g(x)g(x+h)-g(x) even though that difference can be zero for arbitrarily small nonzero hh.

For F(x)=(3x2+1)4F(x)=(3x^2+1)^4, the outer function is u4u^4 and the inner is 3x2+13x^2+1. Hence F′=4(3x2+1)3(6x)F'=4(3x^2+1)^3(6x). For 1+x2\sqrt{1+x^2}, the result is x/1+x2x/\sqrt{1+x^2}. Its denominator is always positive, so the derivative exists for all real xx.

Sine from first principles

The angle-addition identity gives

sin⁡(x+h)−sin⁡xh=sin⁡xcos⁡h−1h+cos⁡xsin⁡hh.\frac{\sin(x+h)-\sin x}h =\sin x\frac{\cos h-1}h+\cos x\frac{\sin h}h.

The two trigonometric limits previously established are 00 and 11. Therefore (sin⁡x)′=cos⁡x(\sin x)'=\cos x. Applying the cosine angle-addition identity similarly yields (cos⁡x)′=−sin⁡x(\cos x)'=-\sin x. These formulas assume radians. They were not used to prove the limits on which their derivation relies.

The quotient rule now gives

(tan⁡x)′=cos⁡2x+sin⁡2xcos⁡2x=sec⁡2x,(\tan x)'=\frac{\cos^2x+\sin^2x}{\cos^2x}=\sec^2x,

where cos⁡x≠0\cos x\ne0. Reciprocals and quotients also give (cot⁡x)′=−csc⁡2x(\cot x)'=-\csc^2x, (sec⁡x)′=sec⁡xtan⁡x(\sec x)'=\sec x\tan x, and (csc⁡x)′=−csc⁡xcot⁡x(\csc x)'=-\csc x\cot x on their domains. Stating a derivative formula does not remove the poles of the original function.

A time-dependent oscillation

For a modeled displacement s(t)=2sin⁡(3t+π/4)s(t)=2\sin(3t+\pi/4) metres, where tt is in seconds and the phase is in radians, the velocity is v(t)=6cos⁡(3t+π/4)v(t)=6\cos(3t+\pi/4) metres per second. Acceleration is a(t)=−18sin⁡(3t+π/4)a(t)=-18\sin(3t+\pi/4) metres per second squared. The factor 33 has units radians per second, with radians dimensionless in the derivative convention. Forgetting the inner derivative loses the frequency factor in every rate.

If an input θ\theta is explicitly measured numerically in degrees, the function is sin⁡(πθ/180)\sin(\pi\theta/180) in radian notation. Its derivative with respect to the degree-valued input is (π/180)cos⁡(πθ/180)(\pi/180)\cos(\pi\theta/180). The familiar cosine-only rule is not unit-independent.

Nested operations

For H(x)=sin⁡((x2+1)3)H(x)=\sin((x^2+1)^3), apply the chain rule at each level:

H′(x)=cos⁡((x2+1)3) 3(x2+1)2 2x.H'(x)=\cos((x^2+1)^3)\,3(x^2+1)^2\,2x.

For J(x)=x2sin⁡(3x)J(x)=x^2\sin(3x), multiplication is outside the composition, so first use the product rule: J′=2xsin⁡(3x)+3x2cos⁡(3x)J'=2x\sin(3x)+3x^2\cos(3x). A derivative can require several rules, each corresponding to a different layer of the expression.

Exercises

  1. Differentiate (1−2x)7(1-2x)^7.
  2. Differentiate cos⁡(x2)\cos(x^2) and sin⁡2x\sin^2x. Explain the different compositions.
  3. Find the derivative of tan⁡(4x)\tan(4x) and state where it exists.
  4. Differentiate xcos⁡xx\cos x twice.
  5. A displacement is 5cos⁡(2t)5\cos(2t) metres. Find velocity and acceleration at t=π/4t=\pi/4 seconds.
  6. Explain why a chain-rule proof must not assume every inner increment is nonzero.
Worked solutions
  1. −14(1−2x)6-14(1-2x)^6.
  2. −2xsin⁡(x2)-2x\sin(x^2) and 2sin⁡xcos⁡x2\sin x\cos x. The first applies cosine to a square; the second squares a sine.
  3. 4sec⁡2(4x)4\sec^2(4x) wherever cos⁡(4x)≠0\cos(4x)\ne0.
  4. First derivative cos⁡x−xsin⁡x\cos x-x\sin x; second −2sin⁡x−xcos⁡x-2\sin x-x\cos x.
  5. v=−10sin⁡(2t)v=-10\sin(2t) and a=−20cos⁡(2t)a=-20\cos(2t), giving −10-10 metres per second and zero acceleration at the stated time.
  6. An inner function may be constant locally or return to its base value at many nearby inputs. The remainder formulation does not divide by that increment.

Open this chapter with navigation and any associated laboratory →

Chapter 10

Exponentials, Logarithms, and Inverses

Differentiate inverse functions and derive logarithmic, exponential, inverse-trigonometric, and real-power rules on their domains.

Prerequisites: Chain rule, inverse functions, logarithm laws, and trigonometric derivatives.

  • Apply the inverse derivative theorem with a nonzero derivative.
  • Differentiate exponential, logarithmic, and inverse-trigonometric functions.
  • Use logarithmic differentiation without losing sign or domain restrictions.

Exponential growth describes a quantity whose instantaneous rate is proportional to its current size. Before using that model, we need the derivative of an exponential and a precise account of the inverse relationship between exponentials and logarithms.

We take the natural exponential and logarithm as elementary functions with their usual continuous inverse relationship and the normalization lim⁡h→0(eh−1)/h=1\lim_{h\to0}(e^h-1)/h=1. This characterization of the base ee is an input here, not a consequence of the differentiation rules already proved. One can instead construct log⁡x\log x by an integral and derive exponentials after the fundamental theorem; that alternative ordering leads to the same functions.

Exponentials and inverse rates

Using ex+h=exehe^{x+h}=e^xe^h in the difference quotient gives (ex)′=ex(e^x)'=e^x. The chain rule gives (eu(x))′=eu(x)u′(x)(e^{u(x)})'=e^{u(x)}u'(x). For a fixed base b>0b>0, write bx=exlog⁡bb^x=e^{x\log b}, obtaining (bx)′=bxlog⁡b(b^x)'=b^x\log b. If b=1b=1, the formula correctly returns zero.

If y=f(x)y=f(x) has a differentiable inverse locally, the identity f(f−1(y))=yf(f^{-1}(y))=y and the chain rule imply

(f−1)′(y)=1f′(f−1(y)),(f^{-1})'(y)=\frac1{f'(f^{-1}(y))},

provided the denominator is nonzero. A standard inverse derivative theorem ensures differentiability of the inverse when ff is continuously differentiable near the point with nonzero derivative. The requirement is substantive: x3x^3 has an inverse, but at zero its derivative is zero and the inverse cube root lacks a finite derivative there.

Applying the theorem to exe^x gives (log⁡x)′=1/x(\log x)'=1/x for x>0x>0. For x≠0x\ne0, (log⁡∣x∣)′=1/x(\log|x|)'=1/x on each half-line. A logarithm of a negative argument remains undefined over the reals; the absolute value supplies a different, well-defined function.

Variable arguments and real powers

For a differentiable nonzero uu, differentiation of log⁡∣u∣\log|u| gives u′/uu'/u wherever uu does not vanish. Thus (log⁡(x2+1))′=2x/(x2+1)(\log(x^2+1))'=2x/(x^2+1) on all real inputs. For log⁡(3x−2)\log(3x-2) we additionally require x>2/3x>2/3.

For a real constant rr and x>0x>0, define xr=erlog⁡xx^r=e^{r\log x}. Differentiation yields rxr−1rx^{r-1}. Some rational powers extend to negative arguments and require separate domain care; this proof covers the positive domain uniformly.

When both base and exponent vary, an ordinary power rule is insufficient. For y=xxy=x^x on x>0x>0, take logarithms:

log⁡y=xlog⁡x,y′y=log⁡x+1,y′=xx(log⁡x+1).\log y=x\log x,\qquad \frac{y'}y=\log x+1, \qquad y'=x^x(\log x+1).

The function is positive on the stated domain, so taking the logarithm is legitimate. For y=(x2+1)sin⁡xy=(x^2+1)^{\sin x}, the same method gives

y′=(x2+1)sin⁡x[cos⁡xlog⁡(x2+1)+2xsin⁡xx2+1].y'=(x^2+1)^{\sin x}\left[\cos x\log(x^2+1)+\frac{2x\sin x}{x^2+1}\right].

Inverse trigonometric functions

Restrict sine to [−π/2,π/2][-\pi/2,\pi/2] to define arcsin⁡\arcsin. At interior inputs ∣x∣<1|x|<1, set y=arcsin⁡xy=\arcsin x. Then sin⁡y=x\sin y=x, so cos⁡y y′=1\cos y\,y'=1. On that branch cosine is positive, giving

(arcsin⁡x)′=11−x2.(\arcsin x)'=\frac1{\sqrt{1-x^2}}.

Likewise (arccos⁡x)′=−1/1−x2(\arccos x)'=-1/\sqrt{1-x^2} for ∣x∣<1|x|<1. For y=arctan⁡xy=\arctan x in (−π/2,π/2)(-\pi/2,\pi/2), sec⁡2y y′=1\sec^2y\,y'=1, hence (arctan⁡x)′=1/(1+x2)(\arctan x)'=1/(1+x^2) for every real xx. The inverse branch determines the sign of a square root; it cannot be ignored.

A growth model

For P(t)=P0ektP(t)=P_0e^{kt} with P0>0P_0>0, P′=kPP'=kP. The parameter kk has units inverse time. Doubling time for k>0k>0 solves ekt=2e^{kt}=2, giving log⁡2/k\log2/k. A proportional growth model assumes a constant rate parameter and no limiting resources or changing conditions. Differentiating the formula does not establish those modeling assumptions.

Exercises

  1. Differentiate e−3x2e^{-3x^2} and 25x2^{5x}.
  2. Differentiate log⁡∣x2−4∣\log|x^2-4| and state the domain.
  3. Differentiate arctan⁡(2x)\arctan(2x) and arcsin⁡(x/2)\arcsin(x/2) on their derivative domains.
  4. Differentiate xsin⁡xx^{\sin x} on x>0x>0.
  5. A positive quantity grows as P0e0.04tP_0e^{0.04t} with tt in years. Find its doubling time and its instantaneous relative rate.
  6. Why can the inverse derivative formula not give a finite derivative for the cube root at zero?
Worked solutions
  1. −6xe−3x2-6xe^{-3x^2} and 5log⁡2 25x5\log2\,2^{5x}.
  2. 2x/(x2−4)2x/(x^2-4), with x≠±2x\ne\pm2.
  3. 2/(1+4x2)2/(1+4x^2) for all real xx; 1/[21−x2/4]1/[2\sqrt{1-x^2/4}] for ∣x∣<2|x|<2.
  4. xsin⁡x[cos⁡xlog⁡x+(sin⁡x)/x]x^{\sin x}[\cos x\log x+(\sin x)/x].
  5. log⁡2/0.04≈17.329\log2/0.04\approx17.329 years, and P′/P=0.04P'/P=0.04 per year.
  6. The original cube function has zero derivative at zero, so the theorem’s nonzero-denominator hypothesis fails. Direct quotients confirm the inverse has unbounded slope.

Open this chapter with navigation and any associated laboratory →

Chapter 11

Implicit Differentiation and Related Rates

Differentiate constraints, identify local branches, and build rate equations before substituting a particular instant.

Prerequisites: Chain rule and elementary plane and solid geometry.

  • Differentiate an equation that defines a local function implicitly.
  • Distinguish derivative failure from failure of a chosen coordinate representation.
  • Solve related-rate models with correct signs and units.

A circle cannot be represented by one real function y=f(x)y=f(x) over its entire horizontal span: most horizontal positions have two vertical coordinates. It can nevertheless have well-defined tangent lines on its local branches. An equation relating the variables can retain information that solving explicitly would split into cases.

For x2+y2=25x^2+y^2=25, treat yy as a differentiable function of xx on a branch. Differentiation gives 2x+2yy′=02x+2yy'=0, hence

y′=−x/y(y≠0).y'=-x/y\qquad(y\ne0).

At (3,4)(3,4), the slope is −3/4-3/4, so the tangent is y−4=−(3/4)(x−3)y-4=-(3/4)(x-3). At (5,0)(5,0), the formula cannot produce a finite slope for yy as a function of xx. The circle still has a tangent, namely the vertical line x=5x=5. A failure of this coordinate representation is not a disappearance of the geometric curve.

For a general differentiable equation F(x,y)=0F(x,y)=0, the chain rule gives Fx+Fyy′=0F_x+F_yy'=0 when the partial derivatives are defined, leading formally to y′=−Fx/Fyy'=-F_x/F_y where Fy≠0F_y\ne0. A local implicit function theorem supplies sufficient conditions for such a branch to exist. In this course, polynomial and elementary examples will also be checked directly when possible; the displayed ratio should not be treated as an unconditional existence theorem.

From the circle’s y′=−x/yy'=-x/y, the quotient rule gives

y′′=−y+xy′y2=−x2+y2y3=−25y3.y''=\frac{-y+xy'}{y^2} =-\frac{x^2+y^2}{y^3}=-\frac{25}{y^3}.

The first equality differentiates yy as a function rather than freezing it as a constant. At (3,4)(3,4), this is −25/64-25/64. The upper branch is concave down there, consistent with its geometry.

In a related-rate problem, several quantities depend on the same time variable. Write the constraint between them, then differentiate with respect to time. For a spherical balloon with radius r(t)r(t) and volume V(t)=4πr(t)3/3V(t)=4\pi r(t)^3/3,

dVdt=4πr2drdt.\frac{dV}{dt}=4\pi r^2\frac{dr}{dt}.

If the radius is 22 centimetres and volume is increasing at 12π12\pi cubic centimetres per second, then r′=12π/(16π)=3/4r'=12\pi/(16\pi)=3/4 centimetres per second. Substituting r=2r=2 before differentiating would turn volume into a constant and destroy the time dependence we need.

A rigid ladder of length 55 metres touches a horizontal floor and vertical wall. Let x(t)x(t) be the foot’s distance from the wall and y(t)>0y(t)>0 the top’s height. While the geometric model applies, x2+y2=25x^2+y^2=25 and

2xx′+2yy′=0,y′=−xyx′.2xx'+2yy'=0,\qquad y'=-\frac{x}{y}x'.

At x=3x=3, y=4y=4, and x′=0.2x'=0.2 metres per second, the top moves at y′=−0.15y'=-0.15 metres per second. The sign describes downward motion. Its downward speed is 0.150.15 metres per second.

As y→0+y\to0^+ with fixed positive x′x', the formula predicts unbounded downward rate. That signals a limitation of extending this idealized contact-and-constant-foot-speed model to the final configuration. Real mechanics, contact forces, and finite energy may invalidate the imposed rate before that limit. Mathematics can expose a model’s consequences without certifying its physical continuation.

Define variables and units; identify the instant; draw or describe the geometry; write the constraint; differentiate while quantities are still functions; substitute the instant’s values; solve for the requested rate; interpret its sign and domain. A clear equation often does more work than a memorized “related rates” formula.

  1. Find the tangent slope to xy=6xy=6 at (2,3)(2,3) by implicit differentiation.
  2. Find y′y' for x2+xy+y2=7x^2+xy+y^2=7, wherever the denominator is nonzero.
  3. A circle’s radius increases at 0.40.4 centimetres per second. Find its area rate when r=5r=5 centimetres.
  4. Repeat the ladder calculation when x=4x=4, x′=0.2x'=0.2.
  5. Explain why a rate measured in cubic metres per second cannot be reported as a radius rate without an additional geometric factor.
  6. For the upper semicircle y=25−x2y=\sqrt{25-x^2}, verify the implicit derivative by explicit differentiation.
Worked solutions
  1. y+xy′=0y+xy'=0 gives y′=−y/x=−3/2y'=-y/x=-3/2.
  2. 2x+y+xy′+2yy′=02x+y+xy'+2yy'=0, so y′=−(2x+y)/(x+2y)y'=-(2x+y)/(x+2y).
  3. A′=2πrr′=4πA'=2\pi rr'=4\pi square centimetres per second.
  4. y=3y=3, hence y′=−(4/3)(0.2)=−4/15y'=-(4/3)(0.2)=-4/15 metres per second.
  5. Volume and length are different quantities. For a sphere the conversion factor is 1/(4πr2)1/(4\pi r^2), with inverse-area units.
  6. The chain rule gives −x/25−x2=−x/y-x/\sqrt{25-x^2}=-x/y for −5<x<5-5<x<5.

Open this chapter with navigation and any associated laboratory →

Chapter 12

The Mean Value Theorem and Local Approximation

Connect interval change to instantaneous rates, prove monotonicity consequences, and bound linearization error.

Prerequisites: Continuity, EVT, and differentiation rules.

  • Apply Rolle's theorem and MVT only on intervals meeting their hypotheses.
  • Deduce monotonicity and the constant-difference theorem.
  • Construct a linearization and bound its remainder using a second derivative.

If a journey covers 100 kilometres in two hours, must the instantaneous speed ever equal 50 kilometres per hour? Under a differentiable position model for one-directional travel, the mean value theorem says yes. Its force comes from linking information about a whole interval to a derivative at some interior point.

An interior extremum has a stationary derivative

Suppose ff has a local maximum at an interior point cc and is differentiable there. For small positive hh, [f(c+h)−f(c)]/h≤0[f(c+h)-f(c)]/h\leq0; for small negative hh, the quotient is nonnegative. If the two-sided limit exists, it must be both nonpositive and nonnegative, hence zero. A local minimum is similar. This is Fermat’s stationary-point theorem. It does not apply to endpoints or points where the derivative fails to exist, and its converse is false: x3x^3 is stationary at zero without an extremum.

Rolle and the mean value theorem

Rolle’s theorem: if ff is continuous on [a,b][a,b], differentiable on (a,b)(a,b), and f(a)=f(b)f(a)=f(b), some c∈(a,b)c\in(a,b) has f′(c)=0f'(c)=0. EVT gives attained extremes. If the function is nonconstant, at least one extreme differs from the equal endpoint value and is attained inside; Fermat’s theorem applies. If it is constant, every interior derivative is zero.

For the mean value theorem, remove the secant line. Define

g(x)=f(x)−f(a)−f(b)−f(a)b−a(x−a).g(x)=f(x)-f(a)-\frac{f(b)-f(a)}{b-a}(x-a).

Its endpoint values are both zero, so Rolle gives an interior cc with

f′(c)=f(b)−f(a)b−a.f'(c)=\frac{f(b)-f(a)}{b-a}.

The continuity and differentiability assumptions remain essential. For f(x)=∣x∣f(x)=|x| on [−1,1][-1,1], the endpoint secant slope is zero but no derivative equals zero: slopes are −1-1 and 11, and the derivative fails at zero.

For f(x)=x2f(x)=x^2 on [1,3][1,3], the secant slope is 44, matched by f′(2)=4f'(2)=4. MVT guarantees at least one matching point, not necessarily a unique one.

Consequences that calculus repeatedly uses

If f′>0f'>0 everywhere on an interval, any two inputs u<vu<v satisfy f(v)−f(u)=f′(c)(v−u)>0f(v)-f(u)=f'(c)(v-u)>0. Thus ff is strictly increasing. If f′=0f'=0 throughout an interval, the same equation shows it is constant. If f′=g′f'=g' on an interval, apply the zero-derivative result to f−gf-g: the functions differ by a constant.

“On an interval” cannot be dropped. The function 11 on negative inputs and 22 on positive inputs has derivative zero throughout its disconnected domain, but is not one constant across that domain. Each connected interval can carry its own constant.

If ∣f′∣≤M|f'|\leq M throughout the interval, MVT yields ∣f(v)−f(u)∣≤M∣v−u∣|f(v)-f(u)|\leq M|v-u|. This turns a derivative bound into a finite-change bound. For example, ∣(sin⁡x)′∣=∣cos⁡x∣≤1|(\sin x)'|=|\cos x|\leq1 implies ∣sin⁡v−sin⁡u∣≤∣v−u∣|\sin v-\sin u|\leq|v-u| for radian inputs.

Linearization with an error estimate

The local linear approximation at aa is L(x)=f(a)+f′(a)(x−a)L(x)=f(a)+f'(a)(x-a). For f(x)=xf(x)=\sqrt{x} near 44, L(x)=2+(x−4)/4L(x)=2+(x-4)/4, so 4.1≈2.025\sqrt{4.1}\approx2.025.

If ff is twice continuously differentiable on an interval joining aa and xx and ∣f′′∣≤M|f''|\leq M there, Taylor’s first-order remainder bound gives

∣f(x)−L(x)∣≤M2∣x−a∣2.|f(x)-L(x)|\leq\frac M2|x-a|^2.

One proof follows from the integral remainder after integration is established; a Cauchy mean value theorem argument also proves it. We state the bound here with its extra smoothness assumption rather than treating it as part of the derivative definition. For the square-root example on [4,4.1][4,4.1], ∣f′′(u)∣=1/(4u3/2)≤1/32|f''(u)|=1/(4u^{3/2})\leq1/32, giving error at most 0.12/64=0.000156250.1^2/64=0.00015625. The negative second derivative also shows the tangent lies above this concave-down curve.

Exercises

  1. Find a point supplied by MVT for x3x^3 on [0,2][0,2].
  2. Explain precisely why Rolle’s theorem fails for ∣x∣|x| on [−1,1][-1,1].
  3. Prove that two antiderivatives of 2x2x on R\mathbb R differ by a constant.
  4. Linearize log⁡x\log x at 11 and approximate log⁡(1.02)\log(1.02).
  5. Bound that approximation error using ∣f′′∣≤1|f''|\leq1 on [1,1.02][1,1.02].
  6. If a differentiable function has ∣f′∣≤3|f'|\leq3 on [2,5][2,5], bound ∣f(5)−f(2)∣|f(5)-f(2)|.
Worked solutions
  1. The secant slope is 44. Solve 3c2=43c^2=4 in (0,2)(0,2): c=2/3c=2/\sqrt3.
  2. It is continuous with equal endpoint values, but not differentiable at the interior point zero.
  3. Their difference has derivative zero throughout the interval R\mathbb R; MVT forces that difference to be constant.
  4. L(x)=x−1L(x)=x-1, giving 0.020.02.
  5. The bound is (1/2)(0.02)2=0.0002(1/2)(0.02)^2=0.0002.
  6. At most 3(5−2)=93(5-2)=9.

Open this chapter with navigation and any associated laboratory →

Chapter 13

Extrema, Concavity, and Curve Analysis

Use sign charts and endpoint checks to classify behavior instead of treating every zero derivative as an extremum.

Prerequisites: MVT, first and second derivatives, and continuity on closed intervals.

  • Identify critical numbers in the domain and test derivative signs.
  • Distinguish local and absolute extrema.
  • Diagnose concavity and genuine inflection points.

Setting a derivative equal to zero produces candidates, not an optimization solution. A stationary point may be a maximum, minimum, or neither. An extreme value can also occur at a nondifferentiable point or an included endpoint. A systematic analysis keeps all of these possibilities visible.

Critical numbers and the first derivative test

A critical number is an interior input in the function’s domain where the derivative is zero or does not exist. An input excluded from the domain is not a critical number, although it can divide intervals of behavior. If ff is continuous at a critical number and f′f' changes from positive to negative across it, ff has a local maximum there. Continuity at the critical number together with a negative-to-positive derivative change gives a local minimum. The conclusion follows from monotonicity on each side, established through MVT.

Continuity at the critical point cannot be omitted. Define q(x)=−∣x∣q(x)=-|x| for x≠0x\ne0 but set q(0)=−1q(0)=-1. The derivative is positive immediately to the left and negative immediately to the right, yet q(0)q(0) is a strict local minimum because nearby nonzero outputs are close to zero and greater than −1-1. Derivative signs control behavior on the neighboring intervals; continuity connects that behavior to the actual value at the join.

For f(x)=x3−3xf(x)=x^3-3x, f′=3(x−1)(x+1)f'=3(x-1)(x+1). It is positive on (−∞,−1)(-\infty,-1), negative on (−1,1)(-1,1), and positive on (1,∞)(1,\infty). Thus f(−1)=2f(-1)=2 is a local maximum and f(1)=−2f(1)=-2 a local minimum. Neither is an absolute extremum on all of R\mathbb R, where the cubic is unbounded in both directions.

On the closed interval [−2,3][-2,3], EVT guarantees absolute extrema. Evaluate all interior critical points and both endpoints: f(−2)=−2f(-2)=-2, f(−1)=2f(-1)=2, f(1)=−2f(1)=-2, and f(3)=18f(3)=18. The absolute minimum is −2-2, attained at both −2-2 and 11; the absolute maximum is 1818 at 33. An algorithm that records only one minimizing input would lose a valid answer.

Concavity and second derivatives

Where f′f' is increasing, the graph is concave up; where f′f' is decreasing, it is concave down. A positive second derivative on an interval is sufficient for concavity up, and a negative second derivative for concavity down. An inflection point is a point on the graph where concavity changes across the point. The function must be continuous at the point under the convention used here. Merely solving f′′=0f''=0 does not establish this change.

For the cubic above, f′′=6xf''=6x changes sign at zero, so (0,0)(0,0) is an inflection point. For x4x^4, f′′=12x2f''=12x^2 is zero at zero but positive on both sides. There is no inflection there. For x1/3x^{1/3}, concavity changes at zero even though the second derivative is undefined there; a point of the graph need not have a second derivative to be an inflection point under this definition.

The second derivative test says: if f′(c)=0f'(c)=0 and f′′(c)>0f''(c)>0, a local minimum occurs; if f′′(c)<0f''(c)<0, a local maximum occurs, under differentiability in a neighborhood. When f′′(c)=0f''(c)=0, the test is inconclusive. The functions x4x^4, −x4-x^4, and x3x^3 demonstrate all three possibilities at zero.

An organized curve investigation

Begin with the domain. Check symmetry and intercepts when useful. Investigate discontinuities and end behavior. Factor the first derivative to make a sign chart. Use the second derivative to inspect concavity. Evaluate important points and distinguish conclusions proved by signs from features merely suggested by a sketch.

For g(x)=x+1/xg(x)=x+1/x on x≠0x\ne0, g′=1−1/x2g'=1-1/x^2 vanishes at ±1\pm1. It is negative for 0<∣x∣<10<|x|<1 and positive for ∣x∣>1|x|>1. Thus g(−1)=−2g(-1)=-2 is a local maximum and g(1)=2g(1)=2 a local minimum. Its disconnected domain prevents treating those as global extrema: near zero the branches diverge, and their end behaviors are unbounded. Domain and asymptote checks come before any global claim.

Exercises

  1. Analyze increasing and decreasing intervals of x3−6x2+9xx^3-6x^2+9x.
  2. Find its local extrema and inflection point.
  3. Find the absolute extrema of x2−4x+1x^2-4x+1 on [0,5][0,5].
  4. Explain why x4x^4 defeats the inference “zero second derivative means inflection.”
  5. Find the absolute minimum of ∣x−2∣|x-2| on [0,3][0,3] and explain which derivative test requires care.
  6. Explain why an excluded vertical-asymptote input is not a critical number.
Worked solutions
  1. f′=3(x−1)(x−3)f'=3(x-1)(x-3): increasing before 11 and after 33, decreasing between them.
  2. Local maximum f(1)=4f(1)=4, local minimum f(3)=0f(3)=0. Since f′′=6x−12f''=6x-12 changes sign at 22, the inflection point is (2,2)(2,2).
  3. The stationary point is 22 with value −3-3. Endpoint values are 11 and 66, giving absolute minimum −3-3 at 22 and maximum 66 at 55.
  4. Its second derivative is zero only at zero but concavity is up on both sides.
  5. The minimum is zero at 22, where the derivative does not exist. Setting only f′=0f'=0 would miss it.
  6. A critical number must belong to the function’s domain. An asymptote still separates intervals to analyze.

Open this chapter with navigation and any associated laboratory →

Chapter 14

Optimization as a Modeling Problem

Translate a constraint into a feasible domain, then justify a global optimum with calculus and boundary analysis.

Prerequisites: Critical numbers, EVT, sign charts, and elementary geometry.

  • Separate objective, constraint, and feasible domain.
  • Reduce a constrained problem to one variable.
  • Establish a global optimum rather than merely a stationary point.

An optimization problem does not begin with a derivative. It begins by deciding what is being optimized, which quantities may change, and what restrictions connect them. Different assumptions create different mathematics even when the objects look similar.

A rectangle beside a wall

A rectangular enclosure needs fencing on three sides because an existing wall supplies the fourth. Let xx be each perpendicular side and yy the parallel fenced side, in metres. With 6060 metres of fencing, 2x+y=602x+y=60. The objective is area A=xyA=xy. Substituting the constraint gives

A(x)=x(60−2x),0≤x≤30.A(x)=x(60-2x),\qquad 0\leq x\leq30.

The closed interval includes degenerate rectangles so that boundary behavior is explicit. The polynomial is continuous, so EVT gives a global maximum. Since A′=60−4xA'=60-4x, the only interior stationary point is x=15x=15, giving y=30y=30 and area 450450 square metres. The endpoint areas are zero. This comparison establishes the absolute maximum. Alternatively, A=450−2(x−15)2A=450-2(x-15)^2 proves the same conclusion without differentiation.

If fencing were required on all four sides, the constraint would be 2x+2y=602x+2y=60 and the maximizing rectangle would be a square of side 1515, area 225225. A correct derivative applied to the wrong constraint solves the wrong problem perfectly.

An open box

Cut squares of side xx from the corners of a 1212 by 1212 centimetre sheet and fold the sides. The volume is

V(x)=x(12−2x)2,0≤x≤6.V(x)=x(12-2x)^2,\qquad 0\leq x\leq6.

Differentiation in factored form gives

V′=(12−2x)(12−6x).V'=(12-2x)(12-6x).

The interior stationary point is x=2x=2, giving dimensions 88 by 88 by 22 and volume 128128 cubic centimetres. Both endpoints give zero volume. The factor 12−2x12-2x also vanishes at 66, but that is an endpoint, not an additional interior stationary box. We have checked the entire feasible interval, including the cases that are geometrically degenerate.

Minimizing on an unbounded interval

A closed cylinder has prescribed volume V0>0V_0>0. With radius r>0r>0 and height hh, πr2h=V0\pi r^2h=V_0 and surface area is

S(r)=2πr2+2V0r.S(r)=2\pi r^2+\frac{2V_0}r.

The derivative vanishes when 4πr−2V0/r2=04\pi r-2V_0/r^2=0, or r3=V0/(2π)r^3=V_0/(2\pi). At this radius, the constraint gives h=2rh=2r. Because S′′=4π+4V0/r3>0S''=4\pi+4V_0/r^3>0 throughout r>0r>0, the derivative is increasing and the stationary point is a global minimum. Also S→∞S\to\infty as r→0+r\to0^+ or r→∞r\to\infty, ruling out an escaped finite boundary optimum. EVT alone cannot be invoked on the open unbounded domain.

The model assumes the same material cost per unit area for the side and both ends, negligible thickness, and no seam allowance. Different manufacturing costs alter the objective and may alter the optimum.

Sensitivity and discrete decisions

A continuous optimum can guide a decision constrained to integer quantities, but rounding must respect the feasible set. Evaluate neighboring feasible integers, not just the nearest decimal in isolation. If an objective is unimodal, those nearby candidates often suffice; justify that structure before excluding distant alternatives.

Dimensions and parameters also matter. For fencing length PP, the wall problem gives x=P/4x=P/4, y=P/2y=P/2, and maximum area P2/8P^2/8. Doubling available fencing quadruples the maximum area in this model. A parameterized result conveys more than one isolated numerical answer.

Exercises

  1. Maximize the area of a four-sided rectangle with perimeter 8080 metres.
  2. Repeat the three-sided wall problem for fencing length P>0P>0.
  3. For a square sheet of side L>0L>0, determine the cut that maximizes the open-box volume.
  4. A cylinder must have volume 16π16\pi cubic centimetres. Find the minimum-area radius and height.
  5. Explain why finding A′(x)=0A'(x)=0 does not by itself prove an absolute maximum.
  6. Identify an assumption that would change the cylinder’s objective without changing its volume constraint.
Worked solutions
  1. Write A=x(40−x)A=x(40-x) on [0,40][0,40]. The maximum is 400400 square metres at x=y=20x=y=20.
  2. A=x(P−2x)A=x(P-2x); its stationary point is P/4P/4, giving y=P/2y=P/2 and area P2/8P^2/8. Endpoint areas vanish.
  3. V=x(L−2x)2V=x(L-2x)^2 on [0,L/2][0,L/2], with derivative (L−2x)(L−6x)(L-2x)(L-6x). The maximum occurs at x=L/6x=L/6, giving volume 2L3/272L^3/27.
  4. r3=16π/(2π)=8r^3=16\pi/(2\pi)=8, so r=2r=2 and h=4h=4 centimetres.
  5. Stationary points can be minima or neither; boundaries, nondifferentiable points, and nonattainment must be checked.
  6. More expensive end material, seam allowance, a missing lid, or thickness-dependent costs would change the surface-cost expression.

Open this chapter with navigation and any associated laboratory →

Chapter 15

Indeterminate Forms and Newton's Method

Use derivative-based limit and root tools while keeping their hypotheses, failure modes, and approximation status visible.

Prerequisites: MVT, derivative rules, limits, and tangent lines.

  • Recognize the forms to which L’Hôpital's rule applies.
  • Derive Newton's update from a tangent approximation.
  • Diagnose cycles and distinguish residual from root error.

Derivatives can help with limits and nonlinear equations, but neither operation becomes unconditional. A limit theorem requires hypotheses; an iterative method requires a starting point and convergence analysis. Getting a plausible decimal is not the same as having either guarantee.

L’Hôpital’s rule

One useful version is this: suppose ff and gg are differentiable on a deleted one-sided neighborhood of aa, g′g' is nonzero there, and both ff and gg tend to zero, or both have infinite magnitude. If f′/g′f'/g' tends to a finite or extended-real limit LL, then f/gf/g tends to LL. Corresponding versions apply to limits at infinity. For a two-sided result, establish the necessary behavior from both sides. The theorem follows from a generalized mean value theorem; it is a statement about a quotient of derivatives, not the derivative of the quotient.

Thus

lim⁡x→0ex−1x=lim⁡x→0ex1=1.\lim_{x\to0}\frac{e^x-1}{x} =\lim_{x\to0}\frac{e^x}{1}=1.

This is a consistency check after exponential differentiation has been established, not an independent proof of the normalization used earlier. Using a theorem to prove one of its own prerequisites would be circular.

For lim⁡x→∞log⁡x/x\lim_{x\to\infty}\log x/x, both numerator and denominator grow without bound and the derivative ratio is 1/x→01/x\to0. In contrast, applying the rule to (x2+1)/x(x^2+1)/x at zero is invalid: the numerator tends to 11, not zero or infinity. The original quotient has opposite unbounded one-sided behavior, while the derivative ratio 2x2x tends to zero. The missing hypothesis changes the answer.

Rewrite other forms first

A product of a vanishing factor and an unbounded factor is not a quotient to which the rule directly applies. For xlog⁡xx\log x as x→0+x\to0^+, write (log⁡x)/(1/x)(\log x)/(1/x). The derivative ratio is (1/x)/(−1/x2)=−x→0(1/x)/(-1/x^2)=-x\to0, so the product tends to zero.

For a variable power, logarithms can expose a quotient. If y=xxy=x^x for x>0x>0, then log⁡y=xlog⁡x→0\log y=x\log x\to0, so y→1y\to1. Do not treat 000^0 as a numerical evaluation; the particular limiting expression has been analyzed.

Newton’s tangent update

To solve f(x)=0f(x)=0, approximate the graph at a current guess xnx_n by its tangent f(xn)+f′(xn)(x−xn)f(x_n)+f'(x_n)(x-x_n). The tangent’s zero is

xn+1=xn−f(xn)f′(xn),x_{n+1}=x_n-\frac{f(x_n)}{f'(x_n)},

provided f′(xn)≠0f'(x_n)\ne0. For f(x)=x2−2f(x)=x^2-2, the update is (xn+2/xn)/2(x_n+2/x_n)/2. Starting at x0=1x_0=1 gives 1.51.5, 1.416666…1.416666\ldots, and 1.414215686…1.414215686\ldots, approaching 2\sqrt2.

Near a simple root rr, a twice continuously differentiable function with f′(r)≠0f'(r)\ne0 has a neighborhood in which suitably close starting guesses converge quadratically. Roughly, the next error is proportional to the square of the present error. This is a local theorem, not a guarantee from every starting value.

For f(x)=x3−2x+2f(x)=x^3-2x+2, starting at zero gives x1=1x_1=1 and x2=0x_2=0, an exact two-cycle. At a point with zero derivative, the displayed update may be undefined. For functions with several roots, a different starting guess can approach a different root or fail to converge.

Certifying an approximation

A small residual ∣f(xn)∣|f(x_n)| need not mean small distance to a root if the graph is flat. If a root rr and estimate xnx_n lie in an interval where ∣f′∣≥m>0|f'|\geq m>0, MVT yields

∣xn−r∣≤∣f(xn)∣/m.|x_n-r|\leq|f(x_n)|/m.

Bisection offers a bracket-based error guarantee without requiring a nonzero derivative. A robust numerical procedure can retain a sign-changing bracket and accept Newton steps only when they stay within it, falling back to bisection otherwise. Root approximation and root existence remain separate pieces of the argument.

Exercises

  1. Evaluate lim⁡x→0(sin⁡x−x)/x3\lim_{x\to0}(\sin x-x)/x^3 by repeated justified applications of L’Hôpital’s rule.
  2. Evaluate lim⁡x→∞x/ex\lim_{x\to\infty}x/e^x.
  3. Explain why the rule cannot be applied directly to sin⁡x/(1+x)\sin x/(1+x) at zero.
  4. Take two Newton steps for x2−3=0x^2-3=0 starting at 22.
  5. Verify the two-cycle for x3−2x+2x^3-2x+2 starting at zero.
  6. If ∣f′∣≥2|f'|\geq2 between an estimate and a known root and the residual is 0.0060.006, bound root error.
Worked solutions
  1. The derivative ratios are (cos⁡x−1)/(3x2)(\cos x-1)/(3x^2), −sin⁡x/(6x)-\sin x/(6x), and −cos⁡x/6-\cos x/6, yielding −1/6-1/6. Each earlier stage is 0/00/0, with nonzero derivative denominator on each deleted side.
  2. The ∞/∞\infty/\infty form gives derivative ratio 1/ex→01/e^x\to0.
  3. The denominator tends to 11, so ordinary quotient limits give zero. The required indeterminate form is absent.
  4. x1=2−(4−3)/4=7/4x_1=2-(4-3)/4=7/4; x2=(7/4+12/7)/2=97/56≈1.732142857x_2=(7/4+12/7)/2=97/56\approx1.732142857.
  5. f(0)=2f(0)=2, f′(0)=−2f'(0)=-2, so the next value is 11. Then f(1)=1f(1)=1, f′(1)=1f'(1)=1, returning zero.
  6. At most 0.0030.003.

Open this chapter with navigation and any associated laboratory →

Chapter 16

Antiderivatives and Initial Values

Reverse differentiation on intervals, retain constants, and connect acceleration, velocity, and position.

Prerequisites: Differentiation rules and the constant-difference consequence of MVT.

  • Find and verify families of elementary antiderivatives.
  • Use initial conditions to determine constants.
  • Distinguish displacement from total distance in a motion model.

An accelerometer can describe how velocity changes. Recovering velocity also requires an initial velocity; recovering position requires another initial condition. This is a general feature of reversing differentiation: a derivative loses information about additive constants.

An antiderivative of ff on an interval is a differentiable function FF satisfying F′=fF'=f there. If FF is one antiderivative, every other one on the same interval has the form F+CF+C, by the mean value theorem. We write

∫f(x) dx=F(x)+C.\int f(x)\,dx=F(x)+C.

This indefinite integral denotes a family of functions, not yet the definite integral of accumulation introduced next. The symbol dxdx records the variable of integration. Changing it without changing the function’s interpretation can change the problem.

Reverse rules and verify by differentiation

For r≠−1r\ne-1, on an interval where the real powers and derivatives are defined,

∫xr dx=xr+1r+1+C.\int x^r\,dx=\frac{x^{r+1}}{r+1}+C.

The exceptional exponent matters: division by r+1r+1 fails at r=−1r=-1. Instead ∫(1/x) dx=log⁡∣x∣+C\int(1/x)\,dx=\log|x|+C on an interval avoiding zero. The disconnected domain permits different constants on the positive and negative half-lines.

Other elementary pairs are exe^x, whose antiderivative is exe^x; cos⁡x\cos x, whose antiderivative is sin⁡x\sin x; and sin⁡x\sin x, whose antiderivative is −cos⁡x-\cos x. Linearity allows termwise integration of finite linear combinations. Products do not integrate by multiplying antiderivatives: (∫f)(∫g)(\int f)(\int g) differentiates with a product-rule extra term.

For example,

∫(6x2−4x+3) dx=2x3−2x2+3x+C.\int(6x^2-4x+3)\,dx=2x^3-2x^2+3x+C.

Differentiate the answer to check every term. A symbolic answer that differentiates to the wrong expression is wrong regardless of how plausible its form looks.

Initial-value information

Suppose F′(x)=3x2−2F'(x)=3x^2-2 and F(1)=5F(1)=5. Antidifferentiation gives F=x3−2x+CF=x^3-2x+C. At 11 this is C−1C-1, so C=6C=6. Thus F=x3−2x+6F=x^3-2x+6. Without the initial value, there would be infinitely many vertical translates satisfying the derivative equation.

For a constant acceleration a(t)=−9.8a(t)=-9.8 metres per second squared, an initial upward velocity v(0)=19.6v(0)=19.6 gives v(t)=19.6−9.8tv(t)=19.6-9.8t. If s(0)=1.5s(0)=1.5 metres, the position is s(t)=1.5+19.6t−4.9t2s(t)=1.5+19.6t-4.9t^2, over the time interval where the free-flight model applies. At the highest point v=0v=0, so t=2t=2 seconds and s=21.1s=21.1 metres. The velocity initial condition and position initial condition determine two separate constants.

Direction changes and distance

For v(t)=2t−2v(t)=2t-2 on [0,3][0,3], a position antiderivative is s(t)=t2−2t+Cs(t)=t^2-2t+C. Displacement is s(3)−s(0)=3s(3)-s(0)=3 units of length. Velocity is negative on [0,1)[0,1) and positive on (1,3](1,3], so the object first moves backward and then forward. Its distance from time zero to one is ∣s(1)−s(0)∣=1|s(1)-s(0)|=1, and from one to three is ∣s(3)−s(1)∣=4|s(3)-s(1)|=4, giving total distance 55.

Taking the absolute value of final displacement would give 33 and miss the reversal. Integration of speed will supply a uniform way to handle such cases. First identify where the sign changes so that the absolute value can be treated correctly.

Exercises

  1. Find all antiderivatives of 4x3−2x+54x^3-2x+5 on R\mathbb R.
  2. Find an antiderivative of x−2+x−1x^{-2}+x^{-1} on (0,∞)(0,\infty).
  3. Solve F′=2x+3F'=2x+3 with F(0)=−4F(0)=-4.
  4. An object has acceleration 6t6t, initial velocity −2-2, and initial position 55. Find velocity and position.
  5. For v(t)=t−2v(t)=t-2 on [0,5][0,5], find displacement and total distance.
  6. Explain why one constant need not describe antiderivatives across a disconnected domain.
Worked solutions
  1. x4−x2+5x+Cx^4-x^2+5x+C.
  2. −1/x+log⁡x+C-1/x+\log x+C; differentiate to verify both terms.
  3. F=x2+3x−4F=x^2+3x-4.
  4. v=3t2−2v=3t^2-2 and s=t3−2t+5s=t^3-2t+5.
  5. A position function is t2/2−2tt^2/2-2t. Displacement is 2.52.5; backward distance before 22 is 22, forward distance after is 4.54.5, giving total distance 6.56.5.
  6. MVT connects two inputs through an interval in the domain. It cannot cross a missing point, so each component interval may have a separate constant.

Open this chapter with navigation and any associated laboratory →

Chapter 17

Accumulation and the Definite Integral

Define accumulation through tagged sums, calculate a sum limit exactly, and state what integrability requires.

Prerequisites: Limits, finite sums, continuity, and elementary area.

  • Construct tagged Riemann sums with units and signs.
  • Define the definite integral independently of antiderivatives.
  • Evaluate the integral of a square by a finite-sum limit.

A pump’s flow rate is measured in litres per minute. Multiplying a constant rate by time gives transferred volume. For a changing rate, dividing the time interval into short pieces gives a sum of approximate transfers. The definite integral is the limit that makes this accumulation independent of the partition and sample choices, when such a limit exists.

Partitions and tagged sums

A partition a=x0<x1<⋯<xn=ba=x_0<x_1<\cdots<x_n=b divides [a,b][a,b] into widths Δxi=xi−xi−1\Delta x_i=x_i-x_{i-1}. Choose a tag ξi\xi_i in each subinterval. The Riemann sum is

∑i=1nf(ξi)Δxi.\sum_{i=1}^n f(\xi_i)\Delta x_i.

The mesh is max⁡iΔxi\max_i\Delta x_i. A bounded function is Riemann integrable when there is a real number II such that every tagged sum is within any requested tolerance of II once the mesh is sufficiently small. We denote that number by ∫abf(x) dx\int_a^b f(x)\,dx. The definition quantifies over all sufficiently fine partitions and tags; convergence of one favored sequence of sums is not by itself the full definition.

Continuous functions on closed bounded intervals are integrable. Bounded monotone functions are also integrable, even if they have jumps. Not every bounded function is: the function equal to one at rational inputs and zero at irrational inputs on [0,1][0,1] admits tagged sums identically one or identically zero at every mesh, because both types of inputs occur in every subinterval. No common limit exists.

A square integrated from sums

For f(x)=x2f(x)=x^2 on [0,1][0,1], use nn equal pieces and right endpoints i/ni/n. The sum is

Rn=1n∑i=1n(in)2=n(n+1)(2n+1)6n3=13+12n+16n2.R_n=\frac1n\sum_{i=1}^n\left(\frac in\right)^2 =\frac{n(n+1)(2n+1)}{6n^3} =\frac13+\frac1{2n}+\frac1{6n^2}.

We used ∑i2=n(n+1)(2n+1)/6\sum i^2=n(n+1)(2n+1)/6, which can be proved by induction. Thus Rn→1/3R_n\to1/3. The left sum is Ln=1/3−1/(2n)+1/(6n2)L_n=1/3-1/(2n)+1/(6n^2) and tends to the same limit. Because x2x^2 is increasing, every tagged sum for this equal partition lies between these endpoint sums; more generally, continuity supplies integrability for arbitrary fine partitions. Hence ∫01x2 dx=1/3\int_0^1x^2\,dx=1/3.

The sums give a finite-error statement too: Rn−Ln=1/nR_n-L_n=1/n. At n=10n=10, the left and right values are 0.2850.285 and 0.3850.385. Both bracket 1/31/3; neither is the exact integral. A picture of ten rectangles can illustrate the approximation without converting it into equality.

Signed accumulation and units

If ff is negative on part of the interval, its rectangular contributions are negative. Therefore the integral represents signed accumulation. Geometric area between a graph and the axis is ∫ab∣f(x)∣ dx\int_a^b|f(x)|\,dx, not necessarily ∫abf(x) dx\int_a^bf(x)\,dx.

For f(x)=xf(x)=x on [−1,1][-1,1], symmetry gives signed integral zero, although geometric area is one. If ff is velocity, the integral gives displacement; integrating ∣f∣|f| gives distance. If ff is flow in litres per minute and xx is minutes, the integral has units litres. An integral’s units multiply the integrand’s units by those of its input variable.

Properties from sums

Integrals are linear, additive across adjacent intervals, and preserve inequalities for integrable functions. Define ∫baf=−∫abf\int_b^af=-\int_a^bf and ∫aaf=0\int_a^af=0. The reversed orientation is a sign convention consistent with additivity, not a negative geometric interval length in the original partition definition.

If m≤f≤Mm\leq f\leq M on [a,b][a,b], then m(b−a)≤∫abf≤M(b−a)m(b-a)\leq\int_a^bf\leq M(b-a). Also ∣∫f∣≤∫∣f∣|\int f|\leq\int|f|. These bounds can detect an impossible answer before an exact calculation is available.

Exercises

  1. Write the right-endpoint sum for f(x)=3xf(x)=3x on [0,2][0,2] and calculate its limit.
  2. Compute left and right sums for x2x^2 on [0,1][0,1] with n=4n=4.
  3. If 2≤f(x)≤52\leq f(x)\leq5 on [1,4][1,4], bound its integral.
  4. Explain why an integral of a velocity can be zero while distance travelled is positive.
  5. Show that changing a bounded integrable function at one point does not alter its integral.
  6. Why does the rational-indicator example not become integrable when a single fine regular grid happens to choose only irrational tags?
Worked solutions
  1. The width is 2/n2/n, tag 2i/n2i/n, giving (12/n2)∑i=6(n+1)/n→6(12/n^2)\sum i=6(n+1)/n\to6.
  2. L4=0.21875L_4=0.21875 and R4=0.46875R_4=0.46875, using the displayed formulas.
  3. Between 66 and 1515.
  4. Opposite-direction displacements can cancel. Distance adds their magnitudes.
  5. If the changed value differs by at most BB, only subintervals touching that point can be affected. Their total contribution difference is bounded by 2B2B times the mesh, tending to zero. Interior tags usually affect one piece; allowing two handles shared endpoints.
  6. Integrability requires all sufficiently fine tagged sums to converge to the same number. Rational tags still give one, while irrational tags give zero.

Open this chapter with navigation and any associated laboratory →

Chapter 18

The Fundamental Theorem of Calculus

Prove that continuous accumulation has the original rate, and derive evaluation by antiderivatives.

Prerequisites: Riemann integrals, continuity, antiderivatives, and the chain rule.

  • Prove the accumulation derivative theorem for continuous integrands.
  • Derive the evaluation theorem and check its hypotheses.
  • Differentiate integrals with variable bounds and interpret net change.

Differentiation began with local rates. Integration began with sums of accumulation. The fundamental theorem explains why these operations are inverse in the appropriate sense. It is the reason an area problem can be solved by finding a function with a prescribed derivative.

Accumulation has the original rate

Let ff be continuous on [a,b][a,b] and define A(x)=∫axf(t) dtA(x)=\int_a^xf(t)\,dt. For an interior xx and small nonzero hh,

A(x+h)−A(x)h=1h∫xx+hf(t) dt.\frac{A(x+h)-A(x)}h=\frac1h\int_x^{x+h}f(t)\,dt.

Subtract f(x)f(x) and use an integral bound. For either sign of hh,

∣A(x+h)−A(x)h−f(x)∣≤sup⁡t between x and x+h∣f(t)−f(x)∣.\left|\frac{A(x+h)-A(x)}h-f(x)\right| \leq\sup_{t\text{ between }x\text{ and }x+h}|f(t)-f(x)|.

Continuity at xx makes the right side tend to zero. Therefore A′(x)=f(x)A'(x)=f(x). At endpoints use the matching one-sided version. The proof compares a short-interval average with the local value. It does not assume ff is constant on the short interval.

Continuity is sufficient, not necessary in every possible case. At a jump, however, the conclusion can fail. If f(t)=0f(t)=0 for t<0t<0 and 11 for t≥0t\geq0, the accumulation from −1-1 is zero for x≤0x\leq0 and xx for x>0x>0. It has no derivative at zero. Substituting f(0)=1f(0)=1 into a theorem whose continuity hypothesis fails would be invalid.

Evaluation by an antiderivative

If F′=fF'=f on the interval and ff is continuous, FF and AA have the same derivative. The constant-difference theorem gives F(x)−A(x)=F(a)F(x)-A(x)=F(a), since A(a)=0A(a)=0. At bb,

∫abf(x) dx=F(b)−F(a).\int_a^bf(x)\,dx=F(b)-F(a).

This is the evaluation form of the fundamental theorem. Different books interchange the labels “first” and “second” for its two parts; the statements, rather than the numbering, are what matter.

For example,

∫02(3x2−2x+1) dx=[x3−x2+x]02=6.\int_0^2(3x^2-2x+1)\,dx =[x^3-x^2+x]_0^2=6.

The constant in an indefinite integral cancels in the endpoint difference. There is no need to write “+C+C” after the numerical value of a definite integral.

Variable bounds use the chain rule

For continuous ff and differentiable gg with values in its integration domain,

ddx∫ag(x)f(t) dt=f(g(x))g′(x).\frac d{dx}\int_a^{g(x)}f(t)\,dt=f(g(x))g'(x).

With two variable bounds, split at a fixed base point:

ddx∫u(x)v(x)f(t) dt=f(v(x))v′(x)−f(u(x))u′(x).\frac d{dx}\int_{u(x)}^{v(x)}f(t)\,dt =f(v(x))v'(x)-f(u(x))u'(x).

Thus d[∫1x2cos⁡t dt]/dx=2xcos⁡(x2)d[\int_1^{x^2}\cos t\,dt]/dx=2x\cos(x^2), and d[∫x2xt2 dt]/dx=2(2x)2−x2=7x2d[\int_x^{2x}t^2\,dt]/dx=2(2x)^2-x^2=7x^2. The integration variable tt is a dummy variable; its name prevents confusing it with the moving bound.

For A(x)=∫0xe−t2 dtA(x)=\int_0^xe^{-t^2}\,dt, A′(x)=e−x2A'(x)=e^{-x^2} even though the integrand has no elementary antiderivative. An accumulation function can be mathematically well-defined and differentiable without having a formula in the elementary functions.

Net change and forbidden crossings

If a quantity QQ has continuous rate Q′Q', then Q(b)−Q(a)=∫abQ′(t) dtQ(b)-Q(a)=\int_a^bQ'(t)\,dt. For v(t)=2t−2v(t)=2t-2 on [0,3][0,3], the displacement is [t2−2t]03=3[t^2-2t]_0^3=3. Distance is −∫01v+∫13v=5-\int_0^1v+\int_1^3v=5.

Do not apply the evaluation formula across a singularity without checking integrability. Writing ∫−111/x2 dx=[−1/x]−11=−2\int_{-1}^11/x^2\,dx=[-1/x]_{-1}^1=-2 is invalid. The integrand is unbounded at zero and does not satisfy the stated theorem; the corresponding improper integral diverges positively. A negative answer for a nonnegative integrand is an immediate sign that a hypothesis has been violated.

Exercises

  1. Evaluate ∫13(2x+1) dx\int_1^3(2x+1)\,dx.
  2. Differentiate ∫0sin⁡x(1+t2) dt\int_0^{\sin x}(1+t^2)\,dt.
  3. Differentiate ∫x2x3et dt\int_{x^2}^{x^3}e^t\,dt.
  4. For A(x)=∫0xcos⁡t dtA(x)=\int_0^x\cos t\,dt, identify intervals where AA is increasing.
  5. Explain why ∫0xe−t2 dt\int_0^xe^{-t^2}\,dt does not need an elementary closed form to have a derivative.
  6. Identify the precise failure in evaluating ∫−11(1/x) dx\int_{-1}^1(1/x)\,dx by the endpoint difference of log⁡∣x∣\log|x|.
Worked solutions
  1. [x2+x]13=12−2=10[x^2+x]_1^3=12-2=10.
  2. (1+sin⁡2x)cos⁡x(1+\sin^2x)\cos x.
  3. 3x2ex3−2xex23x^2e^{x^3}-2xe^{x^2}.
  4. Where cos⁡x>0\cos x>0, namely (−π/2+2kπ,π/2+2kπ)(-\pi/2+2k\pi,\pi/2+2k\pi) for integers kk.
  5. Continuity of the integrand establishes the accumulation function and the derivative theorem applies directly.
  6. The integrand is undefined and unbounded at zero. The two separate one-sided improper integrals do not both converge; endpoint cancellation is not an ordinary integral. A symmetric principal value is a distinct concept.

Open this chapter with navigation and any associated laboratory →

Chapter 19

Substitution and Integral Structure

Reverse the chain rule, transform bounds correctly, and recognize when a proposed substitution has not simplified the integral.

Prerequisites: Chain rule, antiderivatives, and both forms of the fundamental theorem.

  • Derive substitution from the chain rule and FTC.
  • Transform an entire definite integral, including its bounds.
  • Verify antiderivatives and distinguish substitution from unsupported symbol cancellation.

The integral of 2xcos⁡(x2)2x\cos(x^2) contains a composition and the derivative of its inner function. The chain rule says d[sin⁡(x2)]/dx=2xcos⁡(x2)d[\sin(x^2)]/dx=2x\cos(x^2), so reversing it supplies an antiderivative. Substitution organizes that recognition.

The theorem behind the notation

If F′=fF'=f, then the chain rule gives (F∘g)′=f(g)g′(F\circ g)'=f(g)g'. Consequently,

∫f(g(x))g′(x) dx=F(g(x))+C.\int f(g(x))g'(x)\,dx=F(g(x))+C.

Under continuous integrand and continuously differentiable substitution hypotheses on [a,b][a,b], FTC yields

∫abf(g(x))g′(x) dx=∫g(a)g(b)f(u) du.\int_a^bf(g(x))g'(x)\,dx =\int_{g(a)}^{g(b)}f(u)\,du.

The substitution need not be one-to-one for this signed one-dimensional identity when these conditions hold: the chain-rule proof handles reversals automatically. More general change-of-variable theorems require their own conditions. The shorthand u=g(x)u=g(x) and du=g′(x)dxdu=g'(x)dx records this theorem; it is not permission to replace pieces of an integrand inconsistently.

Matching the differential factor

For ∫x1+x2 dx\int x\sqrt{1+x^2}\,dx, let u=1+x2u=1+x^2, giving du=2x dxdu=2x\,dx. Then

∫x1+x2 dx=12∫u1/2 du=13(1+x2)3/2+C.\int x\sqrt{1+x^2}\,dx =\frac12\int u^{1/2}\,du =\frac13(1+x^2)^{3/2}+C.

Differentiate the result: (1/3)(3/2)(1+x2)1/2(2x)=x1+x2(1/3)(3/2)(1+x^2)^{1/2}(2x)=x\sqrt{1+x^2}. This verifies the coefficient, the most common source of a substitution error.

For ∫3x2/(1+x3) dx\int 3x^2/(1+x^3)\,dx, use u=1+x3u=1+x^3, obtaining log⁡∣1+x3∣+C\log|1+x^3|+C on intervals avoiding x=−1x=-1. The absolute value and domain restriction remain necessary.

Definite bounds belong to the new variable

Evaluate ∫012xex2 dx\int_0^1 2x e^{x^2}\,dx. With u=x2u=x^2, the bounds transform from 0,10,1 to 0,10,1, so the result is ∫01eu du=e−1\int_0^1e^u\,du=e-1. For ∫122xex2 dx\int_1^2 2x e^{x^2}\,dx, the transformed bounds are 1,41,4, yielding e4−ee^4-e. The first example happens to preserve the numerical bounds; that coincidence is not a rule.

One may instead find an antiderivative in xx and then use the original bounds. What fails is evaluating a uu-antiderivative at bounds that still represent xx without transformation.

For ∫0π/2sin⁡xcos⁡x dx\int_0^{\pi/2}\sin x\cos x\,dx, let u=cos⁡xu=\cos x, du=−sin⁡x dxdu=-\sin x\,dx. Bounds become 11 and 00, so −∫10u du=1/2-\int_1^0u\,du=1/2. The negative sign and reversed bounds are two coordinated consequences of decreasing cosine; dropping both accidentally may conceal an error in a more complex problem.

Useful patterns and honest limits

The pattern g′/gg'/g integrates to log⁡∣g∣\log|g| on nonzero intervals. The pattern g′egg'e^g integrates to ege^g. The pattern g′/(1+g2)g'/(1+g^2) integrates to arctan⁡g\arctan g. Recognize the derivative factor, including constant multiples, rather than merely spotting a complicated inner expression.

For ∫ex2 dx\int e^{x^2}\,dx, choosing u=x2u=x^2 introduces dx=du/(2u)dx=du/(2\sqrt u) on a positive branch and does not produce the simple exponential integral. The missing factor 2x2x matters. A substitution that makes a formula longer is not progress, and not every elementary integrand has an elementary antiderivative. Definite integration and numerical approximation remain meaningful.

Exercises

  1. Integrate 6x(3x2+1)46x(3x^2+1)^4.
  2. Evaluate ∫02x/(1+x2) dx\int_0^2x/(1+x^2)\,dx.
  3. Evaluate ∫013x21+x3 dx\int_0^1 3x^2\sqrt{1+x^3}\,dx.
  4. Integrate cos⁡(5x)\cos(5x) and verify the factor.
  5. Evaluate ∫0πsin⁡xcos⁡x dx\int_0^\pi\sin x\cos x\,dx; explain why using u=sin⁡xu=\sin x is still consistent despite noninjectivity.
  6. Explain the missing factor in the incorrect assertion ∫ex2 dx=ex2+C\int e^{x^2}\,dx=e^{x^2}+C.
Worked solutions
  1. With u=3x2+1u=3x^2+1, the result is (3x2+1)5/5+C(3x^2+1)^5/5+C.
  2. [log⁡(1+x2)/2]02=(log⁡5)/2[\log(1+x^2)/2]_0^2=(\log5)/2.
  3. u=1+x3u=1+x^3 transforms the bounds to 1,21,2, giving [2u3/2/3]12=2(22−1)/3[2u^{3/2}/3]_1^2=2(2\sqrt2-1)/3.
  4. (1/5)sin⁡(5x)+C(1/5)\sin(5x)+C; the chain-rule factor 55 cancels 1/51/5 on differentiation.
  5. The integral is zero by [sin⁡2x/2]0π[\sin^2x/2]_0^\pi. The substitution identity gives ∫00u du=0\int_0^0u\,du=0 and accounts for the orientation reversal during the path.
  6. Differentiating the proposed answer gives 2xex22xe^{x^2}, which is a different integrand.

Open this chapter with navigation and any associated laboratory →

Chapter 20

Area, Volume, Work, and Average Value

Build integral models from thin pieces with correct geometry, signs, and units.

Prerequisites: Definite integrals, FTC, substitution, and basic geometric formulas.

  • Form area and volume integrals from cross-sectional geometry.
  • Distinguish signed work, net change, and geometric magnitude.
  • Calculate and justify the average value of a continuous function.

The most important step in an integral application is identifying a representative small contribution. Its dimensions and geometry determine the integrand. Once that model is right, the calculation may be routine. A perfectly evaluated integral of the wrong slices is still the wrong volume.

Area between curves

If f(x)≥g(x)f(x)\geq g(x) on [a,b][a,b], a thin vertical strip has approximate area [f(x)−g(x)]Δx[f(x)-g(x)]\Delta x. Summing and taking a limit gives ∫ab(f−g) dx\int_a^b(f-g)\,dx. If their order changes, split at intersection points or integrate ∣f−g∣|f-g|.

Between y=xy=x and y=x2y=x^2 on [0,1][0,1], the line lies above the parabola, so

A=∫01(x−x2) dx=12−13=16.A=\int_0^1(x-x^2)\,dx=\frac12-\frac13=\frac16.

The interval comes from solving x=x2x=x^2, not from choosing visually convenient endpoints. For y=xy=x and the axis on [−1,1][-1,1], signed integration cancels; geometric area requires splitting at zero.

Cross-sectional volume

For a solid with continuous cross-sectional area A(x)A(x) perpendicular to an axis,

V=∫abA(x) dx.V=\int_a^bA(x)\,dx.

Revolving the region under y=xy=x for 0≤x≤20\leq x\leq2 around the horizontal axis gives disks of radius xx, so V=∫02πx2 dx=8π/3V=\int_0^2\pi x^2\,dx=8\pi/3. The square comes from disk area; integrating circumference would have the wrong units and geometry.

Revolving the region between y=xy=x and y=x2y=x^2 on [0,1][0,1] around the horizontal axis gives washers with outer radius xx and inner radius x2x^2:

V=π∫01(x2−x4) dx=2π15.V=\pi\int_0^1(x^2-x^4)\,dx=\frac{2\pi}{15}.

One must subtract disk areas, not square the difference of radii. The expressions R2−r2R^2-r^2 and (R−r)2(R-r)^2 represent different shapes.

For the same region revolved about the vertical axis, cylindrical shells have radius xx, height x−x2x-x^2, and approximate volume 2πx(x−x2)Δx2\pi x(x-x^2)\Delta x. Hence V=2π∫01(x2−x3) dx=π/6V=2\pi\int_0^1(x^2-x^3)\,dx=\pi/6. This shell formula follows by summing thin annular strips; it is not interchangeable with the washer formula without changing the slicing geometry. Equivalent washers in the vertical variable have radii y\sqrt y and yy, producing π∫01(y−y2) dy=π/6\pi\int_0^1(y-y^2)\,dy=\pi/6 as a check.

Work with a varying force

For force F(x)F(x) acting along a one-dimensional displacement coordinate, work is W=∫abF(x) dxW=\int_a^bF(x)\,dx. Its units are force times distance. For an ideal spring with restoring-force magnitude kxkx and quasistatic stretching from zero to extension bb, the external work is ∫0bkx dx=kb2/2\int_0^bkx\,dx=kb^2/2.

If k=100k=100 newtons per metre and b=0.2b=0.2 metres, the work is 22 joules. Using the final force 2020 newtons across the whole extension would give 44 joules, double the correct result because the force started at zero and increased linearly. A restoring force acting on a freely moving object has the opposite sign; define which force’s work is being calculated.

Average value

For an integrable function on [a,b][a,b] with a<ba<b, define

favg=1b−a∫abf(x) dx.f_{\mathrm{avg}}=\frac1{b-a}\int_a^bf(x)\,dx.

The denominator converts accumulated output-times-input units back to output units. For x2x^2 on [0,2][0,2], the average is (1/2)(8/3)=4/3(1/2)(8/3)=4/3. This differs from the average of endpoint values, (0+4)/2=2(0+4)/2=2. Endpoint averaging is exact for a linear function but not generally for a curved one.

If ff is continuous, its minimum mm and maximum MM bound the average. IVT then guarantees some cc with f(c)=favgf(c)=f_{\mathrm{avg}}. This is the mean value theorem for integrals. It asserts existence, not uniqueness or that cc is the interval midpoint.

Exercises

  1. Find the area between y=2xy=2x and y=x2y=x^2 where they intersect.
  2. Revolve the region under x\sqrt x on [0,4][0,4] about the horizontal axis and find its volume.
  3. Find the volume of washers with outer radius 22 and inner radius xx on [0,1][0,1].
  4. Find the work to stretch a spring with k=50k=50 newtons per metre from extension 0.10.1 to 0.30.3 metres.
  5. Find the average of sin⁡x\sin x on [0,π][0,\pi] and an interior input attaining it.
  6. Explain why a volume integral with integrand measured in metres rather than square metres is dimensionally wrong when integrating with respect to metres.
Worked solutions
  1. Intersections are 0,20,2. The area is [x2−x3/3]02=4/3[x^2-x^3/3]_0^2=4/3.
  2. Disk area is π(x)2=πx\pi(\sqrt x)^2=\pi x; volume is 8π8\pi cubic units.
  3. π∫01(4−x2) dx=11π/3\pi\int_0^1(4-x^2)\,dx=11\pi/3 cubic units.
  4. 25(0.32−0.12)=225(0.3^2-0.1^2)=2 joules.
  5. Average 2/π2/\pi; one input is arcsin⁡(2/π)\arcsin(2/\pi), another is π−arcsin⁡(2/π)\pi-\arcsin(2/\pi).
  6. Length times length gives area, not volume. A slice’s cross-sectional area supplies the missing factor of length.

Open this chapter with navigation and any associated laboratory →

Chapter 21

A Calculus Investigation from Rates to Totals

Complete a connected model, defend its assumptions, and review proof, computation, and interpretation across Calculus I.

Prerequisites: The preceding limits, differentiation, and integration chapters.

  • Connect a rate model to cumulative change and average values.
  • Check a model's extrema, units, feasibility, and parameter sensitivity.
  • Complete a cumulative set of computational and proof problems.

A reservoir receives water at modeled rate r(t)=6t−t2r(t)=6t-t^2 cubic metres per hour for 0≤t≤60\leq t\leq6 hours. Water leaves at a constant rate of 55 cubic metres per hour. At time zero the reservoir contains 2020 cubic metres. These are hypothetical quantities chosen for an investigation, not observations from a real reservoir. Assume the reservoir can hold the modeled volume without overflowing during the interval.

Recover the stored quantity

The net rate is V′=6t−t2−5=4−(t−3)2V'=6t-t^2-5=4-(t-3)^2. Antidifferentiation and the initial condition give

V(t)=20+3t2−t33−5t.V(t)=20+3t^2-\frac{t^3}3-5t.

Check by differentiation and by V(0)=20V(0)=20. At six hours, V(6)=26V(6)=26, so net change is 66 cubic metres. The inflow itself accumulates ∫06(6t−t2) dt=36\int_0^6(6t-t^2)\,dt=36 cubic metres, while outflow accumulates 3030. The initial stock, gross inflow, gross outflow, and final stock are distinct quantities whose units agree.

Analyze the whole interval

The rate vanishes at t=1t=1 and 55. It is negative before 11, positive between 11 and 55, and negative after 55. Compare all critical points and endpoints:

V(0)=20,V(1)=53/3,V(5)=85/3,V(6)=26.V(0)=20,\quad V(1)=53/3,\quad V(5)=85/3,\quad V(6)=26.

Thus the minimum stock is 53/353/3 and the maximum 85/385/3 cubic metres. Both lie in the interior. The model never empties, and a capacity of at least 85/385/3 would avoid overflow under these assumptions. A capacity of 2525 would violate the no-overflow model; continuing the same stock formula would no longer describe actual stored volume after overflow begins.

The strongest net inflow occurs at t=3t=3 because V′=4−(t−3)2V'=4-(t-3)^2 is maximized there, with value 44 cubic metres per hour. This is not the time of greatest stock. Maximizing a quantity and maximizing its rate are different tasks.

Total variation and average stock

The stock decreases by 7/37/3 cubic metres before 11, increases by 32/332/3 between 11 and 55, then decreases by 7/37/3 before 66. Its total variation is 46/346/3 cubic metres, compared with net change 66. Total variation sums absolute changes and can be written ∫06∣V′∣ dt\int_0^6|V'|\,dt after splitting at the rate’s zeros.

The average stock is

16∫06V(t) dt=16[20t+t3−t412−5t22]06=23.\frac16\int_0^6V(t)\,dt =\frac16\left[20t+t^3-\frac{t^4}{12}-\frac{5t^2}2\right]_0^6 =23.

In this particular symmetric example, the endpoint average (20+26)/2(20+26)/2 also equals 2323. That coincidence is not a general rule for cubic stock functions; the integral supplies the justification here.

Investigate a changed outflow

Replace the outflow 55 by a parameter c≥0c\geq0. Then Vc′=6t−t2−cV_c'=6t-t^2-c, and its zeros, when 0≤c≤90\leq c\leq9, are 3±9−c3\pm\sqrt{9-c}. For c>9c>9, the net rate is negative throughout the interval. The final stock is 20+36−6c=56−6c20+36-6c=56-6c. To establish physical feasibility, justify the interval minimum rather than relying on an endpoint check without an interval argument. The model’s physical feasibility is an interval question.

Cumulative problems

  1. Prove lim⁡x→1x2=1\lim_{x\to1}x^2=1 using an explicit epsilon–delta construction.
  2. A function equals x2x^2 for x≠2x\ne2 and 99 at 22. Is it continuous or differentiable at 22?
  3. Differentiate log⁡(1+x2)sin⁡x\log(1+x^2)\sin x.
  4. Use MVT to prove that a differentiable function with zero derivative on an interval is constant there.
  5. Find the absolute extrema of x3−3xx^3-3x on [−2,2][-2,2].
  6. Evaluate ∫012x/(1+x2) dx\int_0^1 2x/(1+x^2)\,dx.
  7. Differentiate ∫xx21+t2 dt\int_{x}^{x^2}\sqrt{1+t^2}\,dt.
  8. For the reservoir, calculate the stock at t=3t=3 and explain why it is neither an absolute maximum nor minimum.
  9. Find the area between xx and x3x^3 on [0,1][0,1].
  10. Explain why ∫−111/x2 dx\int_{-1}^11/x^2\,dx cannot be evaluated by a single endpoint antiderivative difference.
Worked cumulative solutions
  1. If ∣x−1∣<1|x-1|<1, then ∣x+1∣<3|x+1|<3. Choose δ=min⁡{1,ε/3}\delta=\min\{1,\varepsilon/3\}. Factoring gives ∣x2−1∣<3∣x−1∣<ε|x^2-1|<3|x-1|<\varepsilon.
  2. Its limit is 44, unequal to its value 99. It is discontinuous and therefore not differentiable.
  3. [2x/(1+x2)]sin⁡x+log⁡(1+x2)cos⁡x[2x/(1+x^2)]\sin x+\log(1+x^2)\cos x.
  4. For any u<vu<v in the interval, MVT gives f(v)−f(u)=f′(c)(v−u)=0f(v)-f(u)=f'(c)(v-u)=0. Hence all values agree.
  5. Endpoint and critical values are −2,2,−2,2-2,2,-2,2 at −2,−1,1,2-2,-1,1,2 respectively. Maximum 22 occurs at −1,2-1,2; minimum −2-2 occurs at −2,1-2,1.
  6. Substitute u=1+x2u=1+x^2 to obtain log⁡2\log2.
  7. 2x1+x4−1+x22x\sqrt{1+x^4}-\sqrt{1+x^2}.
  8. V(3)=23V(3)=23. The rate there is positive and maximal, so stock is still increasing; the stock maximum occurs later at 55.
  9. ∫01(x−x3) dx=1/4\int_0^1(x-x^3)\,dx=1/4.
  10. The integrand is unbounded at zero; the continuous-integrand FTC does not apply and the improper integral diverges.

Independent extension

Choose a new constant outflow and capacity. Determine the earliest time the model would empty or overflow, if either occurs on [0,6][0,6]. Report the stock formula, all candidate crossing intervals, a justified numerical bracket where necessary, a derivative sign chart, and an explicit statement of where your physical model stops applying. Include a computational check, but retain an analytic explanation for each existence or uniqueness claim. This is a mathematical investigation, not an engineering recommendation for an actual water system.

Open this chapter with navigation and any associated laboratory →