1.1: Basic Concepts
- Page ID
- 103463
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)A differential equation is an equation that contains one or more derivatives of an unknown function.
A differential equation is an ordinary differential equation if it involves an unknown function of only one variable, or a partial differential equation if it involves partial derivatives of a function of more than one variable.
Note: Partial differential equations are studied in another class; in this class we’ll consider only ordinary differential equations, and we’ll just call them differential equations.
The order of a differential equation is the order of the highest derivative that it contains.
\begin{array}{rcll} {dy\over dx}-y=x&\mbox{ (first order)}, \\ (x^2-5xy)dx+7xy^2dy=0&\mbox{ (first order)}, \\ {d^2y\over dx^2}+7x {dy\over dx}+3y=5x&\mbox{ (second order)}, \\ x^3y'''+7xy^2=\cos x &\mbox{ (third order)}, \\ 9x^3y^{(n)}+3x^2y'+7y=e^x&\mbox{ (n-th order)}. \nonumber \end{array}
The simplest differential equations are first order equations of the form
\[{dy\over dx}=f(x) \nonumber\]
or equivalently
\[y'=f(x), \nonumber\]
where \(f\) is a known function of \(x\). We already know from calculus how to find functions that satisfy this kind of equation. For example, if
\[y'=x^3, \nonumber\]
then
\[y=\int x^3\, dx={x^4\over4}+c, \nonumber\]
where \(c\) is an arbitrary constant. If \(n>1\) we can find functions \(y\) that satisfy equations of the form
\[\label{eq:1.2.1} y^{(n)}=f(x)\]
by repeated integration. Again, this is a calculus problem.
Except for illustrative purposes in this section, there’s no need to consider differential equations like Equation \ref{eq:1.2.1} as they can usually be solved with calculus. We’ll usually consider differential equations that can be written as
\[\label{eq:1.2.2} y^{(n)}=f(x,y,y', \dots,y^{(n-1)}),\]
where at least one of the functions \(y\), \(y'\), …, \(y^{(n-1)}\) actually appears on the right when written in the form \ref{eq:1.2.2}. Here are some examples:
\begin{array}{rcll} {dy\over dx}-x^2y=0&\mbox{ (first order)}, \\ {dy\over dx}+2xy^2=-2&\mbox{ (first order)}, \\ {d^2y\over dx^2}+2 {dy\over dx}+y=2x&\mbox{ (second order)}, \\ xy'''+y^2=\sin x &\mbox{ (third order)}, \\ y^{(n)}+xy'+3y=x&\mbox{ (n-th order)}. \nonumber \end{array}
Although none of these equations is written as in Equation \ref{eq:1.2.2}, all of them can be written in this form:
\[\begin{array}{rcl} y'=x^2y, \\ y'=-2-2xy^2, \\ y''=2x-2y'-y, \\ y'''= \dfrac{\sin x-y^2}{x}, \\[4pt] y^{(n)}=x-xy'-3y. \end{array}\nonumber \]
Solutions of Differential Equations
A solution of a differential equation is a function that satisfies the differential equation on some open interval; thus, \(y\) is a solution of Equation \ref{eq:1.2.2} if \(y\) is \(n\) times differentiable and
\[y^{(n)}(x)=f(x,y(x),y'(x), \dots,y^{(n-1)}(x)) \nonumber\]
for all \(x\) in some open interval \((a,b)\). In this case, we also say that \(y\) is a solution of Equation \ref{eq:1.2.2} on \((a,b)\).
A solution of a differential equation must be both continuous and differentiable on an open interval.
The graph of a solution of a differential equation is a solution curve. More generally, a curve \(C\) is said to be an integral curve of a differential equation if every function \(y=y(x)\) whose graph is a segment of \(C\) is a solution of the differential equation. Thus, any solution curve of a differential equation is an integral curve, but an integral curve need not be a solution curve (see Example 6).
Verify that
\[\label{eq:1.2.5} y={x^2\over3}+{1\over x}\]
is a solution of
\[\label{eq:1.2.6} xy'+y=x^2\]
on \((0,\infty)\) and on \((-\infty,0)\).
Solution
Substituting Equation \ref{eq:1.2.5} and
\[y'={2x\over3} - {1\over x^2}\nonumber \]
into Equation \ref{eq:1.2.6} yields
\[xy'(x)+y(x)=x \left({2x\over3} - {1\over x^2}\right)+ \left({x^2\over3}+{1\over x}\right)=x^2\nonumber \]
for all \(x\ne0\). Therefore \(y\) is a solution of Equation \ref{eq:1.2.6} on \((-\infty,0)\) and \((0,\infty)\). However, \(y\) isn’t a solution of the differential equation on any open interval that contains \(x=0\), since \(y\) is not defined at \(x=0\).
Figure 1.1.2 shows the graph of Equation \ref{eq:1.2.5}. The part of the graph of Equation \ref{eq:1.2.5} on \((0,\infty)\) is a solution curve of Equation \ref{eq:1.2.6}, as is the part of the graph on \((-\infty,0)\).
Show that if \(c_1\) and \(c_2\) are constants then
\[\label{eq:1.2.7} y=(c_1+c_2x)e^{-x}+2x-4\]
is a solution of \[\label{eq:1.2.8} y''+2y'+y=2x\] on \((-\infty,\infty)\).
Solution
Differentiating Equation \ref{eq:1.2.7} twice yields
\[y'=-(c_1+c_2x)e^{-x}+c_2e^{-x}+2 \nonumber\]
and
\[y''=(c_1+c_2x)e^{-x}-2c_2e^{-x}, \nonumber\]
so
\[\begin{align*} y''+2y'+y&=(c_1+c_2x)e^{-x}-2c_2e^{-x} + 2\left[-(c_1+c_2x)e^{-x}+c_2e^{-x}+2\right] +(c_1+c_2x)e^{-x}+2x-4 \\[4pt] &=(1-2+1)(c_1+c_2x)e^{-x}+(-2+2)c_2e^{-x} + 4+2x-4 \\[4pt] &=2x \end{align*}\]
for all values of \(x\). Therefore \(y\) is a solution of Equation \ref{eq:1.2.8} on \((-\infty,\infty)\).
Note: Functions that satisfy a differential equation at isolated points are neither interesting nor useful, and since they do not satisfy the equation on an open interval we do not consider them solutions.
For example, \(y=x^2\) satisfies
\[xy'+x^2=3x \nonumber\]
if and only if \(x=0\) or \(x=1\), but it is not a solution of this differential equation because it does not satisfy the equation on an open interval.
If \(a\) is any positive constant, the circle
\[\label{eq:1.2.3} x^2+y^2=a^2\]
is an integral curve of
\[\label{eq:1.2.4} y'=-{x\over y}.\]
Solution
To see this, note that the only functions whose graphs are segments of Equation \ref{eq:1.2.3} are
\[y_1=\sqrt{a^2-x^2} \quad \text{and} \quad y_2=-\sqrt{a^2-x^2}.\nonumber \]
We leave it to you to verify that these functions both satisfy Equation \ref{eq:1.2.4} on the open interval \((-a,a)\). However, Equation \ref{eq:1.2.3} is not a solution curve of Equation \ref{eq:1.2.4}, since it is not the graph of a function. However, both \(y_1\) and \(y_2\) are each solution curves.
Initial Value Problems
An nth order differential equation has an infinite family of solutions that depend upon the \(n\) arbitrary constants \(c_1\), \(c_2\), …, \(c_n\). In the absence of additional conditions, there’s no reason to prefer one solution of a differential equation over another. However, we’ll often be interested in finding a particular solution of a differential equation that satisfies one or more specific conditions.
Find a solution of \[y'=x^3 \nonumber\] such that \(y(1)=2\).
Solution
At the beginning of this section we saw that the family of solutions of the differential equation \(y'=x^3\) are
\[y={x^4\over4}+c. \nonumber\]
To determine a value of \(c\) such that \(y(1)=2\), we set \(x=1\) and \(y=2\) here to obtain
\[2=y(1)={1\over4}+c \nonumber\]
so
\[c={7\over4}. \nonumber\]
Therefore the required particular solution of the initial value problem is
\[y={x^4+7\over4}. \nonumber\]
Figure 1.1.2 shows the graph of this solution. Note that imposing the condition \(y(1)=2\) is equivalent to requiring the graph of \(y\) to pass through the point \((1,2)\).
A solution of a differential equation that cannot be obtained from a family of solutions by varying the parameters is called a singular solution.
\(y=({x^2\over4}+c)^2\) is a family of solutions to the differential equation \(y'=xy^{1/2}\) (verify this). But \(y=0\) is also a solution (easily verified) and it cannot be obtained from the family by varying the parameter c. Therefore \(y=0\) is a singular solution.
We can rewrite the initial value problem considered in Example 1.1.7 more briefly as
\[y'=x^3,\quad y(1)=2.\nonumber\]
The requirement \(y(1)=2\) is an initial condition. Initial value problems can also be posed for higher order differential equations. For example,
\[\label{eq:1.2.11} y'' - 2y'+3y=e^x, \quad y(0)=1, \quad y'(0)=2\]
is an initial value problem for a second order differential equation where \(y\) and \(y'\) are required to have specified values at \(x=0\). In general, an initial value problem for an \(n\)-th order differential equation requires \(y\) and its first \(n-1\) derivatives to have specified values at some point \(x_0\). These requirements are the initial conditions. We’ll denote an initial value problem for a differential equation by writing the initial conditions after the equation, as in Equation \ref{eq:1.2.11}. For example, we would write an initial value problem for Equation \ref{eq:1.2.2} as
\[\label{eq:1.2.12} y^{(n)}=f(x,y,y', \dots,y^{(n-1)}),\, y(x_0)=k_0,\, y'(x_0)=k_1,\, \dots,\, y^{(n-1)}=k_{n-1}.\]
Consistent with our earlier definition of a solution of the differential equation in Equation \ref{eq:1.2.12}, we say that \(y\) is a solution of the initial value problem Equation \ref{eq:1.2.12} if \(y\) is \(n\) times differentiable and
\[y^{(n)}(x)=f(x,y(x),y'(x), \dots,y^{(n-1)}(x))\nonumber\]
for all \(x\) in some open interval \((a,b)\) that contains \(x_0\), and \(y\) satisfies the initial conditions in Equation \ref{eq:1.2.12}.
The largest open interval that contains \(x_0\) on which \(y\) is defined and satisfies the differential equation is the interval of validity of \(y\).
In Example 1.1.7 we saw that
\[\label{eq:1.2.13} y={x^4+7\over4}\]
is a solution of the initial value problem
\[y'=x^3,\quad y(1)=2.\nonumber\]
Since the function in Equation \ref{eq:1.2.13} is defined for all \(x\), the interval of validity of this solution is \((-\infty,\infty)\).
In Example 1.1.3 we verified that
\[\label{eq:1.2.14} y={x^2\over3}+{1\over x}\]
is a solution of
\[xy'+y=x^2 \nonumber\]
on \((0,\infty)\) and on \((-\infty,0)\). By evaluating Equation \ref{eq:1.2.14} at \(x=\pm1\), you can see that Equation \ref{eq:1.2.14} is a solution of the initial value problems
\[\label{eq:1.2.15} xy'+y=x^2,\quad y(1)={4\over3}\]
and
\[\label{eq:1.2.16} xy'+y=x^2,\quad y(-1)=-{2\over3}.\]
The interval of validity of Equation \ref{eq:1.2.14} as a solution of Equation \ref{eq:1.2.15} is \((0,\infty)\), since this is the largest interval that contains \(x_0=1\) on which Equation \ref{eq:1.2.14} is defined. Similarly, the interval of validity of Equation \ref{eq:1.2.14} as a solution of Equation \ref{eq:1.2.16} is \((-\infty,0)\), since this is the largest interval that contains \(x_0=-1\) on which Equation \ref{eq:1.2.14} is defined.
Free Fall Under Constant Gravity
The term initial value problem originated in problems of motion where the independent variable is \(t\) (representing elapsed time), and the initial conditions are the position and velocity of an object at the initial (starting) time of an experiment, usually, but not always, 0.
An object falls under the influence of gravity near Earth’s surface, where it can be assumed that the magnitude of the acceleration due to gravity is a constant \(g\).
- Construct a mathematical model for the motion of the object in the form of an initial value problem for a second order differential equation, assuming that the altitude and velocity of the object at time \(t=0\) are known. Assume that gravity is the only force acting on the object.
- Solve the initial value problem derived above to obtain the altitude as a function of time.
Solution a
Let \(y(t)\) be the altitude of the object at time \(t\). Since the acceleration of the object has constant magnitude \(g\) and is in the downward (negative) direction, \(y\) satisfies the second order equation
\[y''=-g, \nonumber\]
where the prime now indicates differentiation with respect to \(t\). If \(y_0\) and \(v_0\) denote the altitude and velocity when \(t=0\), then \(y\) is a solution of the initial value problem
\[\label{eq:1.2.17} y''=-g,\quad y(0)=y_0,\quad y'(0)=v_0.\]
Solution b
Integrating Equation \ref{eq:1.2.17} twice yields the family of solutions to the differential equation \(y''=-g, \nonumber\)
\[\begin{aligned} y'&=-gt+c_1, \\ y&=-{gt^2\over2}+c_1t+c_2.\end{aligned} \]
Imposing the initial conditions \(y(0)=y_0\) and \(y'(0)=v_0\) in these two equations shows that \(c_1=v_0\) and \(c_2=y_0\). Therefore the particular solution of the initial value problem Equation \ref{eq:1.2.17} is
\[y=- {gt^2\over2}+v_0t+y_0. \nonumber\]