Skip to main content
Mathematics LibreTexts

4.4.1: Autonomous Second Order Equations (Exercises)

  • Page ID
    30718
  • \( \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}\)

    Q4.4.1

    In Exercises 4.4.1-4.4.4 find the equations of the trajectories of the given undamped equation. Identify the equilibrium solutions, determine whether they are stable or unstable, and plot some trajectories.

    1. \(y''+y^3=0\)

    2. \(y''+y^2=0\)

    3. \(y''+y|y|=0\)

    4. \(y''+ye^{-y}=0\)

    Q4.4.2

    In Exercises 4.4.5–4.4.8 find the equations of the trajectories of the given undamped equation. Identify the equilibrium solutions, determine whether they are stable or unstable, and find the equations of the separatrices (that is, the curves through the unstable equilibria). Plot the separatrices and some trajectories in each of the regions of Poincaré plane determined by them.

    5. \(y''-y^3+4y=0\)

    6. \(y''+y^3-4y=0\)

    7. \(y''+y(y^2-1)(y^2-4)=0\)

    8. \(y''+y(y-2)(y-1)(y+2)=0\)

    Q4.4.3

    In Exercises 4.4.9–4.4.12 plot some trajectories of the given equation for various values (positive, negative, zero) of the parameter a. Find the equilibria of the equation and classify them as stable or unstable. Explain why the phase plane plots corresponding to positive and negative values of a differ so markedly. Can you think of a reason why zero deserves to be called the critical value of \(a\)?

    9. \(y''+y^2-a=0\)

    10. \(y''+y^3-ay=0\)

    11. \(y''-y^3+ay=0\)

    12. \(y''+y-ay^3=0\)

    Q4.4.4

    In Exercises 4.4.13-4.4.18 plot trajectories of the given equation for \(c = 0\) and small nonzero (positive and negative) values of \(c\) to observe the effects of damping.

    13. \(y''+cy'+y^3=0\)

    14. \(y''+cy'-y=0\)

    15. \(y''+cy'+y^3=0\)

    16. \(y''+cy'+y^2=0\)

    17. \(y''+cy'+y|y|=0\)

    18. \(y''+y(y-1)+cy=0\)

    Q4.4.5

    19. The van der Pol equation

    \[y''-\mu(1-y^2)y'+y=0, \tag{A} \]

    where \(\mu\) is a positive constant and \(y\) is electrical current (Section 6.3), arises in the study of an electrical circuit whose resistive properties depend upon the current. The damping term \(-\mu(1-y^2)y'\) works to reduce \(|y|\) if \(|y|<1\) or to increase \(|y|\) if \(|y|>1\). It can be shown that van der Pol’s equation has exactly one closed trajectory, which is called a limit cycle. Trajectories inside the limit cycle spiral outward to it, while trajectories outside the limit cycle spiral inward to it (Figure [figure:4.4.16}). Use your favorite differential equations software to verify this for \(\mu=.5,1.1.5,2\). Use a grid with \(-4<y<4\) and \(-4<v<4\).

    fig040416.svg
    Figure 4.4E.16 : Trajectories of van der Pol's euqation

    20. Rayleigh’s equation,

    \[y''-\mu(1-(y')^2/3)y'+y=0 \nonumber \]

    also has a limit cycle. Follow the directions of Exercise 4.4.19 for this equation.

    21. In connection with Equation 4.4.16, suppose \(y(0)=0\) and \(y'(0)=v_0\), where \(0<v_0<v_c\).

    1. Let \(T_1\) be the time required for \(y\) to increase from zero to \(y_{\max}=2\sin^{-1}(v_0/v_c)\). Show that \[{dy\over dt}=\sqrt{v_0^2-v_c^2\sin^2y/2},\quad 0\le t<T_1. \tag{A} \]
    2. Separate variables in (A) and show that \[T_1=\int_0^{y_{\max}}{du\over\sqrt{v_0^2-v_c^2\sin^2u/2}} \tag{B} \]
    3. Substitute \(\sin u/2=(v_0/v_c)\sin\theta\) in (B) to obtain \[T_1=2\int_0^{\pi/2}{d\theta\over\sqrt{v_c^2-v_0^2\sin^2\theta}}. \tag{C} \]
    4. Conclude from symmetry that the time required for \((y(t),v(t))\) to traverse the trajectory \[v^2=v_0^2-v_c^2\sin^2y/2 \nonumber \] is \(T=4T_1\), and that consequently \(y(t+T)=y(t)\) and \(v(t+T)=v(t)\); that is, the oscillation is periodic with period \(T\).
    5. Show that if \(v_0=v_c\), the integral in (C) is improper and diverges to \(\infty\). Conclude from this that \(y(t)<\pi\) for all \(t\) and \(\lim_{t\to\infty}y(t)=\pi\).

    22. Give a direct definition of an unstable equilibrium of \(y''+p(y)=0\).

    23. Let \(p\) be continuous for all \(y\) and \(p(0)=0\). Suppose there’s a positive number \(\rho\) such that \(p(y)>0\) if \(0<y\le \rho\) and \(p(y)<0\) if \(-\rho\le y<0\). For \(0<r\le\rho\) let

    \[\alpha(r)=\min\left\{\int_0^r p(x)\,dx,\ \int_{-r}^0 |p(x)|\,dx\right\} \mbox{\quad and \quad} \beta(r)=\max\left\{\int_0^r p(x)\,dx,\ \int_{-r}^0 |p(x)|\,dx\right\}. \nonumber \]

    Let \(y\) be the solution of the initial value problem

    \[y''+p(y)=0,\quad y(0)=v_0,\quad y'(0)=v_0, \nonumber \]

    and define \(c(y_0,v_0)=v_0^2+2\int_0^{y_0}p(x)\,dx\).

    1. Show that \[0<c(y_0,v_0) <v_0^2+2\beta(|y_0|)\mbox{\quad if \quad} 0<|y_0|\le\rho. \nonumber \]
    2. Show that \[v^2+2\int_0^y p(x)\,dx=c(y_0,v_0),\quad t>0. \nonumber \]
    3. Conclude from (b) that if \(c(y_0,v_0)<2\alpha(r)\) then \(|y|<r,\ t>0\).
    4. Given \(\epsilon>0\), let \(\delta>0\) be chosen so that \[\delta^2+2\beta(\delta)<\max\left\{\epsilon^2/2,2\alpha(\epsilon/\sqrt2) \right\}. \nonumber \] Show that if \(\sqrt{y_0^2+v_0^2}<\delta\) then \(\sqrt{y^2+v^2}<\epsilon\) for \(t>0\), which implies that \(\overline y=0\) is a stable equilibrium of \(y''+p(y)=0\).
    5. Now let \(p\) be continuous for all \(y\) and \(p(\overline y)=0\), where \(\overline y\) is not necessarily zero. Suppose there’s a positive number \(\rho\) such that \(p(y)>0\) if \(\overline y<y\le \overline y+\rho\) and \(p(y)<0\) if \(\overline y-\rho\le y<\overline y\). Show that \(\overline y\) is a stable equilibrium of \(y''+p(y)=0\).

    24. Let \(p\) be continuous for all \(y\).

    1. Suppose \(p(0)=0\) and there’s a positive number \(\rho\) such that \(p(y)<0\) if \(0<y\le \rho\). Let \(\epsilon\) be any number such that \(0<\epsilon<\rho\). Show that if \(y\) is the solution of the initial value problem \[y''+p(y)=0,\quad y(0)=y_0,\quad y'(0)=0 \nonumber \] with \(0<y_0<\epsilon\), then \(y(t)\ge\epsilon\) for some \(t>0\). Conclude that \(\overline y=0\) is an unstable equilibrium of \(y''+p(y)=0\).
    2. Now let \(p(\overline y)=0\), where \(\overline y\) isn’t necessarily zero. Suppose there’s a positive number \(\rho\) such that \(p(y)<0\) if \(\overline y<y\le \overline y+\rho\). Show that \(\overline y\) is an unstable equilibrium of \(y''+p(y)=0\).
    3. Modify your proofs of (a) and (b) to show that if there’s a positive number \(\rho\) such that \(p(y)>0\) if \(\overline y-\rho\le y<\overline y\), then \(\overline y\) is an unstable equilibrium of \(y''+p(y)=0\).

    This page titled 4.4.1: Autonomous Second Order Equations (Exercises) is shared under a CC BY-NC-SA 3.0 license and was authored, remixed, and/or curated by William F. Trench.

    • Was this article helpful?