University of Florida/Egm6321/f09.Team2/HW2

THIS IS AN OLD PAGE DO NOT USE.

Go to User:Egm6321.f09.Team2/HW2 for latest information

Problem #1

Problem Statement: Find the Euler Integrating Factor h(y) for a non-linear 1st order ODE for the particular case where we assume h_xN=0

---

Two conditions must be satisfied in order for a nonlinear 1^st order ODE to be exact.

The first condition is that $F(x,y,y^{\prime })=0$ must be in the form $M(x,y)+N(x,y)y^{\prime }=0$

The second condition is that $M_y=N_x$. This example will show how to satisfy the second condition.

Start with equation $M(x,y)+N(x,y)y^{\prime }=0$ which already satisfies the first condition.

If we assume the equation is not exact we can multiply it by an unkown factor $h(x,y)$ (Euler Integrating Factor) to attempt to find a factor that will make the equation exact.

The equation becomes:

$h(x,y)[M(x,y)dx+N(x,y)dy]=0$

$(hM)dx+(hN)dy=0,wherehM=\bar{M},andhN=\bar{N}$

${\bar{M}}_y=(hM{)}_y=h_{y}M+h{M}_y$

${\bar{N}}_x=(hN{)}_y=h_{x}N+h{N}_x$

Set ${\bar{M}}_y={\bar{N}}_x$ and rearrange,

$h_{x}N-h_{y}M+h(N_x-M_y)=0$

To solve for $h_y$ we will let $h_{x}N=0$ which implies that $h_x=0$ since $M\ne 0$

$-h_{y}M+h(N_x-M_y)=0$, rearrage terms

$\frac{h_y}{h}dy=\frac{1}{M}(N_x-M_y)$ take the integral of both sides,

$\int \frac{h_y}{h}dy=\int \frac{1}{M}(N_x-M_y)dy$

If the right hand side of the integral is only a function of y then it can be written as $g(y)$

$\frac{1}{M}(N_x-M_y)=g(y)$

$ln\mid h\mid =\int g(y)dy$

$h(y)=exp\underset{}{\overset{y}{\int }}g(s)ds$

You have lost a negative sign in your expression for $g(y)$--Egm6321.f09.TA 02:38, 28 September 2009 (UTC)

Problem 2

Problem Statement: Given the Non-homogeneous linear 1^st order ODE with VC $y^{\prime }+\frac{1}{x}y=x^2$ show the steps to obtain $y=\frac{x^3}{4}+\frac{C}{x}$ where C=Const.

---

The first step is to check if the given equation is in the proper form:

$P_{(x)}{y}^{\prime }+Q_{(x)}y=R_{(x)}$ ; Yes [pg.(8-1) Eg.(1)]

If $P_{(x)}\ne 0$ for all (x) then:

$y^{\prime }+\frac{Q_{(x)}}{P_{(x)}}y=\frac{R_{(x)}}{P(x)}$

In this form we can use the Integration Factor Method to solve for y

$N_{(x,y)}\frac{dy}{dx}+M_{(x,y)}=0$

Multiply by $h_{(x,y)}$

$h_{(x,y)}[N_{(x,y)}dy+M_{(x,y)}dx]$

$(hM)dx+(hN)dy=0$ where $(hM)=\bar{M};(hN)=\bar{N}$

Satisfy the exactness condition ${\bar{M}}_y={\bar{N}}_x$ where:

${\bar{M}}_y=(hM{)}_y=h_{y}M+h{M}_y$

${\bar{N}}_x=(hN{)}_x=h_{x}N+h{N}_x$

$h_{x}N-h_{y}M+h(Nx-My)=0$

To solve for h(x,y) we must make an assumption that h is a function of (x) only.

You do not have to make this assumption. You have chosen $h(x,y)=h(x)$ instead of $h(x,y)=h(y)$. --Egm6321.f09.TA 02:42, 28 September 2009 (UTC)

Assume h_yM=0, so our equation becomes;

$h_{x}N+h(N_x-M_y)=0$

$\frac{h_x}{h}=-\frac{1}{N_{(x,y)}}(N_x-M_y)$

Let $N_{(x,y)}=1$ and $M_{(x,y)}=\frac{Q_{(x)}}{P_{(x)}}y=\frac{1}{x}y$

$N_x=0$ and $M_y=\frac{1}{x}$

Substitute back into the right hand side of the equationto verify h is a function of (x) only,

$\frac{1}{1}(0-\frac{1}{x})=\frac{1}{x}=f(x)$

Using that result we can now solve for $h$

$\int \frac{h_x}{h}dx=\int \frac{1}{x}dx$

$ln|h|=ln|x|$

$h_x=x$

Now that we have $h_x$ we can multiply through our original equation to obtain,

$(x)y^{\prime }+(x)\frac{1}{x}y=(x)x^2=x{y}^{\prime }+y=x^3$

Integrating yields:

$xy=\frac{x^4}{4}+C$ divide by x,

Show the steps - the terms must be grouped together as the derivative of a product before you can integrate this. It is not directly integrable in this form.--Egm6321.f09.TA 02:42, 28 September 2009 (UTC)

$y=\frac{x^3}{4}+\frac{C}{x}$

Problem #3

Problem Statement: Show that $\frac{1}{2}{x}^2{y}^{\prime }+[x^{4}y+10]=0$ is exact.

---

In this problem, we are asking if: (1) this equation is only a function of x, and (2) it is "exact", meaning exact or exactly integrable by the Integrating Factors Method.

To be exact, it must meet 2 conditions:

(1)it must meet the form $M(x,y)+N(x,y)y^{\prime }=0$
(2)$M_y=N_x$

If $M(x,y)dy=a(x)$ , then $M(x,y)=a(x)y+k$ where k is a constant.

Similarly:
$N(x,y)dx=b(x)$ , then $N(x,y)=\underset{}{\overset{x}{\int }}b(s)ds=c(x)=\bar{b}(x)$. Replacing M and N in our initial form

$M(x,y)dx+N(x,y)dy=[a(x)y+k]dx+c(x)dy=0$

Dividing by $dx$ and re-arranging:
$[a(x)y+k]+\bar{b}(x)y^{\prime }=0$. Our initial problem meets condition 1, where
$\bar{b}(x)=\frac{1}{2}{x}^2=N(x,y)$
$a(x)+k=x^{4}y+10=M(x,y)$ where $k=10$

For condition 2, $M_y=\frac{d}{dy}[a(x)+k]=\frac{d}{dy}(x^{4}y+10)=x^4$ and $N_x=\frac{d}{dx}(\frac{1}{2}{x}^2)=x$

But these two are not equal. Therefore, using the Integrating Factors Method, we must find $h(x,y)$ such that $(hM{)}_y=(hN{)}_x$ and thereby satisfy condition 2.
Differentiating and rearranging terms gives us:
$h_{x}N-h_{y}M+h(N_x-M_y)=0$. If we assume $h_y=0$ since $M_y$ and $N_x$ are both functions of x, and rearrange terms, we get:
$\frac{h_x}{h}=\frac{-1}{N}(N_x-M_y)dx$. Integrating both sides creates $log(h)$. Solving for h:
$h=e$^$(-2)(\underset{}{\overset{x}{\int }}{N}_x-M_y,dx$. We know $N_x$ and $M_y$. Therefore $h=e$^$(-2)(\frac{1}{2}{x}^2-\frac{1}{5}{x}^5)$.
Multiplying $M(x,y)$ and $N(x,y)$ by $h$ causes our function to meet the 2nd condition and therefore be exact.

$(e$^$(\frac{2}{5}{x}^5-\frac{1}{2}{x}^2)(\frac{1}{2}{x}^2{y}^{\prime }+[x^{4}y+10]=0$ is exact.

This equation is incomplete. There are parenthesis missing. Where are the differentials?--Egm6321.f09.TA 02:47, 28 September 2009 (UTC)

You need to show that you have arrived at the correct integrating factor. Specifically: $\frac{\partial \bar{M}}{\partial y}=\frac{\partial \bar{N}}{\partial x}$.--Egm6321.f09.TA 02:47, 28 September 2009 (UTC)

Problem #4

Problem Statement: Show that $(\frac{1}{3}{x}^3)(y^4)y^{\prime }+(5x^3+2)(\frac{1}{5}{y}^5)=0$ is an exact nonlinear, first order ODE.

---

To prove this equation is an exact, non-linear, 1st order ODE, it must meet 2 conditions:
(1) it must fit the form $M(x,y)y^{\prime }+N(x,y)=0$ and
(2) $M_y=N_x$

It is obvious from simple observation that our equation meets condition 1; however it fails to meet condition 2.
$M(x,y)=\frac{1}{3}{x}^3{y}^4$ and $N(x,y)=(5x^3-2)(\frac{1}{5}{y}^5)$ then
$M_y=\frac{4}{3}{x}^3{y}^3$ and $N_x=5x^3{y}^4$. Clearly $M_y\ne N_x$.

To meet condition 2, we must find h(x,y) such that $(hM{)}_y=(hN{)}_x$. As $M_y$ and $N_x$ only differ in $y$, let us assume that $h_x=0$. Differentiating and rearranging terms gives us:
$\frac{h_y}{h}=\frac{1}{M}(N_x-M_y)dx$. Integrating both sides creates $logh$. Solving for h:
$h=e$^$(\frac{1}{M}\underset{}{\overset{y}{\int }}{N}_x-M_y,dy)$. We know $N_x$ and $M_y$. Substituting for $N_x$ and $M_y$ and simplifying:
$h=e$^$(15y-4(log(y))$. This, multiplied by our original equation will cause it to be an exact, nonlinear, first order ODE.

The $\frac{1}{M}$ should be inside the integral.--Egm6321.f09.TA 02:51, 28 September 2009 (UTC)

The simpliest method to prove this equation as an exact, nonlinear, first order ODE is to show that (1) it fits the form for a class of exact, non-linear 1st order ODE's. Specifically:

$\bar{b}(x)c(y)y^{\prime }+a(x)\bar{c}(y)=0$

where $\bar{b}(x)=\underset{}{\overset{x}{\int }}b(x)dx$
and $\bar{c}(y)=\underset{}{\overset{y}{\int }}c(y)dy$

If we assume:
$a(x)=5x^3+2$
$b(x)=x^2$ and
$c(y)=y^4$ then
$\bar{b}(x)=\underset{}{\overset{x}{\int }}{x}^{2}dx=\frac{1}{3}{x}^3$ and $\bar{c}(y)=\underset{}{\overset{y}{\int }}{y}^{4}dy=\frac{1}{5}{y}^5$

Since our equation does fit the form for this class, it is an exact, nonlinear, first order ODE.

Need to verify that ${\bar{N}}_x={\bar{M}}_y$--Egm6321.f09.TA 02:51, 28 September 2009 (UTC)

Problem #5

Problem Statement: Show that the second exactness condition for $xy{y}^{″}+x(y^{\prime }{)}^2+y{y}^{\prime }=0$ is satisfied.

---

The 2nd exactness condition actually consists of 2 parts, and our equation must meet both:
(1) f_xx + (2p) f_xy + (p²) f_yy = g_xp + p(g_yp) - g_y

(2) f_xp + p (f_yp) + 2 (f_y) = g_pp

where p is y^'

In our equation,
f(x,y,p) = $xy$ and g(x,y,p) = $x(y^{\prime }{)}^2+y{y}^{\prime }$

Looking at our equation, for the first half of the condition:
f_xx = $\frac{d}{dx}$[$\frac{d}{dx}$($xy$)] = $\frac{d}{dx}$($y$) = 0

f_xy = $\frac{d}{dy}$[$\frac{d}{dx}$($xy$)] = $\frac{d}{dy}$($y$) = 1

f_yy = $\frac{d}{dy}$[$\frac{d}{dy}$($xy$)] = $\frac{d}{dy}$($x$) = 0

g_xp = $\frac{d}{dp}$[$\frac{d}{dx}$($x(p{)}^2+yp$)] = $\frac{d}{dp}$($p^2$) = 2p

g_yp = $\frac{d}{dp}$[$\frac{d}{dy}$($x(p{)}^2+yp$)] = $\frac{d}{dp}$($p$) = 1

g_y = $\frac{d}{dy}$($x(p{)}^2+yp$) = p

Plugging back into the first half of our condition,
$0+2p(1)+p²(0)=2p+p(1)-p$
$2p=2p$ and therefore the first half of the condition is satisfied.

What is sup? Check to verify that your code has rendered properly.--Egm6321.f09.TA 02:57, 28 September 2009 (UTC)

Looking at the 2nd half of the condition:
f_xp = $\frac{d}{dp}$[$\frac{d}{dy}$($xy$)] = $\frac{d}{dp}$($y$) = 0

f_yp = $\frac{d}{dp}$[$\frac{d}{dx}$($xy$)] = $\frac{d}{dp}$($x$) = 0

f_y = $\frac{d}{dy}$($xy$)] = x

g_pp = $\frac{d}{dp}$[$\frac{d}{dp}$($x(p{)}^2+yp$)] = $\frac{d}{dp}$($2xp+y$) = $2x(p^{\prime })$

Plugging back into the second half of our condition (remembering that p = y'),
$0+p(0)+2(x)p^{\prime }=2x(p^{\prime })$ thereby satisfying the second half of the condition.

What is $p^{\prime }$ ? This would be $y^{″}$ and this should not appear in this problem.--Egm6321.f09.TA 02:57, 28 September 2009 (UTC)

Problem #6

Problem Statement: Derive $f_{xp}+p{f}_{yp}+2f_y=g_{pp}$ using the following relations:

(1) $f(x,y,p):={\varphi }_p(x,y,p)$
(2) $g(x,y,p):={\varphi }_x+{\varphi }_{y}p$
(3) ${\varphi }_{xp}={\varphi }_{px}$
(4) ${\varphi }_{yp}={\varphi }_{py}$

---

Start by taking the partial derivative of g with respect to p twice: $g_{pp}=\frac{{\partial }^{2}g}{\partial p^2}$

$g_p={\varphi }_{xp}+{\varphi }_{yp}p+{\varphi }_y$
$g_{pp}={\varphi }_{xpp}+{\varphi }_{ypp}p+{\varphi }_{yp}+{\varphi }_{yp}$
$=>g_{pp}={\varphi }_{xpp}+{\varphi }_{ypp}p+2{\varphi }_{yp}$

Now use (1), (3), and (4) to substitute ${\varphi }_{xpp}$, ${\varphi }_{ypp}$, and ${\varphi }_{yp}$

$f_x={\varphi }_{px}={\varphi }_{xp}$ by (1) and (3)
$f_{xp}={\varphi }_{xpp}$
$f_y={\varphi }_{py}={\varphi }_{yp}$ by (1) and (4)
$f_{yp}={\varphi }_{ypp}$

Therefore,

$g_{pp}=f_{xp}+f_{yp}p+2f_y=f_{xp}+p{f}_{yp}+2f_y$

Clear and concise.--Egm6321.f09.TA 02:58, 28 September 2009 (UTC)

Problem #7

Problem Statement: Derive $f_{xx}+2p{f}_{xy}+p^2{f}_{yy}=g_{xp}+p{g}_{yp}-g_y$ using the following relations:
(1) $f(x,y,p):={\varphi }_p(x,y,p)$
(2) $g(x,y,p):={\varphi }_x+{\varphi }_{y}p$
(3) ${\varphi }_{px}={\varphi }_{xp}$
(4) ${\varphi }_{py}={\varphi }_{yp}$
(5) ${\varphi }_{xy}={\varphi }_{yx}$

---

Solving (1) for ${\varphi }_p$ and taking the partial with respect to $x$ yields
${\varphi }_{px}=f_x$
Then solving (2) for ${\varphi }_x$ and taking the partial with respect to $p$ yields
${\varphi }_{xp}=g_p-{\varphi }_{yp}p-{\varphi }_y$
Equating these two by (3) and solving for ${\varphi }_y$ yields
${\varphi }_y=-f_x+g_p-{\varphi }_{yp}p$
Performing a similar task, take (1), solve for ${\varphi }_p$ and take the partial with respect to $y$
${\varphi }_{py}=f_y$
Then solving (2) for ${\varphi }_y$ and taking the partial with respect to $p$ yields
${\varphi }_{yp}=\frac{g_p}{p}-\frac{g}{p^2}-\frac{{\varphi }_{xp}}{p}+\frac{{\varphi }_x}{p^2}$
Equating these two by (4) and solving for ${\varphi }_x$ yields
${\varphi }_x=f_y{p}^2-g_{p}p+g+{\varphi }_{xp}p$

Finally, $fx={\varphi }_{px}={\varphi }_{xp}$ and $fy={\varphi }_{py}={\varphi }_{yp}$ as previously shown and so the two equations are
(a) ${\varphi }_y=-f_x+g_p-f_{y}p$
(b) ${\varphi }_x=f_y{p}^2-g_{p}p+g+f_{x}p$
Taking the partial with respect to $x$ for (a) and the partial with respect to $y$ for (b) and equating based on (5) yields
$-f_{xx}+g_{px}-f_{yx}p=f_{yy}{p}^2-g_{py}p+g_y+f_{xy}p$

Noting $f_{yx}=f_{xy}$, $g_{px}=g_{xp}$, $g_{py}=g_{yp}$, and moving the f terms to the left and the g terms to the right

$f_{xx}+2p{f}_{xy}+p^2{f}_{yy}=g_{xp}+p{g}_{yp}-g_y$

You have performed much unnecessary rearranging of $\mathrm{\Phi }$ terms to get here. You need to rearrange until the end. Equate ${\mathrm{\Phi }}_{xy}={\mathrm{\Phi }}_{yx}$ using the appropriate relations and your solution will be must more straight forward.--Egm6321.f09.TA 03:01, 28 September 2009 (UTC)

Problem #8

Problem Statement: Equations 4&5 on p.(10-2).

---

Given:

$(8x^5{y}^{\text{'}})y^{{\text{'}}^{\prime }}+2x^2{y}^{\text{'}}+20x^4(y^{\text{'}}{)}^2+4xy=0$

Show that the Nonlinear 2nd-Order ODE is exact

Solution

Second Condition of Exactness: 10-2 Eq (4),(5)

$F_{xx}+2p{F}_{xy}+p^2{F}_{yy}=G_{xp}+p{G}_{yp}-G_y$

$F_{xp}+p{F}_{yp}+2F_y=G_{pp}$

$F(x,y,p)y^{{\text{'}}^{\prime }}+G(x,y,p)=0$

$F(x,y,p)=8x^{5}pG(x,y,p)=2x^{2}p+20x^4{p}^2+4xy$

The Following partial derivatives are then identified:

$F_x=40x^{4}p{G}_x=4xp+80x^3{p}^2+4y$

$F_{xx}=160x^{3}p{G}_{xp}=4x+160x^{3}p$

$F_{xy}=0G_y=4x$

$F_{xp}=40x^4{G}_{yp}=0$

$F_y=0G_p=2x^2+40x^{4}p$

$F_{yy}=0G_{pp}=40x^4$

$160x^{3}p+2p(0)+p^2(0)=4x+160x^{3}p+(0)-4x160x^{3}p=160x^{3}p$

Page 10-2 Eq(4) is then satisfied.

Applying the results in Page 10-2 Eq(5):

$40x^4+p(0)+2(0)=40x^4$

$40x^4=40x^4$

Page 10-2 Eq(5) is also satisfied.

Because Eq(4) and Eq(5) are both satisfied, then the Secound Exactness Condition is satisfied.

Therefore,

$F(x,y,p)=8x^{5}pG(x,y,p)=2x^{2}p+20x^4{p}^2+4xy$

is Exactness Second ODE Template:Font

Nice work. One caveat: For your final statement, express the equation in terms of $x,y,y^{\prime }$. The equation involves an unknown function $y(x)$, where $p$ was a dummy variable which helped us perform differentiations. --Egm6321.f09.TA 03:05, 28 September 2009 (UTC)

Problem #9

Problem Statement: Verify the exactness of the ODE $\sqrt{x}{y}^{″}+2x{y}^{\prime }+3y=0$.

---

Given

$\sqrt[]{x}{y}^{{\text{'}}^{\prime }}+2x{y}^{\text{'}}+3y=0$

Prove that the equation is not exact.

Solution

First Condition of Exactness Page10-1 Eq(1):

$F(x,y,p)y^{{\text{'}}^{\prime }}+G(x,y,p)=0$

$F(x,y,p)=\sqrt[]{x}G(x,y,p)=2xp+3y$

This satisfy Template:Font the First Condition of Exactness.

The Second Exactness Condition for a second order ODE is as follows:

$F_{xx}+2p{F}_{xy}+p^2{F}_{yy}=G_{xp}+p{G}_{yp}-G_{y}Page10-2Eq(4)$

$F_{xp}+p{F}_{yp}+2F_y=G_{pp}Page10-2Eq(5)$

The following partial derivatives are found:

$F_x=-\frac{1}{2}{x}^{\frac{-1}{2}}{G}_x=2p$

$F_{xx}=-\frac{1}{4}{x}^{\frac{-3}{2}}{G}_{xp}=2$

$F_{xy}=0G_y=3$

$F_{xp}=0G_{yp}=0$

$F_y=0G_p=2x$

$F_{yy}=0G_{pp}=0$

$F_{yp}=0$

Using these values in Eq(4):

$\frac{-1}{4}{x}^{\frac{-3}{2}}+2p(0)+p^2(0)=2+p(0)-3$

$\frac{-1}{4}{x}^{\frac{-3}{2}}=-1$

Page 10-2 Eq(4) is not satisfied therefore the ODE is not exact.

Add equation numbers to your solution here so that you can reference the equations in your solution, rather than a transparency that is located somewhere else on the web. This will make your solution self contained and complete in itself.--Egm6321.f09.TA 03:09, 28 September 2009 (UTC)

Page 10-2 Eq(5) can be populated to find:

$0=0+p(0)+2(0)$

$0=0$

Eq (5) is satisfied.

In order for the ODE to be exact it must satisfy both Eq(4)and Eq(5), since it fails to do this it is concluded that the ODE is not exact.

Contributing Team Members

Egm6321.f09.Team 2.walker 20:58, 20 September 2009 (UTC) (Walker, Matthew)

Joe Gaddone 14:04, 21 September 2009 (UTC)

Egm6321.f09.Team2.sungsik 19:20, 23 September 2009 (UTC)

Kumanchik 19:53, 23 September 2009 (UTC)