Skip to main content

Coordinated Differential Equations

Appendix C Hints and Answers to Selected Exercises

1 A First Look at Differential Equations
1.1 Modeling with Differential Equations
1.1.9 Exercises

1.1.9.29.

Hint.
Rewriting the differential equation as x+axq(t)=0 and using the fact that
x(t)=aCeataeatt0teasq(s)ds+b(t),
we see that
x(t)+ax(t)q(t)=aCeataeatt0teasq(s)ds+q(t)+aCeat+aeatt0teasq(s)dsq(t)=0.

1.1.9.31.

Hint.
Think about the limit of the interaction term as the number of prey becomes very large.

1.4 Analyzing Equations Numerically
1.4.6 Exercises

1.4.6.7.

Hint.
Hints for part (2):
  • For fixed i show that
    ai+1(1+s)ai+t(1+s)[(1+s)ai1+t]+t(1+s){(1+s)[(1+s)ai2+t]+t}+t(1+s)i+1a0+[1+(1+s)+(1+s)2++(1+s)i]t.
  • Now use a geometric sum to show that
    ai+1(1+s)i+1a0+ts[(1+s)i+11]=(1+s)i+1(ts+a0)ts.
  • Apply part (1) to derive
    ai+1e(i+1)s(ts+a0)ts.

1.5 First-Order Linear Equations
1.5.8 Exercises

1.5.8.16.

Hint.
y=3x4+8x3+6x2+126(x+1)2

1.5.8.17.

Hint.
y=(ex+1)/(x+1)2

1.5.8.19.

Hint.
y=exsecx

1.5.8.20.

Hint.
y=sinx/(x+2)

1.5.8.21.

Hint.
If x(t) is the amount of salt in the tank at time t, we know that x(0)=10. The volume of the tank is V=200+5t. We can model the amount of salt in the tank at time t with a differential equation,
dxdt=rate inrate out=10(0.1)5xV=15x200+5t=1x40+t.
The resulting equation
dxdt+140+tx=1
is a first order linear differential equation. An integrating factor for this equation is given by
μ(t)=exp(140+tdt)=40+t.
Multiplying both sides of the differential equation by μ(t), we have
ddt[(40+t)x]=(40+t)dxdt+x=(40+t)(dxdt+140+tx)=40+t.
Integrating both sides of this equation, we obtain
(40+t)x=40t+t22+C.
Using the intial condition x(0)=10, we can determine that C=400 or
x(t)=t2+80t+8002t+80.
The tank is full at time t=400/5=80, and the tank contains x(80)=170/356.67 kilograms of salt when the tank is full.

1.5.8.26.

Hint.
  1. If y=y1+1/v, then y=y1v/v2. Substituting into our original equation, we obtain
    y=y1vv2=p+qy1+ry12vv2.
    On the other hand,
    y=p+q(y1+1v)+r(y1+1v)2=p+qy1+qv+ry12+2ry1v+rv2=y1+qv+2ry1v+rv2.
    Therefore,
    vv2=qv+2ry1v+rv2,
    which is just the first-order linear equation
    v+[q(t)+2r(t)y1(t)]v=r(t).
  2. y=t+1Ct
  3. y(t)=1Ccostsint+sint
  4. y(t)=2+1Cet1

1.6 Existence and Uniqueness of Solutions
1.6.5 Exercises

1.6.5.2.

Hint.
(b) Make sure that the derivative of y(t) exists at t=t0.

2 Second-Order Linear Equations
2.1 Homogeneous Linear Equations
2.1.7 Exercises

2.1.7.31.

Hint.
Pay careful attention to units.

2.1.7.32.

Hint.
  1. Observe that
    ax1+b1+cx1=a(b2a)2ebt/2a+b(b2a)ebt/2a+cebt/2a=ebt/2a(b24ab22a+c)=ebt/2a(b2+4ac4a)=0.
  2. If y=v(t)x1(t)=v(t)ebt/2a is a solution to our differential equation, then
    ay+by+cy=a(vx1+2vx1+vx1)+b(vx1+vx1)+cvx1=avx1+2avx1+bvx1+v(ax1+bx1+cx1)=avebt/2a+[2a(b2a)ebt/2a+bebt/2a]v=avebt/2a=0.
    Since a0, we know that v=0. Hence, v(t)=c1+c2t.

2.1.7.33.

Hint.
  1. x+px+qx=(vx1+2vx1+vx1)+p(vx1+vx1)+q(vx1)=x1v+2vx1+px1v+v(x1+px1+qx1)=x1v+(2x1+px1)v=0.
  2. If u=v, then x1u+(2x1+px1)u=0.
  3. If x1(t)=1/t, then
    2t2x1+3tx1x1=2t2(2t3)+3t(1t2)1t=0.
    If we assume that x=v/t is a second solution, then
    2t2x+3txx=2tvv=0.
    If we let u=v, then a solution of 2tuu=0 is u=t and v=tdt=2t3/2/3. Therefore, the second solution to our equation is
    x=vt=23t.

2.2 Forcing
2.2.7 Exercises

2.2.7.2.

Answer.
The solution to the associated homogeneous equation is c1e2x+c2ex, and we do not see any e2x or ex in the right-hand-side, so we look at Table 2.2.12 and see the solution has the form y(x)=Ax2+Bx+C. We then get y(x)=2Ax+B,y(x)=2A, so
yy2y=2A(2Ax+B)2(Ax2+Bx+C)=
2Ax2+(2A2B)x+(2AB2C)=4x22x9
We then have to equate the coefficients. 2A=4 gives A=2, 2A2B=2 gives B=A+1=3, and 2AB2C=9 gives C=1, so the answer is y(x)=2x2+3x+1

2.2.7.25.

Hint.
Suppose that that f(t) and g(t) are linearly dependent on an interval I=(a,b). Then one function is a multiple of the other, say f(t)=cg(t). Thus, f(t)=cg(t) and
W(f,g)=f(t)g(t)f(t)g(t)=cg(t)g(t)cg(t)g(t)=0.
Conversely, suppose that W(f,g)=0 for all t in (a,b). If g=0, then 0f=g and the two functions are linearly dependent. Assume that g(t0)0 for some t0 in (a,b). Since g is differentiable, it must also be continuous and there is some interval (c,d) contained in (a,b) such that t0(c,d) and g does not vanish on this interval. Therefore,
ddt(fg)=fgfgg2=W(f,g)g2=0,
and f/g is constant on the interval (c,d). Thus, f(t0)=cg(t0) and f(t0)=cg(t0). Since f and cg are both solutions to the differential equation y+py+qy=0 and have the same initial condition, f(t)=cg(t) for all t(a,b) by the existence and uniqueness theorem. Consequently, f and g are linearly dependent.

2.2.7.26.

Hint.
  1. We can rewrite 2t2y+3tyy=0 as
    y+32ty12t2y=0.
    Since p(t)=1/2t, Abel’s Theorem tells us that
    W(y1,y2)=cexp(32tdt)=cexp(32lnt)=ct3/2.
  2. Since y1 and y2 are solutions to our differential equation, we know that
    y1+p(t)y1+q(t)y1=0y2+p(t)y2+q(t)y2=0.
    Multiplying the first equation by y2 and the second equation by y1 and subtracting, we obtain
    (2.2.5)(y1y2y1y2)+p(t)(y1y2y1y2)=0.
    If
    W(t)=W(y1,y2)=y1y2y1y2,
    then
    W=y1y2y1y2,
    and equation (2.2.5) becomes
    W+p(t)W=0.
    This equation is separable with solution
    W(t)=cexp(p(t)dt).

2.2.7.27.

Hint.
  1. If yp=u1y1+u2y2, then
    yp=u1y1+u1y1+u2y2+u2y2=u1y1+u2y2yp=u1y1+u1y1+u2y2+u2y2.
    Substituting these expressions into equation (2.2.6), we have
    yp+pyp+qyp=(u1y1+u1y1+u2y2+u2y2)+p(u1y1+u2y2)+q(u1y1+u2y2)=u1[y1+py1+qy1]+u2[y2+py2+qy2]+u1y1+u2y2=u1y1+u2y2=f(t).
  2. If we solve the system
    u1(t)y1(t)+u2(t)y2(t)=0u1(t)y1(t)+u2(t)y2(t)=f(t).
    for u1 and u2, we obtain
    u1(t)=y2(t)f(t)W[y1,y2](t)u2(t)=y1(t)f(t)W[y1,y2](t).
  3. Integrate the two equations from part (2).
  4. The general solution to the homogeneous equation y+4y=0 is
    yh=c1cos2t+c2sin2t.
    To find a particular solution, assume that the solution has the form
    yp=u1(t)cos2t+u2(t)sin2t.
    By part (2)
    u1(t)=3costu2(t)=32csct3sint.
    Integrating, we obtain
    u1(t)=3sintu2(t)=32ln|csctcott|+3cost.
    Therefore,
    yp(t)=u1(t)y1(t)+u2(t)y2(t)=3sintcos2t+[32ln|csctcott|+3cost]sin2t,
    and the general solution is
    y=yh+yp=c1cos2t+c2sin2t3sintcos2t+[32ln|csctcott|+3cost]sin2t.

2.3 Sinusoidal Forcing
2.3.5 Exercises

2.3.5.2.

Hint.
Assume the complex solution has form yc=Ae3it.

2.3.5.12.

Hint.
Assume the complex solution has form yc=Ae3it.

3 The Laplace Transform
3.1 The Laplace Transform
3.1.7 Exercises

3.1.7.15.

Answer.
If the function in question were simply 1s4, we would observe that L(t3)(s)=6s4, so the answer would be t36. Since we subtract 2 from the argument s, we multiply by e2t, so the answer is e2tt36.

3.1.7.16.

Answer.
We will split this into partial fractions via
1(s+1)2(s24)=A(s+1)+B(s+1)2+Cs+2+Ds2
Clearing denominators gives
1=[A(s+1)+B](s24)+C(s2)(s+1)2+D(s+2)(s+1)2
Plugging in s=1 yields 3B=1, so B=13. Plugging in s=2 yields 36D=1, so D=136. Plugging in s=2 yields 4C=1, so C=14. Plugging in s=0 and using our values for B,C, and D yields 4A+179=1, so A=29 So we are to find the inverse Laplace of
F(s)=29(s+1)13(s+1)21/4s+2+1/36s2
=2/9s+11/3(s+1)21/4s+2+1/36s2
This gives us f(t)=29et13tet14e2t+136e2t

3.2 Solving Initial Value Problems
3.2.6 Exercises

3.2.6.1.

Answer.
Taking the Laplace transformation of both sides yields
s2Y(s)s2sY(s)+23Y(s)=24s+3
so
Y(s)=24(s3)(s+1)(s+3)+s2(s3)(s+1)
A partial fractions decomposition reveals
Y(s)=1s3+3s+1+2s+3+1/4s3+3/4s+1
So the answer is
y(t)=e3t3et+2e3t+14e3t+34et
y(t)=54e3t94et+2e3t

3.4 Convolution
3.4.5 Exercises

3.4.5.13.

Answer.
Taking the Laplace transformation of both sides yields
2s2Y(s)+5sY(s)+2Y(s)=G(s)
Therefore, Y(s)=G(s)2s2+5s+2. We now want to factor 2s2+5s+2. We multiply the s2 and constant coefficients to get 4, and search for two numbers multiplying to 4 and adding to the s coefficient (5). We see that 4 and 1 work, and we can write 2s2+5s+2=2s2+4s+s+2=2s(s+2)+(s+2)=(2s+1)(s+2), so we get
y(t)=L1{G(s)(2s+1)(s+2)}=g(t)L1{1(2s+1)(s+2)}
We could use partial fractions or another convolution to calculate L1{1(2s+1)(s+2)}. If we choose convolution, we see L1{1(2s+1)(s+2)} =12et/2e2t =120teτ/2e2(tτ)dτ =13(et/2e2t), so the final answer is
13g(t)(et/2e2t)=130tg(τ)(eτt2e2(τt))dτ

5 Linear Systems and Linearization
5.1 Linear Algebra in a Nutshell
5.1.6 Exercises

5.1.6.4.

Answer.
det(A)=1(3)4(2)=5

5.1.6.6.

Answer.
det(A)=2(4)(6)(2)=4

5.2 Planar Systems
5.2.6 Exercises

5.2.6.2.

Answer.
This is equivalent to x=Ax, where A=(1222). The characteristic polynomial is (1λ)(2λ)2(2) =λ2λ6, with zeros λ=2,3. We then get
A(2)I=(1+2222+2)=(1224)
So λ=2 has corresponding eigenvector (21). Analogously,
A3I=(132223)=(4221)
So λ=3 has corresponding eigenvector (12). Therefore the set of solutions is
c1e2t(21)+c2e3t(12)

5.2.6.6.

Answer.
From exercise 2, we see the general solution to the system x=x+2y,y=2x+2y is
c1e2t(21)+c2e3t(12)
Plugging in t=0 and setting the top to 19, the bottom to 8, we end up at the system
2c1+c2=19c1+2c2=8
Solving this yields c1=6,c2=7, and we thus have the solution
6e2t(21)+7e3t(12)
which can also be written as
x(t)=12e2t+7e3ty(t)=6e2t+14e3t

5.2.6.10.

Hint.
Assume that your solution must be of the form
xp=(a2t2+a1t+a0b2t2+b1t+b0.)
This is called the method of undetermined coefficients.

5.3 Phase Plane Analysis of Linear Systems
5.3.6 Exercises

5.3.6.2.

Answer.
(a) The characteristic polynomial is (4λ)(3λ)+6, with zeros (and thus eigenvalues) λ=2,3
(b) If the eigenvalues have the same sign, the one with a bigger absolute value is dominant. Since these eigenvalues have different signs, there is no dominant eigenvalue.
(c) For λ=2 we subtract 2 from the main diagonal of the matrix to obtain A+2I=(6231), so v=(13) is an eigenvector.
For λ=3, we subtract 3 from the main diagonal to obtain A3I=(1236), so v=(21)
(d) There are two straight-line solutions, one corresponding to each eigenvalue. For λ=2 we got (13), so the solution set here is the set of multiples of (13), or equivalently, the line 1y=3x.
For λ=3 we got the eigenvalue (21), so the straight-line solution is the line passing through (21) and the origin, or the line 2y=1x.
(e) Since there is one positive and one negative eigenvalue, the equilibrium solution is a saddle.
(f)