4.7 Practice Problems
Problem 4.7.1. Verify that \(y_1 = t\) and \(y_2 = t^2\) solve \(t^2y'' - 2ty' + 2y = 0\) on \(t > 0\). Compute their Wronskian directly and from Abel’s identity.
Show solution
Solution. \(t^2 \cdot 0 - 2t + 2t = 0\) and \(t^2 \cdot 2 - 2t \cdot 2t + 2t^2 = 0\). Directly, \(W = t \cdot 2t - t^2 \cdot 1 = t^2\). In normal form the equation is \(y'' - (2/t)y' + (2/t^2)y = 0\), so \(a_1 = -2/t\), and Abel gives \(W(t) = W(1)\exp \int _1^t2/s\,ds = 1 \cdot t^2\).
Problem 4.7.2. Show that \(\sin t\) and \(t\) cannot both be solutions of an equation \(y'' + p(t)y' + q(t)y = 0\) with \(p, q\) continuous on \((-1, 1)\).
Show solution
Solution. Their Wronskian is \(W = \sin t \cdot 1 - t\cos t\), which is \(0\) at \(t = 0\) but not identically zero (\(W'(t) = t\sin t > 0\) for \(t \neq 0\) near \(0\), so \(W \neq 0\) for small \(t \neq 0\)). By Abel’s identity the Wronskian of two solutions is either never zero or always zero, so they cannot both be solutions.
Problem 4.7.3. Derive the variation of parameters formula for \(y'' + p(t)y' + q(t)y = g(t)\) from Theorem 4.3.1, and use it to solve \(y'' + y = \tan t\) on \((-\pi /2, \pi /2)\).
Show solution
Solution. With \(x = (y, y')\), \(\Phi = \begin {bmatrix} y_1 & y_2\\ y_1' & y_2'\end {bmatrix}\) and forcing \((0, g)^{\mathsf T}\), the condition \(\Phi c' = (0, g)^{\mathsf T}\) gives \(c_1' = -y_2g/W\), \(c_2' = y_1g/W\) by Cramer’s rule, and \(y_p = c_1y_1 + c_2y_2\), which is the stated formula. For \(\tan t\), with \(y_1 = \cos t\), \(y_2 = \sin t\), \(W = 1\): \[ c_1 = -\int \sin t\tan t\,dt = \sin t - \log |\sec t + \tan t|, \qquad c_2 = \int \sin t\,dt = -\cos t , \] so \[ y_p = \cos t\sin t - \cos t\log |\sec t + \tan t| - \sin t\cos t = -\cos t\log |\sec t + \tan t| . \]
Problem 4.7.4. Find a second solution of \((1 - t^2)y'' - 2ty' + 2y = 0\) on \((-1, 1)\), given \(y_1 = t\) (this is Legendre’s equation of order one).
Show solution
Solution. In normal form \(p = -2t/(1 - t^2)\), so \(e^{-\int p} = e^{-\log (1 - t^2)} = 1/(1 - t^2)\). Then \begin {align*} y_2 &= t\int \frac {dt}{t^2(1 - t^2)} = t\int \Bigl (\frac {1}{t^2} + \frac {1}{1 - t^2}\Bigr )dt\\ &= t\Bigl (-\frac 1t + \frac 12\log \frac {1 + t}{1 - t}\Bigr ) = \frac t2\log \frac {1 + t}{1 - t} - 1 , \end {align*}
the Legendre function of the second kind \(Q_1\), which is singular at \(t = \pm 1\).
Problem 4.7.5. Let \(x\) solve \(x' = A(t)x\) and \(y\) solve the adjoint system. Show that if \(y(T) = c\), then \(y(t)^{\mathsf T}x(t) = c^{\mathsf T}x(T)\) for all \(t\), and explain how one backward integration of the adjoint gives the sensitivity of \(c^{\mathsf T}x(T)\) to the initial state \(x(t_0)\).
Show solution
Solution. By Theorem 4.4.1, \(y^{\mathsf T}x\) is constant, so it equals its value \(c^{\mathsf T}x(T)\) at \(t = T\). In particular \(c^{\mathsf T}x(T) = y(t_0)^{\mathsf T}x(t_0)\), so the gradient of \(c^{\mathsf T}x(T)\) with respect to \(x(t_0)\) is \(y(t_0)\). One backward integration of the adjoint from \(y(T) = c\) gives all \(n\) partial derivatives, where finite differences would need \(n + 1\) forward integrations.
Problem 4.7.6. Show that for \(A(t) = a(t)A_0\) with \(A_0\) a constant matrix, \(\Phi (t, t_0) = \exp \bigl (A_0\int _{t_0}^ta(s)\,ds\bigr )\), and use this to solve \(x' = \frac 1t\begin {bmatrix} 0 & 1\\ 1 & 0\end {bmatrix}x\) on \(t > 0\).
Show solution
Solution. \(B(t) = \bigl (\int _{t_0}^ta\bigr )A_0\) is a scalar multiple of \(A_0\), which commutes with \(A(t) = a(t)A_0\), so Proposition 4.5.1 applies. Here \(\int _1^tds/s = \log t\), and since \(A_0^2 = I\), \(e^{\tau A_0} = \cosh \tau \,I + \sinh \tau \,A_0\). With \(\tau = \log t\), \(\cosh \tau = \frac 12(t + t^{-1})\) and \(\sinh \tau = \frac 12(t - t^{-1})\), so \(\Phi (t, 1) = \frac 12\begin {bmatrix} t + t^{-1} & t - t^{-1}\\ t - t^{-1} & t + t^{-1}\end {bmatrix}\). Its determinant is \(1\), as Liouville’s formula requires since \(\tr A = 0\).
Questions on this section
Stuck on something here? Ask below and it stays attached to this topic.