# 2.E: Higher order linear ODEs (Exercises)

These are homework exercises to accompany Libl's "Differential Equations for Engineering" Textmap. This is a textbook targeted for a one semester first course on differential equations, aimed at engineering students. Prerequisite for the course is the basic calculus sequence.

## 2.1: Second order linear ODEs

**Exercise 2.1.2**: Show that \( y = e^x\) and \( y = e^{2x} \) are linearly independent.

**Exercise 2.1.3**: Take \( y'' + 5y = 10x + 5\). Find (guess!) a solution.

**Exercise 2.1.4**: Prove the superposition principle for nonhomogeneous equations. Suppose that \(y_1\) is a solution to \( Ly_1 = f(x) \) and \( y_2\) is a solution to \( Ly_2 = g(x) \) (same linear operator \( L\)). Show that \( y = y_1 + y_2\) solves \( Ly = f(x) + g(x) \).

**Exercise 2.1.5**: For the equation \( x^2y'' - xy' = 0 \), find two solutions, show that they are linearly independent and find the general solution. Hint: Try \(y = x'\).

Note that equations of the form \( ax^2y'' + bxy' + cy = 0 \) are called Euler’s equations or Cauchy-Euler equations. They are solved by trying \( y = x^r\) and solving for \(r\) (we can assume that \( x \geqslant 0\) for simplicity).

**Exercise 2.1.6**: Suppose that \( {(b - a)}^2 - 4ac > 0 \). a) Find a formula for the general solution of \(ax^2y'' + bxy' + cy = 0\). Hint: Try \( y = x^r \) and find a formula for \(r\). b) What happens when \( {(b - a)}^2 -4ac = 0 \) or \( {(b - a)}^2 - 4ac < 0 \)?

We will revisit the case when \( {(b - a)}^2 - 4ac < 0 \) later.

**Exercise 2.1.7**: Same equation as in Exercise 2.1.6. Suppose \( {(b - a)}^2 - 4ac = 0 \). Find a formula for the general solution of \( ax^2y'' + bxy' + cy = 0 \). Hint: Try \( y = x^r \ln x \) for the second solution.

If you have one solution to a second order linear homogeneous equation you can find another one. This is the reduction of order method.

**Exercise 2.1.8** (reduction of order): Suppose \( y_1\) is a solution to \( y'' + p(x)y' + q(x)y = 0 \). Show that

\[ y_2 (x) = y_1 (x) \int \frac { e^{- \int p(x)dx}}{{(y_1(x))}^2} dx \]

is also a solution.

Note: If you wish to come up with the formula for reduction of order yourself, start by trying \( y_2 (x) = y_1(x)v(x)\). Then plug \(y_2\) into the equation, use the fact that \(y_1\) is a solution, substitute \(w = v'\), and you have a first order linear equation in \(w\). Solve for \(w\) and then for \(v\). When solving for \(w\), make sure to include a constant of integration. Let us solve some famous equations using the method.

**Exercise 2.1.9** (Chebyshev’s equation of order 1): Take \( (1 - x^2)y'' - xy' + y = 0 \). a) Show that \( y = x\) is a solution. b) Use reduction of order to find a second linearly independent solution. c) Write down the general solution.

**Exercise 2.1.10** (Hermite’s equation of order 2): Take \( y'' - 2xy' + 4y = 0\). a) Show that \( y = 1 - 2x^2\) is a solution. b) Use reduction of order to find a second linearly independent solution. c) Write down the general solution.

**Exercise 2.1.101**: Are \( \sin (x) \) and \( e^x\) linearly independent? Justify.

**Exercise 2.1.102**: Are \( e^x\) and \( e^{x+2}\) linearly independent? Justify.

**Exercise 2.1.103**: Guess a solution to \( y'' + y' + y = 5 \).

**Exercise 2.1.104**: Find the general solution to \( xy'' + y' = 0 \). Hint: Notice that it is a first order ODE in \(y'\).

## 2.2: Constant coefficient second order linear ODEs

**Exercise 2.2.6:** Find the general solution of \( 2y'' + 2y' - 4y = 0 \).

**Exercise 2.2.7:**Find the general solution of \(y'' + 9y' - 10y = 0 \).

**Exercise 2.2.8:** Solve \( y'' - 8y' + 16y = 0 \) for \(y(0) = 2, y'(0) = 0 \).

**Exercise 2.2.9:** Solve \(y'' + 9y' = 0\) for \( y(0) = 1, y'(0) = 1\).

**Exercise 2.2.10:** Find the general solution of \( 2y'' + 50y = 0 \).

**Exercise 2.2.11:** Find the general solution of \( y'' + 6y' + 13y = 0 \).

**Exercise 2.2.12:** Find the general solution of \( y'' = 0 \) using the methods of this section.

**Exercise 2.2.13:** The method of this section applies to equations of other orders than two. We will see higher orders later. Try to solve the first order equation \( 2y' + 3y = 0 \) using the methods of this section.

**Exercise 2.2.14:** Let us revisit Euler’s equations of Exercise 2.1.6. Suppose now that \( {(b - a)}^2 - 4ac < 0 \). Find a formula for the general solution of \( ax^2y'' + bxy' + cy = 0 \). Hint: Note that \( x^r = e^{r \ln x} \).

**Exercise 2.2.101:** Find the general solution to \(y'' + 4y' + 2y = 0 \).

**Exercise 2.2.102:** Find the general solution to \( y'' - 6y' + 9y = 0\).

**Exercise 2.2.103:** Find the solution to \( 2y'' + y' + y = 0, y(0) = 1, y'(0) = -2 \).

**Exercise 2.2.104: **Find the solution to \( 2y'' + y' - 3y = 0, y(0) = a, y'(0) = b \).

**Exercise 2.2.105:** Find the solution to \( z''(t) = -2z'(t) - 2z (t), z(0) = 2, z'(0) = -2 \).

## 2.3: Higher order linear ODEs

**Exercise 2.3.1:**Find the general solution for \( y''' - y'' + y' - y = 0 \).

**Exercise 2.3.2:** Find the general solution for \(y^{(4)} - 5y''' + 6y'' = 0\).

**Exercise 2.3.3:** Find the general solution for \( y''' + 2y'' + 2y' = 0 \).

**Exercise 2.3.4:** Suppose the characteristic equation for a differential equation is \( {(r - 1)}^2{(r -2)}^2 = 0 \). a) Find such a differential equation. b) Find its general solution.

**Exercise 2.3.5:** Suppose that a fourth order equation has a solution \( y = 2e^{4x} x \cos x \). a) Find such an equation. b) Find the initial conditions that the given solution satisfies.

**Exercise 2.3.6:** Find the general solution for the equation of Exercise 2.3.5.

**Exercise 2.3.7:** Let \( f(x) = e^x - \cos x, g(x) = e^x + \cos x \) and \( h(x) = \cos x \). Are \( f(x), g(x), and h(x) \) and linearly independent? If so, show it, if not, find a linear combination that works.

**Exercise 2.3.8:** Let \(f(x) = 0, g(x) = \cos x \), and \( h(x) = \sin x \). Are \( f(x), g(x), and h(x) \) and linearly independent? If so, show it, if not, find a linear combination that works.

**Exercise 2.3.9:** Are \( x, x^2, and x^4 \) linearly independent? If so, show it, if not, find a linear combination that works.

**Exercise 2.3.10:** Are \( e^x, xe^x, and x^2e^x\) linearly independent? If so, show it, if not, find a linear combination that works.

**Exercise 2.3.101:** Find the general solution of \(y^{(5)} - y^{(4)} = 0 \).

**Exercise 2.3.102:** Suppose that the characteristic equation of a third order differential equation has roots \(3 \pm 2i\). a) What is the characteristic equation? b) Find the corresponding differential equation. c) Find the general solution.

**Exercise 2.3.103:** Solve \( 1001y''' + 3.2y'' + \pi y' - \sqrt {4} y = 0, y(0) = 0, y' (0) = 0, y'' (0) = 0 \).

**Exercise 2.3.104:** Are \(e^x, e^{x+1}, e^{2x}, \sin (x) \) linearly independent? If so, show it, if not find a linear combination that works.

**Exercise 2.3.105:** Are \( \sin (x), x, x \sin (x) \) linearly independent? If so, show it, if not find a linear combination that works.

## 2.4: Mechanical Vibrations

**Exercise 2.4.2**: Consider a mass and spring system with a mass \(m = 2\), spring constant \( k = 3\), and damping constant \(c = 1\). a) Set up and find the general solution of the system. b) Is the system underdamped, overdamped or critically damped? c) If the system is not critically damped, find a \(c\) that makes the system critically damped.

**Exercise 2.4.3**: Do Exercise 2.4.2 for \(m = 3, k = 12, and c = 12 \).

**Exercise 2.4.4**: Using the mks units (meters-kilograms-seconds), suppose you have a spring with spring constant \( 4 \dfrac {N}{m} \). You want to use it to weigh items. Assume no friction. You place the mass on the spring and put it in motion. a) You count and find that the frequency is 0.8 Hz (cycles per second). What is the mass? b) Find a formula for the mass \(m\) given the frequency \(w\) in Hz.

**Exercise 2.4.5**: Suppose we add possible friction to Exercise 2.4.4. Further, suppose you do not know the spring constant, but you have two reference weights 1 kg and 2 kg to calibrate your setup. You put each in motion on your spring and measure the frequency. For the 1 kg weight you measured 1.1 Hz, for the 2 kg weight you measured 0.8 Hz. a) Find \(k\) (spring constant) and \(c\) (damping constant). b) Find a formula for the mass in terms of the frequency in Hz. Note that there may be more than one possible mass for a given frequency. c) For an unknown object you measured 0.2 Hz, what is the mass of the object? Suppose that you know that the mass of the unknown object is more than a kilogram.

**Exercise 2.4.6**: Suppose you wish to measure the friction a mass of 0.1 kg experiences as it slides along a floor (you wish to find \(c\)). You have a spring with spring constant \( k = 5 \dfrac {N}{m} \). You take the spring, you attach it to the mass and fix it to a wall. Then you pull on the spring and let the mass go. You find that the mass oscillates with frequency 1 Hz. What is the friction?

**Exercise 2.4.101**: A mass of \(2\) kilograms is on a spring with spring constant \(k\) newtons per meter with no damping. Suppose the system is at rest and at time \( t = 0 \) the mass is kicked and starts traveling at 2 meters per second. How large does \(k\) have to be to so that the mass does not go further than 3 meters from the rest position?

**Exercise 2.4.102**: Suppose we have an RLC circuit with a resistor of 100 miliohms (0.1 ohms), inductor of inductance of 50 millihenries (0.05 henries), and a capacitor of 5 farads, with constant voltage. a) Set up the ODE equation for the current \(I\). b) Find the general solution. c) Solve for \(I (0) = 10\) and \(I' (0) = 0 \).

**Exercise 2.4.103**: A 5000 kg railcar hits a bumper (a spring) at \( 1 \dfrac {m}{s} \), and the spring compresses by 0.1 m. Assume no damping. a) Find \(k\). b) Find out how far does the spring compress when a 10000 kg railcar hits the spring at the same speed. c) If the spring would break if it compresses further than 0.3 m, what is the maximum mass of a railcar that can hit it at \(1 \dfrac {m}{s} \)? d) What is the maximum mass of a railcar that can hit the spring without breaking at \( 2 \dfrac {m}{s}\)?

## 2.5: Nonhomogeneous Equations

**Exercise 2.5.2:** Find a particular solution of \(y'' - y' - 6y = e^{2x} \).

**Exercise 2.5.3:** Find a particular solution of \(y'' - 4y' + 4y = e^{2x} \).

**Exercise 2.5.4:** Solve the initial value problem \( y'' + 9y = \cos (3x) + \sin (3x) \) for \( y(0) = 2, y'(0) = 1 \).

**Exercise 2.5.5:** Setup the form of the particular solution but do not solve for the coefficients for \( y^{(4)} - 2y''' + y'' = e^x\).

**Exercise 2.5.6:** Setup the form of the particular solution but do not solve for the coefficients for \( y^{(4)} - 2y''' + y'' = e^x + x + \sin x \).

**Exercise 2.5.7:** a) Using variation of parameters find a particular solution of \( y'' - 2y' + y = e^x \). b) Find a particular solution using undetermined coefficients. c) Are the two solutions you found the same? What is going on?

**Exercise 2.5.8:** Find a particular solution of \( y'' - 2y' + y = \sin (x^2) \). It is OK to leave the answer as a definite integral.

**Exercise 2.5.9:** For an arbitrary constant \(c\) find a particular solution to \(y'' - y = e^{cx} \). Hint: Make sure to handle every possible real \(c\).

**Exercise 2.5.101:** Find a particular solution to \( y'' - y' + y = 2 \sin (3x) \)

**Exercise 2.5.102:** a) Find a particular solution to \( y'' + 2y = e^x + x^3 \). b) Find the general solution.

**Exercise 2.5.103:** Solve \( y'' + 2y' + y = x^2, y(0) = 1, y'(0) = 2 \).

**Exercise 2.5.104:** Use variation of parameters to find a particular solution of \(y'' - y = \dfrac {1}{e^x + e^{-x}} \).

**Exercise 2.5.105:** For an arbitrary constant \(c\) find the general solution to \(y'' - 2y = \sin (x + c)\).

## 2.6: Forced Oscillations and Resonance

**Exercise 2.6.1:** Derive a formula for \(x_{sp}\) if the equation is \( mx'' + cx' + kx = F_0 \sin ( \omega t) \). Assume \( c > 0 \).

**Exercise 2.6.2:** Derive a formula for \( x_{sp}\) if the equation is \( mx'' + cx' + kx = F_0 \cos (\omega t) + F_1 \cos (3 \omega t)\). Assume \( c > 0 \).

**Exercise 2.6.3:** Take \( mx'' + cx' + kx = F_0 \cos ( \omega t) \). Fix \(m > 0\) and \( k > 0\). Now think of the function \( C(\omega ) \). For what values of \(c\) (solve in terms of \( m, k,\) and \(F_0\) ) will there be no practical resonance (that is, for what values of \(c\) is there no maximum of \(C(\omega )\) for \(\omega > 0 \) )?

**Exercise 2.6.4:** Take \(mx'' + cx' + kx = F_0 \cos (\omega t) \). Fix \( c > 0\) and \(k > 0 \). Now think of the function \( C(\omega ) \). For what values of \(m\) (solve in terms of \(c, k, \) and \(F_0 \)) will there be no practical resonance (that is, for what values of \(m\) is there no maximum of \(C(\omega )\) for \(\omega > 0\) )?

**Exercise 2.6.5:** Suppose a water tower in an earthquake acts as a mass-spring system. Assume that the container on top is full and the water does not move around. The container then acts as a mass and the support acts as the spring, where the induced vibrations are horizontal. Suppose that the container with water has a mass of \(m = 10, 000 kg \). It takes a force of 1000 newtons to displace the container 1 meter. For simplicity assume no friction. When the earthquake hits the water tower is at rest (it is not moving).

Suppose that an earthquake induces an external force \( F(t) = mA \omega^2 \cos (\omega t). \)

a) What is the natural frequency of the water tower?

b) If \(\omega \) is not the natural frequency, find a formula for the maximal amplitude of the resulting oscillations of the water container (the maximal deviation from the rest position). The motion will be a high frequency wave modulated by a low frequency wave, so simply find the constant in front of the sines.

c) Suppose \(A = 1\) and an earthquake with frequency 0.5 cycles per second comes. What is the amplitude of the oscillations? Suppose that if the water tower moves more than 1.5 meter, the tower collapses. Will the tower collapse?

**Exercise 2.6.101:** A mass of 4 kg on a spring with \(k = 4\) and a damping constant \(c = 1\). Suppose that \(F_0 = 2 \). Using forcing of \( F_0 \cos (\omega t) \). Find the \(\omega \) that causes practical resonance and find the amplitude.

**Exercise 2.6.102:** Derive a formula for \(x_{sp}\) for \( mx'' + cx' + kx = F_0 \cos (\omega t) + A \) where \(A\) is some constant. Assume \( c > 0 \).

**Exercise 2.6.103: **Suppose there is no damping in a mass and spring system with \( m = 5, k = 20, \) and \(F_0 = 5 \). Suppose that \( \omega \) is chosen to be precisely the resonance frequency. a) Find \(\omega \). b) Find the amplitude of the oscillations at time \(t = 100\).