πŸŽ“ BookMCQ
← Back to 9. Mathematical Modelling with Differential Equations

πŸ“ Differential equations terminology and solutions (35 MCQs)

πŸ“– From Calculus β€’ 9. Mathematical Modelling with Differential Equations β€’ 35 questions available

What is Differential equations terminology and solutions?

Definition:
A differential equation involves derivatives of an unknown function, where the order is the highest derivative present, and a solution satisfies the equation identically, such as solving yβ€²+2y=0y' + 2y = 0.

Example:
For yβ€²=3x2y' = 3x^2, integrating gives y=x3+Cy = x^3 + C. If y(0)=1y(0)=1, then C=1C=1 so y=x3+1y = x^3 + 1.

Reason:
Understanding terminology helps classify equations correctly, ensuring appropriate solution methods are applied for accurate results in mathematical modeling.

5
Easy
28
Medium
2
Hard

πŸ“ All Differential equations terminology and solutions MCQs

Q1. A student claims that y=x2+Cy = x^2 + C is the general solution to the differential equation dydx=2x\frac{dy}{dx} = 2x on the interval (βˆ’βˆž,∞)(-\infty, \infty). However, another student argues this is incomplete because it does not account for singular solutions. Evaluating the definition of a general solution for first-order equations, which critique is most mathematically valid regarding the completeness of this family?

A.The solution is incomplete because y=x2+Cy=x^2+C fails at x=0x=0, requiring a piecewise definition.
B.The solution is complete because integrating 2x2x yields exactly one arbitrary constant, covering all possible integral curves for this linear equation. βœ…
C.The solution is incomplete because the integration process assumes x>0x>0, so a separate constant is needed for negative intervals.
D.The solution is incorrect because the derivative of x2+Cx^2+C is 2x+12x+1, not 2x2x.
πŸ’‘ Difficulty: medium | βœ… Correct: B

πŸ“– Explanation: For a first-order ordinary differential equation like dy/dx=2xdy/dx = 2x, the general solution must contain exactly one arbitrary constant corresponding to the order of the equation. Since the function f(x)=2xf(x)=2x is continuous everywhere, the existence and uniqueness theorem guarantees that y=x2+Cy=x^2+C encompasses every possible solution on the entire real line. There are no singular solutions lost during algebraic manipulation here, making the family complete.

Q2. Consider the differential equation dydx=y1/3\frac{dy}{dx} = y^{1/3} with initial condition y(0)=0y(0)=0. Both y=0y=0 and y=(2x/3)3/2y=(2x/3)^{3/2} satisfy this IVP. If a student applies the standard separation of variables method without checking for division by zero, they will only find the non-zero solution. What fundamental theoretical failure explains why uniqueness is violated in this specific scenario?

A.The function f(y)=y1/3f(y)=y^{1/3} is continuous but its partial derivative with respect to yy is undefined at y=0y=0, violating the Lipschitz condition. βœ…
B.The separation of variables method is inherently flawed for fractional exponents and should never be used for such equations.
C.The initial condition y(0)=0y(0)=0 is a singular point where the differential equation ceases to be first-order.
D.The general solution requires two arbitrary constants for fractional powers, and using only one restricts the solution space incorrectly.
πŸ’‘ Difficulty: hard | βœ… Correct: A

πŸ“– Explanation: This problem targets the Existence and Uniqueness Theorem conditions. While f(y)=y1/3f(y)=y^{1/3} is continuous (guaranteeing existence), its derivative 13yβˆ’2/3\frac{1}{3}y^{-2/3} blows up at y=0y=0. This discontinuity in the partial derivative means the Lipschitz condition fails, allowing multiple integral curves to pass through the same initial point. Students often mechanically separate variables, implicitly assuming yβ‰ 0y \neq 0, thereby missing the equilibrium solution and failing to recognize the theoretical boundary of the method's validity.

Q3. An engineer models a system with dydx+p(x)y=q(x)\frac{dy}{dx} + p(x)y = q(x). They derive a general solution y(x)=eβˆ’βˆ«pdx[∫qe∫pdxdx+C]y(x) = e^{-\int p dx}[\int q e^{\int p dx} dx + C]. When applying the initial condition y(x0)=y0y(x_0)=y_0, they obtain a specific value for CC. However, physical observation shows the system diverges at x=xsx=x_s where p(x)p(x) has a vertical asymptote. How should the 'solution' to this initial-value problem be correctly interpreted?

A.The derived formula is globally valid; the physical divergence indicates measurement error in determining p(x)p(x).
B.The solution is valid only on the largest open interval containing x0x_0 where both p(x)p(x) and q(x)q(x) are continuous; xsx_s marks the boundary of this domain. βœ…
C.The constant CC must be re-evaluated as a complex number to bridge the gap across the singularity at xsx_s.
D.The model is invalid because linear differential equations cannot possess vertical asymptotes in their coefficients.
πŸ’‘ Difficulty: medium | βœ… Correct: B

πŸ“– Explanation: A crucial conceptual constraint in differential equations is that solutions to initial-value problems are defined on intervals, not disjoint sets. Even if an algebraic formula exists, the solution to an IVP is restricted to the maximal interval of continuity containing the initial point x0x_0. If p(x)p(x) is discontinuous at xsx_s, the solution cannot cross this barrier. This distinguishes the mathematical object 'solution to an IVP' from a mere algebraic expression, emphasizing domain restrictions dictated by coefficient continuity.

Q4. In analyzing population dynamics, a researcher obtains the implicit relation ln⁑∣y∣=kt+C\ln|y| = kt + C from y'=ky. They conclude y=Aekty = Ae^{kt} where A=±eCA = \pm e^C. A critic notes this derivation excludes y=0y=0. If the physical system allows for extinction (population reaching zero), how does the algebraic step of exponentiation fundamentally alter the solution set compared to the original differential equation?

A.Exponentiation is a bijective map, so no solutions are lost; the critic is mistaken about extinction.
B.The step ∣y∣=ekt+C|y|=e^{kt+C} forces yy to be strictly positive or negative, structurally excluding the equilibrium solution y=0y=0 unless added manually. βœ…
C.The constant CC absorbs the zero case when Cβ†’βˆ’βˆžC \to -\infty, so the formula actually includes extinction naturally.
D.Logarithms are only defined for integers in this context, making the continuous model invalid for discrete populations.
πŸ’‘ Difficulty: medium | βœ… Correct: B

πŸ“– Explanation: This addresses a subtle but critical error analysis issue in solving separable equations. Dividing by yy to separate variables assumes yβ‰ 0y \neq 0. Subsequent exponentiation of ln⁑∣y∣\ln|y| yields ∣y∣>0|y| > 0, reinforcing the exclusion of zero. While y=0y=0 satisfies the original ODE y'=ky, it is not captured by AektAe^{kt} for any finite real AA. In modeling contexts like population biology, this 'lost' equilibrium solution represents extinction, demonstrating that algebraic manipulations can shrink the solution space, requiring explicit verification of singular cases.

Q5. Given the slope field for y' = f(x,y) where vectors are horizontal along y=2y=2 and vertical along x=0x=0. A student sketches an integral curve passing through (0,2)(0,2) that follows the horizontal direction initially but then curves upward. Analyzing the geometric constraints, why is this sketch necessarily incorrect based solely on the properties of slope fields?

A.Slope fields only indicate magnitude, not direction, so curvature cannot be determined visually.
B.At the intersection point (0,2)(0,2), the slope is undefined due to conflicting horizontal and vertical indicators, preventing any smooth curve from passing through.
C.If y'=0 along y=2y=2, then y=2y=2 is an equilibrium solution; by uniqueness, no other integral curve can touch or cross this line. βœ…
D.Vertical slopes imply infinite velocity, which is physically impossible in any realistic biological or mechanical model.
πŸ’‘ Difficulty: medium | βœ… Correct: C

πŸ“– Explanation: This question tests graph-based reasoning and the concept of equilibrium solutions. If the slope field shows horizontal segments (slope = 0) along the line y=2y=2, then the constant function y(x)=2y(x)=2 is a solution. Assuming f(x,y)f(x,y) satisfies conditions for uniqueness (e.g., continuous partial derivative), integral curves cannot intersect. Therefore, no non-constant solution can ever reach, leave, or cross the line y=2y=2. The student’s sketch violates the fundamental geometric property that distinct integral curves of a well-behaved first-order ODE are non-intersecting, making the depicted behavior mathematically impossible.

Q6. A student solves y' = 2xy and gets y=Cex2y = Ce^{x^2}. They are asked to verify if y=0y = 0 is included. They argue 'Yes, because if C=0C=0, y=0y=0.' For the equation y' = y \ln(y) with solution form involving ln⁑∣ln⁑y∣\ln|\ln y|, would setting the arbitrary constant to a specific value similarly recover the equilibrium solution y=1y=1?

A.Yes, because all equilibrium solutions correspond to C=0C=0 in the general exponential form.
B.No, because y=1y=1 makes ln⁑y=0\ln y = 0, causing the logarithmic term in the general solution to be undefined regardless of CC. βœ…
C.Yes, but only if the constant is allowed to be complex rather than real.
D.No, because y=1y=1 is an unstable equilibrium and general solutions only capture stable states.
πŸ’‘ Difficulty: medium | βœ… Correct: B

πŸ“– Explanation: This is a challenging conceptual question comparing different separable structures. In y'=2xy, the integration yields ln⁑∣y∣=x2+C\ln|y|=x^2+C, leading to y=Cex2y=Ce^{x^2} where C=0C=0 gracefully recovers y=0y=0. However, for y'=y\ln y, separation gives ∫dyyln⁑y=ln⁑∣ln⁑y∣=x+C\int \frac{dy}{y \ln y} = \ln|\ln y| = x+C. Here, y=1y=1 implies ln⁑y=0\ln y = 0, making the left side undefined. No finite real CC can satisfy this. Thus, unlike the previous case, the equilibrium y=1y=1 is structurally isolated from the general family and cannot be recovered by parameter choice, highlighting that not all lost solutions behave identically.

Q7. When verifying if y=sin⁑(x)+Cy = \sin(x) + C is a solution to y' = \cos(x) on (βˆ’βˆž,∞)(-\infty, \infty), a student substitutes and confirms equality. However, for y' = \sqrt{1-y^2}, verifying y=sin⁑(x)y=\sin(x) works only on restricted intervals. What distinguishes these two verification scenarios regarding the role of the domain?

A.Trigonometric identities are universally true, whereas square roots require non-negative arguments, restricting the domain of validity.
B.The first equation is linear while the second is nonlinear, and nonlinear solutions always have restricted domains.
C.Verification for the second equation fails because cos⁑(x)\cos(x) becomes negative, violating the non-negativity of the principal square root, thus limiting valid intervals. βœ…
D.Both equations actually have global solutions; the restriction is merely a computational artifact of choosing the principal root.
πŸ’‘ Difficulty: medium | βœ… Correct: C

πŸ“– Explanation: This application question probes understanding of implicit domain constraints in verification. For y&#039;=\sqrt{1-y^2}, the RHS requires 1βˆ’y2β‰₯01-y^2 \geq 0 AND outputs only non-negative values. While y=sin⁑xy=\sin x satisfies ∣yβˆ£β‰€1|y|\leq 1, its derivative cos⁑x\cos x is negative on intervals like (Ο€/2,3Ο€/2)(\pi/2, 3\pi/2). Since β‹…β‰₯0\sqrt{\cdot} \geq 0, the equation cos⁑x=1βˆ’sin⁑2x\cos x = \sqrt{1-\sin^2 x} fails wherever cos⁑x<0\cos x < 0. Thus, y=sin⁑xy=\sin x is only a solution on intervals where cos⁑xβ‰₯0\cos x \geq 0. This contrasts with y&#039;=\cos x, where no sign restriction exists, illustrating how the functional form of the ODE dictates the valid domain of a candidate solution.

Q8. In a mixing problem modeled by dSdt+rVS=rCin\frac{dS}{dt} + \frac{r}{V}S = rC_{in}, the general solution approaches a steady state Sss=VCinS_{ss} = VC_{in} as tβ†’βˆžt \to \infty. If a student incorrectly sets up the equation as dSdt=rCinβˆ’VrS\frac{dS}{dt} = rC_{in} - \frac{V}{r}S, how does this structural error manifest in the long-term behavior of the solution compared to the correct model?

A.The incorrect model predicts oscillatory behavior instead of convergence to a steady state.
B.The incorrect model still converges to VCinVC_{in} but with a different time constant.
C.The incorrect model converges to a steady state of r2CinV\frac{r^2 C_{in}}{V}, which has inconsistent units and wrong magnitude. βœ…
D.The incorrect model predicts unbounded growth because the coefficient of SS becomes positive.
πŸ’‘ Difficulty: medium | βœ… Correct: C

πŸ“– Explanation: This error analysis question tests dimensional consistency and structural understanding of linear ODEs. The correct rate-out term is (S/V)β‹…r(S/V) \cdot r, giving coefficient r/Vr/V with units [time]βˆ’1[time]^{-1}. The student’s erroneous term (V/r)S(V/r)S has coefficient V/rV/r with units [time][time], which is dimensionally invalid for a rate equation. Solving the incorrect ODE S&#039; + (V/r)S = rC_{in} yields steady state S=rCin/(V/r)=r2Cin/VS = rC_{in} / (V/r) = r^2 C_{in}/V. Beyond being numerically wrong, this reveals a deeper misconception: confusing concentration-to-amount conversion with volume-to-flow ratios. Recognizing unit inconsistency in the coefficient is a powerful diagnostic tool for identifying setup errors before solving.

Q9. Consider the IVP y&#039; = 3y^{2/3}, y(0)=0. The functions y=0y=0 and y=x3y=x^3 are both solutions. A numerical solver using Euler’s method with step size h=0.1h=0.1 starting at (0,0)(0,0) produces yn=0y_n = 0 for all steps. Why does this numerical method fail to detect the non-trivial solution y=x3y=x^3 despite it being a valid analytical solution?

A.Euler’s method approximates derivatives using forward differences, which cannot capture solutions that emerge tangentially from an equilibrium. βœ…
B.Numerical solvers always converge to the stable equilibrium; y=x3y=x^3 is unstable and thus invisible to computation.
C.The step size h=0.1h=0.1 is too large; reducing it to h=0.001h=0.001 would reveal the cubic solution.
D.Euler’s method requires the function to be Lipschitz continuous, and since it isn’t, the algorithm defaults to zero.
πŸ’‘ Difficulty: medium | βœ… Correct: A

πŸ“– Explanation: This mixed-concept question bridges numerical methods and theoretical uniqueness. At (0,0)(0,0), f(y)=3y2/3f(y)=3y^{2/3} has zero slope, so Euler’s update yn+1=yn+hf(yn)y_{n+1} = y_n + h f(y_n) yields 0+hβ‹…0=00 + h \cdot 0 = 0 perpetually. The method is β€˜stuck’ at the equilibrium because it uses local linear approximation, and the tangent to y=x3y=x^3 at origin is horizontalβ€”identical to y=0y=0. Without perturbation or higher-order terms, the discrete scheme cannot spontaneously depart from equilibrium. This illustrates a profound limitation: numerical methods reflect local geometry, and when uniqueness fails, they may consistently select one branch (often the trivial one) while remaining blind to others, regardless of step size refinement.

Q10. A physics student models cooling with T&#039; = k(T - T_e). They solve it as T=Te+CektT = T_e + Ce^{kt}. Given T(0)>TeT(0) > T_e and knowing objects cool, they expect C>0C > 0 and k>0k > 0. However, substitution into the ODE reveals a sign contradiction. What is the precise source of the error in interpreting the constants?

A.The constant CC must be negative because temperature decreases, contradicting the initial condition T(0)>TeT(0)>T_e.
B.The proportionality constant kk in Newton’s Law of Cooling must be negative to ensure decay when T>TeT>T_e; using positive kk implies heating.
C.The model should use T&#039; = k(T_e - T) with k>0k>0; the student’s form requires k<0k<0 for cooling, which they overlooked.
D.Both B and C are valid interpretations depending on convention, but the student failed to specify which convention was adopted. βœ…
πŸ’‘ Difficulty: medium | βœ… Correct: D

πŸ“– Explanation: This conceptual understanding question addresses sign conventions in modeling. Newton’s Law states rate of change is proportional to temperature difference. If written as T&#039;=k(T-T_e), then for cooling (T>T_e \Rightarrow T&#039;<0), kk MUST be negative. Alternatively, writing T&#039;=k(T_e-T) with k>0k>0 achieves the same physics. The student’s error isn’t mathematical but interpretive: they assumed k>0k>0 universally while using the (Tβˆ’Te)(T-T_e) form. Both conventions are legitimate, but mixing them causes sign errors. Understanding that the sign of kk is contingent on the algebraic arrangement of the difference term is essential for correct physical interpretation and avoiding contradictions between model structure and observed behavior.

Q11. When solving y&#039; + P(x)y = Q(x), a student computes integrating factor ΞΌ=e∫Pdx\mu = e^{\int P dx} but forgets the absolute value in ∫1xdx\int \frac{1}{x} dx, using ΞΌ=x\mu = x instead of ∣x∣|x|. For an IVP with x0=βˆ’2x_0 = -2, how does this omission affect the solution’s validity?

A.The solution remains valid because the absolute value only affects the constant of integration, not the functional form.
B.The solution is invalid on (βˆ’βˆž,0)(-\infty, 0) because ΞΌ=x\mu = x is negative there, violating the requirement that integrating factors be positive.
C.The solution is algebraically correct but misses solutions defined only for x>0x>0; however, since x0=βˆ’2x_0=-2, the relevant branch is captured correctly by ΞΌ=x\mu=x. βœ…
D.The omission causes the product rule reversal to fail, making the entire integration step mathematically unjustified on the negative interval.
πŸ’‘ Difficulty: medium | βœ… Correct: C

πŸ“– Explanation: This error analysis focuses on domain-specific validity in linear ODEs. The integrating factor ΞΌ=e∫Pdx\mu = e^{\int P dx} is technically eln⁑∣x∣+C=A∣x∣e^{\ln|x| + C} = A|x|. Choosing ΞΌ=x\mu = x (without abs value) works perfectly on intervals where x<0x<0 because d(x y)/dx = x y&#039; + y holds regardless of sign. The absolute value matters when crossing x=0x=0 or combining branches, but for an IVP entirely within x<0x<0, ΞΌ=x\mu=x is a valid integrating factor. The key insight is that integrating factors need not be globally positiveβ€”they must simply satisfy \mu&#039; = P\mu on the interval of interest. This nuanced understanding prevents over-correction while maintaining rigor within specified domains.

Q12. A student observes that for y&#039; = y^2, the solution y=βˆ’1x+Cy = \frac{-1}{x+C} has a vertical asymptote. They claim β€˜the solution ceases to exist beyond this asymptote.’ Another argues β€˜the solution continues on the other side as a separate branch.’ Based on the definition of a solution to an IVP, which perspective is correct and why?

A.The first student is correct; a solution to an IVP must be defined on a single connected interval containing the initial point. βœ…
B.The second student is correct; the algebraic formula defines the solution everywhere except at poles, regardless of connectivity.
C.Both are correct depending on whether the problem is pure mathematics or applied physics.
D.Neither is correct; y=y2y=y^2 has global solutions because the asymptote is removable.
πŸ’‘ Difficulty: medium | βœ… Correct: A

πŸ“– Explanation: This terminology question clarifies the rigorous definition of β€˜solution to an IVP.’ By standard ODE theory, a solution must be a differentiable function on an OPEN INTERVAL containing x0x_0. The maximal interval of existence is bounded by singularities (like vertical asymptotes) or domain boundaries. For y(0)=1y(0)=1, the solution is y=βˆ’1/(xβˆ’1)y=-1/(x-1), valid only on (βˆ’βˆž,1)(-\infty, 1). The branch x>1x>1 is a DIFFERENT solution corresponding to a different initial condition. Treating disconnected branches as one solution violates the interval requirement. This distinction is vital for understanding blow-up phenomena and ensuring that β€˜solution’ refers to a single continuous trajectory, not a collection of algebraic fragments.

Q13. In carbon dating, the model N&#039; = -kN yields N=N0eβˆ’ktN = N_0 e^{-kt}. Archaeologists measure remaining fraction f=N/N0f = N/N_0 to find age t=βˆ’ln⁑(f)/kt = -\ln(f)/k. If measurement error causes f>1f > 1 (e.g., contamination), the formula gives negative age. Beyond physical impossibility, what does this reveal about the inverse mapping from data to model parameters?

A.The logarithmic transformation amplifies small errors near f=1f=1, making age estimates highly sensitive to contamination.
B.The model assumes closed systems; f>1f>1 indicates violation of the conservation assumption underlying the exponential decay law.
C.Negative age is mathematically valid but physically meaningless, showing the model’s domain is restricted to f∈(0,1]f \in (0,1].
D.All of the above are valid interpretations of the anomaly. βœ…
πŸ’‘ Difficulty: medium | βœ… Correct: D

πŸ“– Explanation: This mixed-concept question integrates modeling, error analysis, and domain restrictions. Mathematically, ln⁑(f)\ln(f) for f>1f>1 yields negative tt, exposing the inverse function’s domain as (0,1](0,1]. Physically, f>1f>1 violates the closed-system assumption of radioactive decay, indicating contamination or measurement error. Numerically, near f=1f=1, dt/df=βˆ’1/(kf)dt/df = -1/(kf) shows high sensitivity. Thus, the anomaly simultaneously reveals mathematical domain limits, physical model violations, and numerical instability. Recognizing multi-layered implications of out-of-range data is essential for robust scientific modeling, where anomalies serve as diagnostics for theory, measurement, and computation alike.

Q14. Consider two students solving y&#039; = \sqrt{x^2 + y^2}. Student A uses polar coordinates and finds implicit solutions. Student B claims no elementary solution exists and resorts to numerical methods. If both are partially correct, what does this say about the relationship between β€˜solvable’ and β€˜expressible in elementary functions’?

A.Student B is wrong; polar coordinates always yield elementary solutions for homogeneous equations.
B.β€˜Solvable’ in ODE theory includes quadratures and special functions, not just elementary expressions; Student A’s implicit form may involve non-elementary integrals. βœ…
C.Student A’s solution is invalid because polar transformation introduces spurious solutions not present in Cartesian form.
D.Elementary solvability depends solely on the equation’s order, not its functional form.
πŸ’‘ Difficulty: hard | βœ… Correct: B

πŸ“– Explanation: This Olympiad-style question challenges assumptions about solvability. Many ODEs admit solutions via integration (quadrature) even if those integrals aren’t elementary. Polar substitution for y&#039;=\sqrt{x^2+y^2} leads to dr/dΞΈ=r1+(dr/dΞΈ)2/...dr/d\theta = r \sqrt{1+(dr/d\theta)^2}/..., which may simplify but doesn’t guarantee elementary antiderivatives. β€˜Solving’ an ODE can mean reducing to integrals, expressing via special functions, or characterizing qualitativelyβ€”not necessarily finding y=f(x)y=f(x) with polynomials/exponentials/trig. Confusing β€˜solvable’ with β€˜elementarily expressible’ is a common misconception. Advanced ODE theory embraces broader solution concepts, recognizing that many important equations (e.g., Airy, Bessel) are β€˜solved’ via named non-elementary functions, expanding the notion beyond calculus-level expectations.

Q15. A student verifies that y=x3y = x^3 satisfies y&#039; = 3x^2 and concludes it is THE general solution. Another points out that y=x3+5y = x^3 + 5 also satisfies it. What fundamental misunderstanding about the term β€˜general solution’ does the first student exhibit?

A.Confusing a particular solution with the general solution, which must include all possible solutions via arbitrary constants. βœ…
B.Believing that only monomial solutions qualify as general solutions for polynomial ODEs.
C.Assuming that verification of one solution proves completeness of the solution set.
D.Thinking that the constant of integration is optional for odd-powered functions.
πŸ’‘ Difficulty: easy | βœ… Correct: A

πŸ“– Explanation: This direct recall question tests basic terminology. The general solution of an nth-order ODE must contain n independent arbitrary constants and encompass ALL solutions on the interval. y=x3y=x^3 is merely a particular solution (with C=0). The student’s error is equating β€˜a solution’ with β€˜the general solution,’ ignoring the parametric family structure. This foundational concept underpins all IVP solving: without the arbitrary constant, one cannot satisfy arbitrary initial conditions. Reinforcing this distinction prevents later confusion when encountering singular solutions or multi-parameter families, ensuring students understand that β€˜general’ implies comprehensiveness, not just correctness of a single instance.

Q16. In modeling drug concentration C(t)C(t) with C&#039; = -kC + D(t), a student assumes D(t)D(t) is constant and solves accordingly. Later, they learn D(t)D(t) is actually periodic. How does this change in input function qualitatively alter the nature of the solution compared to the constant-input case?

A.The solution transitions from monotonic approach to steady state to oscillatory behavior around a moving baseline. βœ…
B.The solution remains exponential but with a modified decay constant dependent on the period of D(t)D(t).
C.Periodic input causes resonance, leading to unbounded growth in concentration regardless of decay rate.
D.The qualitative behavior is unchanged; only the numerical values of concentration differ.
πŸ’‘ Difficulty: medium | βœ… Correct: A

πŸ“– Explanation: This application question explores how input structure shapes solution dynamics. For constant DD, the linear ODE yields transient exponential decay plus constant steady state. For periodic D(t)D(t), the particular solution inherits periodicity (via undetermined coefficients or variation of parameters), resulting in sustained oscillations superimposed on transients. The long-term behavior shifts from static equilibrium to dynamic cycling. This illustrates a core principle: the forcing function’s character directly imprints on the solution’s asymptotic structure. Understanding this helps modelers anticipate system responsesβ€”e.g., circadian drug dosing leads to rhythmic plasma levels, not flat concentrationsβ€”and underscores that solving ODEs isn’t just computation but predicting behavioral regimes based on input types.

Q17. When solving y&#039; = f(x)y, a student divides by yy to separate variables, obtaining ln⁑∣y∣=∫fdx+C\ln|y| = \int f dx + C. They then write y=Β±eCe∫fdxy = \pm e^C e^{\int f dx} and rename Β±eC\pm e^C as AA. A peer asks, β€˜Can AA be zero?’ How should this be answered based on the derivation’s logical flow?

A.No, because eCe^C is never zero and the Β± accounts only for sign, so A=0A=0 is excluded by the algebra.
B.Yes, because we must check if y=0y=0 satisfies the original ODE separately, and if so, include A=0A=0 to capture it. βœ…
C.Only if f(x)f(x) is identically zero, otherwise AA must be nonzero.
D.The question is moot because y=0y=0 is never a solution to y&#039;=f(x)y.
πŸ’‘ Difficulty: medium | βœ… Correct: B

πŸ“– Explanation: This conceptual understanding question addresses the recovery of lost solutions. Division by yy assumes yβ‰ 0y \neq 0, so y=0y=0 is excluded from the separated derivation. However, y=0y=0 clearly satisfies y&#039;=f(x)y. The renaming A=Β±eCA = \pm e^C captures all nonzero reals, but not zero. To make y=Ae∫fdxy=A e^{\int f dx} truly general, we must VERIFY y=0y=0 is a solution and THEN allow A=0A=0. This two-step processβ€”derive nonzero family, check equilibrium, merge if validβ€”is essential for completeness. Simply declaring AA arbitrary without verification risks either including invalid solutions or excluding valid ones, highlighting the need for post-derivation validation of singular cases.

Q18. A slope field for y&#039; = x/y shows radial symmetry. A student sketches circular integral curves centered at origin. Another sketches hyperbolas. Given that ydy=xdxy dy = x dx integrates to y2βˆ’x2=Cy^2 - x^2 = C, which sketch is correct and what does this reveal about interpreting slope fields geometrically?

A.Circles are correct because radial symmetry implies rotational invariance of solutions.
B.Hyperbolas are correct; radial symmetry in slope field doesn’t imply circular solutionsβ€”it reflects homogeneity, not solution shape. βœ…
C.Both are correct on different domains; circles for C<0C<0, hyperbolas for C>0C>0.
D.Neither is correct; the actual solutions are parabolas due to the singularity at y=0y=0.
πŸ’‘ Difficulty: medium | βœ… Correct: B

πŸ“– Explanation: This graph-based question tests misinterpretation of visual cues. The slope field y&#039;=x/y has vectors pointing radially outward/inward, suggesting circular symmetry. But integrating ydy=xdxy dy = x dx gives y2βˆ’x2=Cy^2 - x^2 = C, which are hyperbolas (or lines if C=0). The radial appearance arises because slope depends only on ratio x/yx/y (homogeneous), not because solutions are circles. Students often conflate field symmetry with solution geometry. This highlights that slope fields show LOCAL direction, not GLOBAL shape; integration reveals true trajectories. Visual intuition must be verified analytically, as symmetric fields can correspond to diverse curve families, and apparent patterns may reflect equation structure rather than solution topology.

Q19. In verifying solutions to y&#039;&#039; + y = 0, a student checks y=sin⁑xy=\sin x and y=cos⁑xy=\cos x individually. They then assert y=sin⁑x+cos⁑xy=\sin x + \cos x is also a solution β€˜because both parts work.’ For the nonlinear equation y&#039; = y^2, would the same logic apply to y1=βˆ’1/xy_1 = -1/x and y2=βˆ’1/(x+1)y_2 = -1/(x+1)?

A.Yes, because superposition holds for all autonomous differential equations.
B.No, because superposition applies only to linear homogeneous equations; y1+y2y_1 + y_2 does not satisfy y&#039;=y^2. βœ…
C.Yes, but only if the initial conditions for y1y_1 and y2y_2 are compatible.
D.No, because y&#039;=y^2 has no closed-form solutions, making verification impossible.
πŸ’‘ Difficulty: medium | βœ… Correct: B

πŸ“– Explanation: This conceptual question contrasts linear vs. nonlinear solution structures. Linearity enables superposition: if y1,y2y_1, y_2 solve L[y]=0L[y]=0, so does c1y1+c2y2c_1 y_1 + c_2 y_2. Nonlinear equations like y&#039;=y^2 lack this property; (y_1+y_2)&#039; \neq (y_1+y_2)^2 generally. Testing: y1+y2=βˆ’1/xβˆ’1/(x+1)y_1+y_2 = -1/x -1/(x+1), derivative is 1/x2+1/(x+1)21/x^2 + 1/(x+1)^2, but square is [1/x+1/(x+1)]2[1/x + 1/(x+1)]^2, which differs by cross term. The student’s erroneous transfer of linear intuition to nonlinear contexts is a pervasive misconception. Emphasizing that superposition is EXCLUSIVELY a linear phenomenon prevents false generalizations and reinforces the structural divide between linear and nonlinear ODE theory.

Q20. A student solves an IVP and obtains y=x2+4y = \sqrt{x^2 + 4}. They claim the domain is all real numbers. However, the original ODE was y&#039; = x/y with y(0)=2y(0)=2. Is their domain claim correct, and what hidden constraint might they have overlooked?

A.Yes, because x2+4>0x^2+4 > 0 for all real xx, so the square root is always defined and positive. βœ…
B.No, because the ODE requires y≠0y \neq 0, and although x2+4\sqrt{x^2+4} never hits zero, the derivative x/yx/y is undefined at x=0x=0.
C.Yes, but only if we consider complex extensions; over reals, the domain excludes x=0x=0.
D.No, because the solution must satisfy the ODE at every point, and while yy is defined everywhere, the ODE’s right-hand side imposes additional smoothness requirements.
πŸ’‘ Difficulty: easy | βœ… Correct: A

πŸ“– Explanation: This application question verifies domain awareness. The solution y=x2+4y=\sqrt{x^2+4} is indeed defined and positive for all real xx, satisfying yβ‰ 0y \neq 0. At x=0x=0, y=2y=2, and y&#039;=0/2=0, matching derivative of x2+4\sqrt{x^2+4} at 0. No singularity exists. The student is correct. Distractors exploit common fears about division by zero or square roots, but here both are safely avoided. This reinforces that domain determination requires checking BOTH the solution expression AND the ODE’s requirements at every point. Sometimes apparent concerns are unfounded; careful evaluation prevents unnecessary domain restrictions while ensuring genuine constraints aren’t ignored.

Q21. When deriving the logistic model P&#039; = kP(1-P/L), a student separates variables as ∫dPP(1βˆ’P/L)=∫kdt\int \frac{dP}{P(1-P/L)} = \int k dt. After partial fractions, they get ln⁑∣Pβˆ£βˆ’ln⁑∣1βˆ’P/L∣=kt+C\ln|P| - \ln|1-P/L| = kt + C. Exponentiating yields P1βˆ’P/L=Aekt\frac{P}{1-P/L} = Ae^{kt}. If P(0)>LP(0) > L, what sign must AA have, and why is this physically significant?

A.AA must be negative because 1βˆ’P/L<01-P/L < 0, and this reflects population exceeding carrying capacity, leading to decline. βœ…
B.AA must be positive because population is always positive, and the negative denominator is absorbed into the constant.
C.AA can be either sign depending on integration constant choice; physical meaning is irrelevant to algebra.
D.AA must be imaginary because P>LP>L makes the logarithm argument negative, indicating model breakdown.
πŸ’‘ Difficulty: medium | βœ… Correct: A

πŸ“– Explanation: This error analysis connects algebraic signs to physical regimes. If P(0)>LP(0) > L, then 1βˆ’P/L<01-P/L < 0, so P/(1βˆ’P/L)<0P/(1-P/L) < 0. Thus A=[P/(1βˆ’P/L)]eβˆ’kβ‹…0<0A = [P/(1-P/L)] e^{-k \cdot 0} < 0. The negative AA isn’t a flawβ€”it encodes the supra-capacity state. As tt increases, AektAe^{kt} becomes more negative, driving PP down toward LL. Misinterpreting AA’s sign as erroneous would miss this crucial dynamics. The algebra faithfully tracks physical reality: sign changes in intermediate expressions correspond to regime shifts. Understanding this prevents discarding valid solutions and reinforces that constants carry physical information beyond mere placeholders.

Q22. A numerical analyst uses Euler’s method on y&#039; = -100y, y(0)=1 with step size h=0.05h=0.05. The approximation oscillates wildly instead of decaying smoothly. Knowing the exact solution is eβˆ’100te^{-100t}, what does this instability reveal about the interaction between method and equation stiffness?

A.Euler’s method is unconditionally stable for all linear ODEs; the oscillation indicates coding error.
B.The step size exceeds the stability region for this stiff equation; hh must satisfy ∣1+hλ∣<1|1 + h\lambda| < 1 for decay. βœ…
C.Stiffness only affects implicit methods; explicit Euler should handle this fine with any hh.
D.The oscillation is expected and matches the true solution’s rapid decay; it’s a visualization artifact.
πŸ’‘ Difficulty: medium | βœ… Correct: B

πŸ“– Explanation: This challenging question links numerical stability to ODE theory. For y&#039;=\lambda y, Euler gives yn+1=(1+hΞ»)yny_{n+1}=(1+h\lambda)y_n. Stability requires ∣1+hλ∣<1|1+h\lambda|<1. With Ξ»=βˆ’100\lambda=-100, need ∣1βˆ’100h∣<1β‡’h<0.02|1-100h|<1 \Rightarrow h<0.02. Using h=0.05h=0.05 gives multiplier βˆ’4-4, causing alternating sign growth. This isn’t a bugβ€”it’s inherent to explicit methods on stiff systems. The lesson: numerical validity depends on matching step size to equation timescales. Blindly applying algorithms without stability analysis yields garbage. This bridges theory and computation, showing that β€˜solving’ an ODE numerically requires understanding the method’s domain of applicability, not just implementing formulas.

Q23. In solving y&#039; + 2xy = x, a student finds integrating factor ex2e^{x^2} and writes y=eβˆ’x2∫xex2dxy = e^{-x^2} \int x e^{x^2} dx. They evaluate integral as 12ex2+C\frac{1}{2} e^{x^2} + C, getting y=12+Ceβˆ’x2y = \frac{1}{2} + C e^{-x^2}. Another student gets y=12+Ceβˆ’x2+eβˆ’x2y = \frac{1}{2} + C e^{-x^2} + e^{-x^2}. Are these equivalent, and if so, why?

A.No, the second student made an arithmetic error in integration.
B.Yes, because the extra eβˆ’x2e^{-x^2} term is absorbed into the arbitrary constant CC upon renaming. βœ…
C.No, because constants of integration cannot depend on xx, making the second form invalid.
D.Yes, but only if the initial condition specifies C=1C=1 for the second form.
πŸ’‘ Difficulty: easy | βœ… Correct: B

πŸ“– Explanation: This direct recall question tests understanding of arbitrary constants. The integral ∫xex2dx=12ex2+K\int x e^{x^2} dx = \frac{1}{2} e^{x^2} + K. Multiplying by eβˆ’x2e^{-x^2} gives 12+Keβˆ’x2\frac{1}{2} + K e^{-x^2}. If a student mistakenly adds another eβˆ’x2e^{-x^2}, it combines as (K+1)eβˆ’x2(K+1)e^{-x^2}, which is still Ceβˆ’x2C e^{-x^2} for some CC. Arbitrary constants absorb additive terms of the same functional form. This flexibility is why we don’t worry about β€˜correct’ constant placement during integrationβ€”only that ONE arbitrary constant appears in final answer. Recognizing equivalence under constant renaming prevents false distinctions between correct solution forms and builds confidence in solution representation.

Q24. A model for tumor growth uses V&#039; = kV \ln(V_0/V). At V=V0V=V_0, growth stops. A student linearizes near V0V_0 by setting u=V0βˆ’Vu = V_0 - V, obtaining u&#039; \approx -k u. They conclude small deviations decay exponentially. What assumption underlies this local analysis, and when might it fail?

A.Assumes ln⁑(V0/V)β‰ˆ(V0βˆ’V)/V0\ln(V_0/V) \approx (V_0-V)/V_0 via Taylor expansion; fails if deviations are large or VV approaches 0. βœ…
B.Assumes kk is constant; fails if tumor microenvironment changes over time.
C.Assumes exponential decay is universal; fails for non-autonomous systems.
D.No assumptions are needed; linearization is always valid near equilibria.
πŸ’‘ Difficulty: medium | βœ… Correct: A

πŸ“– Explanation: This application question examines local vs. global behavior. Linearization replaces nonlinear ln⁑(V0/V)\ln(V_0/V) with first-order Taylor approx (V0βˆ’V)/V0(V_0-V)/V_0 near V=V0V=V_0. This yields exponential decay of perturbations, valid ONLY locally. If tumor shrinks significantly (Vβ‰ͺV0V \ll V_0) or grows beyond linear regime, approximation breaks. Also, ln⁑(V0/V)\ln(V_0/V) singular at V=0V=0, so linearization invalid near extinction. The key insight: linear stability analysis describes NEIGHBORHOOD behavior, not global dynamics. Confusing local exponential decay with global fate is dangerous in biological modeling. Always validate linear predictions against full nonlinear behavior, especially when trajectories may exit the region where approximation holds.

Q25. When solving y&#039; = \frac{x+y}{x-y}, a student recognizes homogeneity and substitutes v=y/xv=y/x, getting separable equation. After solving, they back-substitute to get implicit solution. A peer suggests using polar coordinates instead. Comparing methods, what advantage does polar offer for THIS specific equation?

A.Polar coordinates eliminate the need for back-substitution entirely.
B.The equation becomes dr/dΞΈ=rdr/d\theta = r, which is simpler than the rational function in vv. βœ…
C.Polar coordinates always yield explicit solutions for homogeneous equations.
D.There is no advantage; both methods require identical computational effort.
πŸ’‘ Difficulty: medium | βœ… Correct: B

πŸ“– Explanation: This mixed-concept question compares solution strategies. The Cartesian homogeneous sub v=y/xv=y/x leads to dv/dx=...dv/dx = ... with messy rational functions. In polar, x=rcos⁑θ,y=rsin⁑θx=r\cos\theta, y=r\sin\theta, the equation simplifies dramatically to dr/dΞΈ=rdr/d\theta = r because numerator/denominator become r(cos⁑+sin⁑)/r(cosβ‘βˆ’sin⁑)r(\cos+\sin)/r(\cos-\sin), and geometric relations reduce complexity. This isn’t universalβ€”it exploits the specific symmetry of (x+y)/(xβˆ’y)(x+y)/(x-y). Choosing coordinate systems aligned with equation structure can transform intractable algebra into simple ODEs. This strategic thinking elevates problem-solving beyond rote method application, encouraging analysis of equation geometry before committing to a technique. Recognizing when polar simplifies homogeneous equations is a valuable heuristic in advanced ODE work.

Q26. A student claims that since y=exy=e^x solves y&#039;=y, then y=ex2y=e^{x^2} must solve y&#039;=2xy by β€˜chain rule analogy.’ While coincidentally correct, their reasoning is flawed. What is the precise logical gap in extending solution forms via analogy?

A.Chain rule applies to differentiation, not to constructing ODE solutions; correctness here is accidental, not deductive. βœ…
B.The analogy fails because ex2e^{x^2} grows faster than exe^x, violating growth rate constraints.
C.Analogies are never valid in mathematics; only formal proofs establish truth.
D.The student confused independent variable scaling with functional composition, which are distinct operations.
πŸ’‘ Difficulty: medium | βœ… Correct: A

πŸ“– Explanation: This error analysis targets faulty generalization. While d/dx(e^{u}) = e^u u&#039;, inferring that y=f(u(x))y=f(u(x)) solves y&#039;=g(x)f(u) whenever ff solves y&#039;=f is unjustified. Here, ex2e^{x^2} works because d/dx(ex2)=2xex2d/dx(e^{x^2})=2x e^{x^2}, matching 2xy2xy. But for y&#039;=3y, analogously guessing e3x2e^{3x^2} fails: derivative is 6xe3x2β‰ 3e3x26x e^{3x^2} \neq 3 e^{3x^2}. Correctness in one case doesn’t validate the method. Solutions must be derived systematically, not guessed by pattern matching. This reinforces that mathematical truth requires verification, not plausibility, and prevents overconfidence in heuristic reasoning that may succeed accidentally but fail catastrophically elsewhere.

Q27. In an RL circuit modeled by LI&#039; + RI = V(t), the homogeneous solution decays as eβˆ’(R/L)te^{-(R/L)t}. A student argues β€˜since resistance dissipates energy, larger RR should slow decay.’ But formula shows larger RR increases decay rate. What resolves this apparent paradox?

A.Energy dissipation rate is I2RI^2 R; larger RR increases instantaneous power loss, accelerating current decay despite higher resistance. βœ…
B.The student confused resistance with inductance; larger LL slows decay, not RR.
C.Decay rate refers to energy, not current; current decays faster but energy decays slower with larger RR.
D.There is no paradox; the student’s intuition is simply wrong and should be discarded.
πŸ’‘ Difficulty: medium | βœ… Correct: A

πŸ“– Explanation: This conceptual question reconciles math with physics. Decay rate R/LR/L governs CURRENT, not energy. Power dissipated is I2RI^2 R; for fixed II, larger RR means MORE energy loss per second, depleting stored magnetic energy faster. Though resistance opposes current flow, in transient decay it enhances dissipation. The student’s intuition conflated steady-state opposition with transient dissipation. Math correctly captures physics: dI/dt=βˆ’(R/L)IdI/dt = -(R/L)I shows stronger damping with larger RR. Resolving such paradoxes requires distinguishing quantities (current vs. energy) and understanding how parameters affect different aspects of system behavior. This deepens modeling competence beyond formula memorization.

Q28. A student solves y&#039; = y \tan x and gets y=Csec⁑xy = C \sec x. They state domain is (βˆ’Ο€/2,Ο€/2)(-\pi/2, \pi/2) due to secant’s singularity. But initial condition is y(Ο€)=2y(\pi) = 2. How should the solution be adjusted, and what principle governs this adjustment?

A.Shift the interval to (Ο€/2,3Ο€/2)(\pi/2, 3\pi/2) containing Ο€\pi; solutions are defined on maximal intervals containing x0x_0. βœ…
B.Use absolute value: y=C∣sec⁑x∣y = C |\sec x| to extend domain across singularities.
C.The solution is invalid because tan⁑x\tan x is discontinuous at Ο€\pi, so no solution exists.
D.Keep original domain; initial conditions outside domain are ignored.
πŸ’‘ Difficulty: medium | βœ… Correct: A

πŸ“– Explanation: This application question emphasizes interval selection for IVPs. tan⁑x\tan x and sec⁑x\sec x have singularities at Ο€/2+nΟ€\pi/2 + n\pi. For x0=Ο€x_0=\pi, the maximal interval of continuity containing Ο€\pi is (Ο€/2,3Ο€/2)(\pi/2, 3\pi/2). On this interval, sec⁑x<0\sec x < 0, so y=Csec⁑xy=C\sec x with C=βˆ’2C=-2 satisfies y(Ο€)=2y(\pi)=2. The solution isn’t β€˜adjusted’ algebraically but CONTEXTUALLY by selecting appropriate branch. Solutions to IVPs are tied to specific intervals of coefficient continuity. Ignoring this leads to domain mismatches. This principle ensures solutions are meaningful at initial points and respects the local nature of ODE theory, where global formulas may require piecewise interpretation across singularities.

Q29. When verifying y=x2sin⁑(1/x)y = x^2 \sin(1/x) (with y(0)=0y(0)=0) as solution to some ODE, a student computes y&#039; = 2x \sin(1/x) - \cos(1/x) for xβ‰ 0x \neq 0 and claims y&#039;(0)=0 by inspection. Why is this derivative evaluation at zero invalid, and what is the correct approach?

A.Derivative at zero must use limit definition; y&#039; formula for xβ‰ 0x \neq 0 doesn’t apply at x=0x=0 due to oscillation. βœ…
B.The function isn’t differentiable at zero because sin⁑(1/x)\sin(1/x) is undefined there.
C.The student’s value is correct; the issue is only with higher derivatives.
D.Oscillatory functions never have derivatives at accumulation points of zeros.
πŸ’‘ Difficulty: medium | βœ… Correct: A

πŸ“– Explanation: This Olympiad-style question tests rigorous differentiation. The expression 2xsin⁑(1/x)βˆ’cos⁑(1/x)2x \sin(1/x) - \cos(1/x) is valid ONLY for xβ‰ 0x \neq 0. At x=0x=0, must use lim⁑hβ†’0[h2sin⁑(1/h)βˆ’0]/h=lim⁑hsin⁑(1/h)=0\lim_{h\to0} [h^2 \sin(1/h) - 0]/h = \lim h \sin(1/h) = 0. So y&#039;(0)=0 IS correct, but NOT from plugging into the xβ‰ 0x \neq 0 formula (which oscillates wildly). The student got right answer via wrong method. This distinction matters: derivative existence at a point requires limit definition when function is defined piecewise or has singular structure. Relying on differentiated formulas outside their domain risks errors in more pathological cases. Rigor demands respecting domain boundaries in calculus operations.

Q30. A population model P&#039; = rP(1-P/K) has equilibria at P=0P=0 and P=KP=K. A student performs phase line analysis and labels P=0P=0 as unstable, P=KP=K as stable. They then claim β€˜all positive initial populations converge to KK.’ Is this always true, and what edge case might invalidate it?

A.Yes, for all P(0)>0P(0)>0, convergence to KK is guaranteed by logistic dynamics.
B.No, if P(0)=0P(0)=0, population stays at zero; also, if model allows P<0P<0, behavior differs, but biologically Pβ‰₯0P \geq 0. βœ…
C.No, convergence only occurs if r>0r>0; if r<0r<0, KK becomes unstable.
D.The claim is true mathematically but false biologically due to Allee effects not captured in logistic model.
πŸ’‘ Difficulty: medium | βœ… Correct: B

πŸ“– Explanation: This conceptual question examines equilibrium classification nuances. Phase line shows P=0P=0 unstable (repeller for P>0P>0), P=KP=K stable (attractor). For ANY P(0)>0P(0)>0, solution β†’ KK. BUT P(0)=0P(0)=0 is equilibrium itself, so no convergence to KK. The student’s β€˜all positive’ excludes zero, so technically correct. However, distractor B highlights the boundary case. In practice, stating β€˜all positive’ is precise, but learners often forget P=0P=0 is a valid (if trivial) solution. Emphasizing that equilibria are solutions themselves prevents overgeneralization. Also, noting biological realism (Allee effects) reminds that mathematical models have limits, but within logistic framework, the statement holds for P(0)>0P(0)>0.

Q31. In solving y&#039; = \sqrt{y}, a student gets 2y=x+C2\sqrt{y} = x + C and solves y=(x+C)2/4y = (x+C)^2/4. They apply y(0)=0y(0)=0 to get C=0C=0, yielding y=x2/4y=x^2/4. But y=0y=0 also satisfies IVP. Why did separation miss y=0y=0, and how to reconcile?

A.Division by y\sqrt{y} assumes y>0y>0; y=0y=0 is equilibrium solution found by inspection, not from general formula. βœ…
B.The constant CC should be infinite to capture y=0y=0, but this isn’t standard.
C.Separation of variables is invalid for square root functions; only linear methods work.
D.The solution y=x2/4y=x^2/4 is actually incorrect; only y=0y=0 satisfies the IVP.
πŸ’‘ Difficulty: medium | βœ… Correct: A

πŸ“– Explanation: This error analysis revisits lost solutions. Separating dy/y=dxdy/\sqrt{y} = dx requires yβ‰ 0y \neq 0. Integration yields 2y=x+C2\sqrt{y}=x+C, implying y>0y>0 (since sqrt β‰₯0). Setting y(0)=0y(0)=0 gives C=0C=0, so y=x2/4y=x^2/4 for xβ‰₯0xβ‰₯0. But y=0y=0 also satisfies ODE and IC. It wasn’t found because division excluded it. Reconciliation: general solution is y=0y=0 OR y=(x+C)2/4y=(x+C)^2/4 for xβ‰₯βˆ’Cxβ‰₯-C. For IVP y(0)=0y(0)=0, BOTH are valid, demonstrating non-uniqueness. This reinforces that separation can miss equilibria, and IVPs with non-Lipschitz RHS may have multiple solutions. Always check y=0y=0 separately when dividing by functions of yy.

Q32. A student uses integrating factor ΞΌ(x)=e∫P(x)dx\mu(x) = e^{\int P(x) dx} for y&#039; + P(x)y = Q(x). They ask, β€˜Why can we ignore the constant of integration when computing ΞΌ\mu?’ What is the most precise justification?

A.Because any constant factor cancels out when multiplying through and dividing by ΞΌ\mu in the final solution formula. βœ…
B.Because the constant of integration is always zero by convention in integrating factor calculations.
C.Because including a constant would make ΞΌ\mu non-differentiable, breaking the product rule.
D.Because the method only works for definite integrals, where constants don’t appear.
πŸ’‘ Difficulty: easy | βœ… Correct: A

πŸ“– Explanation: This direct recall question addresses a common procedural query. Let ΞΌ=e∫Pdx+C=Ke∫Pdx\mu = e^{\int P dx + C} = K e^{\int P dx} where K=eCK=e^C. Multiplying ODE: K e^{\int P} y&#039; + K e^{\int P} P y = K e^{\int P} Q. Left side is d/dx(Ke∫Py)d/dx (K e^{\int P} y). Integrating: Ke∫Py=∫Ke∫PQdx+DK e^{\int P} y = \int K e^{\int P} Q dx + D. Divide by Ke∫PK e^{\int P}: y=eβˆ’βˆ«P[∫Ke∫PQdx+D]/K=eβˆ’βˆ«P[∫e∫PQdx+D/K]y = e^{-\int P} [\int K e^{\int P} Q dx + D]/K = e^{-\int P} [\int e^{\int P} Q dx + D/K]. Constant KK cancels, and D/KD/K becomes new arbitrary constant. Thus, choosing K=1K=1 (i.e., ignoring constant) simplifies without loss of generality. This algebraic cancellation justifies the shortcut.

Q33. In a falling object model with air resistance mv&#039; = mg - cv, terminal velocity is vt=mg/cv_t = mg/c. A student derives v(t)=vt(1βˆ’eβˆ’ct/m)v(t) = v_t (1 - e^{-ct/m}) for v(0)=0v(0)=0. They note v(t)<vtv(t) < v_t for all finite tt. Does this mean terminal velocity is never reached, and how to interpret β€˜terminal’ physically?

A.Mathematically, vtv_t is asymptotic limit, never attained in finite time; physically, it’s approached arbitrarily closely, serving as effective max speed. βœ…
B.Terminal velocity is reached when t=m/ct = m/c, after which velocity is constant.
C.The model is flawed because real objects do reach terminal velocity in finite time.
D.β€˜Terminal’ refers to the end of the fall, not a velocity limit; the student misinterpreted the term.
πŸ’‘ Difficulty: medium | βœ… Correct: A

πŸ“– Explanation: This conceptual question clarifies asymptotic behavior. The solution v(t)=vt(1βˆ’eβˆ’kt)v(t) = v_t(1-e^{-kt}) approaches vtv_t as tβ†’βˆžtβ†’βˆž but equals it only at infinity. Physically, after ~3-5 time constants (m/cm/c), vv is within 5% of vtv_t, indistinguishable in practice. β€˜Terminal’ denotes the limiting value governing long-term dynamics, not a finitely achieved state. Confusing mathematical limit with physical attainment is common. Understanding asymptotes as practical thresholds bridges theory and reality. This also highlights that models describe trends, not exact milestones; engineering tolerances define β€˜reached,’ not pure math. Grasping this distinction prevents misapplication of theoretical results to real-world timing questions.

Q34. A student solves y&#039; = 2x and writes general solution y=x2+Cy = x^2 + C. They then say β€˜for each CC, this is a particular solution.’ Is this terminology correct, and what is the proper distinction between β€˜general’ and β€˜particular’ solutions?

A.Incorrect; y=x2+Cy=x^2+C is the general solution; assigning specific CC yields a particular solution. βœ…
B.Correct; each instance of CC defines a general solution for that parameter value.
C.Terminology varies by textbook; both usages are acceptable in modern mathematics.
D.β€˜Particular solution’ refers only to solutions satisfying initial conditions, not arbitrary constant assignments.
πŸ’‘ Difficulty: easy | βœ… Correct: A

πŸ“– Explanation: This direct recall question reinforces foundational terminology. The GENERAL solution contains arbitrary constants equal to ODE order and represents the COMPLETE family of solutions. A PARTICULAR solution is obtained by specifying values for those constants (often via initial conditions). Saying β€˜y=x2+Cy=x^2+C is a particular solution for each CC’ misuses terms; it’s the general solution PARAMETERIZED by CC. Only when CC is fixed (e.g., C=5C=5) do we have a particular solution. Precise language prevents confusion in communication and assessment. This distinction is universal across ODE theory and essential for correctly describing solution sets, especially when discussing uniqueness, existence, and solution structure in advanced topics.

Q35. When analyzing y&#039; = f(x,y), a student notes that if ff and βˆ‚f/βˆ‚yβˆ‚f/βˆ‚y are continuous near (x0,y0)(x_0,y_0), then unique solution exists. They then encounter y&#039; = y^{1/3} at (0,0)(0,0) and observe non-uniqueness. They conclude β€˜continuity of ff alone is insufficient for uniqueness.’ Is this conclusion valid, and what does it imply about theorem conditions?

A.Valid; continuity of ff guarantees existence but not uniqueness; Lipschitz condition (or continuous βˆ‚f/βˆ‚yβˆ‚f/βˆ‚y) is additionally required. βœ…
B.Invalid; y1/3y^{1/3} is not continuous at 0, so the theorem doesn’t apply.
C.Valid, but only for autonomous equations; non-autonomous equations require different criteria.
D.Invalid; the student miscalculated the partial derivative; it is continuous at 0.
πŸ’‘ Difficulty: medium | βœ… Correct: A

πŸ“– Explanation: This conceptual understanding question clarifies Existence-Uniqueness Theorem hypotheses. Continuity of ff ensures EXISTENCE (Peano’s theorem). UNIQUENESS requires stronger condition: Lipschitz continuity in yy, often verified via continuous βˆ‚f/βˆ‚yβˆ‚f/βˆ‚y. For y1/3y^{1/3}, ff continuous at 0 β‡’ existence, but βˆ‚f/βˆ‚y=(1/3)yβˆ’2/3βˆ‚f/βˆ‚y = (1/3)y^{-2/3} unbounded β‡’ uniqueness fails. Student’s conclusion is correct and important. Many assume continuity suffices for both; this example dispels that myth. Understanding the hierarchy of conditions prevents misapplication of theorems and fosters deeper appreciation of why certain equations behave pathologically. This knowledge is crucial for diagnosing model issues and selecting appropriate solution methods when standard assumptions break down.

πŸ”— Related Topics (MCQs)