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.