G6 - Constructing differential equations

G6 - Constructing differential equations

Translate rate and proportionality language into differential equations, calibrate constants from measurements and state initial conditions. Check signs, units and assumptions, and verify a proposed function by differentiation.

Translating rate and proportionality language

Can an equation describe a curve before we know its formula? Suppose the curve passes through (0,3)(0,3) and its tangent gradient at (x,y)(x,y) is 2x+y2x+y. The statement already tells us how the curve changes: dy/dx=2x+ydy/dx=2x+y.

You already know

A derivative is a rate of change: dy/dxdy/dx is the rate of yy with respect to xx. A tangent gradient need not involve time. The derivative interpretation is recalled, rather than supplied as a formula.

A differential equation involves a derivative of an unknown function. An initial condition specifies the function's value at a starting input, such as y(0)=3y(0)=3. It is different from y′(0)=3y'(0)=3, which specifies an initial gradient.

“Proportional to” introduces a constant multiplier. For example, a mass whose loss rate is proportional to its square root obeys dm/dt=−kmdm/dt=-k\sqrt m, with k>0k>0 and m≥0m\ge0. The minus sign gives a decreasing mass, while kmk\sqrt m is the positive magnitude of the loss. This translation is not supplied in the formula booklet.

Translating proportionality

If a signed rate is proportional to FF, write

rate=kF\text{rate}=kF

State the sign of kk, or use a positive kk and an explicit minus sign for a decrease. Not in the formula booklet: learn it.

Write the changing quantity on top of the derivative and the independent variable below it. Then translate the rate statement before using any numerical data.

Worked example 1: a curve's slope law

Question. A curve passes through (0,3)(0,3). Its gradient at (x,y)(x,y) is 2x+y2x+y. Write a differential equation and initial condition, then find the initial gradient.

Plan. Replace the word gradient by dy/dxdy/dx and substitute the starting point.

dydx=2x+ytranslate the gradienty(0)=3initial pointy′(0)=2(0)+3=3evaluate the rate\begin{aligned} \frac{dy}{dx}&=2x+y && \text{translate the gradient}\\ y(0)&=3 && \text{initial point}\\ y'(0)&=2(0)+3=3 && \text{evaluate the rate} \end{aligned}

Why keep the condition separate? The equation gives a rule for gradients; the condition says which point the curve passes through.

Answer. dy/dx=2x+ydy/dx=2x+y, with y(0)=3y(0)=3; the initial gradient is 33.

Check. Substituting the starting point into the right-hand side gives the stated initial gradient without finding the whole curve.

Change only the sign of the yy-term in the slope law. The initial point stays fixed, but the initial direction changes.

Worked example 2: changing the slope law

Question. A curve passes through (0,3)(0,3). Its gradient at (x,y)(x,y) is 2x−y2x-y. Write the equation and initial condition, then find the initial gradient.

Plan. Translate the new gradient expression and evaluate it at the same point.

dydx=2x−ynew rate lawy(0)=3same initial pointy′(0)=2(0)−3=−3initial rate\begin{aligned} \frac{dy}{dx}&=2x-y && \text{new rate law}\\ y(0)&=3 && \text{same initial point}\\ y'(0)&=2(0)-3=-3 && \text{initial rate} \end{aligned}

Why is the initial gradient negative? At the starting point, 2x=02x=0 and −y=−3-y=-3.

Answer. dy/dx=2x−ydy/dx=2x-y, with y(0)=3y(0)=3; the initial gradient is −3-3.

Check. Changing a sign in the law has not changed the condition y(0)=3y(0)=3.

The same translation works in a physical setting. A tank of height 2.4 m2.4\,\mathrm m, initially containing water to depth 0.3 m0.3\,\mathrm m, has a depth rate proportional to the remaining height if

dhdt=k(2.4−h),k>0,h(0)=0.3\frac{dh}{dt}=k(2.4-h),\quad k>0,\qquad h(0)=0.3

Here tt is in minutes, hh is in metres and kk has units min−1\mathrm{min^{-1}}. For 0≤h<2.40\le h<2.4, the depth increases, more slowly as it approaches the top. The law is a modelling assumption, not a supplied formula.

Watch out

Wrong: “y′(0)=3y'(0)=3” as the translation of “the curve passes through (0,3)(0,3)”. That describes a rate instead of a point. Right: y(0)=3y(0)=3; find any initial gradient from the differential equation.

Change only the coefficient of yy in example 2. Construct the equation; no solution for y(x)y(x) is needed.

Population models and measured rates

A measured quantity and a measured rate can determine a proportionality constant before the differential equation is solved.

If population NN decreases at a rate proportional to NN, write dN/dt=−kNdN/dt=-kN, with k>0k>0. For N=8000N=8000 and a decrease of 120120 fish per year, −120=−8000k-120=-8000k, giving k=0.015 year−1k=0.015\,\mathrm{year^{-1}}. The calibrated law is dN/dt=−0.015NdN/dt=-0.015N, so at N=4000N=4000 it predicts a decrease of 6060 fish per year. A proportional decrease is not a fixed loss at every population.

Growth can depend on more than one factor. “Proportional both to population and to unused capacity” means multiply the factors, rather than add them.

A limited-growth rate law

For constant capacity KK and positive kk, the stated proportionality gives

dNdt=kN(K−N)\frac{dN}{dt}=kN(K-N)

For 0<N<K0<N<K, growth is positive. At N=0N=0 and N=KN=K, the rate is zero. If the model admits N>KN>K, its predicted rate is negative. Not in the formula booklet: learn it.

This is a continuous approximation to a count. Constant capacity and constant proportionality are assumptions; changing resources or migration can make the law unreliable.

Worked example 3: two proportional factors

Question. A lake has capacity 10001000 fish. Its growth rate is proportional to both NN and 1000−N1000-N. Initially N=200N=200, growing at 8080 fish per week. Construct and calibrate the equation, and predict the rate at N=500N=500.

Plan. Multiply the two factors, use the measured rate to find kk, then substitute the second population.

dNdt=kN(1000−N),N(0)=200law and condition80=k(200)(800)initial rate datak=12000divide by 160000dNdt∣N=500=500(500)2000=125predict the rate\begin{aligned} \frac{dN}{dt}&=kN(1000-N),\quad N(0)=200 && \text{law and condition}\\ 80&=k(200)(800) && \text{initial rate data}\\ k&=\frac1{2000} && \text{divide by }160000\\ \left.\frac{dN}{dt}\right|_{N=500}&=\frac{500(500)}{2000}=125 && \text{predict the rate} \end{aligned}

Why multiply NN and 1000−N1000-N? The rate is proportional to both quantities together.

Answer. dN/dt=N(1000−N)/2000dN/dt=N(1000-N)/2000, with N(0)=200N(0)=200; the predicted rate is 125125 fish per week. Here kk has units fish−1 week−1\mathrm{fish^{-1}\,week^{-1}}.

Check. At N=200N=200, the calibrated law gives 200(800)/2000=80200(800)/2000=80, as measured.

Change only the measured growth rate from 8080 to 160160 fish per week. This tests whether the fitted constant responds proportionally to the data.

Worked example 4: recalibrating the same model

Question. Use the same capacity and proportionality law, with N(0)=200N(0)=200, but an initial growth rate of 160160 fish per week. Find the equation and the predicted rate at N=500N=500.

Plan. Recalculate kk from the new rate, leaving both population factors unchanged.

160=k(200)(800)new measured ratek=11000calibratedNdt∣N=500=500(500)1000=250predict\begin{aligned} 160&=k(200)(800) && \text{new measured rate}\\ k&=\frac1{1000} && \text{calibrate}\\ \left.\frac{dN}{dt}\right|_{N=500}&=\frac{500(500)}{1000}=250 && \text{predict} \end{aligned}

Why has kk doubled? The population factors are unchanged, so doubling the measured rate doubles their multiplier.

Answer. dN/dt=N(1000−N)/1000dN/dt=N(1000-N)/1000, with N(0)=200N(0)=200; the predicted rate is 250250 fish per week.

Check. The rate at 500500 has doubled from 125125 to 250250, as has the calibration rate.

There is a useful surprise in the product law. Populations 200200 and 800800 have the same rate, since 200(800)=800(200)200(800)=800(200). More generally,

N(1000−N)=250000−(N−500)2N(1000-N)=250000-(N-500)^2

so the largest rate occurs at half capacity, N=500N=500, without solving for N(t)N(t). This is a completed-square identity, not a supplied formula.

Watch out

Wrong: dN/dt=k[N+(1000−N)]=1000kdN/dt=k[N+(1000-N)]=1000k. Adding the proportional quantities turns the rate into a constant. Right: multiply them, dN/dt=kN(1000−N)dN/dt=kN(1000-N).

Change only the initial measured rate again. Check it back in your final equation.

Continue this lesson.

You’re previewing NovaLearn without an account. Create a free account to read the rest.

Create free account