Series Solutions — Question 7

PDF ↗

Question 7

Consider the nonlinear initial-value problem y′=1+y2,y(0)=0.y'=1+y^2,\qquad y(0)=0. Seek a Maclaurin solution y=∑n≥0anxny=\sum_{n\ge 0}a_nx^n. For the radius question, you may use these facts: sine and cosine are entire; the only complex zeros of cosine are z=π/2+kπz=\pi/2+k\pi, k∈ℤk\in\mathbb Z; and a Taylor series reaches its nearest nonremovable complex singularity.

Tasks

  1. Derive the nonlinear coefficient recurrence using a Cauchy product. Compute all coefficients through degree seven.

  2. Prove that the resulting solution is odd, and explain why the square y2y^2 cannot be handled by squaring individual coefficients separately.

  3. Identify the exact solution by a method independent of the coefficient calculation. Determine its maximal real interval containing 00 and the exact Taylor radius.

  4. For the degree-seven polynomial P7P_7, prove 0<P7(x)<y(x)0<P_7(x)<y(x) when 0<x<π/20<x<\pi/2. Explain why a fixed polynomial cannot approximate this solution uniformly on that whole open interval.

Original worksheet page 1: question and worked solution for 6-3-007
Show solutionHide solution

Question 7 – Solution

Strategy. Nonlinearity couples coefficients by convolution; an exact solution then validates convergence and reveals blow-up.

Step 1: Convolve rather than square termwise. The constant forcing contributes only at n=0n=0: (n+1)an+1=δn0+∑k=0nakan−k,a0=0,\boxed{(n+1)a_{n+1}=\delta_{n0}+\sum_{k=0}^n a_ka_{n-k},\qquad a_0=0,} where δn0=1\delta_{n0}=1 for n=0n=0 and 00 otherwise. This gives P7(x)=x+x33+2x515+17x7315,\boxed{P_7(x)=x+\frac{x^3}{3}+\frac{2x^5}{15}+\frac{17x^7}{315},} with a2=a4=a6=0a_2=a_4=a_6=0.

Step 2: Establish parity and coefficient signs. If yy solves the IVP, −y(−x)-y(-x) does too; local uniqueness gives oddness. The recurrence also proves by induction that all even coefficients vanish and all odd coefficients are positive: each required convolution contains positive odd-index products. For example, the coefficient of x4x^4 in y2y^2 is 2a1a32a_1a_3, not a22a_2^2. Cross terms are essential.

Step 3: Identify the function and its radius. Separation yields arctan⁡y=x\arctan y=x near 00, so y=tan⁡x\boxed{y=\tan x}. It satisfies the equation and initial condition on (−π/2,π/2)(-\pi/2,\pi/2) and blows up at both endpoints, making this the maximal real interval through 00. Since sin⁡(±π/2)≠0\sin(\pm\pi/2)\ne 0, the nearest zeros of cosine are nonremovable poles of the quotient. The supplied facts give R=π/2\boxed{R=\pi/2}. Within this radius, the analytic solution’s coefficients obey the recurrence, which uniquely determines them; hence it is the constructed series. At x=π/2x=\pi/2 the positive-term series diverges: a finite endpoint sum would bound tan⁡x\tan x below that endpoint, contradicting blow-up. Oddness gives divergence at −π/2-\pi/2 as well.

Step 4: Interpret truncation and blow-up. For 0<x<π/20<x<\pi/2, every omitted odd term is positive, so 0<P7(x)<tan⁡x0<P_7(x)<\tan x. But P7P_7 stays bounded as x→(π/2)−x\to(\pi/2)^- while tan⁡x→+∞\tan x\to+\infty. Their difference is unbounded on that interval. Every fixed Taylor truncation has this defect; convergence on smaller compact intervals does not imply uniform approximation up to a blow-up boundary.

See the diagram in the original worksheet below.

Original worksheet page 2: question and worked solution for 6-3-007

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