Variation of Parameters — Question 4

PDF ↗

Question 4

Let gg be continuous on ℝ\mathbb R and consider y‴=g(x)y'''=g(x). Two homogeneous bases are f=(1,x,x2/2),f̃=(1+x,x+x2/2,1+x2/2).f=(1,x,x^2/2),\qquad \widetilde f=(1+x,\ x+x^2/2,\ 1+x^2/2). Regard these as row vectors of functions and write yp=fu=f̃vy_p=fu=\widetilde f\,v, with column vectors of parameters.

Tasks

  1. Find the constant matrix MM such that f̃=fM\widetilde f=fM. Verify that both bases are fundamental.

  2. Derive u′u' by variation of parameters, then obtain v′v' from the change of basis. Explain why the lower-limit-zero constructions produce exactly the same particular solution.

  3. Replace the lower integration limit 00 by an arbitrary real bb. Find the difference between the two resulting particular solutions and prove that it is homogeneous.

  4. Explain why changing basis or choosing separate integration constants does not create more than three free constants in the general solution. Distinguish unique parameter coordinates in a chosen basis from the nonuniqueness of a particular solution.

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

Question 4 – Solution

Strategy. Track the constant coordinate transformation and isolate the effect of the integration base point.

Step 1: Identify the basis transformation. The coefficient columns give M=(101110011),det⁡M=2.M=\begin{pmatrix}1&0&1\\1&1&0\\0&1&1\end{pmatrix},\qquad \det M=2. All six displayed functions are quadratic or lower and solve y‴=0y'''=0. The Wronskian of ff is 11, and that of f̃\widetilde f is det⁡M=2\det M=2, so both triples are fundamental.

Step 2: Transform the coordinates. The standard derivative system gives u′=g(x)(x2/2,−x,1)Tu'=g(x)(x^2/2,-x,1)^T. Since the derivative matrices satisfy F̃=FM\widetilde F=FM, their parameter equations imply Mv′=u′Mv'=u', hence v′=g(x)(x2/4−x/2−1/2−x2/4−x/2+1/2x2/4+x/2+1/2).\boxed{v'=g(x)\begin{pmatrix}x^2/4-x/2-1/2\\-x^2/4-x/2+1/2\\x^2/4+x/2+1/2\end{pmatrix}.} With u(0)=v(0)=0u(0)=v(0)=0, integration gives u=Mvu=Mv at every xx. Thus fu=fMv=f̃vfu=fMv=\widetilde f\,v exactly, not merely up to a homogeneous term. Both equal Y0(x)=12∫0x(x−t)2g(t)dt.Y_0(x)=\frac 12\int_0^x(x-t)^2g(t)\,dt.

Step 3: Change the integration base point. The construction based at bb gives Yb(x)=12∫bx(x−t)2g(t)dtY_b(x)=\frac 12\int_b^x(x-t)^2g(t)\,dt. Therefore Yb−Y0=−x22∫0bg(t)dt+x∫0btg(t)dt−12∫0bt2g(t)dt.\boxed{Y_b-Y_0=-\frac{x^2}{2}\int_0^b g(t)\,dt +x\int_0^b t g(t)\,dt-\frac 12\int_0^b t^2g(t)\,dt.} The three integrals are constants with respect to xx, even when b<0b<0. The difference is a quadratic polynomial and hence homogeneous.

Step 4: Count genuine solution freedom. Once a basis and lower limit are fixed, the invertible derivative matrix uniquely fixes the parameter derivatives, and the specified parameter values fix their integrals. Allowing arbitrary integration constants adds a constant linear combination of the three homogeneous basis functions. Adding a separate homogeneous triple afterward repeats that same freedom rather than creating six independent constants. Particular solutions need not coincide, but their differences are homogeneous; the full solution space of y‴=gy'''=g is one particular solution plus the three-dimensional homogeneous space.

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

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