CP-9 Differential equations — coverage pack

9 specification leaves · notes, questions, answers and worked methods

CP-9.1 · Find and use an integrating factor to solve differential equations of form dy/dx + P(x)y = Q(x) and recognise when it is appropriate to do so.

  • Put a first-order linear equation into the form y+P(x)y=Q(x)y'+P(x)y=Q(x) before choosing the integrating factor I(x)=eP(x)dxI(x)=e^{\int P(x)\,dx}.
  • After multiplication, the left side is the product derivative (Iy)(Iy)', so integrate Iy=IQdx+CIy=\int IQ\,dx+C and then divide by II.
  • Use an integrating factor when yy and yy' occur linearly with coefficients depending only on the independent variable.
  • A common error is to use the coefficient of yy' in the exponent; first divide through so that the coefficient of yy' is 11.

Tier 1 · Easy

  1. 1. Find an integrating factor for dydx+3y=x\dfrac{dy}{dx}+3y=x.[2 marks]

    Answer

    • e3xe^{3x}

    Method: Here P(x)=3P(x)=3, so I(x)=e3dx=e3xI(x)=e^{\int3\,dx}=e^{3x}. A non-zero constant multiple would be equivalent.

Tier 2 · Standard

  1. 1. For x>0x>0, solve dydx+2xy=x2\dfrac{dy}{dx}+\dfrac2x y=x^2.[5 marks]

    Answer

    • y=x35+Cx2y=\dfrac{x^3}{5}+\dfrac{C}{x^2}

    Method: The integrating factor is e2/xdx=x2e^{\int2/x\,dx}=x^2 because x>0x>0. Multiplication gives (x2y)=x4(x^2y)'=x^4. Hence x2y=x5/5+Cx^2y=x^5/5+C, and division by x2x^2 gives y=x3/5+C/x2y=x^3/5+C/x^2.

Tier 3 · Hard

  1. 1. On 0x<π/20\leq x<\pi/2, solve dydx+2tanxy=sinx\dfrac{dy}{dx}+2\tan x\,y=\sin x subject to y(0)=2y(0)=2. Hence find y(π/3)y(\pi/3).[7 marks]

    Answer

    • y=cosx+cos2xy=\cos x+\cos^2x
    • y(π/3)=34y(\pi/3)=\dfrac34

    Method: The integrating factor is e2tanxdx=e2ln(cosx)=sec2xe^{\int2\tan x\,dx}=e^{-2\ln(\cos x)}=\sec^2x. Thus (ysec2x)=sinxsec2x(y\sec^2x)'=\sin x\sec^2x. Integration gives ysec2x=secx+Cy\sec^2x=\sec x+C, so y=cosx+Ccos2xy=\cos x+C\cos^2x. The initial condition gives C=1C=1. At x=π/3x=\pi/3, y=1/2+1/4=3/4y=1/2+1/4=3/4.

CP-9.2 · Find both general and particular solutions to differential equations.

  • A general solution contains the full number of arbitrary constants: one for a first-order equation and two for a second-order equation. Varying the constant generates a family of solution curves, and you are expected to sketch representative members of that family (e.g. y=Cex+x1y=Ce^{-x}+x-1 for several values of CC).
  • A particular solution is obtained when initial or boundary conditions determine all arbitrary constants.
  • Apply conditions only after differentiating the complete general solution as many times as necessary.
  • Do not call a particular integral a particular solution: a particular integral handles the forcing term, while a particular solution also uses the supplied conditions.

Tier 1 · Easy

  1. 1. Find the general solution of dydx=4x3\dfrac{dy}{dx}=4x^3, then find the particular solution for which y(1)=3y(1)=3.[3 marks]

    Answer

    • General: y=x4+Cy=x^4+C
    • Particular: y=x4+2y=x^4+2

    Method: Integrating gives y=x4+Cy=x^4+C. Substituting x=1x=1, y=3y=3 gives 3=1+C3=1+C, so C=2C=2.

Tier 2 · Standard

  1. 1. Solve yy=0y''-y'=0 subject to y(0)=2y(0)=2 and y(0)=3y'(0)=-3.[5 marks]

    Answer

    • y=53exy=5-3e^x

    Method: The auxiliary equation is m2m=0m^2-m=0, with roots 00 and 11. Hence y=A+Bexy=A+Be^x and y=Bexy'=Be^x. From y(0)=3y'(0)=-3, B=3B=-3; then y(0)=A+B=2y(0)=A+B=2 gives A=5A=5.

Tier 3 · Hard

  1. 1. Find the particular solution of y+4y=12xy''+4y=12x for which y(0)=1y(0)=1 and y(0)=0y'(0)=0.[7 marks]

    Answer

    • y=cos(2x)32sin(2x)+3xy=\cos(2x)-\dfrac32\sin(2x)+3x

    Method: The complementary function is Acos2x+Bsin2xA\cos2x+B\sin2x. A linear particular integral px+qpx+q gives 4(px+q)=12x4(px+q)=12x, so p=3p=3 and q=0q=0. Thus y=Acos2x+Bsin2x+3xy=A\cos2x+B\sin2x+3x. The first condition gives A=1A=1. Since y=2Asin2x+2Bcos2x+3y'=-2A\sin2x+2B\cos2x+3, the second gives 2B+3=02B+3=0, so B=3/2B=-3/2.

CP-9.3 · Use differential equations in modelling in kinematics and in other contexts.

  • Define the dependent variable, independent variable, units and positive direction before translating rates into a differential equation.
  • In kinematics, a=dv/dta=dv/dt and v=ds/dtv=ds/dt; resistance acts opposite to motion, so its sign depends on the chosen direction.
  • After solving, interpret equilibrium values, limiting behaviour and constants in the original context rather than leaving a bare formula.
  • A mathematically correct solution can still expose a weak model; check assumptions such as constant coefficients, perfect mixing and whether predicted quantities remain physical.

Tier 1 · Easy

  1. 1. A body at temperature TT is in a room maintained at 18C18^{\circ}\text{C}. State a differential equation expressing that its cooling rate is proportional to its excess temperature, where k>0k>0.[2 marks]

    Answer

    • dTdt=k(T18)\dfrac{dT}{dt}=-k(T-18)

    Method: The excess temperature is T18T-18. Cooling means TT decreases when this excess is positive, so the proportionality constant must appear with a minus sign.

Tier 2 · Standard

  1. 1. A particle falls from rest. Taking downward as positive, its speed satisfies dvdt=100.5v\dfrac{dv}{dt}=10-0.5v. Find v(t)v(t) and the terminal speed predicted by the model.[5 marks]

    Answer

    • v=20(1et/2)v=20(1-e^{-t/2})
    • Terminal speed 2020

    Method: Rewrite as v+0.5v=10v'+0.5v=10. The general solution is v=20+Cet/2v=20+Ce^{-t/2}. Since v(0)=0v(0)=0, C=20C=-20, giving v=20(1et/2)v=20(1-e^{-t/2}). As tt\to\infty, the exponential tends to 00, so the limiting speed is 2020 in the stated units.

Tier 3 · Hard

  1. 1. A well-mixed bioreactor contains 120120 litres of liquid and initially 1212 grams of dissolved nutrient. Solution enters and leaves at 33 litres per minute, keeping the volume constant. The incoming concentration is 0.50.5 grams per litre. Form and solve a differential equation for the nutrient mass MM grams, find when M=36M=36, and state one modelling assumption.[9 marks]

    Answer

    • dMdt=1.5M40\dfrac{dM}{dt}=1.5-\dfrac{M}{40}
    • M=6048et/40M=60-48e^{-t/40}
    • t=40ln2t=40\ln2 minutes
    • For example, the nutrient is mixed uniformly at every instant.

    Method: Nutrient enters at 3(0.5)=1.53(0.5)=1.5 grams per minute. Its concentration in the tank is M/120M/120, so it leaves at 3M/120=M/403M/120=M/40. Thus M+M/40=1.5M'+M/40=1.5, whose solution is M=60+Cet/40M=60+Ce^{-t/40}. The initial mass gives C=48C=-48. Setting M=36M=36 gives 24=48et/4024=48e^{-t/40}, so et/40=1/2e^{-t/40}=1/2 and t=40ln2t=40\ln2. The derivation assumes instantaneous perfect mixing, as well as constant rates and volume.

CP-9.4 · Solve differential equations of form y'' + ay' + by = 0 where a and b are constants by using the auxiliary equation.

  • For y+ay+by=0y''+ay'+by=0, try y=emxy=e^{mx} to obtain the auxiliary equation m2+am+b=0m^2+am+b=0.
  • Convert the roots into the matching real general solution, retaining two independent arbitrary constants.
  • Initial conditions determine the constants only after the full complementary solution and its derivative have been written.
  • Do not leave complex exponentials in a requested real solution; combine conjugate roots into sine and cosine terms.

Tier 1 · Easy

  1. 1. Find the general solution of y+5y+6y=0y''+5y'+6y=0.[3 marks]

    Answer

    • y=Ae2x+Be3xy=Ae^{-2x}+Be^{-3x}

    Method: The auxiliary equation m2+5m+6=0m^2+5m+6=0 factorises as (m+2)(m+3)=0(m+2)(m+3)=0. The distinct roots are 2-2 and 3-3, giving the stated two-exponential solution.

Tier 2 · Standard

  1. 1. Find the general solution of y6y+9y=0y''-6y'+9y=0.[3 marks]

    Answer

    • y=(A+Bx)e3xy=(A+Bx)e^{3x}

    Method: The auxiliary equation is (m3)2=0(m-3)^2=0, so m=3m=3 is a repeated root. The two independent terms are e3xe^{3x} and xe3xxe^{3x}.

Tier 3 · Hard

  1. 1. Solve 2y+4y+10y=02y''+4y'+10y=0 subject to y(0)=2y(0)=2 and y(0)=2y'(0)=2.[6 marks]

    Answer

    • y=2ex(cos2x+sin2x)y=2e^{-x}(\cos2x+\sin2x)

    Method: After division by 22, the auxiliary equation is m2+2m+5=0m^2+2m+5=0, with roots 1±2i-1\pm2i. Thus y=ex(Acos2x+Bsin2x)y=e^{-x}(A\cos2x+B\sin2x). The first condition gives A=2A=2. Differentiating and setting x=0x=0 gives y(0)=A+2B=2y'(0)=-A+2B=2, so B=2B=2.

CP-9.5 · Solve differential equations of form y'' + ay' + by = f(x) by solving the homogeneous case and adding a particular integral to the complementary function (f(x) polynomial, exponential or trigonometric).

  • Write the general solution as complementary function plus particular integral: y=yc+ypy=y_{\mathrm c}+y_{\mathrm p}.
  • Choose a trial particular integral of the same family as f(x)f(x), with enough undetermined coefficients for every term produced by differentiation.
  • If the proposed trial duplicates a term in the complementary function, multiply it by xx repeatedly until it is independent.
  • A common error is to apply initial conditions to the complementary function alone; combine CF and PI first.

Tier 1 · Easy

  1. 1. Find the general solution of yy2y=6y''-y'-2y=6.[4 marks]

    Answer

    • y=Ae2x+Bex3y=Ae^{2x}+Be^{-x}-3

    Method: The auxiliary equation (m2)(m+1)=0(m-2)(m+1)=0 gives yc=Ae2x+Bexy_{\mathrm c}=Ae^{2x}+Be^{-x}. For a constant PI KK, substitution gives 2K=6-2K=6, so K=3K=-3.

Tier 2 · Standard

  1. 1. Find the general solution of y4y=8e2xy''-4y=8e^{2x}.[5 marks]

    Answer

    • y=Ae2x+Be2x+2xe2xy=Ae^{2x}+Be^{-2x}+2xe^{2x}

    Method: The complementary function is Ae2x+Be2xAe^{2x}+Be^{-2x}. Since e2xe^{2x} is already present, try yp=Cxe2xy_{\mathrm p}=Cxe^{2x}. Then yp4yp=4Ce2xy_{\mathrm p}''-4y_{\mathrm p}=4Ce^{2x}, so C=2C=2.

Tier 3 · Hard

  1. 1. Solve y+2y+5y=10cosxy''+2y'+5y=10\cos x subject to y(0)=1y(0)=1 and y(0)=0y'(0)=0.[9 marks]

    Answer

    • y=ex(cos2x+sin2x)+2cosx+sinxy=-e^{-x}(\cos2x+\sin2x)+2\cos x+\sin x

    Method: The auxiliary roots are 1±2i-1\pm2i, so yc=ex(Acos2x+Bsin2x)y_{\mathrm c}=e^{-x}(A\cos2x+B\sin2x). Try yp=acosx+bsinxy_{\mathrm p}=a\cos x+b\sin x. Substitution gives 4a+2b=104a+2b=10 and 2a+4b=0-2a+4b=0, so a=2a=2, b=1b=1. From y(0)=1y(0)=1, A+2=1A+2=1 and A=1A=-1. The derivative at 00 is A+2B+1-A+2B+1; setting this to 00 gives B=1B=-1. Combining the parts yields the stated solution.

CP-9.6 · Understand and use the relationship between the cases when the discriminant of the auxiliary equation is positive, zero and negative and the form of solution of the differential equation.

  • For m2+am+b=0m^2+am+b=0, a positive discriminant gives distinct real roots and two exponential terms.
  • A zero discriminant gives a repeated real root rr and the solution (A+Bx)erx(A+Bx)e^{rx}.
  • A negative discriminant gives roots p±iqp\pm iq and the real solution epx(Acosqx+Bsinqx)e^{px}(A\cos qx+B\sin qx).
  • Do not infer oscillation from a negative real root alone: oscillation comes from a non-zero imaginary part.

Tier 1 · Easy

  1. 1. For y+4y+8y=0y''+4y'+8y=0, state the sign of the auxiliary discriminant and hence write the general solution.[3 marks]

    Answer

    • The discriminant is negative; y=e2x(Acos2x+Bsin2x)y=e^{-2x}(A\cos2x+B\sin2x).

    Method: The discriminant is 424(8)=16<04^2-4(8)=-16<0. The roots are 2±2i-2\pm2i, which give the damped sine-cosine form.

Tier 2 · Standard

  1. 1. Given k>0k>0, find the value of kk for which y+ky+16y=0y''+ky'+16y=0 has a repeated auxiliary root. Write the corresponding general solution.[4 marks]

    Answer

    • k=8k=8; y=(A+Bx)e4xy=(A+Bx)e^{-4x}

    Method: A repeated root requires discriminant k264=0k^2-64=0. Since k>0k>0, k=8k=8. The repeated root is k/2=4-k/2=-4, so the solution is (A+Bx)e4x(A+Bx)e^{-4x}.

Tier 3 · Hard

  1. 1. Classify the auxiliary roots of y+2(p+1)y+(p2+4p+8)y=0y''+2(p+1)y'+(p^2+4p+8)y=0 for all real pp. Also give the general solution at the transition value.[6 marks]

    Answer

    • Distinct real roots for p<7/2p<-7/2, a repeated root for p=7/2p=-7/2, and complex conjugate roots for p>7/2p>-7/2.
    • At p=7/2p=-7/2, y=(A+Bx)e5x/2y=(A+Bx)e^{5x/2}.

    Method: The discriminant is [2(p+1)]24(p2+4p+8)=4(2p+7)[2(p+1)]^2-4(p^2+4p+8)=-4(2p+7). It is positive, zero or negative according as p<7/2p<-7/2, p=7/2p=-7/2 or p>7/2p>-7/2. At equality the repeated root is [2(p+1)]/2=(p+1)=5/2-[2(p+1)]/2=-(p+1)=5/2, giving (A+Bx)e5x/2(A+Bx)e^{5x/2}.

CP-9.7 · Solve the equation for simple harmonic motion x'' = -omega^2 x and relate the solution to the motion.

  • Simple harmonic motion about x=0x=0 satisfies x=ω2xx''=-\omega^2x, so acceleration is proportional to displacement and directed towards equilibrium.
  • The general solution is x=Acosωt+Bsinωtx=A\cos\omega t+B\sin\omega t, equivalently x=Rcos(ωt+ϕ)x=R\cos(\omega t+\phi).
  • The amplitude is R=A2+B2R=\sqrt{A^2+B^2}, the period is 2π/ω2\pi/\omega, and the maximum speed is ωR\omega R.
  • Keep displacement, velocity and acceleration signs distinct; maximum acceleration magnitude occurs at the extreme positions, not at equilibrium.

Tier 1 · Easy

  1. 1. A particle satisfies x=9xx''=-9x. Write its general displacement and state its period.[3 marks]

    Answer

    • x=Acos3t+Bsin3tx=A\cos3t+B\sin3t; period 2π/32\pi/3.

    Method: Comparison with x=ω2xx''=-\omega^2x gives ω=3\omega=3. Substitute this angular frequency into the general SHM solution and use T=2π/ωT=2\pi/\omega.

Tier 2 · Standard

  1. 1. An SHM particle satisfies x=4xx''=-4x, with x(0)=4x(0)=4 and x(0)=6x'(0)=6. Find its displacement, amplitude, period and maximum speed.[6 marks]

    Answer

    • x=4cos2t+3sin2tx=4\cos2t+3\sin2t
    • Amplitude 55, period π\pi, maximum speed 1010.

    Method: Here ω=2\omega=2, so x=Acos2t+Bsin2tx=A\cos2t+B\sin2t. The conditions give A=4A=4 and 2B=62B=6, hence B=3B=3. The amplitude is 42+32=5\sqrt{4^2+3^2}=5, the period is 2π/2=π2\pi/2=\pi, and the maximum speed is ωR=2(5)=10\omega R=2(5)=10.

Tier 3 · Hard

  1. 1. A particle in SHM satisfies x=16xx''=-16x, x(0)=3/2x(0)=\sqrt3/2 and x(0)=2x'(0)=-2. Express xx in the form Rcos(4t+ϕ)R\cos(4t+\phi), find the first positive time at which it crosses equilibrium moving in the negative direction, and find its acceleration when x=1/2x=-1/2.[8 marks]

    Answer

    • x=cos(4t+π/6)x=\cos(4t+\pi/6)
    • t=π/12t=\pi/12
    • x=8x''=8 when x=1/2x=-1/2.

    Method: In x=Acos4t+Bsin4tx=A\cos4t+B\sin4t, the conditions give A=3/2A=\sqrt3/2 and 4B=24B=-2, so B=1/2B=-1/2. Thus R=1R=1 and x=cos(4t+π/6)x=\cos(4t+\pi/6). The first negative-direction equilibrium crossing has phase π/2\pi/2, so 4t+π/6=π/24t+\pi/6=\pi/2 and t=π/12t=\pi/12. Finally x=16x=16(1/2)=8x''=-16x=-16(-1/2)=8.

CP-9.8 · Model damped oscillations using second order differential equations and interpret their solutions.

  • A free damped oscillator is modelled by x+2λx+ω02x=0x''+2\lambda x'+\omega_0^2x=0 with λ>0\lambda>0; the velocity term removes energy.
  • Complex auxiliary roots give underdamped oscillation with a decaying exponential envelope; a repeated root gives critical damping.
  • Two distinct negative real roots give overdamping: the system returns without oscillating, more slowly than at critical damping.
  • A periodic driving force appears on the right side and produces a persistent steady response; do not mistake the decaying complementary function for the complete motion.

Tier 1 · Easy

  1. 1. Classify the motion governed by x+6x+25x=0x''+6x'+25x=0 and write its general solution.[4 marks]

    Answer

    • Underdamped; x=e3t(Acos4t+Bsin4t)x=e^{-3t}(A\cos4t+B\sin4t).

    Method: The auxiliary equation has roots 3±4i-3\pm4i. The non-zero imaginary part produces oscillation, while the negative real part gives the decaying envelope e3te^{-3t}.

Tier 2 · Standard

  1. 1. Find the positive value of cc for which x+cx+9x=0x''+cx'+9x=0 is critically damped. For this value, solve the equation when x(0)=2x(0)=2 and x(0)=4x'(0)=-4, and interpret the long-term motion.[7 marks]

    Answer

    • c=6c=6
    • x=2(1+t)e3tx=2(1+t)e^{-3t}
    • The displacement tends to 00 without oscillation.

    Method: Critical damping requires c236=0c^2-36=0, so positive cc gives c=6c=6. The repeated root is 3-3, hence x=(A+Bt)e3tx=(A+Bt)e^{-3t}. From x(0)=2x(0)=2, A=2A=2; from x(0)=B3A=4x'(0)=B-3A=-4, B=2B=2. Both the exponential and polynomial-exponential term tend to 00, and the repeated real root means no oscillation.

Tier 3 · Hard

  1. 1. A forced damped oscillator satisfies x+4x+20x=10cos(2t)x''+4x'+20x=10\cos(2t), with x(0)=0x(0)=0 and x(0)=0x'(0)=0. Find x(t)x(t) and describe the long-term motion.[10 marks]

    Answer

    • x=e2t(12cos4t38sin4t)+12cos2t+14sin2tx=e^{-2t}(-\dfrac12\cos4t-\dfrac38\sin4t)+\dfrac12\cos2t+\dfrac14\sin2t
    • The transient term decays, leaving a steady oscillation of angular frequency 22 and amplitude 5/4\sqrt5/4.

    Method: The auxiliary roots are 2±4i-2\pm4i, so xc=e2t(Acos4t+Bsin4t)x_{\mathrm c}=e^{-2t}(A\cos4t+B\sin4t). Try xp=acos2t+bsin2tx_{\mathrm p}=a\cos2t+b\sin2t. Substitution gives 16a+8b=1016a+8b=10 and 8a+16b=0-8a+16b=0, hence a=1/2a=1/2 and b=1/4b=1/4. From x(0)=0x(0)=0, A=1/2A=-1/2. At t=0t=0, x=2A+4B+2b=0x'=-2A+4B+2b=0, so B=3/8B=-3/8. The complementary term tends to 00, while the PI remains with amplitude (1/2)2+(1/4)2=5/4\sqrt{(1/2)^2+(1/4)^2}=\sqrt5/4 and angular frequency 22.

CP-9.9 · Analyse and interpret models with one independent variable and two dependent variables as a pair of coupled first order simultaneous equations and solve them, e.g. predator-prey models.

  • For a coupled first-order system, eliminate one dependent variable to obtain a second-order equation for the other, then recover the eliminated variable.
  • Translate each interaction term carefully: a transfer leaving one population or compartment may enter the other with the opposite sign.
  • Use both initial conditions, including any derivative value inferred from an original equation, and verify the final pair in both equations.
  • Interpret signs, equilibrium points and limiting behaviour in context; an algebraic solution that predicts negative populations may invalidate the model after that time.

Tier 1 · Easy

  1. 1. Given u=2u+vu'=2u+v and v=3uv'=3u, eliminate vv to obtain a second-order differential equation for uu.[3 marks]

    Answer

    • u2u3u=0u''-2u'-3u=0

    Method: Differentiate u=2u+vu'=2u+v to get u=2u+vu''=2u'+v'. Substitute v=3uv'=3u, then rearrange to u2u3u=0u''-2u'-3u=0.

Tier 2 · Standard

  1. 1. Solve x=x+2yx'=x+2y, y=2x+yy'=2x+y subject to x(0)=3x(0)=3 and y(0)=1y(0)=1.[7 marks]

    Answer

    • x=2e3t+etx=2e^{3t}+e^{-t}
    • y=2e3tety=2e^{3t}-e^{-t}

    Method: From the first equation, y=12(xx)y=\tfrac12(x'-x), so y=12(xx)y'=\tfrac12(x''-x'). Substituting into y=2x+yy'=2x+y gives 12(xx)=2x+12(xx)\tfrac12(x''-x')=2x+\tfrac12(x'-x), i.e. x2x3x=0x''-2x'-3x=0. The auxiliary equation m22m3=0m^2-2m-3=0 has roots m=3m=3 and m=1m=-1, so x=Ae3t+Betx=Ae^{3t}+Be^{-t}. Then y=12(xx)=Ae3tBety=\tfrac12(x'-x)=Ae^{3t}-Be^{-t}. The conditions x(0)=3x(0)=3, y(0)=1y(0)=1 give A+B=3A+B=3 and AB=1A-B=1, so A=2A=2, B=1B=1: x=2e3t+etx=2e^{3t}+e^{-t} and y=2e3tety=2e^{3t}-e^{-t}.

Tier 3 · Hard

  1. 1. Amounts A(t)A(t) and B(t)B(t) in two connected compartments are modelled by A=2A+BA'=-2A+B and B=2A3BB'=2A-3B, with A(0)=4A(0)=4 and B(0)=1B(0)=1. Find A(t)A(t) and B(t)B(t), determine exactly when BB is greatest, and state the long-term prediction.[10 marks]

    Answer

    • A=3et+e4tA=3e^{-t}+e^{-4t}
    • B=3et2e4tB=3e^{-t}-2e^{-4t}
    • BB is greatest at t=13ln(8/3)t=\dfrac13\ln(8/3), with maximum 94(3/8)1/3\dfrac94(3/8)^{1/3}.
    • Both amounts tend to 00 as tt\to\infty.

    Method: From A=2A+BA'=-2A+B, B=A+2AB=A'+2A. Differentiate and use the second equation to obtain A+5A+4A=0A''+5A'+4A=0. Thus A=Cet+De4tA=Ce^{-t}+De^{-4t}. Now A(0)=4A(0)=4 and A(0)=2(4)+1=7A'(0)=-2(4)+1=-7, giving C=3C=3, D=1D=1. Hence B=A+2A=3et2e4tB=A'+2A=3e^{-t}-2e^{-4t}. For a maximum, B=3et+8e4t=0B'=-3e^{-t}+8e^{-4t}=0, so e3t=8/3e^{3t}=8/3 and t=(1/3)ln(8/3)t=(1/3)\ln(8/3). Writing r=et=(3/8)1/3r=e^{-t}=(3/8)^{1/3} gives e4t=(3/8)re^{-4t}=(3/8)r, hence Bmax=3r2(3/8)r=(9/4)rB_{\max}=3r-2(3/8)r=(9/4)r. Both exponential terms vanish in the long term.