Variation of Parameters — Question 1

PDF ↗

Question 1

Consider the fourth-order IVP (1+x2)y(4)=1,y(0)=y′(0)=y″(0)=y‴(0)=0.(1+x^2)y^{(4)}=1,\qquad y(0)=y'(0)=y''(0)=y'''(0)=0. Use the normalized homogeneous basis 1,x,x2/2,x3/61,x,x^2/2,x^3/6.

Tasks

  1. Normalize the equation and verify that the stated basis is fundamental. Write the four equations for the parameter derivatives in variation of parameters.

  2. Solve those equations. Integrate each parameter from zero and combine the result into one definite integral for y(x)y(x).

  3. Verify the equation and all initial data by differentiating the combined integral. Explain exactly what goes wrong if the leading coefficient is not divided out.

  4. Prove upper and lower bounds for y(x)y(x) when x≥0x\ge 0 using the largest and smallest values of 1/(1+t2)1/(1+t^2) on [0,x][0,x]. State when equality holds.

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

Question 1 – Solution

Strategy. Normalize first, then combine the four parameter integrals into one positive kernel.

Step 1: Form the normalized parameter system. Put g(x)=1/(1+x2)g(x)=1/(1+x^2). The derivative matrix of the given basis is triangular with unit diagonal, so its Wronskian is 11. The parameter derivatives satisfy (1xx2/2x3/601xx2/2001x0001)(u0′u1′u2′u3′)=(000g).\begin{pmatrix} 1&x&x^2/2&x^3/6\\0&1&x&x^2/2\\0&0&1&x\\0&0&0&1 \end{pmatrix} \begin{pmatrix}u_0'\\u_1'\\u_2'\\u_3'\end{pmatrix} =\begin{pmatrix}0\\0\\0\\g\end{pmatrix}. The first three zero rows remove derivative terms introduced by varying the constants; the last row supplies the normalized forcing.

Step 2: Solve and combine. Back substitution gives (u0′,u1′,u2′,u3′)=g(−x3/6,x2/2,−x,1)(u_0',u_1',u_2',u_3')=g(-x^3/6,x^2/2,-x,1). Choose each uj(0)=0u_j(0)=0. In y=∑ujfjy=\sum u_jf_j, the combined numerator is (x−t)3(x-t)^3, yielding y(x)=16∫0x(x−t)31+t2dt.\boxed{y(x)=\frac 16\int_0^x\frac{(x-t)^3}{1+t^2}\,dt.} Oriented integration also defines the solution for negative xx.

Step 3: Verify the initial state and residual. Successive derivatives have kernels (x−t)2/2(x-t)^2/2, x−tx-t, and 11, all divided by 1+t21+t^2. The fourth derivative is 1/(1+x2)1/(1+x^2). The first four state components vanish at zero, and multiplication by 1+x21+x^2 gives the original forcing. If one incorrectly uses g=1g=1, the same construction gives x4/24x^4/24, whose original residual is 1+x21+x^2, not 11.

Step 4: Bound the positive integral. For 0≤t≤x0\le t\le x, 1/(1+x2)≤1/(1+t2)≤11/(1+x^2)\le 1/(1+t^2)\le 1. Hence x424(1+x2)≤y(x)≤x424(x≥0).\boxed{\frac{x^4}{24(1+x^2)}\le y(x)\le\frac{x^4}{24}\quad(x\ge 0).} Both inequalities are strict for x>0x>0, since their integrands differ on an interval of positive length where (x−t)3>0(x-t)^3>0. Both become equalities at x=0x=0. The figure labels the lower and upper bounds LL and UU, respectively.

See the diagram in the original worksheet below.

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

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