University of Florida/Egm4313/s12.team11.perez.gp/R3.2

Problem Statement

Developing the second homogeneous solution for the case of double real root as a limiting case of distinct roots.

Given

Consider two distinct roots of the form:

${\lambda }_1=x$ and ${\lambda }_2=x+ϵ$

(where $ϵ$ is perturbation).

Given

Find the homogeneous L2-ODE-CC having the above distinct roots.

Solution

$(\lambda -{\lambda }_1)(\lambda -({\lambda }_1))=0$

$(\lambda -x)(\lambda -(x+ϵ))=0$

${\lambda }^2-\lambda x-\lambda ϵ-\lambda x+x^2+xϵ=0$

${\lambda }^2-\lambda (2x+ϵ)+x(x+ϵ)=0$

$\therefore y^{″}-y^{\prime }(2\lambda +ϵ)+y\lambda (\lambda +ϵ)=0$ (1)

Given

Show that $\frac{e^{(\lambda +ϵ)x}-e^{\lambda x}}{ϵ}$ is a homogeneous solution. (2)

Solution

Let's find the corresponding derivatives:

$y=\frac{e^{(\lambda +ϵ)x}-e^{\lambda x}}{ϵ}$

$y^{\prime }=\frac{(\lambda +ϵ)e^{\lambda +ϵx}-\lambda e^{\lambda x}}{ϵ}$

$y^{″}=\frac{(\lambda +ϵ{)}^2{e}^{(\lambda +ϵ)x}-{\lambda }^2{e}^{\lambda x}}{ϵ}$

If we now take these three equations and plug them into the homogeneous L2-ODE-CC (1), we get:

$e^{(\lambda +ϵ)x}({\lambda }^2+2\lambda ϵ+{ϵ}^2-2{\lambda }^2-2\lambda ϵ-\lambda ϵ-{ϵ}^2+{\lambda }^2+ϵ\lambda )+e^{\lambda x}(-{\lambda }^2+2{\lambda }^2+\lambda ϵ-{\lambda }^2-\lambda ϵ)=0$

$\therefore 0\equiv 0.$

Since the left and right hand sides of the equation are zero, the solution is in fact a homogeneous equation.

Given

Find the limit of the homogeneous solution in (2) as epsilon approaches zero (think l'Hopital's Rule).

Solution

Using l'Hopital's Rule,

$\underset{ϵ\to 0}{\lim }\frac{e^{(\lambda +ϵ)x}-e^{\lambda x}}{ϵ}=\frac{0}{0}$

(this is an indeterminate form).

L'Hopital's Rule states that we can divide this function into two functions, $f(ϵ)$ and $g(ϵ)$, and then find their derivatives and attempt to find the limit of $\frac{f^{\prime }(ϵ)}{g^{\prime }(ϵ)}$. If a limit exists for this, then a limit exists for our original function.

$\underset{ϵ\to 0}{\lim }\frac{f^{\prime }(ϵ)}{g^{\prime }(ϵ)}=\underset{ϵ\to 0}{\lim }\frac{x{e}^{x(\lambda +ϵ)}}{1}$

$=\frac{x{e}^{(\lambda +0)x}}{1}$

$=x{e}^{x\lambda }.$

Given

Take the derivative of $e^{\lambda x}$ with respect to lambda.

Solution

Taking the derivative with respect to lambda, we find that:

$\frac{d(e^{\lambda x})}{d\lambda }=x{e}^{\lambda x}$.

It is important to remember that we must hold $x$ as a constant when finding this derivative.

Given

Compare the results in parts (3) and (4), and relate to the result by using variation of parameters

Solution

Taking a closer look at Parts 3 and 4 of this problem, we discover that they're in fact equal:

$\frac{d(e^{\lambda x})}{d\lambda }=\underset{e\to 0}{\lim }\frac{e^{(\lambda +ϵ)x}-e^{\lambda x}}{ϵ}=x{e}^{\lambda x}$

Given

Numerical experiment: Compute (2) setting lambda equal to 5 and epsilon equal to 0.001 </math>, and compare to the value obtained from the exact second homogeneous solution.

Solution

After performing these calculations, from (2) we get 148.478.

And from the exact second homogeneous solution, we get 200.05.