Variation of Parameters — Question 7

PDF ↗

Question 7

For the variable-coefficient IVP y‴−xy″=1,y(0)=y′(0)=y″(0)=0,y'''-xy''=1,\qquad y(0)=y'(0)=y''(0)=0, define F(x)=∫0x(x−t)et2/2dtF(x)=\int_0^x(x-t)e^{t^2/2}\,dt and use the homogeneous basis 1,x,F(x)1,x,F(x).

Tasks

  1. Verify the basis and its Wronskian, then derive the three variation-parameter derivatives.

  2. Combine the lower-limit-zero parameter integrals into y(x)=∫0xK(x,t)dty(x)=\int_0^xK(x,t)\,dt. Write KK both using F,F′F,F' and as a definite integral for x≥tx\ge t.

  3. Differentiate to verify the IVP. Prove that yy is strictly convex and y(x)>x3/6y(x)>x^3/6 for every x>0x>0.

  4. Determine whether K(x,t)K(x,t) can depend only on x−tx-t. Use a derivative of the kernel to give a precise obstruction to treating this variable-coefficient response as a translation-invariant convolution.

Original worksheet page 1: question and worked solution for 7-4-007
Show solutionHide solution

Question 7 – Solution

Strategy. Exploit the triangular derivative matrix while retaining both kernel variables.

Step 1: Solve the parameter equations. We have F(0)=F′(0)=0F(0)=F'(0)=0, F″=ex2/2F''=e^{x^2/2} and F‴=xF″F'''=xF''. Thus 1,x,F1,x,F solve the homogeneous equation and have Wronskian ex2/2>0e^{x^2/2}>0. The triangular system gives u2′=e−x2/2,u1′=−F′e−x2/2,u0′=(xF′−F)e−x2/2.u_2'=e^{-x^2/2},\quad u_1'=-F'e^{-x^2/2},\quad u_0'=(xF'-F)e^{-x^2/2}.

Step 2: Combine the parameter integrals. Multiplying by 1,x,F(x)1,x,F(x) and integrating from zero yields K(x,t)=e−t2/2[F(x)−F(t)−(x−t)F′(t)].\boxed{K(x,t)=e^{-t^2/2}\bigl[F(x)-F(t)-(x-t)F'(t)\bigr].} Taylor’s integral identity, obtained by integrating F″F'' twice, gives for x≥tx\ge t K(x,t)=e−t2/2∫tx(x−s)es2/2ds,y(x)=∫0xK(x,t)dt.K(x,t)=e^{-t^2/2}\int_t^x(x-s)e^{s^2/2}\,ds,\qquad y(x)=\int_0^xK(x,t)\,dt. The factor involving tt is essential; it comes from the inverse derivative matrix evaluated at the integration variable.

Step 3: Verify the equation and inequalities. At the diagonal, K(t,t)=Kx(t,t)=0K(t,t)=K_x(t,t)=0, while Kxx(x,t)=e(x2−t2)/2,y″(x)=ex2/2∫0xe−t2/2dt.K_{xx}(x,t)=e^{(x^2-t^2)/2},\qquad y''(x)=e^{x^2/2}\int_0^xe^{-t^2/2}\,dt. Hence y‴=xy″+1y'''=xy''+1 and all three initial data vanish. For x>0x>0, y″>0y''>0, proving strict convexity. Also y‴=1+xy″>1y'''=1+xy''>1 there. Integrating this strict inequality three times from the zero initial vector gives y(x)>x3/6\boxed{y(x)>x^3/6} for every x>0x>0.

Step 4: Test translation invariance. If K(x,t)=H(x−t)K(x,t)=H(x-t), then Kxx=H″(x−t)K_{xx}=H''(x-t) would have the same value at all pairs with the same difference. But for any fixed h>0h>0, Kxx(t+h,t)=eth+h2/2,K_{xx}(t+h,t)=e^{th+h^2/2}, which varies with tt. Therefore this kernel is not a function of x−tx-t alone. The variable coefficient distinguishes the absolute source time from the elapsed time.

Original worksheet page 2: question and worked solution for 7-4-007

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