22 problems found
Two identical particles of unit mass move under gravity in a medium for which the magnitude of the retarding force on a particle is \(k\) times its speed. The first particle is allowed to fall from rest at a point \(A\) whilst, at the same time, the second is projected upwards with speed \(u\) from a point \(B\) a positive distance \(d\) vertically above \(A\). Find their distance apart after a time \(t\) and show that this distance tends to the value \[ d+\frac{u}{k} \] as \(t\rightarrow\infty.\)
Solution: Both particles have equations of motion, \(\ddot{x} = -g-k\dot{x}\), so we can note that the distance between them has the equation of motion: \(\ddot{x} = -k \ddot{x} \Rightarrow x = Ae^{-kt} + B\) \begin{align*} && x(0) &= d \\ \Rightarrow && A+B &= d \\ && x'(0) &= u \\ \Rightarrow && -kA &= u \\ \Rightarrow && A &= -\frac{u}{k} \\ \Rightarrow && B &= d+\frac{u}{k} \\ \Rightarrow && x(t) &= -\frac{u}{k}e^{-kt} + d + \frac{u}{k} \to d + \frac{u}{k} \end{align*} as required.
In the manufacture of Grandma's Home Made Ice-cream, chemicals \(A\) and \(B\) pour at constant rates \(a\) and \(b-a\) litres per second (\(0 < a < b\)) into a mixing vat which mixes the chemicals rapidly and empties at a rate \(b\) litres per second into a second mixing vat. At time \(t=0\) the first vat contains \(K\) litres of chemical \(B\) only. Show that the volume \(V(t)\) (in litres) of the chemical \(A\) in the first vat is governed by the differential equation \[ \dot{V}(t)=-\frac{bV(t)}{K}+a, \] and that \[ V(t)=\frac{aK}{b}(1-\mathrm{e}^{-bt/K}) \] for \(t\geqslant0.\) The second vat also mixes chemicals rapidly and empties at the rate of \(b\) litres per second. If at time \(t=0\) it contains \(L\) litres of chemical \(C\) only (where \(L\neq K\)), how many litres of chemical \(A\) will it contain at a later time \(t\)?
Solution: The total volume in the first vat at time \(t\) is always \(K\), since \(b\) litres per second are coming in and \(b\) litres per second are going out. \begin{align*} &&\frac{\d V}{\d t} &= \underbrace{a}_{\text{incoming chemical }A} - \underbrace{b}_{\text{outgoing volume}} \cdot \underbrace{\frac{V(t)}{K}}_{\text{fraction of outgoing which is }A} \\ &&&= a - b \frac{V}{K} \\ \Rightarrow && \int \frac{1}{a-b\frac{V}{K}}\d V &= \int \d t \\ && - \frac{K}{b} \ln |a - b \frac{V}{K}| &= t +C\\ (t,V) = (0,0): && -\frac{K}{b} \ln a &= C \\ \Rightarrow && 1-\frac{b}{a} \frac{V}{K} &= e^{-bt/K} \\ \Rightarrow && V &= \frac{aK}{b} (1 - e^{-bt/K}) \end{align*} \begin{align*} &&\frac{\d W}{\d t} &= \underbrace{b}_{\text{incoming volume}} \cdot \underbrace{\frac{a}{b} (1 - e^{-bt/K})}_{\text{incoming fraction }A} - \underbrace{b}_{\text{outgoing volume}} \cdot \underbrace{ \frac{W(t)}{L}}_{\text{fraction of outgoing which is }A} \\ &&&= a (1 - e^{-bt/K}) - b \frac{W}{L} \\ \Rightarrow && \frac{\d W}{\d t} + \frac{b}{L} W &= a (1-e^{-bt/K}) \\ && \frac{\d}{\d t} \left ( e^{b/L t} W\right) &= ae^{b/L t}(1-e^{-bt/K}) \\ \Rightarrow && W &= e^{-bt/L} \left ( \frac{aL}{b}e^{b/Lt} - \frac{a}{\frac{b}{L} - \frac{b}{K}}e^{b/L t - b/K} \right) + Ce^{-bt/L} \\ &&&= \frac{aL}{b} \left (1 - \frac{K}{K-L}e^{-b/Kt} \right)+ Ce^{-bt/L} \\ (t,W) = (0,0): && 0 &= \frac{aL}{b} \frac{-L}{K-L} + C \\ \Rightarrow && C &= \frac{aL^2}{b(K-L)} \\ \Rightarrow && W &= \frac{aL}{b} \left (1 - \frac{K}{K-L} e^{-bt/K} + \frac{L}{K-L} e^{-bt/L} \right) \end{align*}
Let \(y,u,v,P\) and \(Q\) all be functions of \(x\). Show that the substitution \(y=uv\) in the differential equation \[ \frac{\mathrm{d}y}{\mathrm{d}x}+Py=Q \] leads to an equation for \(\dfrac{\mathrm{d}v}{\mathrm{d}x}\) in terms of \(x,Q\) and \(u\), provided that \(u\) satisfies a suitable first order differential equation. Hence or otherwise solve \[ \frac{\mathrm{d}y}{\mathrm{d}x}-\frac{2y}{x+1}=\left(x+1\right)^{\frac{5}{2}}, \] given that \(y(1)=0\). For what set of values of \(x\) is the solution valid?
Solution: Suppose \(y = uv\) then and suppose \(\frac{\d u}{\d x} + P u = 0\) then \begin{align*} && \frac{\d y}{\d x} + Py &= Q \\ && uv' + u'v + Puv &= Q \\ && uv' &= Q \\ && \frac{\d v}{\d x} &= \frac{Q}{u} \end{align*} Consider \begin{align*} && 0 &= \frac{\d u}{\d x} - \frac{2u}{x+1} \\ \Rightarrow && \ln u &= 2\ln (1 + x) + C \\ \Rightarrow && u &= A(1+x)^2 \end{align*} and \begin{align*} && \frac{\d v}{\d x} &= \frac1{A}(x+1)^{\frac12} \\ \Rightarrow && v &= \frac2{3A}(x+1)^{\frac32} + k \\ \Rightarrow && y &= \frac23(x+1)^\frac72 + k(x+1)^2 \\ && 0 &= y(1) \\ &&&= \frac23 2^{7/2}+k2^2 \\ \Rightarrow && k &= -\frac{2^{5/2}}{3} \\ \Rightarrow && y &= \frac23 (x+1)^{7/2} - \frac{2^{5/2}}{3}(x+1)^2 \end{align*}
Let \(P,Q\) and \(R\) be functions of \(x\). Prove that, for any function \(y\) of \(x\), the function \[ Py''+Qy'+Ry \] can be written in the form \(\dfrac{\mathrm{d}}{\mathrm{d}x}(py'+qy),\) where \(p\) and \(q\) are functions of \(x\), if and only if \(P''-Q'+R=0.\) Solve the differential equation \[ (x-x^{4})y''+(1-7x^{3})y'-9x^{2}y=(x^{3}+3x^2)\mathrm{e}^{x}, \] given that when \(x=2,y=2\mathrm{e}^{2}\) and \(y'=0.\)
Solution: Suppose \(Py'' + Qy' + Ry = \frac{\d}{\d x}(p y' + qy)\), then \begin{align*} Py'' + Qy' + Ry &= \frac{\d}{\d x}(p y' + qy) \\ &= py'' + p'y' + qy' + q' y \\ &= py'' + (p'+q)y' + q' y \end{align*} Therefore \(P = p, Q = p'+q, R = q'\), Therefore \(q = Q-P'\) and \(R = Q'-P''\) or \(P'' -Q'+R = 0\). \((\Rightarrow)\) Suppose it can be written in that form, then the logic we have applied shows that equation is true. \((\Leftarrow)\) Suppose \(P''-Q'+R = 0\), then letting \(p = P, q = Q-P'\) we find functions of the form which will be expressed correctly. Notice that if \(P = x-x^4, Q = (1-7x^3), R = -9x^2\) then: \begin{align*} P'' - Q' + R &= -12x^2+21x^2-9x^2 \\ &= 0 \end{align*} Therefore we can write our second order ODE as: \begin{align*} && (x^{3}+3x)\mathrm{e}^{x} &= \frac{\d}{\d x} \left ( (x-x^4) y' +(1-7x^3-(1-4x^3))y \right) \\ &&&= \frac{\d}{\d x} \left ( (x-x^4) y' -3x^3y \right) \end{align*} Suppose \(z = (x-x^4)y' - 3x^2y\), then \(z = (2-2^4) \cdot 0 - 3 \cdot 2^2 \cdot 2e^2 = -24e^2\) when \(x = 2\). and we have: \begin{align*} && \frac{\d z}{\d x} &= (x^3+3x^2)e^x \\ \Rightarrow && z &= \int (x^3+3x^2)e^x \d x \\ &&&= x^3 e^x+c \\ \Rightarrow && -48e^2 &= e^2(8) + c \\ \Rightarrow && c &= -56e^2 \\ \Rightarrow && z &= e^x(x^3)-56e^2 \\ \end{align*} So our differential equation is: \begin{align*} && (x-x^4)y'-3x^3 y &= x^3e^x -5 6e^2 \\ \Rightarrow && (1-x^3)y' -3x^2y &= x^2 e^x - \frac{6e^2}{x} \\ \Rightarrow && \frac{\d }{\d x} \left ( (1-x^3)y \right) &= x^2e^x - \frac{56e^2}{x} \\ \Rightarrow && (1-x^3)y &= (x^2-2x+2)e^x - 56e^2 \ln x + k \\ \underbrace{\Rightarrow}_{x=2} && (1-2^3)2e^2 &= (2^2-2\cdot2 + 2)e^2 -56e^2 \ln 2 + k \\ \Rightarrow && k &= -16e^2+56 \ln 2 \cdot e^2 \\ \Rightarrow && y &= \frac{(x^2-2x+2)e^x - 56e^2 \ln x -16e^2+56 \ln 2 \cdot e^2}{(1-x^3)} \end{align*}
The matrix \(\mathbf{F}\) is defined by \[ \mathbf{F}=\mathbf{I}+\sum_{n=1}^{\infty}\frac{1}{n!}t^{n}\mathbf{A}^{n}, \] where $\mathbf{A}=\begin{pmatrix}-3 & -1\\ 8 & 3 \end{pmatrix} \( and \) t \( is a variable scalar. Evaluate \)\mathbf{A}^{2}$, and show that \[ \mathbf{F}=\mathbf{I}\cosh t+\mathbf{A}\sinh t. \] Show also that \(\mathbf{F}^{-1}=\mathbf{I}\cosh t-\mathbf{A}\sinh t\), and that \(\dfrac{\mathrm{d}\mathbf{F}}{\mathrm{d}t}=\mathbf{FA}\). The vector $\mathbf{r}=\begin{pmatrix}x(t)\\ y(t) \end{pmatrix}$ satisfies the differential equation \[ \frac{\mathrm{d}\mathbf{r}}{\mathrm{d}t}+\mathbf{A}\mathbf{r}=\mathbf{0}, \] with \(x=\alpha\) and \(y=\beta\) at \(t=0.\) Solve this equation by means of a suitable matrix integrating factor, and hence show that \begin{alignat*}{1} x(t) & =\alpha\cosh t+(3\alpha+\beta)\sinh t\\ y(t) & =\beta\cosh t-(8\alpha+3\beta)\sinh t. \end{alignat*}
Solution: \begin{align*} \begin{pmatrix} -3 & -1 \\ 8 & 3 \end{pmatrix}^2 &= \begin{pmatrix} 9-8 & 3-3 \\ -24+24 & -8+9 \end{pmatrix} \\ &= \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} \\ &= \mathbf{I} \end{align*} Therefore: \begin{align*} \mathbf{F} &= \mathbf{I}+\sum_{n=1}^{\infty}\frac{1}{n!}t^{n}\mathbf{A}^{n} \\ &= \mathbf{I} + \sum_{n=1}^{\infty} \frac{1}{(2n)!}t^{2n} \mathbf{I} + \sum_{n=0}^{\infty} \frac{1}{(2n+1)!}t^{2n+1} \mathbf{A} \\ &= \cosh t \mathbf{I} + \sinh t \mathbf{A} \end{align*} Notice that \begin{align*} \mathbf{F} (\mathbf{I}\cosh t-\mathbf{A}\sinh t) &= (\mathbf{I}\cosh t+\mathbf{A}\sinh t)(\mathbf{I}\cosh t-\mathbf{A}\sinh t) \\ &= \mathbf{I}^2 \cosh^2 t+\mathbf{A}(\sinh t \cosh t - \cosh t \sinh t) - \mathbf{A}^2\sinh^2 t \\ &= \mathbf{I} \cosh^2 t - \mathbf{I} \sinh^2 t \\ &= \mathbf{I} \end{align*} Therefore \(\mathbf{F}^{-1} = \mathbf{I}\cosh t-\mathbf{A}\sinh t\) \begin{align*} \frac{\d \mathbf{F}}{\d t} &= \frac{\d }{\d t} \left [ \mathbf{I}+\sum_{n=1}^{\infty}\frac{1}{n!}t^{n}\mathbf{A}^{n}\right] \\ &= \sum_{n=1}^{\infty} \frac{1}{(n-1)!}t^{n-1} \mathbf{A}^n \\ &= \left ( \sum_{n=1}^{\infty} \frac{1}{(n-1)!}t^{n-1} \mathbf{A}^{n-1} \right) \mathbf{A} \\ &= \mathbf{F}\mathbf{A} \end{align*} \begin{align*} && \frac{\mathrm{d}\mathbf{r}}{\mathrm{d}t}+\mathbf{A}\mathbf{r}&=\mathbf{0} \\ \Rightarrow && \mathbf{F} \frac{\d \mathbf{r}}{\d t} + \mathbf{FAr} &= \mathbf{0} \\ && \frac{\d }{\d t} \left ( \mathbf{F} \mathbf{r}\right) &= 0 \\ \Rightarrow && \mathbf{Fr} &= \mathbf{c} \\ \Rightarrow && \mathbf{r} &= \mathbf{F}^{-1}\mathbf{c} \\ &&&= ( \mathbf{I}\cosh t-\mathbf{A}\sinh t)\mathbf{c} \\ t = 0: && \binom{\alpha}{\beta} &= \mathbf{c} \\ \Rightarrow && \mathbf{r} &= ( \mathbf{I}\cosh t-\mathbf{A}\sinh t)\binom{\alpha}{\beta} \\ &&&= \binom{\alpha \cosh t}{\beta \cosh t} - \binom{-3\alpha-\beta}{8\alpha + 3\beta}\sinh t \\ &&&= \binom{\alpha \cosh t + (3\alpha + \beta) \sinh t}{\beta \cosh t -(8\alpha + 3\beta)\sinh t} \end{align*} as required
A goalkeeper stands on the goal-line and kicks the football directly into the wind, at an angle \(\alpha\) to the horizontal. The ball has mass \(m\) and is kicked with velocity \(\mathbf{v}_{0}.\) The wind blows horizontally with constant velocity \(\mathbf{w}\) and the air resistance on the ball is \(mk\) times its velocity relative to the wind velocity, where \(k\) is a positive constant. Show that the equation of motion of the ball can be written in the form \[ \frac{\mathrm{d}\mathbf{v}}{\mathrm{d}t}+k\mathbf{v}=\mathbf{g}+k\mathbf{w}, \] where \(\mathbf{v}\) is the ball's velocity relative to the ground, and \(\mathbf{g}\) is the acceleration due to gravity. By writing down horizontal and vertical equations of motion for the ball, or otherwise, find its position at time \(t\) after it was kicked. On the assumption that the goalkeeper moves out of the way, show that if \(\tan\alpha=\left|\mathbf{g}\right|/(k\left|\mathbf{w}\right|),\) then the goalkeeper scores an own goal.
Solution: Applying \(\mathbf{F} = m\mathbf{a} = m \frac{\d \mathbf{v}}{dt}\) we have: \begin{align*} && m \frac{\d \mathbf{v}}{d t} &= m\mathbf{g} - mk(\mathbf{v} - \mathbf{w}) \\ \Rightarrow && \frac{\d \mathbf{v}}{d t} +k \mathbf{v} &= \mathbf{g} + k \mathbf{w} \\ \\ \Rightarrow && e^{k t} \l \frac{\d \mathbf{v}}{d t} +k \mathbf{v} \r &= e^{kt} ( \mathbf{g} + k \mathbf{w}) \\ \Rightarrow && \frac{\d}{\d t} \l e^{kt} \mathbf{v} \r &= e^{kt}( \mathbf{g} + k \mathbf{w}) \\ \Rightarrow && e^{kt} \mathbf{v} &= \frac{1}ke^{kt}( \mathbf{g} + k \mathbf{w}) + c \\ \Rightarrow && \mathbf{v}_0 &= \frac{1}{k} ( \mathbf{g} + k \mathbf{w})+c \\ \Rightarrow && \mathbf{v} &= e^{-kt} \l \mathbf{v_0} - \frac{1}{k}\mathbf{g} - \mathbf{w} \r + \frac{1}{k} \mathbf{g} + \mathbf{w} \\ \Rightarrow && \mathbf{x} &= -\frac{1}{k}e^{-kt} \l \mathbf{v_0} - \frac{1}{k}\mathbf{g} - \mathbf{w} \r + \frac{1}{k} \mathbf{g}t + \mathbf{w}t+C \\ \Rightarrow && \mathbf{0} &= -\frac{1}{k} \l \mathbf{v_0} - \frac{1}{k}\mathbf{g} - \mathbf{w} \r + C \\ \Rightarrow && \mathbf{x} &= \frac1{k}\l 1- e^{-kt} \r\l \mathbf{v_0} - \frac{1}{k}\mathbf{g} - \mathbf{w} \r + \frac{1}{k} \mathbf{g}t + \mathbf{w}t \end{align*} Position at time \(t\) is: \begin{align*} && x_x &= \frac1{k} ( 1-e^{-kt})(u_x - w)+wt \\ && x_y &= \frac1{k} ( 1-e^{-kt})(u_x \frac{g}{kw} - \frac{g}{k})+\frac{1}{k}gt \\ &&&= \frac{g}{kw} \left ( ( 1-e^{-kt})(u_x - w)+wt \right) \\ &&&= \frac{g}{kw} x_x \end{align*} Therefore if \(x_x\) is ever \(0\) then \(x_y\) will also be zero. But the ball must eventually hit the ground, and when it does, it will be in the process of scoring an own goal.
By substituting \(y(x)=xv(x)\) in the differential equation \[ x^{3}\frac{\mathrm{d}v}{\mathrm{d}x}+x^{2}v=\frac{1+x^{2}v^{2}}{\left(1+x^{2}\right)v}, \] or otherwise, find the solution \(v(x)\) that satisfies \(v=1\) when \(x=1\). What value does this solution approach when \(x\) becomes large?
Solution: Let \(y = xv\) then \(y' = v + xv'\) and so \(x^2y' = x^2v + x^3v'\) Our differential equation is now: \begin{align*} && x^2 y' &= \frac{1+y^2}{(1+x^2)\frac{y}{x}} \\ \Rightarrow && xy' &= \frac{(1+y^2)}{(1+x^2)y} \\ \Rightarrow && \frac{y}{1+y^2} \frac{\d y}{\d x} &= \frac{1}{x(1+x^2)} \\ \Rightarrow && \frac{y}{1+y^2} \frac{\d y}{\d x} &= \frac{1}{x} - \frac{x}{1+x^2} \\ \Rightarrow && \frac12 \ln(1+y^2) &= \ln x - \frac12 \ln(1+x^2) + C \\ \Rightarrow && \frac12 \ln (1 + y^2) &= \frac12 \ln \l \frac{x^2}{1+x^2}\r + C \\ \Rightarrow && 1+y^2 &= \frac{Dx^2}{1+x^2} \\ \Rightarrow && D &= 4 \quad \quad: (x = 1, v = 1, y = 1) \\ \Rightarrow && 1 + x^2v^2&= \frac{4x^2}{1+x^2}\\ \Rightarrow && v^2 &= \frac{3x^2-1}{x^2(1+x^2)} \\ \Rightarrow && v &= \sqrt{\frac{3x^2-1}{x^2(1+x^2)}} \\ \end{align*} As \(x \to \infty\), \(v \to 0\)