Math ClubMath Club
v1 · padrão canônico

Lesson 100 — Integration Workshop Trimester 10 (ODEs)

Integrating workshop on ODEs: separable, first-order linear, populations, second-order with constant coefficients, vibrations, RLC, numerical Euler, and Newton cooling.

Used in: AP Calculus BC (USA) · Leistungskurs Mathematik Klasse 12 (Germany) · H2 Mathematics (Singapore) · Spécialité Maths Terminale (France)

y+p(x)y=q(x)1st-order linearay+by+cy=q(x)2nd-order const. coeff.\underbrace{y' + p(x)\,y = q(x)}_{\text{1st-order linear}} \quad \Big| \quad \underbrace{ay'' + by' + cy = q(x)}_{\text{2nd-order const. coeff.}}

Two structural models cover nearly every application from the trimester. The one on the left solves with integrating factor μ=epdx\mu = e^{\int p\,dx}. The one on the right solves via the characteristic equation aλ2+bλ+c=0a\lambda^2 + b\lambda + c = 0. Non-analytic cases receive Euler or RK4.

Choose your door

Rigorous notation, full derivation, hypotheses

Map of techniques and ODE decision criteria

"The existence and uniqueness theorem says: there is exactly one solution curve through each point (x₀, y₀) where f and ∂f/∂y are continuous." — Lebl, Notes on Diffy Qs §1.2

"The idea behind the integrating factor is to write the left-hand side as an exact derivative." — OpenStax Calculus Vol. 2 §4.5

Techniques table for the trimester:

ODECanonical formTechniqueGeneral solution
Separabley=g(x)h(y)y' = g(x)h(y)Separate and integrateImplicit
Linear 1sty+py=qy'+py=qμ=ep\mu=e^{\int p}y=1μμqdx+Cy = \frac{1}{\mu}\int\mu q\,dx + C
MalthusP˙=rP\dot P = rPSeparableP=P0ertP = P_0 e^{rt}
LogisticP˙=rP(1P/K)\dot P = rP(1-P/K)Separable + partial fractionsP=K/(1+Aert)P = K/(1+Ae^{-rt})
CoolingT˙=k(TTa)\dot T = -k(T-T_a)SeparableT=Ta+(T0Ta)ektT = T_a + (T_0-T_a)e^{-kt}
2nd homog.ay+by+cy=0ay''+by'+cy=0Char. aλ2+bλ+c=0a\lambda^2+b\lambda+c=03 cases
2nd forceday+by+cy=qay''+by'+cy=qyh+ypy_h + y_pUndetermined coeff.
Numerical EulerAnyyn+1=yn+hfy_{n+1}=y_n+hfSequence

SVG — ODE decision diagram

ODE arrivesOrder?1st order2nd orderSeparable?→ IntegrateLinear?→ Integr. factorConst. coeff?→ Char.Numericalnot analytic

Worked examples

Exercise list

45 exercises · 11 with worked solution (25%)

Application 19Understanding 7Modeling 12Challenge 5Proof 2
  1. Ex. 100.1Application

    General solution of y=3xyy' = 3xy:

    Select the correct option
    Select an option first
    Show solution
    Separable: dy/y=3xdxdy/y=3x\,dx. Integrate: lny=3x2/2+C1\ln|y|=3x^2/2+C_1. y=Ce3x2/2y=Ce^{3x^2/2}.
    Show step-by-step (with the why)
    1. Step 1. Separable: dy/y=3xdxdy/y = 3x\,dx.
    2. Step 2. Integrate both sides: lny=3x2/2+C1\ln|y| = 3x^2/2 + C_1.
    3. Step 3. Exponentiate: y=eC1e3x2/2|y| = e^{C_1}e^{3x^2/2}. Absorb constant: y=Ce3x2/2y = Ce^{3x^2/2}.
    4. Trick: Separable = isolate dydy from dxdx, integrate each side.
  2. Ex. 100.2Application

    General solution of y=xy2y' = -xy^2:

    Select the correct option
    Select an option first
    Show solution
    Separable: dy/y2=xdxdy/y^2=-x\,dx. Integrate: 1/y=x2/2+C1-1/y=-x^2/2+C_1, so 1/y=x2/2C11/y=x^2/2-C_1. Rename C1=C-C_1=C: y=1/(x2/2+C)y=1/(x^2/2+C).
  3. Ex. 100.3Application

    General solution of y=ycosxy' = y\cos x:

    Select the correct option
    Select an option first
    Show solution
    See the reference indicated in source for detailed solution.
    Show step-by-step (with the why)
    1. Step 1. Separable: dy/y=cosxdxdy/y = \cos x\,dx.
    2. Step 2. Integrate: lny=sinx+C1\ln|y| = \sin x + C_1.
    3. Step 3. Exponentiate: y=Cesinxy = Ce^{\sin x}.
    4. Trick: The exponent is the antiderivative of cosx\cos x, which is sinx\sin x.
  4. Ex. 100.4Application

    General solution of y=x/yy' = x/y (separable ODE):

    Select the correct option
    Select an option first
    Show solution
    Separable: ydy=xdxy\,dy=x\,dx. Integrate: y2/2=x2/2+C1y^2/2=x^2/2+C_1. y2=x2+Cy^2=x^2+C.
  5. Ex. 100.5Application

    General solution of y=1+y2y' = 1 + y^2:

    Select the correct option
    Select an option first
    Show solution
    See the reference indicated in source for detailed solution.
    Show step-by-step (with the why)
    1. Step 1. Separable: dy/(1+y2)=dxdy/(1+y^2) = dx.
    2. Step 2. Integrate: arctany=x+C\arctan y = x + C.
    3. Step 3. Solve: y=tan(x+C)y = \tan(x+C).
    4. Trick: dy/(1+y2)=arctany\int dy/(1+y^2) = \arctan y — standard arctangent integral.
  6. Ex. 100.6Application

    Solution of y=y2exy' = y^2 e^x, y(0)=1y(0) = 1:

    Select the correct option
    Select an option first
    Show solution
    Separable: dy/y2=exdxdy/y^2=e^x\,dx. 1/y=ex+C-1/y=e^x+C. IC y(0)=1y(0)=1: 1=1+CC=2-1=1+C\Rightarrow C=-2. y=1/(2ex)y=1/(2-e^x).
    Show step-by-step (with the why)
    1. dy/y2=exdxdy/y^2=e^x\,dx. Integrate: 1/y=ex+C-1/y=e^x+C.
    2. IC: 1=1+CC=2-1=1+C\Rightarrow C=-2.
    3. y=1/(2ex)y=1/(2-e^x).
  7. Ex. 100.7ApplicationAnswer key

    Solution of yy=xyy' = x, y(0)=1y(0) = 1:

    Select the correct option
    Select an option first
    Show solution
    See the reference indicated in source for detailed solution.
    Show step-by-step (with the why)
    1. Step 1. Rewrite: ydy=xdxy\,dy = x\,dx.
    2. Step 2. Integrate: y2/2=x2/2+Cy^2/2 = x^2/2 + C.
    3. Step 3. IC: y(0)=11/2=0+CC=1/2y(0)=1 \Rightarrow 1/2 = 0 + C \Rightarrow C=1/2.
    4. Step 4. y2=x2+1y^2 = x^2 + 1, so y=x2+1y = \sqrt{x^2+1} (+ sign since y(0)=1>0y(0)=1>0).
    5. Trick: The IC selects the positive branch of the implicit hyperbola.
  8. Ex. 100.8ApplicationAnswer key

    Logistic equation P=0.5P(1P/3)P' = 0.5P(1-P/3), P(0)=1P(0)=1. Solution P(t)P(t):

    Select the correct option
    Select an option first
    Show solution
    Logistic with K=3K=3, r=0.5r=0.5, P0=1P_0=1. A=(31)/1=2A=(3-1)/1=2. P=3/(1+2e0.5t)P=3/(1+2e^{-0.5t}).
  9. Ex. 100.9ApplicationAnswer key

    General solution of y+3y=6y' + 3y = 6:

    Select the correct option
    Select an option first
    Show solution
    Factor e3xe^{3x}. (e3xy)=6e3x(e^{3x}y)'=6e^{3x}. y=2+Ce3xy=2+Ce^{-3x}.
    Show step-by-step (with the why)
    1. Step 1. Standard form: y+3y=6y' + 3y = 6. Integrating factor: μ=e3x\mu = e^{3x}.
    2. Step 2. Multiply: (e3xy)=6e3x(e^{3x}y)' = 6e^{3x}.
    3. Step 3. Integrate: e3xy=2e3x+Ce^{3x}y = 2e^{3x} + C.
    4. Step 4. Divide by e3xe^{3x}: y=2+Ce3xy = 2 + Ce^{-3x}.
    5. Trick: Particular solution = right side divided by coefficient of yy: 6/3=26/3=2.
  10. Ex. 100.10Application

    General solution of y+y2y=0y'' + y' - 2y = 0:

    Select the correct option
    Select an option first
    Show solution
    See the reference indicated in source for detailed solution.
    Show step-by-step (with the why)
    1. Step 1. Characteristic equation: λ2+λ2=0\lambda^2 + \lambda - 2 = 0.
    2. Step 2. Factor: (λ+2)(λ1)=0(\lambda+2)(\lambda-1) = 0. Roots: λ1=2,λ2=1\lambda_1 = -2, \lambda_2 = 1.
    3. Step 3. Distinct real roots: y=C1e2x+C2exy = C_1 e^{-2x} + C_2 e^x.
    4. Trick: Substitute y=eλxy = e^{\lambda x} into the ODE to get the characteristic equation directly.
  11. Ex. 100.11ApplicationAnswer key

    Equilibria and stability of y=2yy2y' = 2y - y^2:

    Select the correct option
    Select an option first
    Show solution
    y=2yy2=y(2y)y'=2y-y^2=y(2-y). Equilibria: y=0y=0 (unstable, since f(0)=2>0f'(0)=2>0) and y=2y=2 (stable, since f(2)=2<0f'(2)=-2<0).
    Show step-by-step (with the why)
    1. Step 1. Equilibria: solve 2yy2=y(2y)=02y - y^2 = y(2-y) = 0. Roots: y=0y^* = 0 and y=2y^* = 2.
    2. Step 2. f(y)=2yy2f(y) = 2y - y^2. Derivative: f(y)=22yf'(y) = 2-2y.
    3. Step 3. At y=0y=0: f(0)=2>0f'(0)=2>0unstable. At y=2y=2: f(2)=2<0f'(2)=-2<0stable.
    4. Trick: Sign of f(y)f'(y^*): positive = unstable (repeller); negative = stable (attractor).
  12. Ex. 100.12Application

    y=yty' = y - t, y(0.2)=0.9y(0.2) = 0.9, h=0.2h = 0.2. Euler: y(0.4)y(0.4) \approx: (Resp: y1=0.9+0.2×0.7=1.04y_1 = 0.9 + 0.2 \times 0.7 = 1.04)

    Select the correct option
    Select an option first
    Show solution
    Euler: f(0.2,0.9)=yt=0.90.2=0.7f(0.2,\,0.9)=y-t=0.9-0.2=0.7. y1=0.9+0.2×0.7=1.04y_1=0.9+0.2\times0.7=1.04.
  13. Ex. 100.13Understanding

    An ordinary differential equation (ODE) is:

    Select the correct option
    Select an option first
    Show solution
    An ordinary differential equation (ODE) relates an unknown function y(t)y(t) with its derivatives with respect to the single independent variable tt.
  14. Ex. 100.14Understanding

    The difference between stable and unstable equilibrium is:

    Select the correct option
    Select an option first
    Show solution
    Stable equilibrium (attractor): every solution nearby converges to it. Unstable (repeller): nearby solutions move away. Determined by the sign of f(y)f'(y^*).
  15. Ex. 100.15UnderstandingAnswer key

    A 1st-order ODE is linear when:

    Select the correct option
    Select an option first
    Show solution
    See the reference indicated in source for detailed solution.
  16. Ex. 100.16Understanding

    An ODE of the type y=f(y/x)y' = f(y/x) is called homogeneous and is solved by:

    Select the correct option
    Select an option first
    Show solution
    Homogeneous ODE: y=f(y/x)y'=f(y/x). Substitution v=y/xv=y/x, y=vxy=vx, y=v+xvy'=v+xv'. Results in separable in vv.
  17. Ex. 100.17ModelingAnswer key

    Solution of the logistic equation P=rP(1P/K)P' = rP(1-P/K), P(0)=P0P(0) = P_0:

    Select the correct option
    Select an option first
    Show solution
    Logistic: P=K/(1+Aert)P=K/(1+Ae^{-rt}) with A=(KP0)/P0A=(K-P_0)/P_0. Rewriting: P=KP0/(P0+(KP0)ert)P=KP_0/(P_0+(K-P_0)e^{-rt}).
  18. Ex. 100.18ModelingAnswer key

    Falling velocity with linear drag (mv˙=mgkvm\dot v = mg - kv, v(0)=0v(0)=0):

    Select the correct option
    Select an option first
    Show solution
    See the reference indicated in source for detailed solution.
  19. Ex. 100.19Modeling

    Coffee at 100°C in a room at 20°C. Temperature T(t)T(t):

    Select the correct option
    Select an option first
    Show solution
    Newton: T=Tamb+(T0Tamb)ekt=20+80ektT=T_{\text{amb}}+(T_0-T_{\text{amb}})e^{-kt}=20+80e^{-kt} with T0=100T_0=100, Tamb=20T_{\text{amb}}=20.
  20. Ex. 100.20Modeling

    ODE of the series RLC circuit:

    Select the correct option
    Select an option first
    Show solution
    KVL in RLC circuit: sum of voltages = source. With I=q˙I=\dot q: Lq¨+Rq˙+q/C=EL\ddot q+R\dot q+q/C=E.
  21. Ex. 100.21Modeling

    The Malthus model P=rPP' = rP is appropriate when:

    Select the correct option
    Select an option first
    Show solution
    Malthus: P=rPP'=rP, exponential growth without limit. Valid for PKP\ll K. For large PP: logistic is more realistic (limited resource).
  22. Ex. 100.22Modeling

    y=xyy' = xy, y(0)=1y(0)=1, h=0.5h=0.5: Euler 1 step, y(0.5)y(0.5) \approx: (Resp: y1=1y_1 = 1, exact 1.133\approx 1.133)

    Select the correct option
    Select an option first
    Show solution
    Euler: f(0,1)=01=0f(0,1)=0\cdot1=0. y1=1+0.5×0=1.0y_1=1+0.5\times0=1.0. Exact: y=ex2/2y=e^{x^2/2}, y(0.5)=e0.1251.133y(0.5)=e^{0.125}\approx1.133. Truncation error large because h=0.5h=0.5 is large.
  23. Ex. 100.23Modeling

    Oscillator with m=1m=1 kg, b=4b=4, k=13k=13, x(0)=3x(0)=3, x˙(0)=0\dot x(0)=0. Solution x(t)x(t):

    Select the correct option
    Select an option first
    Show solution
    See the reference indicated in source for detailed solution.
  24. Ex. 100.24UnderstandingAnswer key

    Comparison between the exponential and logistic models for population growth:

    Select the correct option
    Select an option first
    Show solution
    See the reference indicated in source for detailed solution.
  25. Ex. 100.25Understanding

    A 1st-order ODE is separable when it can be written in the form:

    Select the correct option
    Select an option first
    Show solution
    Separable ODE: y/g(y)=f(x)y'/g(y)=f(x) — separate into different sides. Distractor py+qpy+q: linear, not separable in general. Distractor y2+xy^2+x: does not factor.
  26. Ex. 100.26Understanding

    The maximum growth rate of the logistic equation is:

    Select the correct option
    Select an option first
    Show solution
    P=rP(1P/K)P'=rP(1-P/K). Maximum at: P=K/2P=K/2. Maximum value: r(K/2)(1/2)=rK/4r(K/2)(1/2)=rK/4.
  27. Ex. 100.27Modeling

    Maximum sustainable yield (MSY) in the logistic equation P=rP(1P/K)P' = rP(1-P/K):

    Select the correct option
    Select an option first
    Show solution
    Logistic equation with harvest hh: equilibria where rP(1P/K)=hrP(1-P/K)=h. Maximum sustainable harvest: hmax=rK/4h_{\max}=rK/4. Stock balance is integration of the solution.
  28. Ex. 100.28Application

    General solution of y=tanxy' = \tan x:

    Select the correct option
    Select an option first
    Show solution
    See the reference indicated in source for detailed solution.
  29. Ex. 100.29Application

    Solution of y=x(1y)y' = x(1-y), y(0)=0y(0) = 0:

    Select the correct option
    Select an option first
    Show solution
    See the reference indicated in source for detailed solution.
    Show step-by-step (with the why)
    1. Step 1. Separable: dy/(1y)=xdxdy/(1-y) = x\,dx.
    2. Step 2. Integrate: ln1y=x2/2+C-\ln|1-y| = x^2/2 + C.
    3. Step 3. 1y=eCex2/2|1-y| = e^{-C}e^{-x^2/2}. Absorb: 1y=Aex2/21-y = Ae^{-x^2/2}.
    4. Step 4. IC: y(0)=010=AA=1y(0)=0 \Rightarrow 1-0=A \Rightarrow A=1.
    5. Step 5. y=1ex2/2y = 1 - e^{-x^2/2}.
    6. Trick: Substitute u=1yu = 1-y: u=y=xuu' = -y' = -xu — more direct equation.
  30. Ex. 100.30ApplicationAnswer key

    Solution of y=ysinxy' = y\sin x, y(0)=1y(0) = 1:

    Select the correct option
    Select an option first
    Show solution
    See the reference indicated in source for detailed solution.
    Show step-by-step (with the why)
    1. Step 1. Separable: dy/y=sinxdxdy/y = \sin x\,dx.
    2. Step 2. Integrate: lny=cosx+C\ln|y| = -\cos x + C.
    3. Step 3. y=Aecosxy = Ae^{-\cos x}. IC: y(0)=1A=ecos0=e1=ey(0)=1 \Rightarrow A=e^{\cos 0}=e^1=e.
    4. Step 4. y=eecosx=e1cosxy = e \cdot e^{-\cos x} = e^{1-\cos x}.
    5. Trick: The exponent is cosx+1=1cosx-\cos x + 1 = 1 - \cos x — compact and elegant.
  31. Ex. 100.31Modeling

    Octopus: K=1000K=1000, r=0.3r=0.3/year. Logistic ODE and maximum growth rate:

    Select the correct option
    Select an option first
    Show solution
    See the reference indicated in source for detailed solution.
    Show step-by-step (with the why)
    1. Step 1. Logistic ODE: P˙=rP(1P/K)=0.3P(1P/1000)\dot P = rP(1-P/K) = 0.3P(1-P/1000).
    2. Step 2. Growth rate is a function of PP: g(P)=0.3P(1P/1000)g(P) = 0.3P(1-P/1000).
    3. Step 3. Maximum when g(P)=0g'(P)=0: 0.3(12P/1000)=0P=K/2=5000.3(1-2P/1000)=0 \Rightarrow P=K/2=500.
    4. Step 4. Maximum rate: g(500)=0.3×500×0.5=75g(500) = 0.3 \times 500 \times 0.5 = 75 octopuses/year.
    5. Trick: In logistics, maximum rate always occurs at P=K/2P = K/2.
  32. Ex. 100.32Modeling

    Octopus: K=1000K=1000, r=0.3r=0.3, P(0)=100P(0)=100. P(10)P(10) \approx: (Resp: 691\approx 691 octopuses)

    Select the correct option
    Select an option first
    Show solution
    Logistic: A=(KP0)/P0=(1000100)/100=9A=(K-P_0)/P_0=(1000-100)/100=9. P(10)=1000/(1+9e0.3×10)=1000/(1+9e3)691P(10)=1000/(1+9e^{-0.3\times10})=1000/(1+9e^{-3})\approx691 octopuses.
  33. Ex. 100.33Challenge

    For linear ODE y+p(t)y=0y' + p(t)y = 0, the Wronskian of two solutions satisfies:

    Select the correct option
    Select an option first
    Show solution
    See the reference indicated in source for detailed solution.
  34. Ex. 100.34Challenge

    The solution of yy=1/xy' - y = 1/x involves:

    Select the correct option
    Select an option first
    Show solution
    See the reference indicated in source for detailed solution.
  35. Ex. 100.35Application

    General solution of y=(y21)/yy' = (y^2-1)/y:

    Select the correct option
    Select an option first
    Show solution
    Separate: y/(y21)dy=dxy/(y^2-1)\,dy=dx. Integrate: 12lny21=x+C1\tfrac12\ln|y^2-1|=x+C_1, so lny21=2x+C\ln|y^2-1|=2x+C, i.e., y21=Ce2xy^2-1=Ce^{2x}.
  36. Ex. 100.36ApplicationAnswer key

    General solution of xy=1/yxy' = 1/y (or equivalently yy=1/xyy' = 1/x):

    Select the correct option
    Select an option first
    Show solution
    See the reference indicated in source for detailed solution.
  37. Ex. 100.37Challenge

    For xy=y+xxy' = y + x (homogeneous ODE), the solution is:

    Select the correct option
    Select an option first
    Show solution
    Substitute v=y/xv=y/x, y=vxy=vx, y=v+xvy'=v+xv'. Then xv+v=v+1xv'+v=v+1, so xv=1xv'=1, v=lnx+Cv=\ln x+C. Back: y=x(lnx+C)=Cx+xlnxy=x(\ln x+C)=Cx+x\ln x.
  38. Ex. 100.38Modeling

    Oscillator with m=1m=1 kg and k=10k=10 N/m. Critical damping bcb_c:

    Select the correct option
    Select an option first
    Show solution
    See the reference indicated in source for detailed solution.
  39. Ex. 100.39Challenge

    Bacteria: per capita 3 when P=6P=6 and per capita 2 when P=9P=9. Logistic ODE and carrying capacity:

    Select the correct option
    Select an option first
    Show solution
    Per capita rate is linear: P/P=abPP'/P=a-bP. System: a6b=3a-6b=3 and a9b=2a-9b=2. Subtract: 3b=1b=1/33b=1\Rightarrow b=1/3, a=5a=5. Thus P=P(5P/3)P'=P(5-P/3); carrying capacity K=a/b=15K=a/b=15.
  40. Ex. 100.40Proof

    The mathematical justification for the separation of variables method is:

    Select the correct option
    Select an option first
    Show solution
    See the reference indicated in source for detailed solution.
  41. Ex. 100.41Application

    Doubling time of P=rPP' = rP (r>0r > 0):

    Select the correct option
    Select an option first
    Show solution
    See the reference indicated in source for detailed solution.
  42. Ex. 100.42ApplicationAnswer key

    Empty tank of volume VV, brine cinc_{\text{in}} enters at rate QVQ_V. Amount of salt Q(t)Q(t):

    Select the correct option
    Select an option first
    Show solution
    ODE: Q˙+(QV/V)Q=QVcin\dot Q+(Q_V/V)Q=Q_Vc_{\text{in}}. With Q(0)=0Q(0)=0: Q=cinV(1eQVt/V)Q=c_{\text{in}}V(1-e^{-Q_Vt/V}). Saturates at cinVc_{\text{in}}V.
  43. Ex. 100.43Modeling

    Present value of USD 50 000 to be received in 20 years, continuous rate r=5%r=5\%: (Resp: =50000e1= 50000e^{-1} \approx USD 18 394)

    Select the correct option
    Select an option first
    Show solution
    Continuous present value: PV=50000e0.05×20=50000e118394PV=50000\,e^{-0.05\times20}=50000\,e^{-1}\approx18394 USD.
  44. Ex. 100.44Challenge

    Species with minimum survival threshold TT (Allee effect): ODE and behavior:

    Select the correct option
    Select an option first
    Show solution
    Allee effect (survival threshold TT): for P<TP<T population collapses; for P>TP>T grows to KK. Model: P=rP(1P/K)(1T/P)P'=rP(1-P/K)(1-T/P).
  45. Ex. 100.45Proof

    Prove that P=KP=K is globally stable equilibrium of the logistic equation for P0>0P_0 > 0, r>0r > 0:

    Select the correct option
    Select an option first
    Show solution
    For r>0r>0, P0>0P_0>0: P=rP(1P/K)>0P'=rP(1-P/K)>0 for P<KP<K and P<0P'<0 for P>KP>K. By first-order stability theorem, P=KP=K is globally stable for P0>0P_0>0.

Sources

Updated on 2026-05-06 · Author(s): Clube da Matemática

Found an error? Open an issue on GitHub or submit a PR — open source forever.