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: Systems In Matrix Form

Practice Systems In Matrix Form 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 / 9

0 of 9 answered

The matrix differential equation X′=(120011302)X+(t20sin⁡t)\mathbf{X}' = \begin{pmatrix} 1 & 2 & 0 \\ 0 & 1 & 1 \\ 3 & 0 & 2 \end{pmatrix}\mathbf{X} + \begin{pmatrix} t^2 \\ 0 \\ \sin t \end{pmatrix}X′=​103​210​012​​X+​t20sint​​ represents which system of scalar differential equations?

Select an answer to continue

What this quiz covers

This quiz focuses on Systems In Matrix Form, 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

The matrix differential equation X′=(120011302)X+(t20sin⁡t)\mathbf{X}' = \begin{pmatrix} 1 & 2 & 0 \\ 0 & 1 & 1 \\ 3 & 0 & 2 \end{pmatrix}\mathbf{X} + \begin{pmatrix} t^2 \\ 0 \\ \sin t \end{pmatrix}X′=​103​210​012​​X+​t20sint​​ represents which system of scalar differential equations?

  1. dxdt=x+2y+t2\frac{dx}{dt} = x + 2y + t^2dtdx​=x+2y+t2, dydt=y+z\frac{dy}{dt} = y + zdtdy​=y+z, dzdt=3x+2z\frac{dz}{dt} = 3x + 2zdtdz​=3x+2z
  2. dxdt=x+2y\frac{dx}{dt} = x + 2ydtdx​=x+2y, dydt=y+z+t2\frac{dy}{dt} = y + z + t^2dtdy​=y+z+t2, dzdt=3x+2z+sin⁡t\frac{dz}{dt} = 3x + 2z + \sin tdtdz​=3x+2z+sint
  3. dxdt=x+2y+t2\frac{dx}{dt} = x + 2y + t^2dtdx​=x+2y+t2, dydt=y+z\frac{dy}{dt} = y + zdtdy​=y+z, dzdt=3x+2z+sin⁡t\frac{dz}{dt} = 3x + 2z + \sin tdtdz​=3x+2z+sint (correct answer)
  4. dxdt=x+2y\frac{dx}{dt} = x + 2ydtdx​=x+2y, dydt=y+z\frac{dy}{dt} = y + zdtdy​=y+z, dzdt=3x+2z+sin⁡t\frac{dz}{dt} = 3x + 2z + \sin tdtdz​=3x+2z+sint

Explanation: When you encounter a matrix differential equation like X′=AX+b(t)\mathbf{X}' = A\mathbf{X} + \mathbf{b}(t)X′=AX+b(t), you're looking at a compact way to represent a system of coupled differential equations. To convert it back to scalar form, you need to perform the matrix multiplication and addition component by component. Let's work through this systematically. If X=(xyz)\mathbf{X} = \begin{pmatrix} x \\ y \\ z \end{pmatrix}X=​xyz​​, then X′=(x′y′z′)\mathbf{X}' = \begin{pmatrix} x' \\ y' \\ z' \end{pmatrix}X′=​x′y′z′​​. The right side becomes: (120011302)(xyz)+(t20sin⁡t)=(x+2yy+z3x+2z)+(t20sin⁡t)=(x+2y+t2y+z3x+2z+sin⁡t)\begin{pmatrix} 1 & 2 & 0 \\ 0 & 1 & 1 \\ 3 & 0 & 2 \end{pmatrix}\begin{pmatrix} x \\ y \\ z \end{pmatrix} + \begin{pmatrix} t^2 \\ 0 \\ \sin t \end{pmatrix} = \begin{pmatrix} x + 2y \\ y + z \\ 3x + 2z \end{pmatrix} + \begin{pmatrix} t^2 \\ 0 \\ \sin t \end{pmatrix} = \begin{pmatrix} x + 2y + t^2 \\ y + z \\ 3x + 2z + \sin t \end{pmatrix}​103​210​012​​​xyz​​+​t20sint​​=​x+2yy+z3x+2z​​+​t20sint​​=​x+2y+t2y+z3x+2z+sint​​ This gives us: x′=x+2y+t2x' = x + 2y + t^2x′=x+2y+t2, y′=y+zy' = y + zy′=y+z, z′=3x+2z+sin⁡tz' = 3x + 2z + \sin tz′=3x+2z+sint, which is answer C. Answer A incorrectly omits sin⁡t\sin tsint from the third equation. Answer B misplaces the forcing terms, putting t2t^2t2 in the second equation instead of the first. Answer D omits t2t^2t2 from the first equation entirely. Study tip: Always multiply the matrix and vector first, then add the forcing vector component-wise. Each row of the coefficient matrix corresponds to one scalar equation, and each component of the forcing vector gets added to its corresponding equation position.

Question 2

A coupled oscillator system is described by: m1d2x1dt2=−k1x1+k2(x2−x1)m_1\frac{d^2x_1}{dt^2} = -k_1x_1 + k_2(x_2 - x_1)m1​dt2d2x1​​=−k1​x1​+k2​(x2​−x1​), m2d2x2dt2=−k2(x2−x1)−k3x2m_2\frac{d^2x_2}{dt^2} = -k_2(x_2 - x_1) - k_3x_2m2​dt2d2x2​​=−k2​(x2​−x1​)−k3​x2​. To convert this to matrix form, let y=(x1x1′x2x2′)\mathbf{y} = \begin{pmatrix} x_1 \\ x_1' \\ x_2 \\ x_2' \end{pmatrix}y=​x1​x1′​x2​x2′​​​. What is the correct coefficient matrix A\mathbf{A}A such that y′=Ay\mathbf{y}' = \mathbf{A}\mathbf{y}y′=Ay?

  1. (0100−k1+k2m10k2m100001k2m20−k2+k3m20)\begin{pmatrix} 0 & 1 & 0 & 0 \\ -\frac{k_1+k_2}{m_1} & 0 & \frac{k_2}{m_1} & 0 \\ 0 & 0 & 0 & 1 \\ \frac{k_2}{m_2} & 0 & -\frac{k_2+k_3}{m_2} & 0 \end{pmatrix}​0−m1​k1​+k2​​0m2​k2​​​1000​0m1​k2​​0−m2​k2​+k3​​​0010​​ (correct answer)
  2. (0100−k1−k2m10k2m100001k2m20−k2+k3m20)\begin{pmatrix} 0 & 1 & 0 & 0 \\ -\frac{k_1-k_2}{m_1} & 0 & \frac{k_2}{m_1} & 0 \\ 0 & 0 & 0 & 1 \\ \frac{k_2}{m_2} & 0 & -\frac{k_2+k_3}{m_2} & 0 \end{pmatrix}​0−m1​k1​−k2​​0m2​k2​​​1000​0m1​k2​​0−m2​k2​+k3​​​0010​​
  3. (0100−k1+k2m10k2m100001k2m20−k2−k3m20)\begin{pmatrix} 0 & 1 & 0 & 0 \\ -\frac{k_1+k_2}{m_1} & 0 & \frac{k_2}{m_1} & 0 \\ 0 & 0 & 0 & 1 \\ \frac{k_2}{m_2} & 0 & -\frac{k_2-k_3}{m_2} & 0 \end{pmatrix}​0−m1​k1​+k2​​0m2​k2​​​1000​0m1​k2​​0−m2​k2​−k3​​​0010​​
  4. (0100−k1+k2m10−k2m100001k2m20−k2+k3m20)\begin{pmatrix} 0 & 1 & 0 & 0 \\ -\frac{k_1+k_2}{m_1} & 0 & -\frac{k_2}{m_1} & 0 \\ 0 & 0 & 0 & 1 \\ \frac{k_2}{m_2} & 0 & -\frac{k_2+k_3}{m_2} & 0 \end{pmatrix}​0−m1​k1​+k2​​0m2​k2​​​1000​0−m1​k2​​0−m2​k2​+k3​​​0010​​

Explanation: Converting the second-order system requires x1′=x1′x_1' = x_1'x1′​=x1′​ and x1′′=−k1x1+k2(x2−x1)m1=−(k1+k2)x1+k2x2m1x_1'' = \frac{-k_1x_1 + k_2(x_2-x_1)}{m_1} = \frac{-(k_1+k_2)x_1 + k_2x_2}{m_1}x1′′​=m1​−k1​x1​+k2​(x2​−x1​)​=m1​−(k1​+k2​)x1​+k2​x2​​. Similarly, x2′′=−k2(x2−x1)−k3x2m2=k2x1−(k2+k3)x2m2x_2'' = \frac{-k_2(x_2-x_1) - k_3x_2}{m_2} = \frac{k_2x_1 - (k_2+k_3)x_2}{m_2}x2′′​=m2​−k2​(x2​−x1​)−k3​x2​​=m2​k2​x1​−(k2​+k3​)x2​​. Choice B incorrectly has k1−k2k_1-k_2k1​−k2​ instead of k1+k2k_1+k_2k1​+k2​. Choice C has k2−k3k_2-k_3k2​−k3​ instead of k2+k3k_2+k_3k2​+k3​. Choice D has the wrong sign for the coupling term k2m1\frac{k_2}{m_1}m1​k2​​.

Question 3

Consider the third-order differential equation y′′′−2y′′+3y′−y=e2ty''' - 2y'' + 3y' - y = e^{2t}y′′′−2y′′+3y′−y=e2t. Using the substitution y1=yy_1 = yy1​=y, y2=y′y_2 = y'y2​=y′, y3=y′′y_3 = y''y3​=y′′, which matrix equation represents this system?

  1. ddt(y1y2y3)=(1000101−32)(y1y2y3)+(00e2t)\frac{d}{dt}\begin{pmatrix} y_1 \\ y_2 \\ y_3 \end{pmatrix} = \begin{pmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 1 & -3 & 2 \end{pmatrix}\begin{pmatrix} y_1 \\ y_2 \\ y_3 \end{pmatrix} + \begin{pmatrix} 0 \\ 0 \\ e^{2t} \end{pmatrix}dtd​​y1​y2​y3​​​=​101​01−3​002​​​y1​y2​y3​​​+​00e2t​​
  2. ddt(y1y2y3)=(010001−13−2)(y1y2y3)+(00e2t)\frac{d}{dt}\begin{pmatrix} y_1 \\ y_2 \\ y_3 \end{pmatrix} = \begin{pmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ -1 & 3 & -2 \end{pmatrix}\begin{pmatrix} y_1 \\ y_2 \\ y_3 \end{pmatrix} + \begin{pmatrix} 0 \\ 0 \\ e^{2t} \end{pmatrix}dtd​​y1​y2​y3​​​=​00−1​103​01−2​​​y1​y2​y3​​​+​00e2t​​
  3. ddt(y1y2y3)=(0100011−32)(y1y2y3)+(0e2t0)\frac{d}{dt}\begin{pmatrix} y_1 \\ y_2 \\ y_3 \end{pmatrix} = \begin{pmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ 1 & -3 & 2 \end{pmatrix}\begin{pmatrix} y_1 \\ y_2 \\ y_3 \end{pmatrix} + \begin{pmatrix} 0 \\ e^{2t} \\ 0 \end{pmatrix}dtd​​y1​y2​y3​​​=​001​10−3​012​​​y1​y2​y3​​​+​0e2t0​​
  4. ddt(y1y2y3)=(0100011−32)(y1y2y3)+(00e2t)\frac{d}{dt}\begin{pmatrix} y_1 \\ y_2 \\ y_3 \end{pmatrix} = \begin{pmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ 1 & -3 & 2 \end{pmatrix}\begin{pmatrix} y_1 \\ y_2 \\ y_3 \end{pmatrix} + \begin{pmatrix} 0 \\ 0 \\ e^{2t} \end{pmatrix}dtd​​y1​y2​y3​​​=​001​10−3​012​​​y1​y2​y3​​​+​00e2t​​ (correct answer)

Explanation: When you encounter a higher-order differential equation like this third-order equation, converting it to a first-order system using substitution is a fundamental technique. The key is systematically defining new variables for each derivative and then expressing their relationships. Starting with the substitution y1=yy_1 = yy1​=y, y2=y′y_2 = y'y2​=y′, y3=y′′y_3 = y''y3​=y′′, you need to find the derivatives of each new variable. From the definitions: y1′=y′=y2y_1' = y' = y_2y1′​=y′=y2​ and y2′=y′′=y3y_2' = y'' = y_3y2′​=y′′=y3​. For y3′=y′′′y_3' = y'''y3′​=y′′′, you use the original equation. Since y′′′−2y′′+3y′−y=e2ty''' - 2y'' + 3y' - y = e^{2t}y′′′−2y′′+3y′−y=e2t, solving for y′′′y'''y′′′ gives y′′′=2y′′−3y′+y+e2t=2y3−3y2+y1+e2ty''' = 2y'' - 3y' + y + e^{2t} = 2y_3 - 3y_2 + y_1 + e^{2t}y′′′=2y′′−3y′+y+e2t=2y3​−3y2​+y1​+e2t. This creates the system: y1′=y2y_1' = y_2y1′​=y2​, y2′=y3y_2' = y_3y2′​=y3​, y3′=y1−3y2+2y3+e2ty_3' = y_1 - 3y_2 + 2y_3 + e^{2t}y3′​=y1​−3y2​+2y3​+e2t. In matrix form, this becomes answer choice D. Choice A is wrong because the coefficient matrix has the identity structure in the wrong positions—the first two rows should have the "shift" pattern [0,1,0][0,1,0][0,1,0] and [0,0,1][0,0,1][0,0,1], not [1,0,0][1,0,0][1,0,0] and [0,1,0][0,1,0][0,1,0]. Choice B has the correct matrix structure but incorrect signs in the bottom row: it shows [−1,3,−2][-1,3,-2][−1,3,−2] instead of [1,−3,2][1,-3,2][1,−3,2]. Choice C places the forcing term e2te^{2t}e2t in the wrong position—it should be in the third component, not the second. Strategy tip: Always work systematically through the substitution, carefully tracking signs when solving for the highest derivative. The coefficient matrix follows a predictable pattern with the challenging part being the bottom row.

Question 4

The system dxdt=Ax\frac{d\mathbf{x}}{dt} = \mathbf{A}\mathbf{x}dtdx​=Ax where A=(abcd)\mathbf{A} = \begin{pmatrix} a & b \\ c & d \end{pmatrix}A=(ac​bd​) represents the scalar system dxdt=3x−2y\frac{dx}{dt} = 3x - 2ydtdx​=3x−2y, dydt=5x+y\frac{dy}{dt} = 5x + ydtdy​=5x+y. What is det⁡(A−λI)\det(\mathbf{A} - \lambda \mathbf{I})det(A−λI)?

  1. λ2+4λ+13\lambda^2 + 4\lambda + 13λ2+4λ+13
  2. λ2−4λ+13\lambda^2 - 4\lambda + 13λ2−4λ+13 (correct answer)
  3. λ2−4λ+3\lambda^2 - 4\lambda + 3λ2−4λ+3
  4. λ2−4λ−13\lambda^2 - 4\lambda - 13λ2−4λ−13

Explanation: When you encounter a linear system of differential equations, you're typically working toward finding eigenvalues, which requires computing the characteristic polynomial det⁡(A−λI)\det(\mathbf{A} - \lambda \mathbf{I})det(A−λI). First, identify the coefficient matrix A\mathbf{A}A from the given system. The equations dxdt=3x−2y\frac{dx}{dt} = 3x - 2ydtdx​=3x−2y and dydt=5x+y\frac{dy}{dt} = 5x + ydtdy​=5x+y give us: A=(3−251)\mathbf{A} = \begin{pmatrix} 3 & -2 \\ 5 & 1 \end{pmatrix}A=(35​−21​) Next, form A−λI\mathbf{A} - \lambda \mathbf{I}A−λI: A−λI=(3−λ−251−λ)\mathbf{A} - \lambda \mathbf{I} = \begin{pmatrix} 3-\lambda & -2 \\ 5 & 1-\lambda \end{pmatrix}A−λI=(3−λ5​−21−λ​) Calculate the determinant: det⁡(A−λI)=(3−λ)(1−λ)−(−2)(5)=(3−λ)(1−λ)+10\det(\mathbf{A} - \lambda \mathbf{I}) = (3-\lambda)(1-\lambda) - (-2)(5) = (3-\lambda)(1-\lambda) + 10det(A−λI)=(3−λ)(1−λ)−(−2)(5)=(3−λ)(1−λ)+10 Expanding: 3−3λ−λ+λ2+10=λ2−4λ+133 - 3\lambda - \lambda + \lambda^2 + 10 = \lambda^2 - 4\lambda + 133−3λ−λ+λ2+10=λ2−4λ+13 This matches answer choice B. Now for the wrong answers: Choice A (λ2+4λ+13\lambda^2 + 4\lambda + 13λ2+4λ+13) has the wrong sign on the middle term—this happens if you incorrectly expand (3−λ)(1−λ)(3-\lambda)(1-\lambda)(3−λ)(1−λ) or mix up the trace calculation. Choice C (λ2−4λ+3\lambda^2 - 4\lambda + 3λ2−4λ+3) gets the linear term right but miscalculates the constant term, likely forgetting to add the off-diagonal product (−2)(5)=−10(-2)(5) = -10(−2)(5)=−10. Choice D (λ2−4λ−13\lambda^2 - 4\lambda - 13λ2−4λ−13) makes sign errors in both the expansion and the off-diagonal calculation. Study tip: Remember that for a 2×2 matrix, the characteristic polynomial is always λ2−trace(A)λ+det⁡(A)\lambda^2 - \text{trace}(\mathbf{A})\lambda + \det(\mathbf{A})λ2−trace(A)λ+det(A). Here, trace = 3+1=43+1=43+1=4 and determinant = 3(1)−(−2)(5)=133(1)-(-2)(5)=133(1)−(−2)(5)=13, giving λ2−4λ+13\lambda^2-4\lambda+13λ2−4λ+13.

Question 5

Consider the system of differential equations: dxdt=3x−2y+t\frac{dx}{dt} = 3x - 2y + tdtdx​=3x−2y+t, dydt=x+4y−et\frac{dy}{dt} = x + 4y - e^tdtdy​=x+4y−et. Which matrix form correctly represents this system?

  1. ddt(xy)=(3−214)(xy)+(t−et)\frac{d}{dt}\begin{pmatrix} x \\ y \end{pmatrix} = \begin{pmatrix} 3 & -2 \\ 1 & 4 \end{pmatrix}\begin{pmatrix} x \\ y \end{pmatrix} + \begin{pmatrix} t \\ -e^t \end{pmatrix}dtd​(xy​)=(31​−24​)(xy​)+(t−et​) (correct answer)
  2. ddt(xy)=(3214)(xy)+(t−et)\frac{d}{dt}\begin{pmatrix} x \\ y \end{pmatrix} = \begin{pmatrix} 3 & 2 \\ 1 & 4 \end{pmatrix}\begin{pmatrix} x \\ y \end{pmatrix} + \begin{pmatrix} t \\ -e^t \end{pmatrix}dtd​(xy​)=(31​24​)(xy​)+(t−et​)
  3. ddt(xy)=(3−214)(xy)+(tet)\frac{d}{dt}\begin{pmatrix} x \\ y \end{pmatrix} = \begin{pmatrix} 3 & -2 \\ 1 & 4 \end{pmatrix}\begin{pmatrix} x \\ y \end{pmatrix} + \begin{pmatrix} t \\ e^t \end{pmatrix}dtd​(xy​)=(31​−24​)(xy​)+(tet​)
  4. ddt(xy)=(3−2−14)(xy)+(t−et)\frac{d}{dt}\begin{pmatrix} x \\ y \end{pmatrix} = \begin{pmatrix} 3 & -2 \\ -1 & 4 \end{pmatrix}\begin{pmatrix} x \\ y \end{pmatrix} + \begin{pmatrix} t \\ -e^t \end{pmatrix}dtd​(xy​)=(3−1​−24​)(xy​)+(t−et​)

Explanation: The correct matrix form is A. The coefficient matrix has entries from the coefficients of x and y: first row [3, -2] from 3x - 2y, second row [1, 4] from x + 4y. The forcing vector contains the non-homogeneous terms [t, -e^t]. Choice B incorrectly has +2 instead of -2. Choice C has +e^t instead of -e^t. Choice D has -1 instead of +1 in the coefficient matrix.

Question 6

A system of differential equations has the form dydt=Py+q(t)\frac{d\mathbf{y}}{dt} = \mathbf{P}\mathbf{y} + \mathbf{q}(t)dtdy​=Py+q(t) where P\mathbf{P}P is a 3×33 \times 33×3 matrix and q(t)\mathbf{q}(t)q(t) is a forcing vector. If the original system consists of one third-order equation y′′′+ay′′+by′+cy=f(t)y''' + ay'' + by' + cy = f(t)y′′′+ay′′+by′+cy=f(t), which statement about the structure of P\mathbf{P}P is correct?

  1. P\mathbf{P}P has the form (001100−c−b−a)\begin{pmatrix} 0 & 0 & 1 \\ 1 & 0 & 0 \\ -c & -b & -a \end{pmatrix}​01−c​00−b​10−a​​ and q(t)=(00f(t))\mathbf{q}(t) = \begin{pmatrix} 0 \\ 0 \\ f(t) \end{pmatrix}q(t)=​00f(t)​​
  2. P\mathbf{P}P has the form (010001−c−b−a)\begin{pmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ -c & -b & -a \end{pmatrix}​00−c​10−b​01−a​​ and q(t)=(00f(t))\mathbf{q}(t) = \begin{pmatrix} 0 \\ 0 \\ f(t) \end{pmatrix}q(t)=​00f(t)​​ (correct answer)
  3. P\mathbf{P}P has the form (010001−c−b−a)\begin{pmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ -c & -b & -a \end{pmatrix}​00−c​10−b​01−a​​ and q(t)=(f(t)00)\mathbf{q}(t) = \begin{pmatrix} f(t) \\ 0 \\ 0 \end{pmatrix}q(t)=​f(t)00​​
  4. P\mathbf{P}P has the form (010001cba)\begin{pmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ c & b & a \end{pmatrix}​00c​10b​01a​​ and q(t)=(00f(t))\mathbf{q}(t) = \begin{pmatrix} 0 \\ 0 \\ f(t) \end{pmatrix}q(t)=​00f(t)​​

Explanation: When you encounter a higher-order differential equation that needs to be converted to a system, you're using a standard technique that transforms one nth-order equation into n first-order equations by introducing new variables. To convert y′′′+ay′′+by′+cy=f(t)y''' + ay'' + by' + cy = f(t)y′′′+ay′′+by′+cy=f(t) into a first-order system, define y1=yy_1 = yy1​=y, y2=y′y_2 = y'y2​=y′, and y3=y′′y_3 = y''y3​=y′′. This gives us the relationships: y1′=y2y_1' = y_2y1′​=y2​, y2′=y3y_2' = y_3y2′​=y3​, and from the original equation, y3′=y′′′=−ay′′−by′−cy+f(t)=−ay3−by2−cy1+f(t)y_3' = y''' = -ay'' - by' - cy + f(t) = -ay_3 - by_2 - cy_1 + f(t)y3′​=y′′′=−ay′′−by′−cy+f(t)=−ay3​−by2​−cy1​+f(t). Writing this as dydt=Py+q(t)\frac{d\mathbf{y}}{dt} = \mathbf{P}\mathbf{y} + \mathbf{q}(t)dtdy​=Py+q(t) where y=(y1y2y3)\mathbf{y} = \begin{pmatrix} y_1 \\ y_2 \\ y_3 \end{pmatrix}y=​y1​y2​y3​​​, we get the system with matrix P=(010001−c−b−a)\mathbf{P} = \begin{pmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ -c & -b & -a \end{pmatrix}P=​00−c​10−b​01−a​​ and forcing vector q(t)=(00f(t))\mathbf{q}(t) = \begin{pmatrix} 0 \\ 0 \\ f(t) \end{pmatrix}q(t)=​00f(t)​​. This confirms answer B is correct. Answer A has the wrong structure in the first two rows—the 1's are misplaced, disrupting the derivative relationships. Answer C places the forcing term f(t)f(t)f(t) in the wrong position; it should appear in the equation for y3′y_3'y3′​, not y1′y_1'y1′​. Answer D has the correct matrix structure but wrong signs—the coefficients in the bottom row should be negative, not positive. Remember: the companion matrix always has 1's on the superdiagonal, zeros elsewhere in the first n−1n-1n−1 rows, and the negative coefficients in the last row. The forcing term appears in the position corresponding to the highest derivative.

Question 7

The system dxdt=ax+by+f(t)\frac{dx}{dt} = ax + by + f(t)dtdx​=ax+by+f(t), dydt=cx+dy+g(t)\frac{dy}{dt} = cx + dy + g(t)dtdy​=cx+dy+g(t) has the matrix representation X′=AX+F(t)\mathbf{X}' = \mathbf{A}\mathbf{X} + \mathbf{F}(t)X′=AX+F(t) where A=(2−341)\mathbf{A} = \begin{pmatrix} 2 & -3 \\ 4 & 1 \end{pmatrix}A=(24​−31​) and F(t)=(sin⁡tcos⁡t)\mathbf{F}(t) = \begin{pmatrix} \sin t \\ \cos t \end{pmatrix}F(t)=(sintcost​). What are the values of aaa, ccc, and g(t)g(t)g(t)?

  1. a=2a = 2a=2, c=4c = 4c=4, g(t)=cos⁡tg(t) = \cos tg(t)=cost (correct answer)
  2. a=2a = 2a=2, c=−3c = -3c=−3, g(t)=cos⁡tg(t) = \cos tg(t)=cost
  3. a=4a = 4a=4, c=2c = 2c=2, g(t)=sin⁡tg(t) = \sin tg(t)=sint
  4. a=2a = 2a=2, c=4c = 4c=4, g(t)=sin⁡tg(t) = \sin tg(t)=sint

Explanation: From the matrix form A=(abcd)=(2−341)\mathbf{A} = \begin{pmatrix} a & b \\ c & d \end{pmatrix} = \begin{pmatrix} 2 & -3 \\ 4 & 1 \end{pmatrix}A=(ac​bd​)=(24​−31​), we read a=2a = 2a=2 and c=4c = 4c=4. The forcing vector F(t)=(f(t)g(t))=(sin⁡tcos⁡t)\mathbf{F}(t) = \begin{pmatrix} f(t) \\ g(t) \end{pmatrix} = \begin{pmatrix} \sin t \\ \cos t \end{pmatrix}F(t)=(f(t)g(t)​)=(sintcost​) gives g(t)=cos⁡tg(t) = \cos tg(t)=cost. Choice B confuses ccc with b=−3b = -3b=−3. Choice C swaps aaa and ccc and confuses g(t)g(t)g(t) with f(t)f(t)f(t). Choice D has the correct aaa and ccc but wrong g(t)g(t)g(t).

Question 8

A population model tracks three species with interactions described by: dxdt=2x−xy+z\frac{dx}{dt} = 2x - xy + zdtdx​=2x−xy+z, dydt=xy−3y\frac{dy}{dt} = xy - 3ydtdy​=xy−3y, dzdt=−z+2y\frac{dz}{dt} = -z + 2ydtdz​=−z+2y. The linearization about the equilibrium point (1,2,4)(1, 2, 4)(1,2,4) gives which matrix form?

  1. ddt(uvw)=(0−112−100−2−1)(uvw)\frac{d}{dt}\begin{pmatrix} u \\ v \\ w \end{pmatrix} = \begin{pmatrix} 0 & -1 & 1 \\ 2 & -1 & 0 \\ 0 & -2 & -1 \end{pmatrix}\begin{pmatrix} u \\ v \\ w \end{pmatrix}dtd​​uvw​​=​020​−1−1−2​10−1​​​uvw​​
  2. ddt(uvw)=(2−112−3002−1)(uvw)\frac{d}{dt}\begin{pmatrix} u \\ v \\ w \end{pmatrix} = \begin{pmatrix} 2 & -1 & 1 \\ 2 & -3 & 0 \\ 0 & 2 & -1 \end{pmatrix}\begin{pmatrix} u \\ v \\ w \end{pmatrix}dtd​​uvw​​=​220​−1−32​10−1​​​uvw​​
  3. ddt(uvw)=(0−112−3002−1)(uvw)\frac{d}{dt}\begin{pmatrix} u \\ v \\ w \end{pmatrix} = \begin{pmatrix} 0 & -1 & 1 \\ 2 & -3 & 0 \\ 0 & 2 & -1 \end{pmatrix}\begin{pmatrix} u \\ v \\ w \end{pmatrix}dtd​​uvw​​=​020​−1−32​10−1​​​uvw​​
  4. ddt(uvw)=(0−112−1002−1)(uvw)\frac{d}{dt}\begin{pmatrix} u \\ v \\ w \end{pmatrix} = \begin{pmatrix} 0 & -1 & 1 \\ 2 & -1 & 0 \\ 0 & 2 & -1 \end{pmatrix}\begin{pmatrix} u \\ v \\ w \end{pmatrix}dtd​​uvw​​=​020​−1−12​10−1​​​uvw​​ (correct answer)

Explanation: When analyzing nonlinear systems of differential equations, linearization helps you understand local behavior near equilibrium points. You linearize by computing the Jacobian matrix of partial derivatives and evaluating it at the equilibrium point. For this system, you need to find the Jacobian matrix of the right-hand sides: J=(∂∂x(2x−xy+z)∂∂y(2x−xy+z)∂∂z(2x−xy+z)∂∂x(xy−3y)∂∂y(xy−3y)∂∂z(xy−3y)∂∂x(−z+2y)∂∂y(−z+2y)∂∂z(−z+2y))J = \begin{pmatrix} \frac{\partial}{\partial x}(2x - xy + z) & \frac{\partial}{\partial y}(2x - xy + z) & \frac{\partial}{\partial z}(2x - xy + z) \\ \frac{\partial}{\partial x}(xy - 3y) & \frac{\partial}{\partial y}(xy - 3y) & \frac{\partial}{\partial z}(xy - 3y) \\ \frac{\partial}{\partial x}(-z + 2y) & \frac{\partial}{\partial y}(-z + 2y) & \frac{\partial}{\partial z}(-z + 2y) \end{pmatrix}J=​∂x∂​(2x−xy+z)∂x∂​(xy−3y)∂x∂​(−z+2y)​∂y∂​(2x−xy+z)∂y∂​(xy−3y)∂y∂​(−z+2y)​∂z∂​(2x−xy+z)∂z∂​(xy−3y)∂z∂​(−z+2y)​​ Computing the partial derivatives: J=(2−y−x1yx−3002−1)J = \begin{pmatrix} 2 - y & -x & 1 \\ y & x - 3 & 0 \\ 0 & 2 & -1 \end{pmatrix}J=​2−yy0​−xx−32​10−1​​ Evaluating at the equilibrium point (1,2,4)(1, 2, 4)(1,2,4): J(1,2,4)=(2−2−1121−3002−1)=(0−112−2002−1)J(1,2,4) = \begin{pmatrix} 2 - 2 & -1 & 1 \\ 2 & 1 - 3 & 0 \\ 0 & 2 & -1 \end{pmatrix} = \begin{pmatrix} 0 & -1 & 1 \\ 2 & -2 & 0 \\ 0 & 2 & -1 \end{pmatrix}J(1,2,4)=​2−220​−11−32​10−1​​=​020​−1−22​10−1​​ Answer D matches this calculation with the correct entry −2-2−2 in position (2,2)(2,2)(2,2). Answer A has −1-1−1 instead of −2-2−2 in this position. Answer B incorrectly keeps the original constant terms like 222 and −3-3−3 that should be removed during linearization. Answer C also has −3-3−3 instead of −2-2−2 in the (2,2)(2,2)(2,2) position. Remember: linearization requires computing partial derivatives of each equation with respect to each variable, then substituting the equilibrium values. Don't confuse this with simply copying coefficients from the original equations.

Question 9

A mechanical system with damping is governed by mx¨+cx˙+kx=F(t)m\ddot{x} + c\dot{x} + kx = F(t)mx¨+cx˙+kx=F(t) where m=2m = 2m=2, c=6c = 6c=6, k=8k = 8k=8, and F(t)=4cos⁡(3t)F(t) = 4\cos(3t)F(t)=4cos(3t). Using the state variables x1=xx_1 = xx1​=x and x2=x˙x_2 = \dot{x}x2​=x˙, which matrix form is correct?

  1. ddt(x1x2)=(01−4−3)(x1x2)+(04cos⁡(3t))\frac{d}{dt}\begin{pmatrix} x_1 \\ x_2 \end{pmatrix} = \begin{pmatrix} 0 & 1 \\ -4 & -3 \end{pmatrix}\begin{pmatrix} x_1 \\ x_2 \end{pmatrix} + \begin{pmatrix} 0 \\ 4\cos(3t) \end{pmatrix}dtd​(x1​x2​​)=(0−4​1−3​)(x1​x2​​)+(04cos(3t)​)
  2. ddt(x1x2)=(01−8−6)(x1x2)+(04cos⁡(3t))\frac{d}{dt}\begin{pmatrix} x_1 \\ x_2 \end{pmatrix} = \begin{pmatrix} 0 & 1 \\ -8 & -6 \end{pmatrix}\begin{pmatrix} x_1 \\ x_2 \end{pmatrix} + \begin{pmatrix} 0 \\ 4\cos(3t) \end{pmatrix}dtd​(x1​x2​​)=(0−8​1−6​)(x1​x2​​)+(04cos(3t)​)
  3. ddt(x1x2)=(01−4−3)(x1x2)+(02cos⁡(3t))\frac{d}{dt}\begin{pmatrix} x_1 \\ x_2 \end{pmatrix} = \begin{pmatrix} 0 & 1 \\ -4 & -3 \end{pmatrix}\begin{pmatrix} x_1 \\ x_2 \end{pmatrix} + \begin{pmatrix} 0 \\ 2\cos(3t) \end{pmatrix}dtd​(x1​x2​​)=(0−4​1−3​)(x1​x2​​)+(02cos(3t)​) (correct answer)
  4. ddt(x1x2)=(10−4−3)(x1x2)+(02cos⁡(3t))\frac{d}{dt}\begin{pmatrix} x_1 \\ x_2 \end{pmatrix} = \begin{pmatrix} 1 & 0 \\ -4 & -3 \end{pmatrix}\begin{pmatrix} x_1 \\ x_2 \end{pmatrix} + \begin{pmatrix} 0 \\ 2\cos(3t) \end{pmatrix}dtd​(x1​x2​​)=(1−4​0−3​)(x1​x2​​)+(02cos(3t)​)

Explanation: When converting a second-order differential equation to a first-order system, you're transforming one equation with higher derivatives into multiple equations with only first derivatives. This is essential for numerical methods and state-space analysis. Start by defining your state variables: x1=xx_1 = xx1​=x and x2=x˙x_2 = \dot{x}x2​=x˙. This immediately gives you the first equation: x1˙=x2\dot{x_1} = x_2x1​˙​=x2​. For the second equation, you need x2˙=x¨\dot{x_2} = \ddot{x}x2​˙​=x¨, so solve the original equation for x¨\ddot{x}x¨: 2x¨+6x˙+8x=4cos⁡(3t)2\ddot{x} + 6\dot{x} + 8x = 4\cos(3t)2x¨+6x˙+8x=4cos(3t) x¨=−3x˙−4x+2cos⁡(3t)\ddot{x} = -3\dot{x} - 4x + 2\cos(3t)x¨=−3x˙−4x+2cos(3t) x2˙=−4x1−3x2+2cos⁡(3t)\dot{x_2} = -4x_1 - 3x_2 + 2\cos(3t)x2​˙​=−4x1​−3x2​+2cos(3t) This gives you the system in matrix form: ddt(x1x2)=(01−4−3)(x1x2)+(02cos⁡(3t))\frac{d}{dt}\begin{pmatrix} x_1 \\ x_2 \end{pmatrix} = \begin{pmatrix} 0 & 1 \\ -4 & -3 \end{pmatrix}\begin{pmatrix} x_1 \\ x_2 \end{pmatrix} + \begin{pmatrix} 0 \\ 2\cos(3t) \end{pmatrix}dtd​(x1​x2​​)=(0−4​1−3​)(x1​x2​​)+(02cos(3t)​) Option A has the correct matrix but wrong forcing term (forgot to divide F(t)F(t)F(t) by mmm). Option B forgot to divide all coefficients by m=2m = 2m=2, keeping −8-8−8 and −6-6−6 instead of −4-4−4 and −3-3−3. Option D has the wrong matrix structure with 111 in the wrong position, breaking the fundamental relationship x1˙=x2\dot{x_1} = x_2x1​˙​=x2​. Study tip: Always divide the entire equation by the coefficient of the highest derivative first, then carefully track how each term transforms. The top row of your matrix should always be [0,1,0,...][0, 1, 0, ...][0,1,0,...] for these standard conversions.