Undetermined Coefficients — Question 4

Section 26-3 · Calculus pages 4627–4628

Show solutionHide solution

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.

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.