Equilibrium Solutions — Question 3

PDF ↗

Question 3

Two proposed relaxation laws are A: y′=−y3,B: y′=y3.\text{A: }y'=-y^3,\qquad \text{B: }y'=y^3. Both have an equilibrium at 00 and derivative of the right-hand side equal to zero there.

Tasks

  1. Classify 00 for each law using the sign of the right-hand side. Explain why their identical derivative-test data do not imply identical stability.

  2. Solve both IVPs with y(0)=ay(0)=a, including a=0a=0, and specify their maximal forward time ranges.

  3. For law A, prove stability directly from the solution and determine its large-time decay for a≠0a\ne 0.

  4. For law B, use the solution to prove instability and find the blow-up time for a≠0a\ne 0. Distinguish leaving a small neighborhood from blowing up.

Original worksheet page 1: question and worked solution for 2-8-003
Show solutionHide solution

Question 3 – Solution

Strategy. When the linear term vanishes, retain the nonlinear term. Separation gives both a stability proof and the actual time scales.

Step 1: Classify from signs. For A, −y3-y^3 is positive below 00 and negative above, so trajectories move toward 00. For B the signs reverse, so trajectories move away. Thus 00 is asymptotically stable for A and unstable for B. The shared value f′(0)=0f'(0)=0 is an inconclusive test, not evidence of neutral stability.

Step 2: Solve with the initial sign preserved. For a≠0a\ne 0, differentiating y−2y^{-2} gives 22 for A and −2-2 for B. Applying the initial value yields yA(t)=a1+2a2t,t≥0,\boxed{y_A(t)=\frac{a}{\sqrt{1+2a^2t}},\qquad t\ge 0}, yB(t)=a1−2a2t,0≤t<12a2.\boxed{y_B(t)=\frac{a}{\sqrt{1-2a^2t}},\qquad 0\le t<\frac 1{2a^2}}. The full maximal intervals containing 00 are (−1/(2a2),∞)(-1/(2a^2),\infty) and (−∞,1/(2a2))(-\infty,1/(2a^2)), respectively. Direct differentiation verifies the equations. If a=0a=0, the unique solution of either smooth equation is zero for all real time.

Step 3: Prove stability and attraction for A. For t≥0t\ge 0, |yA(t)|≤|a||y_A(t)|\le|a|. Given ε>0\varepsilon>0, choosing δ=ε\delta=\varepsilon ensures |a|<δ|a|<\delta implies |yA(t)|<ε|y_A(t)|<\varepsilon for all future time. Also yA(t)→0y_A(t)\to 0 for every aa, so the attraction is global. More precisely, 2t|yA(t)|→1(a≠0).\sqrt{2t}\,|y_A(t)|\longrightarrow 1\quad(a\ne 0). The decay is algebraic rather than exponential.

Step 4: Prove instability for B. Choose any fixed ε>0\varepsilon>0 and any nonzero |a|<ε|a|<\varepsilon. The magnitude first reaches ε\varepsilon at tε=12a2−12ε2>0,t_\varepsilon=\frac 1{2a^2}-\frac 1{2\varepsilon^2}>0, then exceeds it. Arbitrarily small nonzero perturbations therefore leave that neighborhood. Blow-up occurs later, at 1/(2a2)1/(2a^2), with the sign of aa. Instability only requires departure from a prescribed neighborhood; it does not require an infinite solution value.

Original worksheet page 2: question and worked solution for 2-8-003

Original worksheet layout. Use Enlarge or open the PDF for a closer view.