Home

Tutoring

Subjects

Live Classes

Study Coach

Essay Review

On-Demand Courses

Colleges

Games


Sign up

Log in

Opening subject page...

Loading your content

Practice

  • All Subjects
  • Algebra Flashcards
  • SAT Math Practice Tests
  • Math Question of the Day
  • Live Classes
  • On-Demand Courses

Varsity Tutors

  • Find a Tutor
  • Test Prep
  • Online Classes
  • K-12 Learning
  • College Search
  • VarsityTutors.com

© 2026 Varsity Tutors. All rights reserved.

← Back to quizzes

Differential Equations Quiz

Differential Equations Quiz: Setting Up Des From Word Problems

Practice Setting Up Des From Word Problems in Differential Equations with focused quiz questions that help you check what you know, review explanations, and build confidence with test-style prompts.

Question 1 / 18

0 of 18 answered

A tank initially contains 100 gallons of pure water. Brine containing 2 pounds of salt per gallon flows into the tank at a rate of 3 gallons per minute, while the well-mixed solution flows out at a rate of 2 gallons per minute. If S(t)S(t)S(t) represents the amount of salt in the tank at time ttt minutes, which differential equation correctly models this situation?

Select an answer to continue

What this quiz covers

This quiz focuses on Setting Up Des From Word Problems, giving you a quick way to practice the rules, question types, and explanations that matter most for Differential Equations.

How to use this quiz

Try each quiz question before looking at the correct answer. Use the explanations to review missed ideas, then come back to similar questions until the pattern feels familiar.

All questions

Question 1

A tank initially contains 100 gallons of pure water. Brine containing 2 pounds of salt per gallon flows into the tank at a rate of 3 gallons per minute, while the well-mixed solution flows out at a rate of 2 gallons per minute. If S(t)S(t)S(t) represents the amount of salt in the tank at time ttt minutes, which differential equation correctly models this situation?

  1. dSdt=6−2S100+t\frac{dS}{dt} = 6 - \frac{2S}{100 + t}dtdS​=6−100+t2S​ (correct answer)
  2. dSdt=6−2S100+2t\frac{dS}{dt} = 6 - \frac{2S}{100 + 2t}dtdS​=6−100+2t2S​
  3. dSdt=6−3S100+t\frac{dS}{dt} = 6 - \frac{3S}{100 + t}dtdS​=6−100+t3S​
  4. dSdt=3−2S100+t\frac{dS}{dt} = 3 - \frac{2S}{100 + t}dtdS​=3−100+t2S​

Explanation: The rate of change of salt equals salt in minus salt out. Salt flows in at (2 lb/gal)(3 gal/min) = 6 lb/min. The tank volume at time t is 100 + (3-2)t = 100 + t gallons. Salt flows out at rate (S/(100+t) lb/gal)(2 gal/min) = 2S/(100+t) lb/min. Thus dS/dt = 6 - 2S/(100+t). Choice B incorrectly uses 100+2t for volume (using outflow rate instead of net rate). Choice C uses wrong outflow rate coefficient. Choice D uses inflow rate instead of salt inflow rate.

Question 2

According to Newton's law of cooling, the rate of change of an object's temperature is proportional to the difference between its temperature and the ambient temperature. An object is placed in an environment where the ambient temperature, M(t)M(t)M(t), is not constant but increases linearly from an initial temperature of M0M_0M0​ at a rate of α\alphaα degrees per hour. Let T(t)T(t)T(t) be the object's temperature. The object cools when it is warmer than its surroundings. Let k>0k > 0k>0 be the constant of proportionality. Which differential equation models the temperature T(t)T(t)T(t)?

  1. dTdt=−k(T−(M0+αt))\frac{dT}{dt} = -k(T - (M_0 + \alpha t))dtdT​=−k(T−(M0​+αt)) (correct answer)
  2. dTdt=−k(T−M0)+α\frac{dT}{dt} = -k(T - M_0) + \alphadtdT​=−k(T−M0​)+α
  3. dTdt=−k(T−M0)\frac{dT}{dt} = -k(T - M_0)dtdT​=−k(T−M0​)
  4. dTdt=−k(T−(M0+α))\frac{dT}{dt} = -k(T - (M_0 + \alpha))dtdT​=−k(T−(M0​+α))

Explanation: Newton's law of cooling states that dT/dt=−k(T−M)dT/dt = -k(T - M)dT/dt=−k(T−M), where MMM is the ambient temperature and k>0k > 0k>0. The problem states that the ambient temperature M(t)M(t)M(t) increases linearly with time from an initial value of M0M_0M0​ at a rate of α\alphaα. This means M(t)=M0+αtM(t) = M_0 + \alpha tM(t)=M0​+αt. Substituting this into the cooling law gives dT/dt=−k(T(t)−M(t))dT/dt = -k(T(t) - M(t))dT/dt=−k(T(t)−M(t)), which becomes dT/dt=−k(T−(M0+αt))dT/dt = -k(T - (M_0 + \alpha t))dT/dt=−k(T−(M0​+αt)).

Question 3

The population of a species, P(t)P(t)P(t), is modeled by the logistic equation with a growth rate rrr and carrying capacity KKK. The species is subject to harvesting. The harvesting rate is not constant; it is proportional to the current population, with a proportionality constant hhh that oscillates seasonally according to h(t)=h0(1+cos⁡(ωt))h(t) = h_0(1 + \cos(\omega t))h(t)=h0​(1+cos(ωt)). Which differential equation describes the population P(t)P(t)P(t)?

  1. dPdt=rP(1−PK)−h0(1+cos⁡(ωt))\frac{dP}{dt} = rP\left(1 - \frac{P}{K}\right) - h_0(1+\cos(\omega t))dtdP​=rP(1−KP​)−h0​(1+cos(ωt))
  2. dPdt=rP(1−PK)−h0(1+cos⁡(ωt))P\frac{dP}{dt} = rP\left(1 - \frac{P}{K}\right) - h_0(1+\cos(\omega t))PdtdP​=rP(1−KP​)−h0​(1+cos(ωt))P (correct answer)
  3. dPdt=rP(1−PK−h0(1+cos⁡(ωt)))\frac{dP}{dt} = rP\left(1 - \frac{P}{K - h_0(1+\cos(\omega t))}\right)dtdP​=rP(1−K−h0​(1+cos(ωt))P​)
  4. dPdt=rP(1−PK)+h0(1+cos⁡(ωt))P\frac{dP}{dt} = rP\left(1 - \frac{P}{K}\right) + h_0(1+\cos(\omega t))PdtdP​=rP(1−KP​)+h0​(1+cos(ωt))P

Explanation: The logistic growth model without harvesting is dP/dt=rP(1−P/K)dP/dt = rP(1 - P/K)dP/dt=rP(1−P/K). Harvesting reduces the population, so a negative term must be added. The problem states the harvesting rate is proportional to the current population, with the proportionality constant being h(t)=h0(1+cos⁡(ωt))h(t) = h_0(1 + \cos(\omega t))h(t)=h0​(1+cos(ωt)). Therefore, the harvesting term is h(t)P=h0(1+cos⁡(ωt))Ph(t)P = h_0(1 + \cos(\omega t))Ph(t)P=h0​(1+cos(ωt))P. Combining these gives the full model: dP/dt=rP(1−P/K)−h0(1+cos⁡(ωt))PdP/dt = rP(1 - P/K) - h_0(1 + \cos(\omega t))PdP/dt=rP(1−P/K)−h0​(1+cos(ωt))P.

Question 4

In a population of NNN individuals, the number of people who have a certain piece of information is I(t)I(t)I(t). The information spreads at a rate proportional to the product of those who have the information and those who do not. Additionally, those who have the information can forget it and become uninformed again, at a rate proportional to the square root of the number of informed people. Let kkk be the transmission rate constant and γ\gammaγ be the forgetting rate constant. Which differential equation models I(t)I(t)I(t)?

  1. dIdt=kI(N−I)−γI\frac{dI}{dt} = kI(N-I) - \gamma IdtdI​=kI(N−I)−γI
  2. dIdt=kI(N−I)−γI\frac{dI}{dt} = kI(N-I) - \gamma \sqrt{I}dtdI​=kI(N−I)−γI​ (correct answer)
  3. dIdt=kI(N−I)+γI\frac{dI}{dt} = kI(N-I) + \gamma \sqrt{I}dtdI​=kI(N−I)+γI​
  4. dIdt=kI(N−I)−γI2\frac{dI}{dt} = kI(N-I) - \gamma I^2dtdI​=kI(N−I)−γI2

Explanation: The rate of change, dI/dtdI/dtdI/dt, is the rate of spreading minus the rate of forgetting. The rate of spreading is proportional to the product of the number of informed people, III, and the number of uninformed people, S=N−IS = N-IS=N−I. So, the spreading term is kI(N−I)kI(N-I)kI(N−I). The rate of forgetting is stated to be proportional to the square root of the number of informed people, which is γI\gamma\sqrt{I}γI​. Since forgetting decreases the number of informed people, this term is subtracted. Combining these gives the equation: dI/dt=kI(N−I)−γIdI/dt = kI(N-I) - \gamma\sqrt{I}dI/dt=kI(N−I)−γI​.

Question 5

The value of an investment, V(t)V(t)V(t), grows due to interest at a rate proportional to its current value, with a constant of proportionality rrr. In addition, funds are continuously withdrawn at a rate proportional to the square of the investment's value, with a constant kkk. A constant amount DDD is also deposited continuously into the account per year. Which differential equation models the value of the investment V(t)V(t)V(t)?

  1. dVdt=rV−kV+D\frac{dV}{dt} = rV - k\sqrt{V} + DdtdV​=rV−kV​+D
  2. dVdt=rV+kV2−D\frac{dV}{dt} = rV + kV^2 - DdtdV​=rV+kV2−D
  3. dVdt=r−kV2+D\frac{dV}{dt} = r - kV^2 + DdtdV​=r−kV2+D
  4. dVdt=rV−kV2+D\frac{dV}{dt} = rV - kV^2 + DdtdV​=rV−kV2+D (correct answer)

Explanation: When modeling real-world scenarios with differential equations, you need to carefully translate each described rate into mathematical terms, paying close attention to whether each effect increases or decreases the quantity. Let's build this equation piece by piece. The investment grows due to interest "at a rate proportional to its current value" with proportionality constant rrr. This gives us +rV+rV+rV (positive because it's growth). Funds are "continuously withdrawn at a rate proportional to the square of the investment's value" with constant kkk. Since withdrawal decreases the investment, this contributes −kV2-kV^2−kV2 (negative because it's a loss). Finally, a constant amount DDD is deposited continuously, contributing +D+D+D. Combining these three effects: dVdt=rV−kV2+D\frac{dV}{dt} = rV - kV^2 + DdtdV​=rV−kV2+D Choice A incorrectly uses −kV-k\sqrt{V}−kV​ instead of −kV2-kV^2−kV2. The problem specifically states the withdrawal rate is proportional to the square of the value, not the square root. Choice B has +kV2+kV^2+kV2 instead of −kV2-kV^2−kV2, treating withdrawals as deposits, and −D-D−D instead of +D+D+D, treating deposits as withdrawals—both signs are backwards. Choice C uses just rrr instead of rVrVrV, missing that interest growth is proportional to the current value, not constant. Study tip: When translating word problems into differential equations, identify each rate of change separately, determine whether it increases (+) or decreases (-) your quantity, then combine all terms. Words like "withdrawn," "lost," or "decays" signal negative terms, while "deposited," "grows," or "increases" signal positive terms.

Question 6

An object of mass mmm falls from rest. The force of gravity is mgmgmg. Air resistance is proportional to the object's velocity v(t)v(t)v(t), with a proportionality constant k>0k>0k>0. The object also has a small rocket attached which provides a constant upward thrust TTT. Let the downward direction be positive. Which differential equation describes the velocity v(t)v(t)v(t) of the object?

  1. mdvdt=mg−kv+Tm\frac{dv}{dt} = mg - kv + Tmdtdv​=mg−kv+T
  2. mdvdt=mg+kv+Tm\frac{dv}{dt} = mg + kv + Tmdtdv​=mg+kv+T
  3. mdvdt=mg−kv−Tm\frac{dv}{dt} = mg - kv - Tmdtdv​=mg−kv−T (correct answer)
  4. mdvdt=−mg−kv−Tm\frac{dv}{dt} = -mg - kv - Tmdtdv​=−mg−kv−T

Explanation: When solving physics problems involving forces and motion, you need to carefully apply Newton's second law: Fnet=maF_{net} = maFnet​=ma, where the net force equals mass times acceleration. Since acceleration is dvdt\frac{dv}{dt}dtdv​, you're looking for an equation of the form mdvdt=Fnetm\frac{dv}{dt} = F_{net}mdtdv​=Fnet​. To find the net force, identify each force and its direction. Since downward is positive, gravity acts as +mg+mg+mg (positive). Air resistance always opposes motion, so when the object moves downward (positive velocity), air resistance acts upward (negative direction), giving −kv-kv−kv. The rocket thrust TTT acts upward, so it contributes −T-T−T to the net force. Therefore: mdvdt=mg−kv−Tm\frac{dv}{dt} = mg - kv - Tmdtdv​=mg−kv−T, which is answer C. Here's why the other options are wrong: Option A incorrectly shows thrust as +T+T+T, treating it as if it acts downward when it actually provides upward force. Option B makes two errors: it shows air resistance as +kv+kv+kv (as if it helps the object fall faster rather than opposing motion) and thrust as +T+T+T (downward instead of upward). Option D gets the signs completely backwards, showing gravity as negative when downward should be positive according to the given coordinate system. Study tip: Always establish your coordinate system first, then systematically go through each force, asking "Which direction does this force point?" and "What sign does that direction get in my coordinate system?" This methodical approach prevents sign errors, which are the most common mistakes in force problems.

Question 7

A cylindrical tank of height HHH and cross-sectional area AAA is initially full of water. At time t=0t=0t=0, a valve at the bottom is opened. Water flows out at a rate proportional to the current water pressure at the bottom of the tank. The pressure at the bottom is proportional to the depth of the water. Let h(t)h(t)h(t) be the height of the water at time ttt. Which differential equation describes h(t)h(t)h(t)? Let k>0k>0k>0 be the final proportionality constant.

  1. dhdt=−kh2\frac{dh}{dt} = -kh^2dtdh​=−kh2
  2. dhdt=−kh\frac{dh}{dt} = -k\sqrt{h}dtdh​=−kh​
  3. dhdt=−kh\frac{dh}{dt} = -khdtdh​=−kh (correct answer)
  4. dhdt=−k\frac{dh}{dt} = -kdtdh​=−k

Explanation: When modeling fluid flow from tanks, you need to connect physical relationships to create the differential equation. The key insight is translating "rate proportional to pressure" and "pressure proportional to depth" into mathematical form. Start with the physical setup: water flows out at a rate proportional to pressure at the bottom, and this pressure is proportional to the water depth h(t)h(t)h(t). So the outflow rate is proportional to h(t)h(t)h(t): outflow rate =kh= kh=kh for some constant k>0k > 0k>0. Now convert to the differential equation. The rate of change of water height equals the rate water flows in minus the rate it flows out. Since no water flows in and water flows out at rate khkhkh, you get: dhdt=0−kh=−kh\frac{dh}{dt} = 0 - kh = -khdtdh​=0−kh=−kh This gives us answer C. Looking at the wrong answers: Answer A gives dhdt=−kh2\frac{dh}{dt} = -kh^2dtdh​=−kh2, which would result from outflow rate being proportional to the square of the depth - not what the problem states. Answer B gives dhdt=−kh\frac{dh}{dt} = -k\sqrt{h}dtdh​=−kh​, which actually describes Torricelli's law where outflow depends on h\sqrt{h}h​ due to gravitational effects, but that's not the relationship described here. Answer D gives dhdt=−k\frac{dh}{dt} = -kdtdh​=−k, representing constant outflow rate regardless of water level - this ignores the pressure dependence entirely. Study tip: In tank problems, always identify what the outflow rate depends on first, then translate that physical relationship directly into the differential equation. The rate of change of height equals inflow minus outflow.

Question 8

A student is learning a new skill. Let P(t)P(t)P(t) be their proficiency level at time ttt, with P=0P=0P=0 representing no skill and P=1P=1P=1 representing mastery. The rate of learning is proportional to the product of their current proficiency and the amount of proficiency they have yet to gain. Let the constant of proportionality be k>0k>0k>0. Which differential equation models the student's proficiency P(t)P(t)P(t)?

  1. dPdt=k(1−P)\frac{dP}{dt} = k(1-P)dtdP​=k(1−P)
  2. dPdt=kP(P−1)\frac{dP}{dt} = kP(P-1)dtdP​=kP(P−1)
  3. dPdt=kP\frac{dP}{dt} = kPdtdP​=kP
  4. dPdt=kP(1−P)\frac{dP}{dt} = kP(1-P)dtdP​=kP(1−P) (correct answer)

Explanation: When you encounter word problems asking you to translate descriptions into differential equations, your key task is converting each phrase into mathematical language. Pay special attention to how rates of change depend on the current state of the system. The problem states that the learning rate is proportional to the product of two quantities: current proficiency PPP and "the amount of proficiency they have yet to gain." If mastery is P=1P = 1P=1 and current proficiency is PPP, then the remaining proficiency to gain is (1−P)(1-P)(1−P). Since the rate is proportional to their product with constant kkk, we get dPdt=kP(1−P)\frac{dP}{dt} = kP(1-P)dtdP​=kP(1−P) Looking at the wrong answers: Choice A gives dPdt=k(1−P)\frac{dP}{dt} = k(1-P)dtdP​=k(1−P), which makes learning rate depend only on remaining proficiency, ignoring current skill level entirely. Choice B has dPdt=kP(P−1)\frac{dP}{dt} = kP(P-1)dtdP​=kP(P−1), where (P−1)(P-1)(P−1) is negative since P<1P < 1P<1, making dPdt<0\frac{dP}{dt} < 0dtdP​<0 and causing proficiency to decrease over time—the opposite of learning. Choice C shows dPdt=kP\frac{dP}{dt} = kPdtdP​=kP, representing exponential growth that ignores the constraint that proficiency can't exceed mastery level. The correct answer D captures the logistic growth model, where learning accelerates as you gain skill but slows as you approach mastery. Notice how dPdt=0\frac{dP}{dt} = 0dtdP​=0 when P=0P = 0P=0 (no learning without some skill) or P=1P = 1P=1 (no learning once mastery is achieved). Study tip: When translating rate problems, identify each factor mentioned in the description and ensure your equation includes all of them mathematically.

Question 9

The concentration of a pollutant in a lake is being reduced by a natural cleaning process. The rate of reduction is proportional to the current concentration. However, a nearby factory discharges more pollutant into the lake at a rate D(t)=D0e−αtD(t) = D_0e^{-\alpha t}D(t)=D0​e−αt. Let k>0k > 0k>0 be the proportionality constant relating the removal rate to the amount of pollutant. Which differential equation models the amount of pollutant A(t)A(t)A(t) in the lake?

  1. dAdt=−kA−D0e−αt\frac{dA}{dt} = -k A - D_0 e^{-\alpha t}dtdA​=−kA−D0​e−αt
  2. dAdt=−kAV+D0e−αt\frac{dA}{dt} = -k \frac{A}{V} + D_0 e^{-\alpha t}dtdA​=−kVA​+D0​e−αt
  3. dAdt=−kA+D0e−αt\frac{dA}{dt} = -k A + D_0 e^{-\alpha t}dtdA​=−kA+D0​e−αt (correct answer)
  4. dAdt=kA+D0e−αt\frac{dA}{dt} = k A + D_0 e^{-\alpha t}dtdA​=kA+D0​e−αt

Explanation: When modeling population dynamics or concentration changes, you need to carefully account for all inflows and outflows affecting the system. The key is translating each physical process into its mathematical representation. The rate of change dAdt\frac{dA}{dt}dtdA​ equals the difference between what's being added and what's being removed. Here, two processes affect the pollutant amount: natural removal (proportional to current concentration) and factory discharge. The natural cleaning removes pollutant at a rate proportional to the current amount, so this contributes −kA-kA−kA to the rate equation (negative because it decreases the total). The factory adds pollutant at rate D0e−αtD_0e^{-\alpha t}D0​e−αt, contributing positively. Therefore: dAdt=−kA+D0e−αt\frac{dA}{dt} = -kA + D_0e^{-\alpha t}dtdA​=−kA+D0​e−αt, which is answer C. Answer A incorrectly makes the factory discharge negative, suggesting it removes rather than adds pollutant. This contradicts the problem statement that the factory "discharges more pollutant into the lake." Answer B introduces an unnecessary volume term VVV in the denominator. While concentration equals amount divided by volume, the problem states that removal is proportional to concentration, which means it's proportional to A/VA/VA/V. But since we're modeling the amount AAA, not concentration, and kkk already incorporates any volume effects, this extra division is incorrect. Answer D has the wrong sign on the removal term, making it +kA+kA+kA instead of −kA-kA−kA. This would mean more pollutant causes faster accumulation rather than faster removal. Study tip: In rate problems, carefully identify each process as adding (+) or removing (-) from your quantity of interest, then sum all contributions.

Question 10

In a reversible chemical reaction, substance A converts to substance B at a rate k1[A]k_1[A]k1​[A], and substance B converts back to A at a rate k2[B]2k_2[B]^2k2​[B]2, where [X][X][X] denotes the concentration of substance X. The total concentration [A]+[B][A] + [B][A]+[B] is a constant CCC. Let x(t)=[A]x(t) = [A]x(t)=[A]. Which differential equation describes the concentration of substance A?

  1. dxdt=−k1x+k2(C−x)2\frac{dx}{dt} = -k_1 x + k_2 (C-x)^2dtdx​=−k1​x+k2​(C−x)2 (correct answer)
  2. dxdt=k1x−k2(C−x)2\frac{dx}{dt} = k_1 x - k_2 (C-x)^2dtdx​=k1​x−k2​(C−x)2
  3. dxdt=−k1x+k2(C−x)\frac{dx}{dt} = -k_1 x + k_2 (C-x)dtdx​=−k1​x+k2​(C−x)
  4. dxdt=−k1(C−x)+k2x2\frac{dx}{dt} = -k_1 (C-x) + k_2 x^2dtdx​=−k1​(C−x)+k2​x2

Explanation: The rate of change of the concentration of A, dx/dtdx/dtdx/dt, is the rate at which A is formed minus the rate at which A is consumed. A is consumed in the forward reaction A→BA \rightarrow BA→B at a rate of k1[A]=k1xk_1[A] = k_1xk1​[A]=k1​x. A is formed in the reverse reaction B→AB \rightarrow AB→A at a rate of k2[B]2k_2[B]^2k2​[B]2. Since [A]+[B]=C[A] + [B] = C[A]+[B]=C, we have [B]=C−[A]=C−x[B] = C - [A] = C-x[B]=C−[A]=C−x. So, the rate of formation of A is k2(C−x)2k_2(C-x)^2k2​(C−x)2. Combining these gives the net rate of change for A: dx/dt=(rate of formation)−(rate of consumption)=k2(C−x)2−k1xdx/dt = (\text{rate of formation}) - (\text{rate of consumption}) = k_2(C-x)^2 - k_1xdx/dt=(rate of formation)−(rate of consumption)=k2​(C−x)2−k1​x.

Question 11

An electric circuit consists of a resistor of resistance RRR, an inductor of inductance LLL, and a time-varying voltage source E(t)=V0sin⁡(ωt)E(t) = V_0 \sin(\omega t)E(t)=V0​sin(ωt), all connected in series. According to Kirchhoff's voltage law, the sum of voltage drops across the inductor (L(dI/dt)L(dI/dt)L(dI/dt)) and the resistor (IRIRIR) equals the source voltage. A student wants to model the total magnetic flux Φ(t)\Phi(t)Φ(t) in the inductor, which is related to the current by Φ(t)=LI(t)\Phi(t) = LI(t)Φ(t)=LI(t). Which of the following is a correct differential equation for the magnetic flux Φ(t)\Phi(t)Φ(t)?

  1. dΦdt+RLΦ=V0sin⁡(ωt)\frac{d\Phi}{dt} + \frac{R}{L}\Phi = V_0 \sin(\omega t)dtdΦ​+LR​Φ=V0​sin(ωt) (correct answer)
  2. LdΦdt+RΦ=V0sin⁡(ωt)L\frac{d\Phi}{dt} + R\Phi = V_0 \sin(\omega t)LdtdΦ​+RΦ=V0​sin(ωt)
  3. 1LdΦdt+RΦ=V0sin⁡(ωt)\frac{1}{L}\frac{d\Phi}{dt} + R\Phi = V_0 \sin(\omega t)L1​dtdΦ​+RΦ=V0​sin(ωt)
  4. dΦdt+RΦ=V0sin⁡(ωt)\frac{d\Phi}{dt} + R\Phi = V_0 \sin(\omega t)dtdΦ​+RΦ=V0​sin(ωt)

Explanation: Kirchhoff's voltage law gives the equation for the current I(t)I(t)I(t): LdIdt+RI=V0sin⁡(ωt)L \frac{dI}{dt} + RI = V_0 \sin(\omega t)LdtdI​+RI=V0​sin(ωt). We are given the relation between magnetic flux Φ\PhiΦ and current III: Φ=LI\Phi = LIΦ=LI. From this, we can express III as I=Φ/LI = \Phi/LI=Φ/L. We can also find the derivative of the current: dIdt=1LdΦdt\frac{dI}{dt} = \frac{1}{L}\frac{d\Phi}{dt}dtdI​=L1​dtdΦ​. Now, substitute these expressions for III and dI/dtdI/dtdI/dt into Kirchhoff's law: L(1LdΦdt)+R(ΦL)=V0sin⁡(ωt)L \left(\frac{1}{L}\frac{d\Phi}{dt}\right) + R\left(\frac{\Phi}{L}\right) = V_0 \sin(\omega t)L(L1​dtdΦ​)+R(LΦ​)=V0​sin(ωt). Simplifying this gives the differential equation for the magnetic flux Φ(t)\Phi(t)Φ(t): dΦdt+RLΦ=V0sin⁡(ωt)\frac{d\Phi}{dt} + \frac{R}{L}\Phi = V_0 \sin(\omega t)dtdΦ​+LR​Φ=V0​sin(ωt).

Question 12

An inverted conical tank with a height of 4 meters and a top radius of 2 meters is filled with water. Water leaks out through a small hole of area aaa at the vertex. According to Torricelli's law, the rate of change of the volume of water, VVV, is given by dV/dt=−a2ghdV/dt = -a\sqrt{2gh}dV/dt=−a2gh​, where hhh is the height of the water and ggg is the acceleration due to gravity. Which differential equation describes the water height h(t)h(t)h(t)?

  1. dhdt=−a2gπh1/2\frac{dh}{dt} = -\frac{a\sqrt{2g}}{\pi}h^{1/2}dtdh​=−πa2g​​h1/2
  2. dhdt=−4a2gπh−3/2\frac{dh}{dt} = -\frac{4a\sqrt{2g}}{\pi}h^{-3/2}dtdh​=−π4a2g​​h−3/2 (correct answer)
  3. dhdt=−a2g4πh1/2\frac{dh}{dt} = -\frac{a\sqrt{2g}}{4\pi}h^{1/2}dtdh​=−4πa2g​​h1/2
  4. dhdt=−a2gπh−1/2\frac{dh}{dt} = -\frac{a\sqrt{2g}}{\pi}h^{-1/2}dtdh​=−πa2g​​h−1/2

Explanation: First, relate the volume VVV of water in the cone to its height hhh. Let rrr be the radius of the water surface at height hhh. By similar triangles, r/h=R/H=2/4=1/2r/h = R/H = 2/4 = 1/2r/h=R/H=2/4=1/2, so r=h/2r = h/2r=h/2. The volume is V=13πr2h=13π(h/2)2h=π12h3V = \frac{1}{3}\pi r^2 h = \frac{1}{3}\pi (h/2)^2 h = \frac{\pi}{12}h^3V=31​πr2h=31​π(h/2)2h=12π​h3. Differentiate with respect to time ttt using the chain rule: dV/dt=ddh(π12h3)dhdt=3π12h2dhdt=π4h2dhdtdV/dt = \frac{d}{dh}(\frac{\pi}{12}h^3) \frac{dh}{dt} = \frac{3\pi}{12}h^2 \frac{dh}{dt} = \frac{\pi}{4}h^2 \frac{dh}{dt}dV/dt=dhd​(12π​h3)dtdh​=123π​h2dtdh​=4π​h2dtdh​. Now, set this equal to the expression from Torricelli's law: π4h2dhdt=−a2gh\frac{\pi}{4}h^2 \frac{dh}{dt} = -a\sqrt{2gh}4π​h2dtdh​=−a2gh​. Solving for dh/dtdh/dtdh/dt gives dhdt=−4a2gπhh2=−4a2gπh−3/2\frac{dh}{dt} = -\frac{4a\sqrt{2g}}{\pi} \frac{\sqrt{h}}{h^2} = -\frac{4a\sqrt{2g}}{\pi}h^{-3/2}dtdh​=−π4a2g​​h2h​​=−π4a2g​​h−3/2.

Question 13

A predator-prey system involves foxes (F) and rabbits (R). Rabbits grow exponentially at rate a=0.8a = 0.8a=0.8 per month in the absence of foxes, but are consumed by foxes at a rate proportional to the product FRFRFR with constant b=0.03b = 0.03b=0.03. Foxes die exponentially at rate c=0.6c = 0.6c=0.6 per month without rabbits, but their population increases due to consuming rabbits at a rate proportional to FRFRFR with efficiency constant d=0.02d = 0.02d=0.02. However, environmental carrying capacity limits rabbit growth when R>100R > 100R>100, reducing growth by a factor of (1−R/200)(1 - R/200)(1−R/200). Which system correctly models this situation?

  1. dRdt=0.8R−0.03FR,dFdt=−0.6F+0.02FR(1−R/200)\frac{dR}{dt} = 0.8R - 0.03FR, \quad \frac{dF}{dt} = -0.6F + 0.02FR(1 - R/200)dtdR​=0.8R−0.03FR,dtdF​=−0.6F+0.02FR(1−R/200)
  2. dRdt=0.8R−0.03FR−R2/200,dFdt=−0.6F+0.02FR\frac{dR}{dt} = 0.8R - 0.03FR - R^2/200, \quad \frac{dF}{dt} = -0.6F + 0.02FRdtdR​=0.8R−0.03FR−R2/200,dtdF​=−0.6F+0.02FR
  3. dRdt=0.8R(1−R/200)−0.03FR,dFdt=−0.6F+0.02FR\frac{dR}{dt} = 0.8R(1 - R/200) - 0.03FR, \quad \frac{dF}{dt} = -0.6F + 0.02FRdtdR​=0.8R(1−R/200)−0.03FR,dtdF​=−0.6F+0.02FR (correct answer)
  4. dRdt=0.8R(1−R/100)−0.03FR,dFdt=−0.6F+0.02FR\frac{dR}{dt} = 0.8R(1 - R/100) - 0.03FR, \quad \frac{dF}{dt} = -0.6F + 0.02FRdtdR​=0.8R(1−R/100)−0.03FR,dtdF​=−0.6F+0.02FR

Explanation: When you encounter predator-prey problems, you need to carefully translate each biological process into mathematical terms, paying close attention to what affects each population's growth rate. Let's build the correct model step by step. For rabbits: they grow exponentially at rate 0.8 (giving us 0.8R0.8R0.8R), lose individuals through predation at rate 0.03FR0.03FR0.03FR, and face carrying capacity limits that reduce growth by factor (1−R/200)(1 - R/200)(1−R/200). The carrying capacity affects the natural growth term, not predation, so we get: dRdt=0.8R(1−R/200)−0.03FR\frac{dR}{dt} = 0.8R(1 - R/200) - 0.03FRdtdR​=0.8R(1−R/200)−0.03FR. For foxes: they die at rate 0.6 without food (giving us −0.6F-0.6F−0.6F) and gain population through successful predation at rate 0.02FR0.02FR0.02FR. No carrying capacity affects foxes, so: dFdt=−0.6F+0.02FR\frac{dF}{dt} = -0.6F + 0.02FRdtdF​=−0.6F+0.02FR. Now let's examine why the other options fail: Option A incorrectly applies the carrying capacity factor to the fox equation rather than the rabbit growth term, and doesn't properly model logistic growth for rabbits. Option B adds an incorrect −R2/200-R^2/200−R2/200 term to rabbits (this isn't how carrying capacity was described) and omits the carrying capacity factor from the natural growth term. Option D uses the wrong carrying capacity parameter (R/100R/100R/100 instead of R/200R/200R/200), which would make the effective carrying capacity 100 instead of 200. Study tip: In predator-prey problems, always identify which population each constraint affects. Carrying capacity typically modifies natural growth rates through logistic factors like (1−N/K)(1 - N/K)(1−N/K), where KKK is the carrying capacity.

Question 14

A chemical reaction converts substance A to substance B at a rate proportional to the amount of A present. Simultaneously, substance B decomposes back to A at a rate proportional to the amount of B present. If [A](t)[A](t)[A](t) and [B](t)[B](t)[B](t) represent concentrations at time ttt, the forward reaction rate constant is kf=0.3k_f = 0.3kf​=0.3, the reverse reaction rate constant is kr=0.1k_r = 0.1kr​=0.1, and the total concentration [A]+[B]=C[A] + [B] = C[A]+[B]=C (constant), which single differential equation in [A][A][A] only captures this system?

  1. d[A]dt=0.3[A]−0.1(C−[A])\frac{d[A]}{dt} = 0.3[A] - 0.1(C - [A])dtd[A]​=0.3[A]−0.1(C−[A])
  2. d[A]dt=−0.3[A]+0.1[B]\frac{d[A]}{dt} = -0.3[A] + 0.1[B]dtd[A]​=−0.3[A]+0.1[B]
  3. d[A]dt=−0.2[A]+0.1C\frac{d[A]}{dt} = -0.2[A] + 0.1Cdtd[A]​=−0.2[A]+0.1C
  4. d[A]dt=−0.3[A]+0.1(C−[A])\frac{d[A]}{dt} = -0.3[A] + 0.1(C - [A])dtd[A]​=−0.3[A]+0.1(C−[A]) (correct answer)

Explanation: When you encounter a reversible chemical reaction problem, you need to carefully track the direction and magnitude of each reaction step, then combine them to find the net rate of change. Let's build the differential equation systematically. Substance A decreases due to the forward reaction (A → B) at rate kf[A]=0.3[A]k_f[A] = 0.3[A]kf​[A]=0.3[A], so this contributes −0.3[A]-0.3[A]−0.3[A] to d[A]dt\frac{d[A]}{dt}dtd[A]​. Simultaneously, A increases due to the reverse reaction (B → A) at rate kr[B]=0.1[B]k_r[B] = 0.1[B]kr​[B]=0.1[B], contributing +0.1[B]+0.1[B]+0.1[B] to d[A]dt\frac{d[A]}{dt}dtd[A]​. Since the total concentration is constant ([A]+[B]=C[A] + [B] = C[A]+[B]=C), we can substitute [B]=C−[A][B] = C - [A][B]=C−[A] to eliminate B from our equation. This gives us: d[A]dt=−0.3[A]+0.1(C−[A])\frac{d[A]}{dt} = -0.3[A] + 0.1(C - [A])dtd[A]​=−0.3[A]+0.1(C−[A]) This matches answer choice D exactly. Option A has the wrong signs—it suggests A increases from the forward reaction, which contradicts the problem statement. Option B correctly identifies the individual reaction terms but doesn't eliminate [B] as required for a single equation in [A] only. Option C appears to combine terms prematurely: −0.3[A]+0.1(C−[A])=−0.3[A]+0.1C−0.1[A]=−0.4[A]+0.1C-0.3[A] + 0.1(C - [A]) = -0.3[A] + 0.1C - 0.1[A] = -0.4[A] + 0.1C−0.3[A]+0.1(C−[A])=−0.3[A]+0.1C−0.1[A]=−0.4[A]+0.1C, not −0.2[A]+0.1C-0.2[A] + 0.1C−0.2[A]+0.1C. Remember: in coupled reaction problems, always write the rate equation for each individual process first, then use conservation laws (like constant total concentration) to eliminate unwanted variables. Don't try to simplify too early—work systematically through each step.

Question 15

A pendulum's angular displacement θ(t)\theta(t)θ(t) from vertical satisfies the nonlinear equation d2θdt2+gLsin⁡θ=0\frac{d^2\theta}{dt^2} + \frac{g}{L}\sin\theta = 0dt2d2θ​+Lg​sinθ=0 for small oscillations, where g=9.8g = 9.8g=9.8 m/s² and L=2L = 2L=2 m. However, the pendulum also experiences air resistance proportional to angular velocity with damping coefficient c=0.5c = 0.5c=0.5, and a small periodic driving force F0cos⁡(ωt)F_0\cos(\omega t)F0​cos(ωt) with F0=0.1F_0 = 0.1F0​=0.1 and ω=2\omega = 2ω=2. Which equation correctly models this driven, damped pendulum using the small-angle approximation sin⁡θ≈θ\sin\theta \approx \thetasinθ≈θ?

  1. d2θdt2+0.5dθdt+9.8θ=0.1cos⁡(2t)\frac{d^2\theta}{dt^2} + 0.5\frac{d\theta}{dt} + 9.8\theta = 0.1\cos(2t)dt2d2θ​+0.5dtdθ​+9.8θ=0.1cos(2t)
  2. d2θdt2+0.5dθdt+4.9sin⁡θ=0.1cos⁡(2t)\frac{d^2\theta}{dt^2} + 0.5\frac{d\theta}{dt} + 4.9\sin\theta = 0.1\cos(2t)dt2d2θ​+0.5dtdθ​+4.9sinθ=0.1cos(2t)
  3. d2θdt2+0.5dθdt+4.9θ=0.1cos⁡(2t)\frac{d^2\theta}{dt^2} + 0.5\frac{d\theta}{dt} + 4.9\theta = 0.1\cos(2t)dt2d2θ​+0.5dtdθ​+4.9θ=0.1cos(2t) (correct answer)
  4. d2θdt2+0.5θ+4.9dθdt=0.1cos⁡(2t)\frac{d^2\theta}{dt^2} + 0.5\theta + 4.9\frac{d\theta}{dt} = 0.1\cos(2t)dt2d2θ​+0.5θ+4.9dtdθ​=0.1cos(2t)

Explanation: When you encounter pendulum problems with multiple forces, you need to systematically build the differential equation by identifying each physical effect and its mathematical representation. Start with the basic pendulum equation d2θdt2+gLsin⁡θ=0\frac{d^2\theta}{dt^2} + \frac{g}{L}\sin\theta = 0dt2d2θ​+Lg​sinθ=0. The restoring force coefficient is gL=9.82=4.9\frac{g}{L} = \frac{9.8}{2} = 4.9Lg​=29.8​=4.9. Since you're told to use the small-angle approximation sin⁡θ≈θ\sin\theta \approx \thetasinθ≈θ, this becomes d2θdt2+4.9θ=0\frac{d^2\theta}{dt^2} + 4.9\theta = 0dt2d2θ​+4.9θ=0. Next, add the damping term. Air resistance proportional to angular velocity contributes +cdθdt=+0.5dθdt+c\frac{d\theta}{dt} = +0.5\frac{d\theta}{dt}+cdtdθ​=+0.5dtdθ​ to the left side. Finally, the periodic driving force F0cos⁡(ωt)=0.1cos⁡(2t)F_0\cos(\omega t) = 0.1\cos(2t)F0​cos(ωt)=0.1cos(2t) appears on the right side as an external forcing term. Combining everything: d2θdt2+0.5dθdt+4.9θ=0.1cos⁡(2t)\frac{d^2\theta}{dt^2} + 0.5\frac{d\theta}{dt} + 4.9\theta = 0.1\cos(2t)dt2d2θ​+0.5dtdθ​+4.9θ=0.1cos(2t), which is answer C. Let's examine why the other options fail: Answer A incorrectly uses g=9.8g = 9.8g=9.8 instead of calculating gL=4.9\frac{g}{L} = 4.9Lg​=4.9. Answer B keeps sin⁡θ\sin\thetasinθ instead of applying the required small-angle approximation sin⁡θ≈θ\sin\theta \approx \thetasinθ≈θ. Answer D scrambles the damping and restoring force terms, placing the velocity term with θ\thetaθ and the displacement term with dθdt\frac{d\theta}{dt}dtdθ​. Study tip: For driven, damped oscillator problems, always follow the pattern: acceleration + damping×velocity + restoring×displacement = driving force. Double-check that coefficients match the physical setup and any approximations are properly applied.

Question 16

A drug is administered intravenously at a constant rate of 5 mg/hour into a patient's bloodstream. The drug is eliminated from the body according to first-order kinetics with elimination constant k=0.2k = 0.2k=0.2 per hour. The drug also binds reversibly to proteins in the blood: free drug becomes bound at rate kb[Df]k_b[D_f]kb​[Df​] and bound drug becomes free at rate ku[Db]k_u[D_b]ku​[Db​], where [Df][D_f][Df​] and [Db][D_b][Db​] are concentrations of free and bound drug, respectively. If kb=0.8k_b = 0.8kb​=0.8 per hour and ku=0.3k_u = 0.3ku​=0.3 per hour, and only free drug is eliminated, which equation correctly models the free drug concentration?

  1. d[Df]dt=5−0.2[Df]−0.8[Df]+0.3[Db]−0.2[Db]\frac{d[D_f]}{dt} = 5 - 0.2[D_f] - 0.8[D_f] + 0.3[D_b] - 0.2[D_b]dtd[Df​]​=5−0.2[Df​]−0.8[Df​]+0.3[Db​]−0.2[Db​]
  2. d[Df]dt=5−0.2[Df]−0.8[Df]+0.3[Db]\frac{d[D_f]}{dt} = 5 - 0.2[D_f] - 0.8[D_f] + 0.3[D_b]dtd[Df​]​=5−0.2[Df​]−0.8[Df​]+0.3[Db​] (correct answer)
  3. d[Df]dt=5−(0.2+0.8)[Df]+0.3[Db]\frac{d[D_f]}{dt} = 5 - (0.2 + 0.8)[D_f] + 0.3[D_b]dtd[Df​]​=5−(0.2+0.8)[Df​]+0.3[Db​]
  4. d[Df]dt=5−0.2([Df]+[Db])−0.8[Df]+0.3[Db]\frac{d[D_f]}{dt} = 5 - 0.2([D_f] + [D_b]) - 0.8[D_f] + 0.3[D_b]dtd[Df​]​=5−0.2([Df​]+[Db​])−0.8[Df​]+0.3[Db​]

Explanation: When modeling drug concentrations with binding and elimination, you need to carefully track all the processes that add to or remove free drug from the bloodstream. Think systematically about each rate term and its sign. For free drug concentration [Df][D_f][Df​], four processes are occurring simultaneously:

  • Input: Constant IV infusion adds 5 mg/hour (positive term: +5)
  • Elimination: Only free drug is eliminated at rate k[Df]=0.2[Df]k[D_f] = 0.2[D_f]k[Df​]=0.2[Df​] (negative term: -0.2[D_f])
  • Binding: Free drug becomes bound at rate kb[Df]=0.8[Df]k_b[D_f] = 0.8[D_f]kb​[Df​]=0.8[Df​] (negative term: -0.8[D_f])
  • Unbinding: Bound drug becomes free at rate ku[Db]=0.3[Db]k_u[D_b] = 0.3[D_b]ku​[Db​]=0.3[Db​] (positive term: +0.3[D_b])
This gives: d[Df]dt=5−0.2[Df]−0.8[Df]+0.3[Db]\frac{d[D_f]}{dt} = 5 - 0.2[D_f] - 0.8[D_f] + 0.3[D_b]dtd[Df​]​=5−0.2[Df​]−0.8[Df​]+0.3[Db​] Answer B is correct because it includes exactly these four terms with proper signs. Answer A incorrectly subtracts 0.2[Db]0.2[D_b]0.2[Db​], suggesting bound drug is also eliminated, which contradicts the problem statement that only free drug is eliminated. Answer C factors the free drug loss terms but is mathematically equivalent to B, making B the more standard form. Answer D incorrectly applies elimination to both free and bound drug by writing 0.2([Df]+[Db])0.2([D_f] + [D_b])0.2([Df​]+[Db​]), again violating the condition that only free drug is eliminated. Study tip: In compartmental models, systematically identify every process affecting your variable of interest, assign the correct sign (input = positive, output = negative), and double-check that your terms match the problem's constraints about which processes affect which compartments.

Question 17

A hot object with initial temperature 200°F is placed in a room maintained at 70°F. The object cools according to Newton's Law of Cooling with a cooling constant of k=0.1k = 0.1k=0.1 per minute. Simultaneously, an external heat source adds thermal energy at a rate equivalent to raising the object's temperature by 15°F per minute if no other thermal effects were present. What differential equation models the temperature T(t)T(t)T(t) of the object?

  1. dTdt=−0.1(T−70)+15\frac{dT}{dt} = -0.1(T - 70) + 15dtdT​=−0.1(T−70)+15 (correct answer)
  2. dTdt=−0.1(T−85)+15\frac{dT}{dt} = -0.1(T - 85) + 15dtdT​=−0.1(T−85)+15
  3. dTdt=−0.1T+7+15\frac{dT}{dt} = -0.1T + 7 + 15dtdT​=−0.1T+7+15
  4. dTdt=−0.1(T−70+15)\frac{dT}{dt} = -0.1(T - 70 + 15)dtdT​=−0.1(T−70+15)

Explanation: Newton's Law gives cooling rate -k(T - T_ambient) = -0.1(T - 70). The external heating adds +15°F/min. The total rate is dT/dt = -0.1(T - 70) + 15. Choice B incorrectly shifts the ambient temperature by the heating rate. Choice C incorrectly expands Newton's Law as -0.1T + 0.1(70) = -0.1T + 7. Choice D incorrectly includes the heating term inside the temperature difference for Newton's Law.

Question 18

An epidemic spreads through a population of 10,000 people. The rate of new infections is proportional to the product of infected individuals I(t)I(t)I(t) and susceptible individuals S(t)S(t)S(t), with transmission rate β=0.0001\beta = 0.0001β=0.0001 per person per day. Infected individuals recover at rate γ=0.05\gamma = 0.05γ=0.05 per day, becoming immune. Additionally, susceptible individuals can become immune through vaccination at a constant rate of 50 people per day, regardless of infection status. If R(t)R(t)R(t) represents recovered/immune individuals, which system models this SIR epidemic with vaccination?

  1. dSdt=−0.0001SI,dIdt=0.0001SI−0.05I−50,dRdt=0.05I+50\frac{dS}{dt} = -0.0001SI, \quad \frac{dI}{dt} = 0.0001SI - 0.05I - 50, \quad \frac{dR}{dt} = 0.05I + 50dtdS​=−0.0001SI,dtdI​=0.0001SI−0.05I−50,dtdR​=0.05I+50
  2. dSdt=−0.0001SI−50,dIdt=0.0001SI−0.05I,dRdt=0.05I+50\frac{dS}{dt} = -0.0001SI - 50, \quad \frac{dI}{dt} = 0.0001SI - 0.05I, \quad \frac{dR}{dt} = 0.05I + 50dtdS​=−0.0001SI−50,dtdI​=0.0001SI−0.05I,dtdR​=0.05I+50 (correct answer)
  3. dSdt=−0.0001SI−50,dIdt=0.0001SI−0.05I+50,dRdt=0.05I\frac{dS}{dt} = -0.0001SI - 50, \quad \frac{dI}{dt} = 0.0001SI - 0.05I + 50, \quad \frac{dR}{dt} = 0.05IdtdS​=−0.0001SI−50,dtdI​=0.0001SI−0.05I+50,dtdR​=0.05I
  4. dSdt=−0.0001SI−50,dIdt=0.0001SI−0.05I,dRdt=0.05I\frac{dS}{dt} = -0.0001SI - 50, \quad \frac{dI}{dt} = 0.0001SI - 0.05I, \quad \frac{dR}{dt} = 0.05IdtdS​=−0.0001SI−50,dtdI​=0.0001SI−0.05I,dtdR​=0.05I

Explanation: When modeling epidemic dynamics with differential equations, you need to carefully track how individuals flow between compartments (Susceptible → Infected → Recovered) and account for all processes affecting each population. Let's build the system step by step. The susceptible population S(t)S(t)S(t) decreases through two processes: infection (rate βSI=0.0001SI\beta SI = 0.0001SIβSI=0.0001SI) and vaccination (constant 50 people/day). So dSdt=−0.0001SI−50\frac{dS}{dt} = -0.0001SI - 50dtdS​=−0.0001SI−50. The infected population I(t)I(t)I(t) increases from new infections (+0.0001SI+0.0001SI+0.0001SI) and decreases from recovery (−0.05I-0.05I−0.05I). Vaccination doesn't directly affect the infected group since those 50 people come from the susceptible population. Therefore: dIdt=0.0001SI−0.05I\frac{dI}{dt} = 0.0001SI - 0.05IdtdI​=0.0001SI−0.05I. The recovered/immune population R(t)R(t)R(t) gains people through natural recovery (+0.05I+0.05I+0.05I) and vaccination (+50+50+50). So dRdt=0.05I+50\frac{dR}{dt} = 0.05I + 50dtdR​=0.05I+50. Choice B correctly captures this flow: susceptible individuals leave through infection and vaccination, infected individuals are created by infection and removed by recovery, and recovered individuals are created by both recovery and vaccination. Choice A incorrectly subtracts vaccination from the infection rate rather than the susceptible population. Choice C mistakenly adds vaccination to the infected population, which makes no biological sense. Choice D omits vaccination from the recovered equation entirely. Study tip: Always verify that your SIR system conserves population (the derivatives should sum to zero when there's no birth/death) and that each process appears with correct signs in the appropriate equations.