Type a command or search the notes. Tab cycles completions. Up and down recall history when blank. Escape dismisses, then leaves the prompt.
MATH
[1251]

Seen Questions

Worked solutions on exact equations, integrating factors, characteristic roots and first- and second-order models.

These past exam questions test exact equations, integrating factors, characteristic equations and forcing terms.

The first question also touches the level-set and gradient geometry and multivariable linear approximation of a later chapter. Its solution derives the tangent-plane and normal-vector facts it needs, so it can be read with the ODE material alone.

[DEFINITION]

Exact differential equation

An equation in differential form

M(x,y) dx+N(x,y) dy=0M(x,y)\,dx+N(x,y)\,dy=0

is exact on a domain when there is a potential function F(x,y)F(x,y) such that

Fx=M,Fy=N.F_x=M,\qquad F_y=N.

Here FxF_x means differentiate with respect to xx while holding yy fixed, and FyF_y means differentiate with respect to yy while holding xx fixed. The later multivariable-calculus chapter develops this notation and the level-curve geometry in more detail.

[METHOD]

Approximate a function linearly

For a differentiable f:R2→Rf:\mathbb{R}^2\to\mathbb{R} and a base point a=(x0,y0)\mathbf{a}=(x_0,y_0), approximate nearby values by

f(x)≈f(a)+∇f(a)⋅(x−a).f(\mathbf{x}) \approx f(\mathbf{a}) +\nabla f(\mathbf{a})\cdot(\mathbf{x}-\mathbf{a}).

To use the approximation, evaluate ff and its first partial derivatives at a\mathbf{a}, then substitute the displacement x−a\mathbf{x}-\mathbf{a}.

[THEOREM]

Gradient normal to a level set

Let F:D→RF:D\to\mathbb{R} be C1C^1 on an open neighbourhood of x0\mathbf{x}_0 and let

S={x∈D:F(x)=c}.S=\{\mathbf{x}\in D:F(\mathbf{x})=c\}.

If x0∈S\mathbf{x}_0\in S and ∇F(x0)≠0\nabla F(\mathbf{x}_0)\ne\mathbf{0}, then ∇F(x0)\nabla F(\mathbf{x}_0) is normal to SS at x0\mathbf{x}_0.

[EXAMPLE]

Exact equations and level-set geometry

Answer the following questions about exact differential equations.

(a)

First classify the differential equation

(3x2y2+7cos⁡(6x)) dx+(32x4y+7ey) dy=0.\left(3x^2y^2+7\cos(6x)\right)\,dx +\left(\frac32x^4y+7e^y\right)\,dy=0.

Select all options which correctly describe this equation:

  • non-linear;
  • first order;
  • second order;
  • initial value problem;
  • boundary value problem;
  • linear;
  • constant coefficient;
  • fifth order.

Now consider a differential equation of the form

M(x,y) dx+N(x,y) dy=0.(1)M(x,y)\,dx+N(x,y)\,dy=0. \tag{1}

Here M,N:R2→RM,N:\mathbb R^2\to\mathbb R are functions which are continuous everywhere, differentiable everywhere and whose first-order partial derivatives are also continuous everywhere.

(b)

Suppose equation (1) is not exact. To use the integrating-factor method, determine when multiplying by a differentiable, nonvanishing function u(x,y)u(x,y) produces an equivalent exact equation. State a necessary condition.

(c)

Now suppose equation (1) is exact. For a constant cc, its solution may be written in the form

F(x,y)=c,F(x,y)=c,

where F(x,y)F(x,y) is to be constructed from MM and NN.

(d)

Suppose (1) is an exact differential equation with solution F(x,y)=cF(x,y)=c, where cc is a constant. This solution can be interpreted as a level curve of the surface

z=F(x,y).z=F(x,y).

To understand the behaviour of the solution near a given point (x0,y0)(x_0,y_0), Liu wishes to use Taylor polynomials to find an approximation. Give the linear approximation to the surface near (x0,y0)(x_0,y_0).

(e)

For the geometry of solutions, consider the following example.

Let F:R2→RF:\mathbb R^2\to\mathbb R be a potential for the exact differential form below, with continuous second-order derivatives. Consider its level curve

F(x,y)=20F(x,y)=20

which passes through the point (5,4)(5,4) and for which locally the function y(x)y(x) is implicitly defined by the exact differential equation

6x2y6 dx+12x3y5 dy=0.6x^2y^6\,dx+12x^3y^5\,dy=0.

Compute a normal vector to the surface z=F(x,y)z=F(x,y) at the point (5,4,20)(5,4,20).

Solution

The first equation is in differential form, so it contains dxdx and dydy but no derivative of order greater than one. It is first order. It is non-linear because, after writing it as an equation for dy/dxdy/dx, the terms y2y^2 and eye^y occur non-linearly. It is neither an initial-value nor a boundary-value problem because no data select a particular solution and it is not a constant-coefficient linear equation. The correct selections in part (a) are

non-linear and first order.{\text{non-linear and first order}.}

For equation (1), multiplying by uu gives

uM dx+uN dy=0.uM\,dx+uN\,dy=0.

Exactness requires the mixed partials of the new coefficients to agree. Therefore, a necessary condition is

∂∂y(uM)=∂∂x(uN).{ \frac{\partial}{\partial y}\bigl(uM\bigr) =\frac{\partial}{\partial x}\bigl(uN\bigr) }.

Expanding both sides gives the equivalent first-order equation for the integrating factor:

M uy+u My=N ux+u Nx.M\,u_y+u\,M_y=N\,u_x+u\,N_x.

The unknown function uu must satisfy this equation: multiplying by an arbitrary non-zero function is not enough, because the new coefficients must pass the exactness test. On a simply connected domain, with the required continuous first partial derivatives, this equality is also sufficient for a potential; nonvanishing uu is what preserves the original solution curves.

If the original equation is exact, there is a potential function FF such that

Fx=M,Fy=N.F_x=M,\qquad F_y=N.

Integrating Fx=MF_x=M with respect to xx, while treating yy as constant, gives

F(x,y)=∫M(x,y) dx+h(y),F(x,y)=\int M(x,y)\,dx+h(y),

where hh is an arbitrary function of yy. Differentiating this expression with respect to yy gives

Fy(x,y)=∂∂y(∫M(x,y) dx)+h′(y).F_y(x,y) =\frac{\partial}{\partial y}\left(\int M(x,y)\,dx\right)+h'(y).

Matching this with N(x,y)N(x,y) determines h′(y)h'(y) and hence h(y)h(y). Therefore the requested construction is

F(x,y)=∫M(x,y) dx+h(y),Fy=N.{ F(x,y)=\int M(x,y)\,dx+h(y), \qquad F_y=N. }

The function FF is determined only up to an additive constant, which is absorbed into the constant cc in F(x,y)=cF(x,y)=c.

For small changes Δx\Delta x and Δy\Delta y, the partial derivatives give the first-order change in height:

F(x0+Δx,y0+Δy)≈F(x0,y0)+Fx(x0,y0)Δx+Fy(x0,y0)Δy.F(x_0+\Delta x,y_0+\Delta y) \approx F(x_0,y_0)+F_x(x_0,y_0)\Delta x+F_y(x_0,y_0)\Delta y.

Replacing Δx\Delta x and Δy\Delta y by x−x0x-x_0 and y−y0y-y_0 gives the tangent plane to z=F(x,y)z=F(x,y) at (x0,y0,F(x0,y0))(x_0,y_0,F(x_0,y_0)):

z=F(x0,y0)+Fx(x0,y0)(x−x0)+Fy(x0,y0)(y−y0).{ z=F(x_0,y_0) +F_x(x_0,y_0)(x-x_0) +F_y(x_0,y_0)(y-y_0) }.

Since Fx=MF_x=M and Fy=NF_y=N, the same plane can be written as

z=F(x0,y0)+M(x0,y0)(x−x0)+N(x0,y0)(y−y0).z=F(x_0,y_0) +M(x_0,y_0)(x-x_0) +N(x_0,y_0)(y-y_0).

For part (e), the differential equation identifies the first partial derivatives of a potential:

Fx=6x2y6,Fy=12x3y5.F_x=6x^2y^6, \qquad F_y=12x^3y^5.

Integrating the first expression with respect to xx gives

F(x,y)=2x3y6+h(y).F(x,y)=2x^3y^6+h(y).

Differentiating with respect to yy gives

Fy=12x3y5+h′(y).F_y=12x^3y^5+h'(y).

Comparison with the coefficient of dydy gives h′(y)=0h'(y)=0, so hh is a constant. That constant can be chosen so that F(5,4)=20F(5,4)=20.

For the graph z=F(x,y)z=F(x,y), at (x,y,F(x,y))(x,y,F(x,y)) changing xx while holding yy fixed gives the tangent direction (1,0,Fx(x,y))(1,0,F_x(x,y)); changing yy while holding xx fixed gives (0,1,Fy(x,y))(0,1,F_y(x,y)). A normal vector must be perpendicular to both. The vector (−Fx(x,y),−Fy(x,y),1)(-F_x(x,y),-F_y(x,y),1) has zero dot product with each direction, since −Fx+Fx=0-F_x+F_x=0 and −Fy+Fy=0-F_y+F_y=0. Thus one normal vector is

(−Fx(x,y)−Fy(x,y)1).\begin{pmatrix} -F_x(x,y)\\ -F_y(x,y)\\ 1 \end{pmatrix}.

At (5,4)(5,4),

Fx(5,4)=6(5)2(4)6=150⋅4096=614400,Fy(5,4)=12(5)3(4)5=1500⋅1024=1536000.\begin{aligned} F_x(5,4) &=6(5)^2(4)^6 =150\cdot4096 =614400,\\ F_y(5,4) &=12(5)^3(4)^5 =1500\cdot1024 =1536000. \end{aligned}

Consequently, one normal vector at (5,4,20)(5,4,20) is

n=(−614400−15360001).{ \mathbf n= \begin{pmatrix} -614400\\ -1536000\\ 1 \end{pmatrix}. }

An exact differential form gives a level curve of a potential. The graph of the same potential is a surface in three dimensions and its partial derivatives supply the tangent-plane and normal information. The linear approximation is the tangent plane, while the original solution curve is obtained by intersecting that surface with the horizontal plane z=cz=c.

[METHOD]

Integrating factor

For y′+p(x)y=q(x)y'+p(x)y=q(x) on an interval where pp and qq are continuous:

  1. Choose μ(x)=e∫p(x) dx.\mu(x)=e^{\int p(x)\,dx}.
  2. Since μ′=pμ\mu'=p\mu, multiply the equation by μ\mu: μy′+μpy=(μy)′=μq.\mu y'+\mu p y=(\mu y)'=\mu q.
  3. Integrate: μy=∫μq dx+C.\mu y=\int\mu q\,dx+C.
  4. Divide by μ\mu: y=1μ(∫μq dx+C).y=\frac1{\mu}\left(\int\mu q\,dx+C\right).
[DEFINITION]

Characteristic equation

For the homogeneous equation

ay′′+by′+cy=0,ay''+by'+cy=0,

the exponential trial function y=erxy=e^{rx} gives the characteristic equation

ar2+br+c=0.ar^2+br+c=0.
[METHOD]

Undetermined coefficients

For suitable forcing functions:

  • a polynomial of degree nn suggests a general polynomial of degree nn;
  • eβxe^{\beta x} suggests CeβxCe^{\beta x};
  • pn(x)eβxp_n(x)e^{\beta x} suggests qn(x)eβxq_n(x)e^{\beta x}, where qnq_n has degree nn;
  • sin⁡(βx)\sin(\beta x) or cos⁡(βx)\cos(\beta x) suggests C1cos⁡(βx)+C2sin⁡(βx);C_1\cos(\beta x)+C_2\sin(\beta x);
  • eαxsin⁡(βx)e^{\alpha x}\sin(\beta x) or eαxcos⁡(βx)e^{\alpha x}\cos(\beta x) suggests eαx(C1cos⁡(βx)+C2sin⁡(βx)).e^{\alpha x} \bigl(C_1\cos(\beta x)+C_2\sin(\beta x)\bigr).

Substitute the trial into the equation and match coefficients. If the trial overlaps with a homogeneous mode whose characteristic root has multiplicity ss, multiply the whole trial by xsx^s (so a simple overlap requires one factor of xx, while a repeated root requires x2x^2).

[EXAMPLE]

First- and second-order price models

Suppose a price B(t)B(t) is modelled by

dBdt=−a(t)(B(t)−G(t))+P(t).(1)\frac{dB}{dt}=-a(t)\bigl(B(t)-G(t)\bigr)+P(t). \tag{1}

where G(t)G(t) is a target price, a(t)a(t) is a relaxation rate and P(t)P(t) is an additional price-change rate. These are prescribed functions of tt, not of BB. Their particular values are not specified in part (a).

(a)

Which classes must Equation (1) belong to? Select all that apply:

  • linear;
  • non-linear;
  • homogeneous;
  • inhomogeneous;
  • an exact differential equation;
  • a separable differential equation.

(b)

At the launch of the coin, measurements gave

a(t)=t9,G(t)=22t4+4t+55,P(t)=3cos⁡(t).a(t)=t^9, \qquad G(t)=2\sqrt[5]{2t^4+4t+5}, \qquad P(t)=3\cos(t).

Rewrite the equation as

ddt(ϕ(t)B(t))=ϕ(t)(a(t)G(t)+P(t)),\frac{d}{dt}\bigl(\phi(t)B(t)\bigr) =\phi(t)\bigl(a(t)G(t)+P(t)\bigr),

so that

B(t)=1ϕ(t)(∫ϕ(t)(a(t)G(t)+P(t)) dt+C).B(t)=\frac1{\phi(t)} \left(\int\phi(t)\bigl(a(t)G(t)+P(t)\bigr)\,dt+C\right).

Find ϕ(t)\phi(t).

(c)

A second illustrative price N(t)N(t) satisfies

7N′′(t)=−70N′(t)−175N(t)+8e−5t.7N''(t)=-70N'(t)-175N(t)+8e^{-5t}.
(i)

To find the homogeneous equation, solve the associated characteristic equation, written in terms of a variable ν\nu:

0= ?0=\ ?
(ii)

Find a particular solution NpN_p that contains no undetermined coefficients.

Solution

Expand the right-hand side of (1) and collect the terms involving BB:

dBdt+a(t)B(t)=a(t)G(t)+P(t).\frac{dB}{dt}+a(t)B(t)=a(t)G(t)+P(t).

This is linear in the unknown function BB. It is not necessarily inhomogeneous: aG+PaG+P could vanish identically (for example, a=G=1a=G=1 and P=−1P=-1). Nor must it be separable, although some choices of the prescribed functions make it so. In differential form,

(a(t)B−a(t)G(t)−P(t)) dt+dB=0,\bigl(a(t)B-a(t)G(t)-P(t)\bigr)\,dt+dB=0,

the partial derivative of the dtdt coefficient with respect to BB is a(t)a(t), whereas the partial derivative of the dBdB coefficient with respect to tt is zero. Exactness is therefore not guaranteed either. The only classification that must hold in part (a) is

linear.{\text{linear}.}

For part (b), the coefficient of BB in standard linear form is p(t)=a(t)=t9p(t)=a(t)=t^9. The integrating factor is

ϕ(t)=e∫t9 dt=et10/10.\phi(t)=e^{\int t^9\,dt} =e^{t^{10}/10}.

Since ϕ′(t)=t9ϕ(t)\phi'(t)=t^9\phi(t), multiplying the standard equation by ϕ\phi gives

ϕB′+t9ϕB=ϕB′+ϕ′B=ddt(ϕB)=ϕ(t9G(t)+3cos⁡t).\begin{aligned} \phi B'+t^9\phi B &=\phi B'+\phi'B\\ &=\frac{d}{dt}(\phi B)\\ &=\phi\bigl(t^9G(t)+3\cos t\bigr). \end{aligned}

Thus,

ϕ(t)=exp⁡(t10/10).{\phi(t)=\exp(t^{10}/10).}

The general solution includes an integration constant:

B(t)=e−t10/10(∫et10/10(t9G(t)+3cos⁡t) dt+C).B(t) =e^{-t^{10}/10} \left( \int e^{t^{10}/10} \bigl(t^9G(t)+3\cos t\bigr)\,dt+C \right).

For part (c), first move all terms involving NN to the left:

7N′′+70N′+175N=8e−5t.7N''+70N'+175N=8e^{-5t}.

The associated homogeneous equation is

7N′′+70N′+175N=0.7N''+70N'+175N=0.

Using the trial solution N=eνtN=e^{\nu t} gives

7ν2eνt+70νeνt+175eνt=0.7\nu^2e^{\nu t}+70\nu e^{\nu t}+175e^{\nu t}=0.

Since eνt≠0e^{\nu t}\ne0, the characteristic equation is

0=7ν2+70ν+175=7(ν+5)2.{0=7\nu^2+70\nu+175 =7(\nu+5)^2.}

The repeated root is ν=−5\nu=-5, so the homogeneous solution is

Nh=(C1+C2t)e−5t.N_h=(C_1+C_2t)e^{-5t}.

The forcing term 8e−5t8e^{-5t} has the same exponential as the repeated characteristic root. Both the usual trial Ce−5tCe^{-5t} and the once-corrected trial Cte−5tCte^{-5t} belong to the homogeneous solution space, so use

Np=At2e−5t.N_p=At^2e^{-5t}.

Write the differential operator as

7D2+70D+175=7(D+5)2,D=ddt.7D^2+70D+175=7(D+5)^2, \qquad D=\frac{d}{dt}.

For Np=At2e−5tN_p=At^2e^{-5t},

(D+5)(At2e−5t)=2At e−5t,(D+5)\bigl(At^2e^{-5t}\bigr)=2At\,e^{-5t},

and hence

(D+5)2(At2e−5t)=2Ae−5t.(D+5)^2\bigl(At^2e^{-5t}\bigr) =2A e^{-5t}.

Substitution into the left-hand side gives

7(D+5)2Np=14Ae−5t.7(D+5)^2N_p=14A e^{-5t}.

Matching this with 8e−5t8e^{-5t} gives

14A=8⟹A=47.14A=8 \quad\Longrightarrow\quad A=\frac47.

Therefore

Np(t)=47t2e−5t.{N_p(t)=\frac47t^2e^{-5t}.}

The complete general solution is

N(t)=(C1+C2t+47t2)e−5t.{ N(t)=\left(C_1+C_2t+\frac47t^2\right)e^{-5t}. }

The first model is a first-order linear equation: the integrating factor creates the product rule. The second uses the constant-coefficient method for a second-order equation. Its repeated characteristic root explains the factors 1,t1,t in the homogeneous solution and the extra factor t2t^2 in the particular trial.

On this page