Bernoulli Differential Equations — Question 10

PDF ↗

Question 10

Consider positive solutions, near x=0x=0, of the exponent-dependent IVP y′+y=yn,y(0)=2,n∈ℝ.y'+y=y^n,\qquad y(0)=2,\qquad n\in\mathbb R. For y>0y>0, interpret yn=exp⁡(nln⁡y)y^n=\exp(n\ln y).

Tasks

  1. Solve the cases n=0n=0 and n=1n=1 directly. Explain why they do not require the usual nonlinear Bernoulli reduction.

  2. For n≠1n\ne 1, derive the solution formula and state the positivity condition needed to interpret its real power.

  3. Check that the general formula agrees with the direct answer at n=0n=0 and verifies the initial condition for every admissible nn.

  4. For each fixed real xx, prove that the formula tends to the n=1n=1 solution as n→1n\to 1. Do not substitute n=1n=1 into an expression with exponent 1/(1−n)1/(1-n).

Original worksheet page 1: question and worked solution for 2-4-010
Show solutionHide solution

Question 10 – Solution

Strategy. Handle the exceptional equations directly, then analyze the apparent singularity of the general expression through its logarithm.

Step 1: Solve the elementary exceptional cases. If n=0n=0, the equation is y′+y=1y'+y=1, giving y=1+e−x\boxed{y=1+e^{-x}}. If n=1n=1, the two yy terms cancel, giving y=2\boxed{y=2}. Both are positive on ℝ\mathbb R; the first is linear and the second has zero derivative.

Step 2: Derive the parameterized formula. For n≠1n\ne 1, set h=1−n≠0h=1-n\ne 0 and v=yhv=y^h. Then v′+hv=h,v(0)=2h,v'+hv=h,\qquad v(0)=2^h, so v=1+(2h−1)e−hxv=1+(2^h-1)e^{-hx}. Therefore yn(x)=[1+(21−n−1)e−(1−n)x]1/(1−n).\boxed{y_n(x)=\bigl[1+(2^{1-n}-1)e^{-(1-n)x}\bigr]^{1/(1-n)}}. Use the connected interval through 00 on which the bracket is strictly positive. This is necessary because the original positive solution has v=yh>0v=y^h>0; a negative base cannot be accepted by choosing a special rational exponent.

Step 3: Check the consistency conditions. At n=0n=0, the formula becomes 1+e−x1+e^{-x}. At x=0x=0, its bracket is 21−n2^{1-n}, and the positive-real power gives 22. Conversely, differentiating v=yhv=y^h on the stated interval recovers the original equation, so the inversion is valid there.

Step 4: Take the limit without dividing by zero. For fixed xx, let B(h)=1+(ehln⁡2−1)e−hxB(h)=1+(e^{h\ln 2}-1)e^{-hx}. Since B(0)=1B(0)=1, it is positive for all sufficiently small |h||h|. Now ln⁡yn=ln⁡B(h)h,B′(0)=ln⁡2.\ln y_n=\frac{\ln B(h)}h,\qquad B'(0)=\ln 2. The derivative definition gives lim⁡h→0ln⁡B(h)/h=ln⁡2\lim_{h\to 0}\ln B(h)/h=\ln 2. Exponentiating yields lim⁡n→1yn(x)=2\boxed{\lim_{n\to 1}y_n(x)=2}, exactly the direct n=1n=1 solution. Joint continuity on the compact segment between 00 and the fixed xx keeps the bracket positive there for small |h||h|, so these values belong to the IVP branch.

Original worksheet page 2: question and worked solution for 2-4-010

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