Processing math: 100%

2.2 HW2

  2.2.1 Section 2.2 problem 1
  2.2.2 Section 2.2 problem 2
  2.2.3 Section 2.2 problem 3
  2.2.4 Section 2.2 problem 4
  2.2.5 Section 2.2 problem 5
  2.2.6 Section 2.2 problem 6
  2.2.7 Section 2.2 problem 7
  2.2.8 Section 2.2 problem 8
  2.2.9 Section 2.3 problem 1
  2.2.10 Section 2.3 problem 2
  2.2.11 Section 2.3 problem 3
  2.2.12 Section 2.3 problem 4
  2.2.13 Section 2.4 problem 1
  2.2.14 Section 2.4 problem 2
  2.2.15 Section 2.4 problem 3
  2.2.16 Section 2.4 problem 4
  2.2.17 Section 2.4 problem 5
  2.2.18 Section 2.4 problem 6
PDF (letter size)
PDF (legal size)

2.2.1 Section 2.2 problem 1

Solve y^{\prime }=\frac{x^{2}}{y}

This is first order non-linear ODE. In the form y^{\prime }=f\left ( x,y\right ) . The function f\left ( x,y\right ) is continuous everywhere except at the line y=0.  Now the ODE is solved by separation\begin{align*} y\frac{dy}{dx} & =x^{2}\\ ydy & =x^{2}dx \end{align*}

Integrating\begin{align*} \int ydy & =\int x^{2}dx\\ \frac{y^{2}}{2} & =\frac{x^{3}}{3}+c \end{align*}

Since initial conditions is not gives, the solution is left in implicit form (as mentioned in discussion class, Thursday Sept. 29, 2016)

\fbox{$y^2=\frac{2}{3}x^3+c_0\qquad y\neq 0$}

2.2.2 Section 2.2 problem 2

Solve y^{\prime }=\frac{x^{2}}{y\left ( 1+x^{3}\right ) }

This is first order non-linear ODE. In the form y^{\prime }=f\left ( x,y\right ) . The function f\left ( x,y\right ) is continuous everywhere except at line y=0 and at line x=-1. Now the ODE is solved by separation\begin{align*} y\frac{dy}{dx} & =\frac{x^{2}}{\left ( 1+x^{3}\right ) }\\ ydy & =\frac{x^{2}}{\left ( 1+x^{3}\right ) }dx \end{align*}

Integrating \int ydy=\int \frac{x^{2}}{\left ( 1+x^{3}\right ) }dx

To integrate \int \frac{x^{2}}{\left ( 1+x^{3}\right ) }dx let u=1+x^{3}, hence \frac{du}{dx}=3x^{2}. Therefore the integral becomes \int \frac{x^{2}}{u}\frac{du}{3x^{2}}=\frac{1}{3}\int \frac{du}{u}=\frac{1}{3}\ln \left \vert u\right \vert =\frac{1}{3}\ln \left \vert 1+x^{3}\right \vert . Hence the above becomes\begin{align*} \frac{y^{2}}{2} & =\frac{1}{3}\ln \left \vert 1+x^{3}\right \vert +c\\ y^{2} & =\frac{2}{3}\ln \left \vert 1+x^{3}\right \vert +c_{1} \end{align*}

Since initial condition is not gives, the solution is left in implicit form

\fbox{$y^2=\frac{2}{3}\ln \left \vert 1+x^3\right \vert +c_1\qquad y\neq 0,x\neq -1$}

2.2.3 Section 2.2 problem 3

Solve y^{\prime }=-y^{2}\sin x

This is first order non-linear ODE. In the form y^{\prime }=f\left ( x,y\right ) . The function f\left ( x,y\right ) is continuous everywhere and \frac{\partial f}{\partial y}=-2y\sin \left ( x\right ) is also continuous everywhere but unbounded at y=-\infty . This is separable, assuming y\neq 0 and dividing by y^{2} the ODE becomes \begin{align*} \frac{1}{y^{2}}\frac{dy}{dx} & =-\sin \left ( x\right ) \\ \frac{dy}{y^{2}} & =-\sin \left ({}\right ) xdx \end{align*}

Integrating\begin{align*} \int \frac{dy}{y^{2}} & =-\int \sin \left ( x\right ) dx\\ -\frac{1}{y} & =\cos \left ( x\right ) +c\\ \frac{1}{y} & =-\cos \left ( x\right ) +c_{1} \end{align*}

Therefore the solution is \fbox{$y\left ( x\right ) =\frac{1}{c_1-\cos \left ( x\right ) }\qquad y\neq 0$}

The reason for y\neq 0 was the assumption to divide by y^{2} above. Another solution is \fbox{$y\left ( x\right ) =0$}

2.2.4 Section 2.2 problem 4

Solve y^{\prime }=\frac{3x^{2}-1}{3+2y}

This is first order non-linear ODE. In the form y^{\prime }=f\left ( x,y\right ) . The function f\left ( x,y\right ) is continuous everywhere except at 3+2y=0 or y=-\frac{3}{2}. Now the ODE is solved by separation.\begin{align*} \left ( 3+2y\right ) \frac{dy}{dx} & =3x^{2}-1\\ \left ( 3+2y\right ) dy & =\left ( 3x^{2}-1\right ) dx \end{align*}

Integrating\begin{align*} \int \left ( 3+2y\right ) dy & =\int \left ( 3x^{2}-1\right ) dx\\ y^{2}+3y & =x^{3}-x+c \end{align*}

Complete the square

y^{2}+3y+\left ( \frac{3}{2}\right ) ^{2}=x^{3}-x+c+\left ( \frac{3}{2}\right ) ^{2}

Since initial condition is not gives, the solution is left in implicit form.

\fbox{$\left ( y+\frac{3}{2}\right ) ^2=x^3-x+c_0\qquad y\neq -\frac{3}{2}$}

2.2.5 Section 2.2 problem 5

y^{\prime }=\cos ^{2}\left ( x\right ) \cos ^{2}\left ( 2y\right )

This is first order non-linear ODE. In the form y^{\prime }=f\left ( x,y\right ) . The function f\left ( x,y\right ) is continuous. \begin{align*} \frac{\partial f}{\partial y} & =\cos ^{2}\left ( x\right ) 2\cos \left ( 2y\right ) \left ( -2\sin \left ( 2y\right ) \right ) \\ & =-4\cos ^{2}\left ( x\right ) \cos \left ( 2y\right ) \sin \left ( 2y\right ) \end{align*}

Which is continuous everywhere and bounded. Hence a solution exist and is unique. Now the ODE is solved by separation.

Case \cos ^{2}\left ( 2y\right ) \neq 0

To divide by \cos ^{2}\left ( 2y\right ) , then for \cos ^{2}\left ( 2y\right ) \neq 0 or \cos \left ( 2y\right ) \neq 0 or 2y\neq \left ( n+\frac{1}{2}\right ) \pi \ or y\neq \left ( n+\frac{1}{2}\right ) \frac{\pi }{4} for all integers.\begin{align} \frac{1}{\cos ^{2}\left ( 2y\right ) }\frac{dy}{dx} & =\cos ^{2}\left ( x\right ) \nonumber \\ \int \frac{dy}{\cos ^{2}\left ( 2y\right ) } & =\int \cos ^{2}\left ( x\right ) dx\tag{1} \end{align}

Now \int \frac{dy}{\cos ^{2}\left ( 2y\right ) }=\frac{1}{2}\tan \left ( 2y\right ) and \begin{align*} \int \cos ^{2}\left ( x\right ) dx & =\int \frac{1+\cos \left ( 2x\right ) }{2}dx\\ & =\int \left ( \frac{1}{2}+\frac{\cos \left ( 2x\right ) }{2}\right ) dx\\ & =\frac{1}{2}x+\frac{1}{2}\frac{\sin \left ( 2x\right ) }{2}+c_{1}\\ & =\frac{x}{2}+\frac{\sin \left ( 2x\right ) }{4}+c_{1} \end{align*}

Hence (1) becomes\begin{align*} \frac{1}{2}\tan \left ( 2y\right ) & =\frac{x}{2}+\frac{\sin \left ( 2x\right ) }{4}+c_{1}\\ \tan \left ( 2y\right ) & =x+\frac{1}{2}\sin \left ( 2x\right ) +c \end{align*}

Since initial condition is not gives, the solution is left in implicit form.

Case \cos ^{2}\left ( 2y\right ) =0

This is when \cos \left ( 2y\right ) =0 or 2y=\left ( n+\frac{1}{2}\right ) \pi or y=\left ( n+\frac{1}{2}\right ) \frac{\pi }{2}\ for all integers. In this case the solution is y=\left ( n+\frac{1}{2}\right ) \frac{\pi }{2}

Summary of solution y\left ( x\right ) \left \{ \begin{array} [c]{lll}\tan \left ( 2y\right ) =x+\frac{1}{2}\sin \left ( 2x\right ) +c & & \cos ^{2}\left ( 2y\right ) \neq 0\\ \left ( n+\frac{1}{2}\right ) \frac{\pi }{2} & & \cos ^{2}\left ( 2y\right ) =0 \end{array} \right .

2.2.6 Section 2.2 problem 6

Solve xy^{\prime }=\left ( 1-y^{2}\right ) ^{\frac{1}{2}}

This is nonlinear first order of the form y^{\prime }=f\left ( x,y\right ) where f\left ( x,y\right ) =\frac{\left ( 1-y^{2}\right ) ^{\frac{1}{2}}}{x}. This is continuous everywhere except at x=0. ODE is solved by separation.

Case 1-y^{2}\neq 0

Or y^{2}\neq 1 or y\neq \pm 1, then dividing by \left ( 1-y^{2}\right ) ^{\frac{1}{2}} and integrating\begin{align*} \int \frac{dy}{\left ( 1-y^{2}\right ) ^{\frac{1}{2}}} & =\int \frac{dx}{x}\\ \arcsin \left ( y\right ) & =\ln \left \vert x\right \vert +c\\ y\left ( x\right ) & =\sin \left ( \ln \left \vert x\right \vert +c\right ) \end{align*}

Hence the solution is y\left ( x\right ) =\sin \left ( \ln \left \vert x\right \vert +c\right ) \qquad y\neq \pm 1,x\neq 0

Case 1-y^{2}=0

Then y\left ( x\right ) =\pm 1

Summary of solutions

\left \{ \begin{array} [c]{lll}y\left ( x\right ) =\sin \left ( \ln \left \vert x\right \vert +c\right ) & & y\neq \pm 1,x\neq 0\\ y\left ( x\right ) =\pm 1 & & x\neq 0 \end{array} \right .

2.2.7 Section 2.2 problem 7

Solve \frac{dy}{dx}=\frac{x-e^{-x}}{y+e^{y}}

This is non-linear first order ODE of the form y^{\prime }=f\left ( x,y\right ) . The function f\left ( x,y\right ) is continuous everywhere except at y which is the solution of e^{y}+y=0. Using a computer, this is y_{c}=-0.567143\cdots . The ODE is solved by separation\begin{align*} \int \left ( y+e^{y}\right ) dy & =\int \left ( x-e^{-x}\right ) dx\\ \frac{y^{2}}{2}+e^{y} & =\frac{x^{2}}{2}+e^{-x}+c \end{align*}

Hence the solution is given by y^{2}+2e^{y}-x^{2}-2e^{-x}=c_{1}\qquad y\neq y_{c}

2.2.8 Section 2.2 problem 8

Solve \frac{dy}{dx}=\frac{x^{2}}{1+y^{2}}

This is non-linear first order ODE of the form y^{\prime }=f\left ( x,y\right ) where f\left ( x,y\right ) is continuous everywhere except at y=\pm 1. The ODE is solved by separation\begin{align*} \left ( 1+y^{2}\right ) dy & =x^{2}dx\\ y+\frac{y^{3}}{3} & =\frac{x^{3}}{3}+c_{1} \end{align*}

Hence the solution is given by y^{3}+3y-x^{3}=c\qquad y=\pm 1

Since initial condition is not gives, the solution is left in implicit form.

2.2.9 Section 2.3 problem 1

pict

To reduce confusion, let x be the substance which causes the concentration in the die. Let Q\left ( t\right ) be the mass (normally called the amount, but saying mass is more clear than saying amount) of x at time t. Hence Q\left ( 0\right ) =200g since initial concentration was 1[g/L] and the volume is 200[L].

The goal is to find an ODE that describes how Q\left ( t\right ) changes in time. That is, how the mass of x in the tank changes in time. Knowing the mass of x at any time in the tank, gives the concentration also, since the tank volume is fixed at 200[L]. So the concentration can always be found using \frac{Q\left ( t\right ) }{200}. Using

\begin{equation} \frac{dQ}{dt}=R_{in}-R_{out}\tag{1} \end{equation}

Where R_{in} is rate of x moving into the tank, i.e. how many grams of x is being poured in per minute, which is zero, since fresh water is moving in. R_{out} is rate of x moving out, i.e. how many grams of x is leaving the tank per minute. This is found as follows\begin{align*} R_{out} & =\frac{Q\left ( t\right ) }{200}\frac{\left [ \text{gram}\right ] }{\left [ \text{L}\right ] }\times 2\left [ \frac{\text{L}}{\text{min}}\right ] \\ & =\frac{2}{200}Q\left ( t\right ) \frac{\left [ \text{gram}\right ] }{\left [ \text{min}\right ] } \end{align*}

Hence (1) becomes\begin{align*} \frac{dQ}{dt} & =0-100Q\left ( t\right ) \\ & =-100Q\left ( t\right ) \end{align*}

Solving the ODE, for Q\left ( t\right ) \neq 0\begin{align*} \frac{dQ}{Q} & =-100dt\\ \ln \left \vert Q\right \vert & =-100t+c \end{align*}

Since Q represent mass, it can not be negative, then there is no need to use \left \vert Q\right \vert .

\begin{align*} \ln Q & =-100t+c\\ Q\left ( t\right ) & =Ae^{-100t} \end{align*}

At t=0,Q\left ( 0\right ) =200[g], hence A=200 from the above. The solution becomes Q\left ( t\right ) =200e^{-100t}

Since initial Q was 200[g] then 1\% of that is 2. Solving for time gives\begin{align*} 2 & =200e^{-100t_{0}}\\ 0.01 & =e^{-100t_{0}}\\ \ln \left ( 0.01\right ) & =-100t_{0} \end{align*}

Solving on the computer gives \fbox{$t_0=460.517$[$\min $]}

Hence it takes 460.517 minutes for the mass of x to reach 1\% of its original amount of 200 gram. This is also the same amount of time for the concentration of x to reach 1\% of its original amount of 1\left [ g/L\right ] . It is easier to work with mass in the ODE, and then convert to concentration when needed.

2.2.10 Section 2.3 problem 2

pict

Let y\left ( t\right ) be the mass of salt at time t in the tank in grams. Hence y\left ( 0\right ) =0 since tank initially contains pure water. The goal is to find an ODE that describes how y\left ( t\right ) changes in time. That is, how the mass of salt in the tank changes in time. Using

\begin{equation} \frac{dy}{dt}=R_{in}-R_{out} \tag{1} \end{equation}

Where R_{in} is rate of mass of salt moving into the tank, i.e. how many grams of salt is being poured in per minute, which is \begin{align*} R_{in} & =\gamma \left [ \frac{\text{gram}}{\text{L}}\right ] \times 2\left [ \frac{\text{L}}{\min }\right ] \\ & =2\gamma \left [ \frac{\text{gram}}{\min }\right ] \end{align*}

And R_{out} is rate of salt moving out, i.e. how many grams of salt is leaving the tank per minute. This is found as follows\begin{align*} R_{out} & =\frac{y\left ( t\right ) }{120}\frac{\left [ \text{gram}\right ] }{\left [ \text{L}\right ] }\times 2\left [ \frac{\text{L}}{\min }\right ] \\ & =\frac{1}{60}y\left ( t\right ) \left [ \frac{\text{gram}}{\min }\right ] \end{align*}

Hence (1) becomes \frac{dy\left ( t\right ) }{dt}=2\gamma -\frac{1}{60}y\left ( t\right )

With y\left ( 0\right ) =0. The ODE is linear and first order, of the form y^{\prime }+p\left ( t\right ) y=g\left ( t\right ) with p\left ( t\right ) =\frac{1}{60} and g\left ( t\right ) =2\gamma . Since both p\left ( t\right ) ,g\left ( t\right ) are continuous then a solution exist and is unique. y^{\prime }+\frac{1}{60}y=2\gamma
Integrating factor is e^{\int \frac{1}{60}dt}=e^{\frac{1}{60}t}, therefore \frac{d}{dt}\left ( e^{\frac{1}{60}t}y\right ) =2\gamma e^{\frac{1}{60}t}
Integrating\begin{align*} e^{\frac{1}{60}t}y & =2\gamma \int e^{\frac{1}{60}t}dt\\ & =2\gamma \frac{e^{\frac{1}{60}t}}{\frac{1}{60}}+c\\ & =120\gamma e^{\frac{1}{60}t}+c \end{align*}

Hence y\left ( t\right ) =120\gamma +c^{\frac{-t}{60}}

In the above, y\left ( t\right ) is the mass of salt in grams in the tank at time t. Hence the concentration of salt in the tank at time t can always be found by dividing y\left ( t\right ) by the volume of the tank. In the limit, as t\rightarrow \infty then from above \lim _{t\rightarrow \infty }y\left ( t\right ) =120\gamma

2.2.11 Section 2.3 problem 3

pict

This problem is solved in two stages. The first ODE is used to find what the amount of salt in the tank will be after 10 minutes. Then a new ODE is set up, with this value as its initial conditions, in order to find the amount of salt in the tank after an additional 10 minutes.

First 10 minutes

Let y_{1}\left ( t\right ) be the mass of salt at time t in the tank in lbs. Hence y_{1}\left ( 0\right ) =0 since tank initially contains pure water. The goal is to find an ODE that describes how y_{1}\left ( t\right ) changes in time. That is, how the mass of salt in the tank changes in time. Using

\begin{equation} \frac{dy_{1}}{dt}=R_{in}-R_{out}\tag{1} \end{equation}

Where R_{in} is rate of mass of salt moving into the tank, i.e. how many lbs of salt is being poured in per minute, which is \begin{align*} R_{in} & =\frac{1}{2}\left [ \frac{\text{lb}}{\text{gallon}}\right ] \times 2\left [ \frac{\text{gallon}}{\min }\right ] \\ & =1\left [ \frac{\text{lb}}{\min }\right ] \end{align*}

And R_{out} is rate of salt moving out, i.e. how many grams of salt is leaving the tank per minute. This is found as follows\begin{align*} R_{out} & =\frac{y_{1}\left ( t\right ) }{100}\frac{\left [ \text{lb}\right ] }{\left [ \text{gallon}\right ] }\times 2\left [ \frac{\text{gallon}}{\min }\right ] \\ & =\frac{1}{50}y_{1}\left ( t\right ) \left [ \frac{\text{lb}}{\min }\right ] \end{align*}

Hence (1) becomes \frac{dy_{1}\left ( t\right ) }{dt}=1-\frac{1}{50}y_{1}\left ( t\right )

With y_{1}\left ( 0\right ) =0. The ODE is linear and first order, of the form y^{\prime }+p\left ( t\right ) y=g\left ( t\right ) with p\left ( t\right ) =\frac{1}{50} and g\left ( t\right ) =1. Since both p\left ( t\right ) ,g\left ( t\right ) are continuous then a solution exist and is unique. y_{1}^{\prime }+\frac{1}{50}y_{1}=1
Integrating factor is e^{\int \frac{1}{50}dt}=e^{\frac{1}{50}t}, therefore \frac{d}{dt}\left ( e^{\frac{1}{50}t}y_{1}\right ) =e^{\frac{1}{50}t}
Integrating\begin{align*} e^{\frac{1}{50}t}y & =\int e^{\frac{1}{50}t}dt\\ & =50e^{\frac{1}{50}t}+c \end{align*}

Hence y_{1}\left ( t\right ) =50+c^{\frac{-t}{50}}

To find c, from initial conditions\begin{align*} 0 & =y_{1}\left ( 0\right ) \\ & =50+c\\ c & =-50 \end{align*}

Hence the solution to the first phase is\begin{align*} y_{1}\left ( t\right ) & =50-50^{\frac{-t}{50}}\\ & =50\left ( 1-e^{\frac{-t}{50}}\right ) \end{align*}

After t=10 minutes y_{1}\left ( 10\right ) =50\left ( 1-e^{\frac{-1}{5}}\right )

The above value is now used as initial conditions for new problem. The new problem will use t=0 as initial time for simplicity, but it is understood that 10 minutes has already elapsed in global scale.

Second phase

Let y_{2}\left ( t\right ) be the mass of salt at time t in the tank in grams. Hence \begin{align*} y_{2}\left ( 0\right ) & =y_{1}\left ( 10\right ) \\ & =50\left ( 1-e^{\frac{-1}{5}}\right ) \end{align*}

From phase one above, this is the amount of salt in lbs in the tank at this moment. The goal is to find an ODE that describes how y_{2}\left ( t\right ) changes in time. That is, how the mass of salt in the tank changes in time. Using\begin{equation} \frac{dy_{2}}{dt}=R_{in}-R_{out}\tag{2} \end{equation}

Where R_{in} is rate of mass of salt moving into the tank, i.e. how many lbs of salt is being poured in per minute. But now R_{in}=0 since fresh water is poured in. And R_{out} is rate of salt moving out, i.e. how many lbs of salt is leaving the tank per minute. This is found as follows\begin{align*} R_{out} & =\frac{y_{2}\left ( t\right ) }{100}\frac{\left [ \text{lb}\right ] }{\left [ \text{gallon}\right ] }\times 2\left [ \frac{\text{gallon}}{\min }\right ] \\ & =\frac{1}{50}y_{2}\left ( t\right ) \left [ \frac{\text{lb}}{\min }\right ] \end{align*}

Hence (2) becomes\begin{align*} \frac{dy_{2}\left ( t\right ) }{dt} & =0-\frac{1}{50}y_{2}\left ( t\right ) \\ & =-\frac{1}{50}y_{2}\left ( t\right ) \end{align*}

The ODE is linear and first order, of the form y^{\prime }+p\left ( t\right ) y=g\left ( t\right ) with p\left ( t\right ) =\frac{1}{50} and g\left ( t\right ) =0. Since both p\left ( t\right ) ,g\left ( t\right ) are continuous then a solution exist and is unique. This is separable.\begin{align} \frac{dy_{2}}{y_{2}} & =-\frac{1}{50}dt\nonumber \\ \ln \left \vert y_{2}\right \vert & =-\frac{t}{50}+c_{1}\nonumber \\ y_{2}\left ( t\right ) & =ce^{\frac{-t}{50}}\tag{3} \end{align}

To find c, from initial conditions y_{2}\left ( 0\right ) =50\left ( 1-e^{\frac{-1}{5}}\right ) , hence 50\left ( 1-e^{\frac{-1}{5}}\right ) =c

Hence the solution (3) to the second phase is y_{2}\left ( t\right ) =50\left ( 1-e^{\frac{-1}{5}}\right ) e^{\frac{-t}{50}}
After t=10 minutes (which will be 20 in global scale)\begin{align*} y_{2}\left ( 10\right ) & =50\left ( 1-e^{\frac{-1}{5}}\right ) e^{\frac{-1}{5}}\\ & =7.4205\text{ lbs} \end{align*}

Therefore after 20 minutes from the global initial time (or 10 minutes from the start of the second phase), the mass of salt in tank is 7.4205 lbs. Therefore the concentration at the same moment, if needed, will be \frac{7.420\,5}{100}\left [ \frac{\text{lbs}}{\text{gallon}}\right ] = 0.074 \left [ \frac{\text{lbs}}{\text{gallon}}\right ] .

2.2.12 Section 2.3 problem 4

pict

Let y\left ( t\right ) be the mass of salt at time t in the tank in lbs. Hence y\left ( 0\right ) =100 since tank initially contains that much salt. The goal is to find an ODE that describes how y\left ( t\right ) changes in time. That is, how the mass of salt in the tank changes in time. Using\begin{equation} \frac{dy}{dt}=R_{in}-R_{out} \tag{1} \end{equation}

Where R_{in} is rate of mass of salt moving into the tank, i.e. how many lbs of salt is being poured in per minute, which is \begin{align*} R_{in} & =1\left [ \frac{\text{lb}}{\text{gallon}}\right ] \times 3\left [ \frac{\text{gallon}}{\min }\right ] \\ & =3\left [ \frac{\text{lb}}{\min }\right ] \end{align*}

And R_{out} is rate of salt moving out, i.e. how many lbs of salt is leaving the tank per minute. This is found as follows\begin{equation} R_{out}=\frac{y\left ( t\right ) }{V\left ( t\right ) }\frac{\left [ \text{lb}\right ] }{\left [ \text{gallon}\right ] }\times 2\left [ \frac{\text{gallon}}{\min }\right ] \tag{2} \end{equation}

Where V\left ( t\right ) is the volume of the whole mixture at time t. This is different from earlier problems where volume was constant. This is because in this problem the rate of pouring into the tank is larger than the rate of flow out of the tank. The volume at time t can easily be found as\begin{align*} V\left ( t\right ) & =200\left [ \text{gallon}\right ] +3\left [ \frac{\text{gallon}}{\min }\right ] t\left [ \min \right ] -2\left [ \frac{\text{gallon}}{\min }\right ] t[\min \\ & =\left ( 200+t\right ) \left [ \text{gallon}\right ] \end{align*}

This means at any time t, there will be 200+t gallons of mixture in the tank. This value is now used in (2) above to complete the solution. Note that the tank will overflow when 200+t=500 since 500 is the maximum size of the tank. Going back to (2) now it becomes\begin{align*} R_{out} & =\frac{y\left ( t\right ) }{200+t}\frac{\left [ \text{lb}\right ] }{\left [ \text{gallon}\right ] }\times 2\left [ \frac{\text{gallon}}{\min }\right ] \\ & =\frac{2y\left ( t\right ) }{200+t} \end{align*}

Therefore (1) becomes y^{\prime }=3-\frac{2y}{200+t}

This is linear ODE of first order of the form y^{\prime }+p\left ( t\right ) y=g\left ( t\right ) where p\left ( t\right ) =\frac{2}{200+t} and g\left ( t\right ) . Both are continuous for t\geq 0 hence there will be a unique solution for t\geq 0. Now the ODE is solved using an integration factor \frac{dy}{dt}+\frac{2y}{200+t}=3
The integrating factor is e^{2\int \frac{1}{200+t}dt}. To evaluate \int \frac{1}{200+t}dt, let u=200+t, hence \frac{du}{dt}=1 and the integral becomes \int \frac{1}{u}du=\ln \left \vert u\right \vert Therefore \int \frac{1}{200+t}dt=\ln \left \vert 200+t\right \vert and the integrating factor is e^{2\ln \left \vert 200+t\right \vert }=\left \vert 200+t\right \vert ^{2}=\left ( 200+t\right ) ^{2}. Therefore now that the integrating is found, the solution can be written as \frac{d}{dt}\left ( y\left ( 200+t\right ) ^{2}\right ) =3\left ( 200+t\right ) ^{2}
Integrating both sides gives y\left ( 200+t\right ) ^{2}=3\int \left ( 200+t\right ) ^{2}dt
Let u=200+t, then \frac{du}{dt}=t, hence \int \left ( 200+t\right ) ^{2}dt=\int u^{2}du=\frac{u^{3}}{3}+c_{1}=\frac{\left ( 200+t\right ) ^{3}}{3}+c_{1}. Therefore the above becomes\begin{align*} y\left ( 200+t\right ) ^{2} & =3\left ( \frac{\left ( 200+t\right ) ^{3}}{3}+c_{1}\right ) \\ & =\left ( 200+t\right ) ^{3}+c \end{align*}

Solving for y\left ( t\right ) gives\begin{equation} \fbox{$y\left ( t\right ) =\left ( 200+t\right ) +\frac{c}{\left ( 200+t\right ) ^2}$} \tag{3} \end{equation}

Now c is found from initial conditions. Given that y\left ( 0\right ) =100, then from the above\begin{align*} 100 & =200+\frac{c}{\left ( 200\right ) ^{2}}\\ & =200+\frac{c}{40000}\\ c & =\left ( -100\right ) \left ( 40000\right ) \\ & =-4\times 10^{6} \end{align*}

Therefore the solution (3) becomes\begin{equation} y\left ( t\right ) =\left ( 200+t\right ) -\frac{4\times 10^{6}}{\left ( 200+t\right ) ^{2}} \tag{4} \end{equation}

Now the above ODE is only valid until the tank overflows. This value of time is found by solving 200+t=500 for t, which gives t=300. Hence (4) becomes\begin{equation} \fbox{$y\left ( t\right ) =\left ( 200+t\right ) -\frac{4\times 10^6}{\left ( 200+t\right ) ^2}\qquad 0\leq t\leq 300$} \tag{5} \end{equation}
At t=300 minutes, the mass of salt in lbs is therefore y\left ( 300\right ) which is\begin{align*} y\left ( 300\right ) & =\left ( 200+300\right ) -\frac{4\times 10^{6}}{\left ( 200+300\right ) ^{2}}\\ & =484\text{ }\left [ \text{lbs}\right ] \end{align*}

And since the volume now is 500 gallons, then the concentration at time t=300 minutes is \frac{484}{500}\left [ \frac{\text{lbs}}{\text{gallon}}\right ] =0.968\,\left [ \frac{\text{lbs}}{\text{gallon}}\right ]

If the tank had infinite capacity, then using the solution found in (5) and dividing by current volume, which was found before to be 200+t and then taking the limit t\rightarrow \infty gives the answer. Let \rho \left ( t\right ) be now the concentration in \left [ \frac{\text{lbs}}{\text{gallon}}\right ] at any time t. Then\begin{align*} \rho \left ( t\right ) & =\frac{\left ( 200+t\right ) -\frac{4\times 10^{6}}{\left ( 200+t\right ) ^{2}}}{V\left ( t\right ) }=\frac{\left ( 200+t\right ) -\frac{4\times 10^{6}}{\left ( 200+t\right ) ^{2}}}{200+t}\\ & =1-\frac{4\times 10^{6}}{200+t} \end{align*}

As t\rightarrow \infty then \rho \left ( t\right ) \rightarrow 1. Therefore at 300 minutes the concentration is 96.8\% of the theoretical limit. The following is a plot of \rho \left ( t\right ) as function of time. At t=0 the concentration is 0.5 since this is the initial condition.

pict

2.2.13 Section 2.4 problem 1

Determine an interval which the given initial value problem is valid. \left ( t-3\right ) y^{\prime }+\ln \left ( t\right ) y=2t with y\left ( 1\right ) =2.

This is linear first order ODE. In standard form it becomes y^{\prime }+\frac{\ln \left ( t\right ) }{t-3}y=\frac{2t}{t-3}, and comparing to y^{\prime }+p\left ( t\right ) y=g\left ( t\right ) then \begin{align*} p\left ( t\right ) & =\frac{\ln \left ( t\right ) }{t-3}\\ g\left ( t\right ) & =\frac{2t}{t-3} \end{align*}

p\left ( t\right ) is not not continuous at t=3 and also at t=0 since \ln \left ( 0\right ) =-\infty . g\left ( t\right ) is not continuous at t=3. Therefore the region must include initial point, which is t=1 but not include t=3 nor t=0. Hence 0<t<3

And for forward only ODE the region is 1\leq t<3

2.2.14 Section 2.4 problem 2

Determine an interval which the given initial value problem is valid. t\left ( t-4\right ) y^{\prime }+y=2t with y\left ( 2\right ) =1.

This is linear first order ODE. In standard form it becomes y^{\prime }+\frac{1}{t\left ( t-4\right ) }y=\frac{2}{\left ( t-4\right ) }, and comparing to y^{\prime }+p\left ( t\right ) y=g\left ( t\right ) then \begin{align*} p\left ( t\right ) & =\frac{1}{t\left ( t-4\right ) }\\ g\left ( t\right ) & =\frac{2}{\left ( t-4\right ) } \end{align*}

p\left ( t\right ) is not continuous at t=0 and t=4 while g\left ( t\right ) is not continuous at t=4. Therefore the region must include initial point, which is t=2 but not include t=4 nor t=0. Hence 0<t<4

And for forward only ODE the region is 2\leq t<4

2.2.15 Section 2.4 problem 3

Determine an interval which the given initial value problem is valid. y^{\prime }+\tan \left ( t\right ) y=\sin \left ( t\right ) with y\left ( \pi \right ) =0.

This is linear first order ODE. Comparing to y^{\prime }+p\left ( t\right ) y=g\left ( t\right ) then \begin{align*} p\left ( t\right ) & =\tan \left ( t\right ) \\ g\left ( t\right ) & =\sin \left ( t\right ) \end{align*}

g\left ( t\right ) is continuous everywhere but p\left ( t\right ) is not continuous at \left \{ \cdots ,-\frac{\pi }{2},\frac{\pi }{2},\frac{3\pi }{2},\cdots \right \}  therefore the region must be between \frac{\pi }{2} and \frac{3\pi }{2} since the initial point \pi is inside this region. Hence \frac{\pi }{2}<t<1.5\pi

2.2.16 Section 2.4 problem 4

Determine an interval which the given initial value problem is valid. \left ( 4-t^{2}\right ) y^{\prime }+2ty=3t^{2} with y\left ( -3\right ) =1.

This is linear first order ODE. In standard form it becomes y^{\prime }+\frac{2t}{\left ( 4-t^{2}\right ) }y=\frac{3t^{2}}{\left ( 4-t^{2}\right ) }, and comparing to y^{\prime }+p\left ( t\right ) y=g\left ( t\right ) then\begin{align*} p\left ( t\right ) & =\frac{2t}{\left ( 4-t^{2}\right ) }\\ g\left ( t\right ) & =\frac{3t^{2}}{\left ( 4-t^{2}\right ) } \end{align*}

p\left ( t\right ) is not not continuous at t^{2}=4 or t=\pm 2 and the same for g\left ( t\right ) . Therefore the region must include initial point, which is t=-3 but not include t=\pm 2. Hence -\infty <t<-2

And for forward only ODE the region is -3\leq t<-2

2.2.17 Section 2.4 problem 5

Determine an interval which the given initial value problem is valid. \left ( 4-t^{2}\right ) y^{\prime }+2ty=3t^{2} with y\left ( 1\right ) =-3.

This is linear first order ODE. In standard form it becomes y^{\prime }+\frac{2t}{\left ( 4-t^{2}\right ) }y=\frac{3t^{2}}{\left ( 4-t^{2}\right ) }, and comparing to y^{\prime }+p\left ( t\right ) y=g\left ( t\right ) then \begin{align*} p\left ( t\right ) & =\frac{2t}{\left ( 4-t^{2}\right ) }\\ g\left ( t\right ) & =\frac{3t^{2}}{\left ( 4-t^{2}\right ) } \end{align*}

p\left ( t\right ) is not not continuous at t^{2}=4 or t=\pm 2 and the same for g\left ( t\right ) . Therefore the region must include initial point, which is t=1 but not include t=\pm 2. Hence -2<t<2

And for forward only ODE the region is 1\leq t<2

2.2.18 Section 2.4 problem 6

Determine an interval which the given initial value problem is valid. \ln \left ( t\right ) y^{\prime }+y=\frac{1}{\tan \left ( t\right ) } with y\left ( 2\right ) =3.

This is linear first order ODE. In standard form it becomes y^{\prime }+\frac{1}{\ln \left ( t\right ) }y=\frac{1}{\tan \left ( t\right ) \ln \left ( t\right ) }, and comparing to y^{\prime }+p\left ( t\right ) y=g\left ( t\right ) then \begin{align*} p\left ( t\right ) & =\frac{1}{\ln \left ( t\right ) }\\ g\left ( t\right ) & =\frac{1}{\tan \left ( t\right ) \ln \left ( t\right ) } \end{align*}

When t=1 then \ln \left ( t\right ) =0 and p\left ( t\right ) becomes unbounded. And since for real t then t must remain positive, else \ln \left ( t\right ) becomes complex. Then p\left ( t\right ) says that t\geq 0 and t\neq 1. Looking at g\left ( t\right ) then \tan \left ( t\right ) =0 when t=\left \{ \cdots ,-\pi ,\pi ,\cdots \right \} hence the region that includes initial point t_{0}=2 must be inside these. Therefore the singular points are t=1,-\pi ,\pi and t\geq 0. Putting all these together, the region is 1<t<\pi

And for forward only ODE the region is 2\leq t<\pi