Variation of Parameters — Question 10

PDF ↗

Question 10

For t≥0t\ge 0, solve y″+y=e−t2,y(0)=y′(0)=0.y''+y=e^{-t^2},\qquad y(0)=y'(0)=0. Define the convergent integrals C=∫0∞e−s2cos⁡sdsC=\int_0^\infty e^{-s^2}\cos s\,ds and S=∫0∞e−s2sin⁡sdsS=\int_0^\infty e^{-s^2}\sin s\,ds. No special-function evaluations of CC or SS are required.

Tasks

  1. Use variation of parameters to express the exact solution through two finite definite integrals, and verify its initial data and residual.

  2. Show that the solution approaches the periodic function P(t)=Csin⁡t−Scos⁡tP(t)=C\sin t-S\cos t in the sense that y(t)−P(t)→0y(t)-P(t)\to 0.

  3. Prove the quantitative estimate |y(t)−P(t)|≤e−t2/(2t)|y(t)-P(t)|\le e^{-t^2}/(2t) for every t>0t>0.

  4. Prove S>0S>0 by pairing successive positive and negative half-waves of sine. Use this to show that yy does not tend to zero, even though the forcing does.

Original worksheet page 1: question and worked solution for 3-10-010
Show solutionHide solution

Question 10 – Solution

Strategy. Keep the variation-of-parameters integrals exact, then bound their tails to identify the surviving oscillation.

Step 1: Write and check the exact response. With W=1W=1, the parameter derivatives are u1′=−e−t2sin⁡tu_1'=-e^{-t^2}\sin t and u2′=e−t2cos⁡tu_2'=e^{-t^2}\cos t. Hence y(t)=sin⁡t∫0te−s2cos⁡sds−cos⁡t∫0te−s2sin⁡sds.\boxed{y(t)=\sin t\int_0^t e^{-s^2}\cos s\,ds -\cos t\int_0^t e^{-s^2}\sin s\,ds.} Equivalently, y=∫0tsin⁡(t−s)e−s2dsy=\int_0^t\sin(t-s)e^{-s^2}\,ds. Differentiating gives y′=∫0tcos⁡(t−s)e−s2dsy'=\int_0^t\cos(t-s)e^{-s^2}\,ds and y″=e−t2−yy''=e^{-t^2}-y; both initial data vanish.

Step 2: Isolate the tails. Absolute convergence follows from the integrability of e−s2e^{-s^2} (for s≥1s\ge 1, it is at most e−se^{-s}). Subtracting the full integrals gives y(t)−P(t)=−∫t∞sin⁡(t−s)e−s2ds.y(t)-P(t)=-\int_t^\infty\sin(t-s)e^{-s^2}\,ds. Thus |y−P|≤∫t∞e−s2ds→0|y-P|\le\int_t^\infty e^{-s^2}\,ds\to 0.

Step 3: Obtain an explicit error bound. For s≥t>0s\ge t>0, 1≤s/t1\le s/t. Consequently |y−P|≤1t∫t∞se−s2ds=e−t22t.|y-P|\le\frac 1t\int_t^\infty s e^{-s^2}\,ds =\boxed{\frac{e^{-t^2}}{2t}}. Combining the tails before estimating avoids an unnecessary sum of two bounds.

Step 4: Prove a nonzero surviving oscillation. Absolute convergence permits pairing intervals. Their kkth pair contributes ∫0πsin⁡u[e−(2kπ+u)2−e−((2k+1)π+u)2],du>0,k≥0.\int_0^\pi\sin u[e^{-(2k\pi+u)^2} -e^{-((2k+1)\pi+u)^2}],du>0,\qquad k\ge 0. Therefore S>0S>0. Along t=2nπt=2n\pi, y(t)→−Sy(t)\to-S; along t=(2n+1)πt=(2n+1)\pi, y(t)→Sy(t)\to S. The solution has no limit and does not decay to zero. A decaying forcing can leave a persistent homogeneous oscillation.

See the diagram in the original worksheet below.

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

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