Differential equations

$\newenvironment {prompt}{}{} \newcommand {\todo }[0]{} \newcommand {\mooculus }[0]{\textsf {\textbf {MOOC}\textnormal {\textsf {ULUS}}}} \newcommand {\npnoround }[0]{\nprounddigits {-1}} \newcommand {\npnoroundexp }[0]{\nproundexpdigits {-1}} \newcommand {\npunitcommand }[1]{\ensuremath {\mathrm {#1}}} \newcommand {\CancelColor }[0]{} \newcommand {\cancelto }[0]{1} \newcommand {\RR }[0]{\mathbb R} \newcommand {\R }[0]{\mathbb R} \newcommand {\N }[0]{\mathbb N} \newcommand {\Z }[0]{\mathbb Z} \newcommand {\sagemath }[0]{\textsf {SageMath}} \newcommand {\d }[0]{\mathop {}\!d} \newcommand {\l }[0]{\ell } \newcommand {\ddx }[0]{\frac {d}{\d x}} \newcommand {\zeroOverZero }[0]{\ensuremath {\boldsymbol {\tfrac {0}{0}}}} \newcommand {\inftyOverInfty }[0]{\ensuremath {\boldsymbol {\tfrac {\infty }{\infty }}}} \newcommand {\zeroOverInfty }[0]{\ensuremath {\boldsymbol {\tfrac {0}{\infty }}}} \newcommand {\zeroTimesInfty }[0]{\ensuremath {\small \boldsymbol {0\cdot \infty }}} \newcommand {\inftyMinusInfty }[0]{\ensuremath {\small \boldsymbol {\infty -\infty }}} \newcommand {\oneToInfty }[0]{\ensuremath {\boldsymbol {1^\infty }}} \newcommand {\zeroToZero }[0]{\ensuremath {\boldsymbol {0^0}}} \newcommand {\inftyToZero }[0]{\ensuremath {\boldsymbol {\infty ^0}}} \newcommand {\numOverZero }[0]{\ensuremath {\boldsymbol {\tfrac {\#}{0}}}} \newcommand {\dfn }[0]{\textbf } \newcommand {\unit }[0]{\mathop {}\!\mathrm } \newcommand {\eval }[1]{\bigg [ #1 \bigg ]} \newcommand {\seq }[1]{\left ( #1 \right )} \newcommand {\epsilon }[0]{\varepsilon } \newcommand {\phi }[0]{\varphi } \newcommand {\iff }[0]{\Leftrightarrow } \DeclareMathOperator {\arccot }{arccot} \DeclareMathOperator {\arcsec }{arcsec} \DeclareMathOperator {\arccsc }{arccsc} \DeclareMathOperator {\si }{Si} \DeclareMathOperator {\proj }{\vec {proj}} \DeclareMathOperator {\scal }{scal} \DeclareMathOperator {\sign }{sign} \newcommand {\arrowvec }[0]{\overrightarrow } \newcommand {\vec }[0]{\mathbf } \newcommand {\veci }[0]{{\boldsymbol {\hat {\imath }}}} \newcommand {\vecj }[0]{{\boldsymbol {\hat {\jmath }}}} \newcommand {\veck }[0]{{\boldsymbol {\hat {k}}}} \newcommand {\vecl }[0]{\boldsymbol {\l }} \newcommand {\uvec }[1]{\mathbf {\hat {#1}}} \newcommand {\utan }[0]{\mathbf {\hat {t}}} \newcommand {\unormal }[0]{\mathbf {\hat {n}}} \newcommand {\ubinormal }[0]{\mathbf {\hat {b}}} \newcommand {\dotp }[0]{\bullet } \newcommand {\cross }[0]{\boldsymbol \times } \newcommand {\grad }[0]{\boldsymbol \nabla } \newcommand {\divergence }[0]{\grad \dotp } \newcommand {\curl }[0]{\grad \cross } \newcommand {\lto }[0]{\mathop {\longrightarrow \,}\limits } \newcommand {\bar }[0]{\overline } \newcommand {\surfaceColor }[0]{violet} \newcommand {\surfaceColorTwo }[0]{redyellow} \newcommand {\sliceColor }[0]{greenyellow} \newcommand {\vector }[1]{\left \langle #1\right \rangle } \newcommand {\sectionOutcomes }[0]{} \newcommand {\HyperFirstAtBeginDocument }[0]{\AtBeginDocument } \newcommand {\epstopdfsetup }[0]{\setkeys {ETE}} changed the theorem header of amsthm failed\MessageBreak }}} \newcommand {\GetTitleStringSetup }[0]{\setkeys {gettitlestring}} \newcommand {\Sectionformat }[2]{#1}$ $% This file was *autogenerated* from digInDifferentialEquations.sagetex.sage with % sagetex.py version 2015/08/26 v3.0-92d9f7a %42b10ab9e218196906c11244c82c00ab% md5sum of corresponding .sage file (minus "goboom","current_tex_line",and pause/unpause lines)$

Differential equations show you relationships between rates of functions.

A differential equation is simply an equation with a derivative in it. Here is an example: $a\cdot f''(x) + b\cdot f'(x) + c\cdot f(x) = g(x).$

What is a differential equation?

When a mathematician solves a differential equation, they are finding functions satisfying the equation.

Which of the following functions solve the differential equation $f^{(4)}(x) = f(x)?$

Remember, $f^{(4)}$ is the fourth derivative of $f$ .

It turns out that the complete solution to this differential equation is $c_1\sin (x)+c_2\cos (x)+c_3e^x+c_4e^{-x}$ . In other words, every solution to this differential equation can be written in this form. You should check that these are all solutions (for example $f(x) = \sin (x)+3\cos (x)-7e^x+\pi e^{-x}$ is a solution). Proving that these are all of the solutions is beyond the scope of this course.

The differential equation above is an example of a fourth order differential equation, because the highest derivative in the equation is a ‘‘fourth’’ derivative. In general the highest derivative in a differential equation is the order.

Differential equations are one of the most practical objects of mathematical study. They appear constantly in every field of science and engineering. They are a powerful way to model many diverse situations.

Modeling with differential equations

Setting up differential equations is a skill to be acquired. However, you can try your hand with our next question.

Imagine that a glass of water has initial temperature $5^\circ \unit {C}$ , and that the ambient temperature is $22^\circ \unit {C}$ . The water will warm up over time. Assume that the rate of change in the temperature of the water is directly proportional to the difference between the current water temperature and the ambient temperature. Which of the following differential equations must be satisfied by the function $H(t)$ which measures the temperature of the water with respect to time?

This is just a straight translation job. ‘‘The rate of change in the temperature of the water’’ is $y'$ . ‘‘Directly proportional to’’ means that it is equal to some constant (say $k$ ) times whatever it is proportional to. ‘‘The difference between the current water temperature and the ambient temperature” is either $22-y$ or $y-22$ , since $y$ is the temperature of the water and $22$ is the ambient temperature. Think about which we should choose before looking at the next hint. Will it be $y'=k(22-y)$ or $y'=k(y-22)$ where $k>0$ ?

Since the temperature of the water is increasing over time, we want $y'>0$ . Since the temperature will be increasing (and it is reasonable to assume it never surpasses the ambient temperature!) $22-y$ is positive. So we can conclude that $y' = k(22-y)$ for some $k>0$ .

The differential equation does not involve the number $5$ . If we wanted to incorporate that piece of data into our model we could ask ‘‘Which solution(s) to this differential equation satisfy $H(0) = 5$ ?’’ This is known as an initial value problem.

Sometimes the rates in question are constant.

One can approximate the force of gravity as constant near the Earth. So the acceleration of a falling object is a constant $g>0$ . If $h(t)$ is the height of an object at time $t$ , which differential equation must $h$ satisfy?

The acceleration of an object is the second derivative of its position, so the differential equation should say the second derivative, $y''$ , is constant. Should it be a positive of negative constant?

A falling object will fall quicker and quicker, so the second derivative of its height should be negative. Thus $y''=-g$ is the correct answer.

Initial value problems

We have already seen, and solved, a particular kind of differential equation in this course. Namely a solution $F$ to the differential equation $y' = f(x)$ is just an antiderivative of $f$ ! We know the ‘‘general solution’’ of this differential equation is just $F(x)+C$ , as long as the domain of $f$ is an interval. We can use this idea to solve differential equations of the form $y^{(n)} = f(x)$ , by just repeatedly integrating and solving for ‘‘ $C$ ’’ when we can.

Find the general solution to the differential equation $y''(x) = x$ .

Since $y''(x) = x$ , we know that $\begin{align*} y'(x) &= \int x \d x\\ &= \answer[given]{\frac{x^2}{2}} + C_1. \end{align*}$

This further implies that $y(x) = \int \frac {x^2}{2} + C_1 \d x,$ so we must have that $y(x) = \answer [given]{\frac {x^3}{6}}+C_1x+C_2,$ for some constant $C_2$ . This is the general solution of the differential equation, in fact, every solution of this differential equation is of the form $y(x) = \frac {1}{6}x^3+C_1x+C_2$ .

We can use the general solution to give specific solutions.

Find the solution to the differential equation $y''(x) = x$ that passes through the points $(0,1)$ and $(1,2)$ .

From our work above we know that $y(x) = \frac {x^3}{6}+C_1x+C_2.$ To find the particular solution we are interested in, we solve the system of equations $\begin{align*} \answer[given]{1}= y(\answer[given]{0}) &= \frac{1}{6}\cdot 0^3+C_1\cdot 0+C_2,\\ \answer[given]{2}= y(\answer[given]{1}) &= \frac{1}{6}\cdot 1^3+C_1\cdot 1+C_2. \end{align*}$

The first line tells us that $C_2 = \answer [given]{1}$ , and now the second line is $2= \frac {1}{6}+C_1+1.$ and so $C_1 = \answer [given]{\frac {5}{6}}$ . Our solution is now $y(x) = \answer [given]{\frac {x^3}{6}+\frac {5 x}{6}+1}.$

Press...	...to do
left/right arrows	Move cursor
shift+left/right arrows	Select region
ctrl+a	Select all
ctrl+x/c/v	Cut/copy/paste
ctrl+z/y	Undo/redo
ctrl+left/right	Add entry to list or column to matrix
shift+ctrl+left/right	Add copy of current entry/column to to list/matrix
ctrl+up/down	Add row to matrix
shift+ctrl+up/down	Add copy of current row to matrix
ctrl+backspace	Delete current entry in list or column in matrix
ctrl+shift+backspace	Delete current row in matrix

Type...	...to get
norm	$\|\|\blue{[?]}\|\|$
text	$\text{\blue{[?]}}$
sym_name	$\backslash\texttt{\blue{[?]}}$
abs	$\left\|\blue{[?]}\right\|$
sqrt	$\sqrt{\blue{[?]}}$
paren	$\left(\blue{[?]}\right)$
floor	$\lfloor \blue{[?]} \rfloor$
factorial	$\blue{[?]}!$
exp	${\blue{[?]}}^{\blue{[?]}}$
sub	${\blue{[?]}}_{\blue{[?]}}$
frac	$\dfrac{\blue{[?]}}{\blue{[?]}}$
int	$\displaystyle\int{\blue{[?]}}d\blue{[?]}$
defi	$\displaystyle\int_{\blue{[?]}}^{\blue{[?]}}\blue{[?]}d\blue{[?]}$
deriv	$\displaystyle\frac{d}{d\blue{[?]}}\blue{[?]}$
sum	$\displaystyle\sum_{\blue{[?]}}^{\blue{[?]}}\blue{[?]}$
prod	$\displaystyle\prod_{\blue{[?]}}^{\blue{[?]}}\blue{[?]}$
root	$\sqrt[\blue{[?]}]{\blue{[?]}}$
vec	$\left\langle \blue{[?]} \right\rangle$
mat	$\left(\begin{matrix} \blue{[?]} \end{matrix}\right)$
*	$\cdot$
infinity	$\infty$
arcsin	$\arcsin\left(\blue{[?]}\right)$
arccos	$\arccos\left(\blue{[?]}\right)$
arctan	$\arctan\left(\blue{[?]}\right)$
sin	$\sin\left(\blue{[?]}\right)$
cos	$\cos\left(\blue{[?]}\right)$
tan	$\tan\left(\blue{[?]}\right)$
sec	$\sec\left(\blue{[?]}\right)$
csc	$\csc\left(\blue{[?]}\right)$
cot	$\cot\left(\blue{[?]}\right)$
log	$\log\left(\blue{[?]}\right)$
ln	$\ln\left(\blue{[?]}\right)$
alpha	$\alpha$
beta	$\beta$
gamma	$\gamma$
delta	$\delta$
epsilon	$\epsilon$
zeta	$\zeta$
eta	$\eta$
theta	$\theta$
iota	$\iota$
kappa	$\kappa$
lambda	$\lambda$
mu	$\mu$
nu	$\nu$
xi	$\xi$
omicron	$\omicron$
pi	$\pi$
rho	$\rho$
sigma	$\sigma$
tau	$\tau$
upsilon	$\upsilon$
phi	$\phi$
chi	$\chi$
psi	$\psi$
omega	$\omega$
Gamma	$\Gamma$
Delta	$\Delta$
Theta	$\Theta$
Lambda	$\Lambda$
Xi	$\Xi$
Pi	$\Pi$
Sigma	$\Sigma$
Phi	$\Phi$
Psi	$\Psi$
Omega	$\Omega$

Modeling with differential equations

Initial value problems

Controls

Symbols

Settings