Skip to content

Common Pitfalls

:::caution Common Pitfall When using undetermined coefficients, always check whether your guess Overlaps with the homogeneous solution. For y"4y=e2xy"' - 4y = e^{2x}Guessing yp=Ae2xy_p = Ae^{2x} fails Because e2xe^{2x} satisfies the homogeneous equation. You must use yp=Axe2xy_p = Axe^{2x} instead. :::

:::caution Common Pitfall The Laplace transform of y(t)y'(t) is sY(s)y(0)sY(s) - y(0)Not sY(s)sY(s). The Initial conditions are built into the transform. Forgetting them leads to incorrect solutions. :::

:::caution Common Pitfall Separation of variables can miss solutions. When dividing by h(y)h(y) to Separate, check whether h(y)=0h(y) = 0 gives any valid solutions. For example, solving y=y2y' = y^2 by Separating gives y=1/(x+C)y = -1/(x + C)But misses the solution y=0y = 0. :::

:::caution Common Pitfall The Fourier series of a function converges to the function only at points Of continuity. At jump discontinuities, it converges to the midpoint of the jump. The Gibbs Phenomenon causes overshoots near jumps that do not vanish as more terms are added. :::

:::caution Common Pitfall When solving PDEs by separation of variables, the boundary conditions Determine the eigenvalues. Using the wrong boundary conditions (e.g., Neumann instead of Dirichlet) Leads to a completely different set of eigenfunctions and eigenvalues. :::

:::caution Common Pitfall Not every first-order ODE can be solved analytically. Equations like y=x2+y2y' = x^2 + y^2 have no closed-form solution in terms of elementary functions. Numerical methods (Euler, Runge-Kutta) may be necessary. :::

:::caution Common Pitfall The linearization of a nonlinear system near a critical point is only Valid for hyperbolic critical points (no eigenvalues on the imaginary axis). If eigenvalues lie on The imaginary axis, the nonlinear system can behave very differently from its linearization. :::

:::caution Common Pitfall When computing the inverse Laplace transform, always check that the Partial fraction decomposition is correct before inverting term-by-term. A common error is Forgetting to include all powers of irreducible quadratic factors. :::

:::caution Common Pitfall For the Euler-Cauchy equation x2y+axy+by=0x^2 y'' + axy' + by = 0Remember that xrx^r with complex r=α±iβr = \alpha \pm i\beta gives solutions involving cos(βlnx)\cos(\beta \ln x) and sin(βlnx)\sin(\beta \ln x)Not cos(βx)\cos(\beta x) and sin(βx)\sin(\beta x). :::

:::caution Common Pitfall The Wronskian W(y1,y2)W(y_1, y_2) being zero at a single point does not Necessarily mean the solutions are linearly dependent. For linear ODEs with continuous coefficients, W0W \equiv 0 everywhere or W0W \neq 0 everywhere. Check Abel’s identity. :::

9.9 Worked Example: Lyapunov Function

Problem. Show that (0,0)(0, 0) is asymptotically stable for the system x=xy3x' = -x - y^3, y=x3yy' = x^3 - y.

Solution

Solution. Critical point: xy3=0-x - y^3 = 0 and x3y=0x^3 - y = 0 gives x=0,y=0x = 0, y = 0.

Choose the Lyapunov function V(x,y)=x22+y22V(x, y) = \frac{x^2}{2} + \frac{y^2}{2}. V(0,0)=0V(0,0) = 0 and V>0V > 0 for (x,y)(0,0)(x, y) \neq (0, 0).

V˙=Vxx+Vyy=x(xy3)+y(x3y)\dot{V} = \frac{\partial V}{\partial x}x' + \frac{\partial V}{\partial y}y' = x(-x - y^3) + y(x^3 - y)

=x2xy3+x3yy2=x2y2+xy(x2y2)= -x^2 - xy^3 + x^3 y - y^2 = -x^2 - y^2 + xy(x^2 - y^2)

This is not negative definite. Let us try V(x,y)=x44+y44V(x, y) = \frac{x^4}{4} + \frac{y^4}{4}:

V˙=x3(xy3)+y3(x3y)=x4x3y3+x3y3y4=(x4+y4)\dot{V} = x^3(-x - y^3) + y^3(x^3 - y) = -x^4 - x^3 y^3 + x^3 y^3 - y^4 = -(x^4 + y^4).

Since V˙=(x4+y4)<0\dot{V} = -(x^4 + y^4) \lt 0 for (x,y)(0,0)(x, y) \neq (0, 0)The origin is asymptotically Stable by Lyapunov’s theorem. \blacksquare

9.10 Worked Example: Van der Pol Oscillator

Problem. Show that the Van der Pol equation x+μ(x21)x+x=0x'' + \mu(x^2 - 1)x' + x = 0 with μ>0\mu > 0 Has a unique limit cycle.

Solution

_Solution.* Write as a system: x=yx' = y, y=xμ(x21)yy' = -x - \mu(x^2 - 1)y.

The only critical point is (0,0)(0, 0). The Jacobian at the origin is J=(011μ)J = \begin{pmatrix} 0 & 1 \\ -1 & \mu \end{pmatrix}.

tr(J)=μ>0\mathrm{tr}(J) = \mu > 0, det(J)=1>0\det(J) = 1 > 0, τ24Δ=μ24\tau^2 - 4\Delta = \mu^2 - 4.

For 0<μ<20 \lt \mu \lt 2: complex eigenvalues with positive real part (unstable spiral). For μ2\mu \geq 2: two positive real eigenvalues (unstable node).

The origin is always unstable. To show existence of a limit cycle, we use a Liénard’s theorem Or construct a trapping region. Define

L(x)=x3/3xL(x) = x^3/3 - x and write the system as x=yx' = y, y=xμL(x)yy' = -x - \mu L'(x) y.

The function L(x)L(x) has zeros at x=±1x = \pm 1. For μ>0\mu > 0The damping is negative for x<1\lvert x \rvert \lt 1 (energy input) and positive for x>1\lvert x \rvert > 1 (energy Dissipation). This creates a unique stable limit cycle that passes through x=±1x = \pm 1.

A detailed …/1-number-and-algebra/3_proof-and-logic requires constructing an annular region and applying the Poincaré-Bendixson Theorem (the inner boundary encloses the unstable origin; the outer boundary is chosen so that Trajectories point inward). \blacksquare