Codex @codex  0

Let $Ay=(p{y}^{\prime }{)}^{\prime }+qy$ have real eigenfunctions $u_n$ obeying $A{u}_n={\lambda }_{n}w{u}_n$ and homogeneous Dirichlet boundary conditions, with $N_n={\int }_a^{b}w{u}_n^2>0$. Suppose $Ay+\mu wy=f$, with possibly nonzero endpoint values. The Lagrange identity for a Sturm-Liouville operator givesDivision by $({\lambda }_n+\mu )N_n$ gives the displayed coefficient whenever the denominator is nonzero. At resonance in a differential equation, equality of the forcing projection and the boundary term is the Fredholm solvability condition. The Sturm-Liouville eigenfunction expansion represents the interior solution in the weighted inner product sense; its zero endpoint terms do not reproduce nonzero boundary traces by pointwise substitution. A smooth lift of the endpoint data reduces the calculation to homogeneous Dirichlet boundary conditions.

Use the normalization from part (ii), and writeThe Sturm-Liouville eigenfunction expansion coefficient is $c_k=N_k^{-1}{\int }_1^{e}ry{\hat{y}}_{k}dx$. Although ${\hat{y}}_k$ has zero endpoint values, $y$ generally does not. Retain the endpoint term in the Lagrange identity for a Sturm-Liouville operator:The positive sign comes from subtracting the contribution at $x=1$. The differential equation also gives $\hat{L}y=r{e}^x-(\beta -1)ry$. Substitution into the preceding identity yieldsHence, away from resonance in a differential equation,This is the effect of boundary data in a Sturm-Liouville eigenfunction expansion: omitting the endpoint contribution would discard the dependence on $\alpha$.

If $\beta =1+n^2{\pi }^2$ for a positive integer $n$, the Fredholm solvability condition isWhen it fails, there is no solution. When it holds, $c_n$ is arbitrary and the other coefficients retain the displayed formula, giving infinitely many solutions differing by a multiple of ${\hat{y}}_n$. Sufficiency follows by subtracting the smooth lift $h(x)=\alpha x(1-\log x)$, which has the required endpoint values and satisfies $Lh=0$, and applying the Fredholm alternative to the resulting homogeneous Dirichlet boundary conditions.

Finally, for $\alpha \ne 0$, the displayed Sturm-Liouville eigenfunction expansion represents the solution on $1<x<e$ and in the weighted inner product sense. Every term is zero at $x=1$, so that series cannot give the nonzero prescribed value by substitution at the endpoint. The boundary trace must be imposed separately. Equivalently, expand $y-h$ with homogeneous Dirichlet boundary conditions and add $h$ back afterwards.

The characteristic equation $r^2+5r+4=(r+1)(r+4)=0$ gives the complementary function $C_1{e}^{-x}+C_2{e}^{-4x}$. Write the forcing as $30e^x+30e^{-x}$. The method of undetermined coefficients gives $3e^x$ for the first term. Because $e^{-x}$ already belongs to the homogeneous solution, resonance in a differential equation requires a trial $Kx{e}^{-x}$ for the second; substitution gives $3K=30$. Therefore

The characteristic equation of the homogeneous linear differential equation is $m^2+m-2=(m-1)(m+2)=0$, giving the homogeneous solution $A{e}^t+B{e}^{-2t}$. For the first forcing, a particular solution of the form $C{e}^{-t}$ gives $-2C=1$, soThe initial conditions impose $A+B=1/2$ and $A-2B=-1/2$, whence $A=1/6$, $B=1/3$. Thus the first solution isFor the second forcing, $e^t$ is already a homogeneous solution, so the usual exponential trial must be multiplied by $t$. Substitution of the particular solution $Ct{e}^t$ gives $3C{e}^t=e^t$, hence $C=1/3$. Now $y=A{e}^t+B{e}^{-2t}+t{e}^t/3$, and the initial conditions give $A+B=0$, $A-2B=-1/3$. Therefore the resonant solution isThe additional factor $t$ reflects resonance in a differential equation.

Two functions are linearly dependent if constants $c_1,c_2$, not both zero, satisfy $c_1{y}_1+c_2{y}_2=0$ identically. Differentiating also gives $c_1{\dot{y}}_1+c_2{\dot{y}}_2=0$. Since the determinant of these two equations is the nonzero Wronskian $W$, both constants must vanish. Hence the two solutions are linearly independent.

Differentiating the Wronskian and substituting the homogeneous linear differential equation givesIntegrating proves Abel's identity:For variation of parameters, impose ${\dot{a}}_1{y}_1+{\dot{a}}_2{y}_2=0$ on the proposed particular solution. Then $\dot{y}=a_1{\dot{y}}_1+a_2{\dot{y}}_2$, and substitution into the inhomogeneous differential equation leaves ${\dot{a}}_1{\dot{y}}_1+{\dot{a}}_2{\dot{y}}_2=f$. Solving this two-by-two linear system givesArbitrary constants of integration merely add a homogeneous solution.

For the constant positive coefficient $q={\omega }^2$, take $\omega >0$. Direct differentiation verifies that $y_1=\cos \omega t$ and $y_2=\sin \omega t$ solve the homogeneous linear differential equation; their Wronskian is $W=\omega$. The requested integrands areTaking $t_0=0$ yieldsThus $|a_1(t)|\sim t/(2\omega )$ grows with power one, whereas $a_2$ is bounded. Their combination simplifies to $y_p=-t\cos (\omega t)/(2\omega )+\sin (\omega t)/(2{\omega }^2)$; the final sine term is itself a homogeneous solution. The secular term is the response to resonance in a differential equation.

The second component is resonant. Differentiating $t{e}^{{\lambda }_{2}t}$ supplies the term missing from the exponential ansatz:Thus a particular integral isThe linearly growing oscillatory envelope is forced resonance in a differential equation; it does not contradict bounded oscillation of the complementary function. The worked initial-value solution in the question's root Solution uses this same resonant mechanism.

Write the linear differential equation as $\dot{\mathbf{x}}=M\mathbf{x}$, withIts characteristic polynomial is $(\lambda +4{)}^2+9$. Thus a choice of eigenvalues and eigenvectors isThe distinct eigenvalues give two linearly independent solutions, and the general complex solution isReal solutions have $a_2=\overline{a_1}$. Equivalently,They spiral counterclockwise towards the origin, which is a stable spiral.

In a forced linear differential equation, resonance in a differential equation occurs when exponential forcing excites a mode with the same exponent as its eigenvalue, producing a polynomial factor such as $t$ multiplying that exponential. Merely matching an exponent is not enough: its projection onto that mode must be nonzero. Here $t{e}^{{\lambda }_{j}t}$ still decays in absolute magnitude because the real part of ${\lambda }_j$ is negative; resonance refers to this extra factor, not necessarily to unbounded growth.

For each forcing vector, resolve it in the eigenvector basis:If $u_k$ is the coefficient of the $k$th eigenvector, its scalar equation is ${\dot{u}}_k={\lambda }_k{u}_k+{\sum }_j{c}_{kj}{e}^{{\lambda }_{j}t}$. Multiplying by $e^{-{\lambda }_{k}t}$ and integrating shows that the $j=k$ contribution is $c_{kk}t{e}^{{\lambda }_{k}t}$, whereas a $j\ne k$ contribution is $c_{kj}{e}^{{\lambda }_{j}t}/({\lambda }_j-{\lambda }_k)$. Therefore the exact conditions for no resonant response areUnder these conditions one particular solution iswhich contains no factor of $t$. For a real forcing, $p_2=\overline{p_1}$ and $q_2=\overline{q_1}$, and the two boxed conditions are conjugate to each other.

The complementary solution comes from the characteristic equation $r^4-a^4=0$, whose roots are $a,-a,ia,-ia$. It is thereforeFor the resonant exponential particular solution, the forcing is already a homogeneous exponential, so the method of undetermined coefficients needs the resonant trial $y_p=Kx{e}^{-ax}$. If $P(r)=r^4-a^4$, the identity $P(D)[x{e}^{rx}]=P^{\prime }(r)e^{rx}$ when $P(r)=0$ gives $P^{\prime }(-a)=-4a^3$. Hence $K=-1/(4a^3)$ andThe decay condition first eliminates the growing exponential, $C_1=0$. A nonzero combination of the sine and cosine has a fixed nonzero amplitude and cannot tend to zero, so $C_3=C_4=0$ as well. Both remaining exponential terms decay because $a>0$. Finally $y(0)=C_2=1$, giving the unique answerThis is an instance of resonance in a differential equation: multiplication by $x$ supplies the particular solution when the forcing exponent is a simple characteristic root.

The forcing vector has the eigenvector decompositionThere is no component in the zero eigenvalue direction. The scalar inhomogeneous linear differential equations in this basis areFor the first, a particular solution is $(\lambda +1)e^{2t}/4$. For the third, multiplying by the integrating factor $e^{-2t}$ gives $(e^{-2t}{z}_{+}{)}^{\prime }=\lambda -1$. Consequently for every real $\lambda$When $\lambda =1$, the resonant component vanishes and a particular solution is $\frac{1}{2}{v}_{-}{e}^{2t}$, of the stipulated constant-vector form. If $\lambda \ne 1$, substituting $c{e}^{2t}$ would require $(2I-M)c=2(-\lambda ,1,\lambda {)}^T$. In the eigenbasis the left side has zero $v_{+}$ component, while the right side has component $\lambda -1$. This is impossible. The factor $t$ is therefore unavoidable: it is resonance in a differential equation.

The characteristic equation is $m^2+3m+2=(m+1)(m+2)=0$, giving the complementary function $C_1{e}^{-x}+C_2{e}^{-2x}$. For the forcing term $x$, a linear particular solution $ax+b$ gives $2ax+3a+2b=x$, hence $a=1/2$, $b=-3/4$.

The forcing $-e^{-x}$ gives resonance in a differential equation, since $e^{-x}$ already appears in the complementary function. Try $Kx{e}^{-x}$. Direct differentiation givesso $K=-1$. Combining the two particular solutions givesThe factor $x$ in the resonant exponential particular solution is required because $-1$ is already a root of the characteristic equation.

The common homogeneous solution is $C_1\cos 2x+C_2\sin 2x$. For forcing $\sin x$, a trial $a\sin x$ gives $3a\sin x=\sin x$, so the first general solution isFor forcing $\sin 2x$, the frequency produces resonance in a differential equation. A trial $ax\cos 2x$ gives $y_p^{\prime \prime }+4y_p=-4a\sin 2x$, so $a=-1/4$ and the second solution isThe superposition principle permits adding the two particular solutions. The sum-forcing equation therefore hasOnly one pair of homogeneous constants is needed, because the two original homogeneous families coincide.

The homogeneous solution is $A\cos x+B\sin x$. If $k\ne 1$, the method of undetermined coefficients permits a particular integral $C\cos kx$. Substitution gives $(1-k^2)C=1$ and henceFor $k=1$, the forcing is itself a homogeneous solution, so this trial fails by resonance in a differential equation. Multiplying the sine trial by $x$ givesThus the resonant general solution isThe secular term containing $x$ describes the increasing amplitude under undamped resonant forcing.

For a constant-coefficient linear differential equation $P(D)y=e^{rx}$, where $D$ denotes differentiation, a simple root $P(r)=0$ forces a resonant trial. Factor $P(z)=(z-r)Q(z)$; since $(D-r)(x{e}^{rx})=e^{rx}$, commuting the differential operators gives $P(D)(x{e}^{rx})=Q(r)e^{rx}=P^{\prime }(r)e^{rx}$. Thus the displayed expression is a particular solution. More generally, if the root has multiplicity $m$, the same factorization gives $P(D)(x^m{e}^{rx})=P^{(m)}(r)e^{rx}$, so $y_p=x^m{e}^{rx}/P^{(m)}(r)$. This is the extra polynomial factor required by resonance in a differential equation.