Processing math: 100%

2.2 HW 1

  2.2.1 Problem 8, section 1.2
  2.2.2 Problem 17, section 1.2
  2.2.3 Problem 10, section 1.4
  2.2.4 Problem 18, section 1.4
  2.2.5 Problem 4, section 1.5
  2.2.6 Problem 6(a), section 1.5
  2.2.7 Key solution for HW 1
PDF (letter size)
PDF (legal size)

2.2.1 Problem 8, section 1.2

Solve\frac{dy}{dt}+\sqrt{1+t^{2}}y=0,y\left ( 0\right ) =\sqrt{5}

Solution

This is separable first order ODE. Therefore\begin{equation} \int \frac{dy}{y}=-\int \sqrt{1+t^{2}}dt\tag{1} \end{equation} The LHS becomes\begin{equation} \int \frac{dy}{y}=\ln \left \vert y\right \vert \tag{2} \end{equation} For the RHS of (1), the integral \int \sqrt{1+t^{2}}dt can be evaluated as follows. Let t=\sinh \left ( \theta \right ) . Hence \frac{dt}{d\theta }=\cosh \left ( \theta \right ) . Therefore\begin{align*} \int \sqrt{1+t^{2}}dt & =\int \sqrt{1+\sinh ^{2}\left ( \theta \right ) }\cosh \left ( \theta \right ) d\theta \\ & =\int \cosh ^{2}\left ( \theta \right ) d\theta \\ & =\int \frac{1}{2}\left ( 1+\cosh \left ( 2\theta \right ) \right ) d\theta \\ & =\frac{1}{2}\left ( \int d\theta +\int \cosh \left ( 2\theta \right ) d\theta \right ) \\ & =\frac{1}{2}\left ( \theta +\frac{\sinh \left ( 2\theta \right ) }{2}\right ) \\ & =\frac{1}{2}\theta +\frac{\sinh \left ( 2\theta \right ) }{4} \end{align*}

Since \sinh \left ( 2\theta \right ) =2\sinh \theta \cosh \theta , the above becomes \int \sqrt{1+t^{2}}dt=\frac{1}{2}\theta +\frac{\sinh \theta \cosh \theta }{2} Since \cosh ^{2}\left ( \theta \right ) -\sinh ^{2}\left ( \theta \right ) =1 then \cosh ^{2}\theta =1+\sinh ^{2}\left ( \theta \right ) and the above becomes \int \sqrt{1+t^{2}}dt=\frac{1}{2}\left ( \theta +\sinh \theta \sqrt{1+\sinh ^{2}\left ( \theta \right ) }\right ) But t=\sinh \left ( \theta \right ) and \theta =\operatorname{arcsinh}\left ( t\right ) . Therefore the above becomes\begin{equation} \int \sqrt{1+t^{2}}dt=\frac{1}{2}\left ( \operatorname{arcsinh}\left ( t\right ) +t\sqrt{1+t^{2}}\right ) \tag{3} \end{equation} Using (2,3) in (1) gives\begin{equation} \ln \left \vert y\right \vert =-\frac{1}{2}\left ( \operatorname{arcsinh}\left ( t\right ) +t\sqrt{1+t^{2}}\right ) +C\tag{4} \end{equation} Where C is arbitrary constant of integration. Writing \operatorname{arcsinh}\left ( t\right ) using known identity as \ln \left \vert t+\sqrt{1+t^{2}}\right \vert . And since \sqrt{1+t^{2}} is always larger than t, then the absolute sign is not needed.  Eq. (4) becomes\begin{align*} \ln \left \vert y\right \vert & =-\frac{1}{2}\left ( \ln \left ( t+\sqrt{1+t^{2}}\right ) +t\sqrt{1+t^{2}}\right ) +C\\ \left \vert y\right \vert & =e^{-\frac{1}{2}\left ( \ln \left ( t+\sqrt{1+t^{2}}\right ) +t\sqrt{1+t^{2}}\right ) }e^{C}\\ y & =C_{1}e^{-\frac{1}{2}\left ( \ln \left ( t+\sqrt{1+t^{2}}\right ) +t\sqrt{1+t^{2}}\right ) }\\ & =C_{1}e^{-\frac{1}{2}\ln \left ( t+\sqrt{1+t^{2}}\right ) }e^{t\sqrt{1+t^{2}}} \end{align*}

Therefore the general solution is y\left ( t\right ) =C_{1}\frac{e^{t\sqrt{1+t^{2}}}}{\left ( t+\sqrt{1+t^{2}}\right ) ^{\frac{1}{2}}} Now initial conditions are used to determine C_{1}.  From y\left ( 0\right ) =\sqrt{5} then the above gives \sqrt{5}=C_{1} Therefore the particular solution is y\left ( t\right ) =\sqrt{5}\frac{e^{t\sqrt{1+t^{2}}}}{\left ( t+\sqrt{1+t^{2}}\right ) ^{\frac{1}{2}}}

2.2.2 Problem 17, section 1.2

Find a continuous solution of the IVP y+y^{\prime }=g\left ( t\right ) ,y\left ( 0\right ) =0\ where g\left ( t\right ) =\left \{ \begin{array} [c]{ccc}2 & & 0\leq t\leq 1\\ 0 & & t>1 \end{array} \right . Solution

This is linear first order ODE. The integrating factor is \mu =e^{\int dt}=e^{t}. Hence the ODE becomes\begin{align*} \frac{d}{dt}\left ( y\mu \right ) & =\mu g\left ( t\right ) \\ \frac{d}{dt}\left ( ye^{t}\right ) & =e^{t}g\left ( t\right ) \end{align*}

Integrating gives\begin{equation} ye^{t}=\int e^{t}g\left ( t\right ) dt+C\tag{1} \end{equation} Breaking the problem into two phases, and solving the above for 0\leq t\leq 1 gives\begin{align*} ye^{t} & =\int 2e^{t}dt+C\\ & =2e^{t}+C\\ y\left ( t\right ) & =2+Ce^{-t} \end{align*}

Applying initial conditions gives 0=2+C, or C=-2 and the above becomes\begin{equation} y\left ( t\right ) =2-2e^{-t}\qquad 0\leq t\leq 1\tag{2} \end{equation} The above solution is valid for 0\leq t\leq 1.

To solve for t>1, initial conditions are first found for t=1. At t=1 the above gives y\left ( 1\right ) =2-\frac{2}{e} Hence for t>1, initial conditions are y\left ( 1\right ) =2-\frac{2}{e}. Now the second phase is solved. From (1) ye^{t}=\int e^{t}g\left ( t\right ) dt+C But now g\left ( t\right ) =0. The above simplifies to\begin{align} ye^{t} & =C\tag{3}\\ y & =Ce^{-t}\nonumber \end{align}

But at t=1,y=2-\frac{2}{e}. Therefore\begin{align*} 2-\frac{2}{e} & =Ce^{-1}\\ C & =2e-2\\ & =2\left ( e-1\right ) \end{align*}

Substituting the above C into (3) gives\begin{equation} y=2\left ( e-1\right ) e^{-t}\qquad t>1\tag{4} \end{equation} Using (2,4) the final solution is therefore y\left ( t\right ) =\left \{ \begin{array} [c]{ccc}2-2e^{-t} & & 0\leq t\leq 1\\ 2\left ( e-1\right ) e^{-t} & & t>1 \end{array} \right .

pict
Figure 2.1:Plot of the solution y(t)

2.2.3 Problem 10, section 1.4

Solve \cos y\frac{dy}{dt}=\frac{-t\sin y}{1+t^{2}},y\left ( 1\right ) =\frac{\pi }{2}

Solution

This is separable first order ODE \int \frac{\cos y}{\sin y}dy=-\int \frac{t}{1+t^{2}}dt But \int \frac{\cos y}{\sin y}dy=\int \frac{\frac{d}{dy}\sin y}{\sin y}dy=\ln \left \vert \sin \left ( y\right ) \right \vert and \int \frac{t}{1+t^{2}}dt=\frac{1}{2}\ln \left \vert 1+t^{2}\right \vert =\frac{1}{2}\ln \left ( 1+t^{2}\right ) since 1+t^{2} is positive. Hence the above becomes \ln \left \vert \sin \left ( y\right ) \right \vert =-\frac{1}{2}\ln \left ( 1+t^{2}\right ) +C Where C is the integration constant. Hence\begin{align*} \left \vert \sin \left ( y\right ) \right \vert & =e^{-\frac{1}{2}\ln \left ( 1+t^{2}\right ) +C}\\ & =e^{-\frac{1}{2}\ln \left ( 1+t^{2}\right ) }e^{C} \end{align*}

Therefore\begin{align} \sin \left ( y\right ) & =C_{1}e^{-\frac{1}{2}\ln \left ( 1+t^{2}\right ) }\nonumber \\ & =C_{1}\frac{1}{\sqrt{1+t^{2}}} \tag{1} \end{align}

From initial conditions y\left ( 1\right ) =\frac{\pi }{2} the above becomes\begin{align*} \sin \left ( \frac{\pi }{2}\right ) & =C_{1}\frac{1}{\sqrt{2}}\\ C_{1} & =\sqrt{2} \end{align*}

Hence (1) becomes\begin{align*} \sin \left ( y\right ) & =\sqrt{2}\frac{1}{\sqrt{1+t^{2}}}\\ y\left ( t\right ) & =\arcsin \left ( \frac{\sqrt{2}}{\sqrt{1+t^{2}}}\right ) \end{align*}

2.2.4 Problem 18, section 1.4

Solve \frac{dy}{dt}=\frac{t+y}{t-y}

Solution

Let u=\frac{y}{t} or y=ut. Hence \frac{dy}{dt}=u+t\frac{du}{dt}. Therefore the ODE becomes\begin{align*} u+t\frac{du}{dt} & =\frac{t+ut}{t-ut}\\ u+t\frac{du}{dt} & =\frac{t\left ( 1+u\right ) }{t\left ( 1-u\right ) }\\ t\frac{du}{dt} & =\frac{\left ( 1+u\right ) }{\left ( 1-u\right ) }-u\\ & =-\frac{u^{2}+1}{u-1}\\ & =\frac{1+u^{2}}{1-u} \end{align*}

This is now separable ODE. Therefore\begin{align} \frac{1-u}{1+u^{2}}\frac{du}{dt} & =\frac{1}{t}\nonumber \\ \int \frac{1-u}{1+u^{2}}du & =\int \frac{1}{t}dt \tag{1} \end{align}

But \begin{align*} \int \frac{1-u}{1+u^{2}}du & =\int \frac{1}{1+u^{2}}du-\int \frac{u}{1+u^{2}}du\\ & =\arctan \left ( u\right ) -\frac{1}{2}\ln \left \vert 1+u^{2}\right \vert \end{align*}

but 1+u^{2} is positive. Hence \int \frac{1-u}{1+u^{2}}du=\arctan \left ( u\right ) -\frac{1}{2}\ln \left ( 1+u^{2}\right ) And \int \frac{1}{t}dt=\ln \left \vert t\right \vert . Hence (1) becomes \arctan \left ( u\right ) -\frac{1}{2}\ln \left ( 1+u^{2}\right ) =\ln \left \vert t\right \vert +C But \frac{y}{t}, and the above becomes \arctan \left ( \frac{y}{t}\right ) -\frac{1}{2}\ln \left ( 1+\left ( \frac{y}{t}\right ) ^{2}\right ) =\ln \left \vert t\right \vert +C The above solution is implicit in y\left ( t\right ) .

2.2.5 Problem 4, section 1.5

Suppose that a population doubles its original size in 100 years, and triples it in 200 years. Show that this population cannot satisfy the Malthusian law of population growth.

Solution

In Malthusian law of population growth, the rate at which population changes is fixed in the model. It is given by a below \frac{dp}{dt}=ap\left ( t\right ) Where a is constant. But the problem says the population is doubled in first 100 years. So if p_{0} was initial population, then after 100 years the population now has become 2p_{0}. There one will expect that after another 100 years the population will double again to become 4p_{0}.

But the problem says that the population triples in 200 years, becoming 3p_{0} and not 4p_{0}. This shows that the rate of growth is not constant. Hence this do not satisfy Malthusian law of population growth.

2.2.6 Problem 6(a), section 1.5

A population grows according to the logistic law, with a limiting population of 5\times 10^{8} individuals. When the population is low it doubles every 40 minutes. What will the population be after two hours if initially it was (a) 10^{8} ?

Solution

In the logistic law, the population model is given by \frac{dp}{dt}=ap-bp^{2} Where p\left ( t\right ) is population at time t and a is the growth rate (constant) and b is the competition rate (also constant). In this model \lim _{t\rightarrow \infty }p\left ( t\right ) =\frac{a}{b} Therefore \begin{equation} \frac{a}{b}=5\times 10^{8}\tag{1} \end{equation} The problem says that a=100\% (per 40 minute) or a=1 (per 40 minute). Therefore a=\frac{1}{40} per minute. And p_{0}=10^{8}. Using the solution of this model, given in the textbook at page 30 as\begin{equation} p\left ( t\right ) =\frac{ap_{0}}{bp_{0}+\left ( a-bp_{0}\right ) e^{-a\left ( t-t_{0}\right ) }}\tag{3} \end{equation} And using t_{0}=0, then the population size at t is now be calculated.  From (1), b=\frac{\frac{1}{40}}{5\times 10^{8}}=\frac{1}{2}\times 10^{-10}=\allowbreak 5\times 10^{-11}. Eq. (3) now becomes\begin{align*} p\left ( t\right ) & =\frac{\frac{1}{40}\left ( 10^{8}\right ) }{\left ( 5\times 10^{-11}\right ) \left ( 10^{8}\right ) +\left ( \frac{1}{40}-\left ( 5\times 10^{-11}\right ) \left ( 10^{8}\right ) \right ) e^{-\frac{1}{40}t}}\\ & =\frac{\frac{1}{40}\left ( 10^{8}\right ) }{\left ( 5\times \frac{1}{1000}\right ) +\left ( \frac{1}{40}-5\times \frac{1}{1000}\right ) e^{-\frac{1}{40}t}} \end{align*}

For t=120 (minutes) the above becomes\begin{align*} p\left ( 120\right ) & =\frac{\frac{1}{40}\left ( 10^{8}\right ) }{\left ( 5\times \frac{1}{1000}\right ) +\left ( \frac{1}{40}-5\times \frac{1}{1000}\right ) e^{-\frac{1}{40}120}}\\ & =\frac{\frac{1}{40}\left ( 10^{8}\right ) }{\left ( 5\times \frac{1}{1000}\right ) +\left ( \frac{1}{40}-5\times \frac{1}{1000}\right ) e^{-3}}\\ & =4.1696\times 10^{8} \end{align*}

Hence \fbox{$p\left ( 120\right ) =4.1696\times 10^8$} The inflection point is \begin{align*} \frac{a}{2b} & =\frac{\frac{1}{40}}{\left ( 2\right ) \left ( 5\times 10^{-11}\right ) }\\ & =2.5\times 10^{8} \end{align*}

The following plot was generated to compare the population p\left ( t\right ) between case (a) and case (b). It shows that when starting with initial population of p_{0}=10^{8} which is case (a) and when starting with p_{0}=10^{9} which is case (b), both populations will eventually reach the limiting population of 5\times 10^{8}. The S curve shows up only when starting with population below the limiting population.

pict
Figure 2.2:Population p(t) change depends on p_0

2.2.7 Key solution for HW 1

PDF