# University of Florida/Egm4313/s12.team11.imponenti/R3.1

## Report 3, Problem 1

### Problem Statement

Find the solution to the following L2-ODE-CC: ${\displaystyle y''-10y'+25y=r(x)\!}$

With the following excitation: ${\displaystyle r(x)=7e^{5x}-2x^{2}\!}$

And the following initial conditions: ${\displaystyle y(0)=4,y'(0)=-5\!}$

Plot this solution and the solution in the example on p.7-3

### Homogeneous Solution

To find the homogeneous solution we need to find the roots of our equation

${\displaystyle \lambda ^{2}-10\lambda +25=0\!}$ ${\displaystyle (\lambda -5)(\lambda -5)=0\!}$ ${\displaystyle \lambda =5\!}$

We know the homogeneous solution for the case of a real double root with ${\displaystyle \lambda =5\!}$ to be

${\displaystyle y_{h}=c_{1}e^{5x}+c_{2}xe^{5x}\!}$

### Particular Solution

For the given excitation we must use the Sum Rule to the particular solution as follows

${\displaystyle y_{p}=y_{p1}+y_{p2}\!}$ where ${\displaystyle y_{p1}\!}$ and ${\displaystyle y_{p2}\!}$ are the solutions to ${\displaystyle r_{1}(x)=7e^{5x}\!}$ and ${\displaystyle r_{2}(x)=-2x^{2}\!}$, respectively

#### First Particular Solution

${\displaystyle r_{1}(x)=7e^{5x}\!}$,

from table 2.1, K 2011, pg. 82 we have

${\displaystyle y_{p1}=Ce^{5x}\!}$

but this corresponds to one of our homogeneous solutions so we must use the modification rule to get

${\displaystyle y_{p1}=Cx^{2}e^{5x}\!}$

Plugging this into the original L2-ODE-CC then substituting;

${\displaystyle y_{p1}''-10y_{p1}'+25y_{p1}=r_{1}(x)\!}$

${\displaystyle (Cx^{2}e^{5x})''-10(Cx^{2}e^{5x})'+25(Cx^{2}e^{5x})=r_{1}(x)\!}$

${\displaystyle 25Cx^{2}e^{5x}+10Cxe^{5x}+10Cxe^{5x}+2Ce^{5x}-10(5Cx^{2}e^{5x}+2Cxe^{5x})+25Cx^{2}e^{5x}=7e^{5x}\!}$

${\displaystyle e^{5x}[25Cx^{2}+10Cx+10Cx+2C-50Cx^{2}-20Cx+25Cx^{2}]=7e^{5x}\!}$

${\displaystyle e^{5x}2C=7e^{5x}\!}$

so ${\displaystyle C={\frac {7}{2}}\!}$ and the first particular solution is,

${\displaystyle y_{p1}={\frac {7}{2}}x^{2}e^{5x}\!}$

#### Second Particular Solution

${\displaystyle r_{2}(x)=-2x^{2}\!}$,

from table 2.1, K 2011, pg. 82 we have

${\displaystyle y_{p2}=a_{2}x^{2}+a_{1}x+a_{0}\!}$

Plugging this into the original L2-ODE-CC then substituting;

${\displaystyle y_{p2}''-10y_{p2}'+25y_{p2}=r_{2}(x)\!}$

${\displaystyle (a_{2}x^{2}+a_{1}x+a_{0})''-10(a_{2}x^{2}+a_{1}x+a_{0})'+25(a_{2}x^{2}+a_{1}x+a_{0})=-2x^{2}\!}$

${\displaystyle 2a_{2}-10(2a_{2}x+a_{1})+25(a_{2}x^{2}+a_{1}x+a_{0})=-2x^{2}\!}$

grouping like terms we get three equations to solve for the three unknowns, these are written in matrix form

${\displaystyle 25a_{2}x^{2}+(25a_{1}-20a_{2})x+(25a_{0}-10a_{1}+2a_{2})=-2x^{2}\!}$

${\displaystyle {\begin{bmatrix}2&-10&25\\-20&25&0\\25&0&0\end{bmatrix}}{\begin{bmatrix}a_{2}\\a_{1}\\a_{0}\end{bmatrix}}={\begin{bmatrix}0\\0\\-2\end{bmatrix}}\!}$

solving by back subsitution leads to ${\displaystyle a_{2}=-{\frac {2}{25}},a_{1}={\frac {8}{125}},a_{0}={\frac {4}{125}}\!}$

so the second particular solution is,

${\displaystyle y_{p2}=-{\frac {2}{25}}x^{2}+{\frac {8}{125}}x+{\frac {4}{125}}\!}$

### General Solution

The general solution is the summation of the homogeneous and particular solutions

${\displaystyle y=y_{h}+y_{p1}+y_{p2}\!}$

${\displaystyle y=c_{1}e^{5x}+c_{2}xe^{5x}+{\frac {7}{2}}x^{2}e^{5x}-{\frac {2}{25}}x^{2}+{\frac {8}{125}}x+{\frac {4}{125}}\!}$

${\displaystyle y=e^{5x}(c_{1}+c_{2}x+{\frac {7}{2}}x^{2})-{\frac {2}{25}}x^{2}+{\frac {8}{125}}x+{\frac {4}{125}}\!}$

Applying the first initial condition ${\displaystyle y(0)=4\!}$

${\displaystyle y(0)=c_{1}+{\frac {4}{125}}=4\!}$

${\displaystyle c_{1}={\frac {496}{125}}\!}$

Second initial condition ${\displaystyle y'(0)=-5\!}$

${\displaystyle y'={\frac {d}{dx}}y=e^{5x}[5(c_{1}+c_{2}x+{\frac {7}{2}}x^{2})+c_{2}+7x]-{\frac {4}{25}}x+{\frac {8}{125}}\!}$

${\displaystyle y'=e^{5x}[{\frac {35}{2}}x^{2}+(5c_{2}+7)x+5c_{1}+c_{2}]-{\frac {4}{25}}x+{\frac {8}{125}}\!}$

${\displaystyle y'(0)=5c_{1}+c_{2}+{\frac {8}{125}}=-5\!}$

${\displaystyle c_{2}=-{\frac {3113}{125}}\!}$

The general solution to the differential equation is therefore

                    ${\displaystyle y(x)=e^{5x}({\frac {496}{125}}-{\frac {3113}{125}}x+{\frac {7}{2}}x^{2})-{\frac {2}{25}}x^{2}+{\frac {8}{125}}x+{\frac {4}{125}}\!}$


### Plot

Below is a plot of this solution and the solution to in the example on p.7-3

our solution ${\displaystyle y(x)=e^{5x}({\frac {496}{125}}-{\frac {3113}{125}}x+{\frac {7}{2}}x^{2})-{\frac {2}{25}}x^{2}+{\frac {8}{125}}x+{\frac {4}{125}}\!}$ (shown in red)

example on p.7-3 ${\displaystyle y(x)=e^{5x}(4-25*x+*{\frac {7}{2}}x^{2})\!}$ (shown in blue)

Egm4313.s12.team11.imponenti 22:31, 20 February 2012 (UTC)