Euler Equations — Question 6

PDF ↗

Question 6

For x>0x>0, consider x2y″+xy′=x21+x2.x^2y''+xy'=\frac{x^2}{1+x^2}. Instead of prescribing data at a positive point, seek the solution that is bounded near 00 and extends continuously with y(0)=0y(0)=0.

Tasks

  1. Integrate the equation to obtain its general solution in definite-integral form on (0,∞)(0,\infty), displaying both arbitrary constants.

  2. Impose boundedness and the specified value at zero. Prove that they select exactly one solution on the positive half-line.

  3. Derive its power series about zero, give its first three nonzero terms, and determine the radius and both endpoint outcomes. Justify termwise integration.

  4. Prove that the selected solution has a C2C^2 extension to zero, compute y′(0)y'(0) and y″(0)y''(0), and certify a two-term approximation on [0,1][0,1] using a rigorous remainder bound.

Original worksheet page 1: question and worked solution for 6-4-006
Show solutionHide solution

Question 6 – Solution

Strategy. Integrate the weighted derivative first, then use regularity at zero to remove the logarithmic homogeneous mode.

Step 1: Integrate the Euler operator. Dividing by xx gives (xy′)′=x/(1+x2)(xy')'=x/(1+x^2). Hence xy′=12ln⁡(1+x2)+C,y=12∫1xln⁡(1+s2)sds+Cln⁡x+D.xy'=\tfrac 12\ln(1+x^2)+C,\qquad y=\frac 12\int_1^x\frac{\ln(1+s^2)}s\,ds+C\ln x+D. These formulas display both constants and solve the equation for every x>0x>0.

Step 2: Apply the singular-endpoint conditions. Since ln⁡(1+s2)/s∼s\ln(1+s^2)/s\sim s at zero, its integral there converges. Boundedness therefore forces C=0C=0; the value y(0)=0y(0)=0 fixes the remaining constant. The unique selected solution is y(x)=12∫0xln⁡(1+s2)sds.\boxed{y(x)=\frac 12\int_0^x\frac{\ln(1+s^2)}s\,ds.} Any two bounded solutions differ by a constant, and their common value at zero makes that constant zero.

Step 3: Expand the integral with justification. The divided logarithmic series ∑n≥1(−1)n−1s2n−1/n\sum_{n\ge 1}(-1)^{n-1}s^{2n-1}/n is uniformly convergent on |s|≤ρ<1|s|\le\rho<1, by comparison with ∑n≥1ρ2n−1/n\sum_{n\ge 1}\rho^{2n-1}/n. Termwise integration therefore yields, for |x|<1|x|<1, y(x)=14∑n=1∞(−1)n−1x2nn2=x24−x416+x636−⋯.\boxed{y(x)=\frac 14\sum_{n=1}^{\infty}\frac{(-1)^{n-1}x^{2n}}{n^2} =\frac{x^2}{4}-\frac{x^4}{16}+\frac{x^6}{36}-\cdots.} The ratio of successive nonzero terms tends to x2x^2, so the radius is 11. At x=±1x=\pm 1 the series converges absolutely by comparison with ∑1/n2\sum 1/n^2. The series and the even integral extension agree there by continuity: the series is uniformly convergent on [−1,1][-1,1] by the same comparison.

Step 4: Recover derivatives and a certificate. The local power series supplies a smooth even extension with y′(0)=0y'(0)=0 and y″(0)=1/2y''(0)=1/2. With the forcing defined continuously as zero at 00, the original equation holds there as well. On [0,1][0,1], the alternating term magnitudes decrease, so for P4=x2/4−x4/16P_4=x^2/4-x^4/16, 0≤y(x)−P4(x)≤x636.\boxed{0\le y(x)-P_4(x)\le\frac{x^6}{36}.} The bound includes the endpoint x=1x=1; it is a certificate, not a claim that the two-term approximation is highly accurate there.

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

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