Undetermined Coefficients — Question 4

PDF ↗

Question 4

On ℝ\mathbb R, consider y(4)−y=6x2+8e−x+10cos⁡x.y^{(4)}-y=6x^2+8e^{-x}+10\cos x. Its three forcing terms require different coefficient calculations.

Tasks

  1. Factor the characteristic polynomial, write the full real homogeneous solution, and select a separate particular-solution trial for each forcing term.

  2. Determine all coefficients and assemble one particular solution.

  3. Verify each forcing contribution directly using L=D4−1L=D^4-1, rather than relying only on the trial-selection rule. Write the complete general solution.

  4. Decide whether any solution is bounded on [0,∞)[0,\infty). Prove your answer even when homogeneous constants are chosen to cancel all removable growing terms.

Original worksheet page 1: question and worked solution for 7-3-004
Show solutionHide solution

Question 4 – Solution

Strategy. Use linearity to separate polynomial forcing from two distinct resonant frequencies.

Step 1: Identify roots and trials. The characteristic roots are 1,−1,i,−i1,-1,i,-i, all simple, so yh=Aex+Be−x+Ccos⁡x+Dsin⁡x.y_h=Ae^x+Be^{-x}+C\cos x+D\sin x. Use ax2+bx+cax^2+bx+c for 6x26x^2, kxe−xkxe^{-x} for 8e−x8e^{-x}, and x(mcos⁡x+nsin⁡x)x(m\cos x+n\sin x) for 10cos⁡x10\cos x. Only the last two blocks are resonant.

Step 2: Compute the coefficients. The fourth derivative of a quadratic is zero, giving a=−6,b=c=0a=-6,b=c=0. For a simple root λ\lambda of P(r)=r4−1P(r)=r^4-1, the shift identity gives P(D)(xeλx)=P′(λ)eλxP(D)(xe^{\lambda x})=P'(\lambda)e^{\lambda x}. Since P′(−1)=−4P'(-1)=-4, we obtain k=−2k=-2.

For the trigonometric terms, L(xcos⁡x)=4sin⁡xL(x\cos x)=4\sin x and L(xsin⁡x)=−4cos⁡xL(x\sin x)=-4\cos x. Therefore m=0,n=−5/2m=0,n=-5/2, and yp=−6x2−2xe−x−52xsin⁡x.\boxed{y_p=-6x^2-2xe^{-x}-\frac 52x\sin x.}

Step 3: Verify all three residuals. Direct differentiation gives L(−6x2)=6x2,D4(xe−x)=e−x(x−4),L(-6x^2)=6x^2,\qquad D^4(xe^{-x})=e^{-x}(x-4), D4(xsin⁡x)=xsin⁡x−4cos⁡x.D^4(x\sin x)=x\sin x-4\cos x. Thus the second and third terms produce 8e−x8e^{-x} and 10cos⁡x10\cos x exactly. The full solution is y=yh+yp\boxed{y=y_h+y_p}; no homogeneous constants are determined without additional data.

Step 4: Rule out bounded solutions. If A≠0A\ne 0, the AexAe^x term dominates every other term, so boundedness is impossible. If A=0A=0, division by x2x^2 gives limx→∞y(x)x2=−6,\lim_{x\to\infty}\frac{y(x)}{x^2}=-6, because the remaining oscillatory growth is only O(x)O(x) and the other exponential terms decay. Hence even after removing the positive exponential, quadratic forcing remains. No solution is bounded on the positive half-line.

Original worksheet page 2: question and worked solution for 7-3-004

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