THIS IS AN OLD PAGE DO NOT USE.
Go to User:Egm6321.f09.Team2/HW2 for latest information
Problem Statement: Find the Euler Integrating Factor h(y) for a non-linear 1st order ODE for the particular case where we assume hxN=0
Two conditions must be satisfied in order for a nonlinear 1st order ODE to be exact.
The first condition is that
must be in the form
The second condition is that
. This example will show how to satisfy the second condition.
Start with equation
which already satisfies the first condition.
If we assume the equation is not exact we can multiply it by an unkown factor
(Euler Integrating Factor) to attempt to find a factor that will make the equation exact.
The equation becomes:
Set
and rearrange,
To solve for
we will let
which implies that
since
, rearrage terms
take the integral of both sides,
If the right hand side of the integral is only a function of y then it can be written as
You have lost a negative sign in your expression for
--Egm6321.f09.TA 02:38, 28 September 2009 (UTC)
Problem Statement: Given the Non-homogeneous linear 1st order ODE with VC
show the steps to obtain
where C=Const.
The first step is to check if the given equation is in the proper form:
; Yes [pg.(8-1) Eg.(1)]
If
for all (x) then:
In this form we can use the Integration Factor Method to solve for y
Multiply by
where
Satisfy the exactness condition
where:
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
instead of
. --Egm6321.f09.TA 02:42, 28 September 2009 (UTC)
Assume hyM=0, so our equation becomes;
Let
and
and
Substitute back into the right hand side of the equationto verify h is a function of (x) only,
Using that result we can now solve for
Now that we have
we can multiply through our original equation to obtain,
Integrating yields:
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)
Problem Statement: Show that
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 ![{\displaystyle M(x,y)+N(x,y)y'=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c016c53d543695e82d466f48466412a1895a88a1)
(2)![{\displaystyle M_{y}=N_{x}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9be8caccbdf6e11e4c3665e179bae54e6b1ce2fa)
If
, then
where k is a constant.
Similarly:
, then
. Replacing M and N in our initial form
![{\displaystyle M(x,y)dx+N(x,y)dy=[a(x)y+k]dx+c(x)dy=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/149ab9f33bbbcadbde73128d82e92556e44a2148)
Dividing by
and re-arranging:
. Our initial problem meets condition 1, where
![{\displaystyle {\overline {b}}(x)={\frac {1}{2}}x^{2}=N(x,y)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/67a8feb0066ea1f518bf66ba313689eeda4f0461)
where ![{\displaystyle k=10}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b698dab3ec76554ed1b958de53897071b95f5bdb)
For condition 2,
and ![{\displaystyle N_{x}={\frac {d}{dx}}({\frac {1}{2}}x^{2})=x}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0cad0a2340130044f97c943f5e5435a1e14ca20c)
But these two are not equal. Therefore, using the Integrating Factors Method, we must find
such that
and thereby satisfy condition 2.
Differentiating and rearranging terms gives us:
. If we assume
since
and
are both functions of x, and rearrange terms, we get:
. Integrating both sides creates
. Solving for h:
^
. We know
and
. Therefore
^
.
Multiplying
and
by
causes our function to meet the 2nd condition and therefore be exact.
^
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:
.--Egm6321.f09.TA 02:47, 28 September 2009 (UTC)
Problem Statement: Show that
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
and
(2) ![{\displaystyle M_{y}=N_{x}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9be8caccbdf6e11e4c3665e179bae54e6b1ce2fa)
It is obvious from simple observation that our equation meets condition 1; however it fails to meet condition 2.
and
then
and
. Clearly
.
To meet condition 2, we must find h(x,y) such that
. As
and
only differ in
, let us assume that
. Differentiating and rearranging terms gives us:
. Integrating both sides creates
. Solving for h:
^
. We know
and
. Substituting for
and
and simplifying:
^
. This, multiplied by our original equation will cause it to be an exact, nonlinear, first order ODE.
The
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:
![{\displaystyle {\overline {b}}(x)c(y)y'+a(x){\overline {c}}(y)=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0e03f8214f883b84d5970f07551533d18a8a4ad6)
where ![{\displaystyle {\overline {b}}(x)=\int _{}^{x}b(x)\,dx}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fc12b66d313b915ba75edf5d0c0eb5a54084874d)
and ![{\displaystyle {\overline {c}}(y)=\int _{}^{y}c(y)\,dy}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8b357c65a106476365ad250ae31e66417d64b10f)
If we assume:
![{\displaystyle a(x)=5x^{3}+2}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7b5093c25804564cdf94e0d366283fe0d90b9148)
and
then
and
![{\displaystyle {\overline {c}}(y)=\int _{}^{y}y^{4}\,dy={\frac {1}{5}}y^{5}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/00fb1e445f3eb6e83d17c44308c543b5545c4a5e)
Since our equation does fit the form for this class, it is an exact, nonlinear, first order ODE.
Need to verify that
--Egm6321.f09.TA 02:51, 28 September 2009 (UTC)
Problem Statement: Show that the second exactness condition for
is satisfied.
The 2nd exactness condition actually consists of 2 parts, and our equation must meet both:
(1) fxx + (2p) fxy + (p2) fyy = gxp + p(gyp) - gy
(2) fxp + p (fyp) + 2 (fy) = gpp
where p is y'
In our equation,
f(x,y,p) =
and g(x,y,p) =
Looking at our equation, for the first half of the condition:
fxx =
[
(
)] =
(
) = 0
fxy =
[
(
)] =
(
) = 1
fyy =
[
(
)] =
(
) = 0
gxp =
[
(
)] =
(
) = 2p
gyp =
[
(
)] =
(
) = 1
gy =
(
) = p
Plugging back into the first half of our condition,
![{\displaystyle 0+2p(1)+p<sup>2</sup>(0)=2p+p(1)-p}](https://wikimedia.org/api/rest_v1/media/math/render/svg/16e1fe6b4bbbbe3bf447716f42746ee7eea731ac)
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:
fxp =
[
(
)] =
(
) = 0
fyp =
[
(
)] =
(
) = 0
fy =
(
)] = x
gpp =
[
(
)] =
(
) =
Plugging back into the second half of our condition (remembering that p = y'),
thereby satisfying the second half of the condition.
What is
? This would be
and this should not appear in this problem.--Egm6321.f09.TA 02:57, 28 September 2009 (UTC)
Problem Statement: Derive
using the following relations:
(1) ![{\displaystyle f(x,y,p):=\phi _{p}(x,y,p)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d16be4178358dcfc2fd7356459e8474a10e45a3e)
(2) ![{\displaystyle g(x,y,p):=\phi _{x}+\phi _{y}p}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3a2274837019e6f234a452df657e6513fccd5a75)
(3) ![{\displaystyle \phi _{xp}=\phi _{px}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e09a974ab64dd2b45791a7a8da65537698cc40b8)
(4) ![{\displaystyle \phi _{yp}=\phi _{py}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2fc73cc942aa0b8592118c629d2578eba8408ad2)
Start by taking the partial derivative of g with respect to p twice:
![{\displaystyle g_{p}=\phi _{xp}+\phi _{yp}p+\phi _{y}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/459f2e71d4890959a34e2338b684a444a5729837)
![{\displaystyle g_{pp}=\phi _{xpp}+\phi _{ypp}p+\phi _{yp}+\phi _{yp}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/575ce58214368cd3343c4a1fe77926afb7265e7c)
Now use (1), (3), and (4) to substitute
,
, and
by (1) and (3)
![{\displaystyle f_{xp}=\phi _{xpp}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4a09be709b0dba5312db0b79f91c5a7eeb14e3dd)
by (1) and (4)
![{\displaystyle f_{yp}=\phi _{ypp}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/df5d88ffad1775cc532bbf70794e27d961797757)
Therefore,
Problem Statement: Derive
using the following relations:
(1) ![{\displaystyle f(x,y,p):=\phi _{p}(x,y,p)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d16be4178358dcfc2fd7356459e8474a10e45a3e)
(2) ![{\displaystyle g(x,y,p):=\phi _{x}+\phi _{y}p}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3a2274837019e6f234a452df657e6513fccd5a75)
(3) ![{\displaystyle \phi _{px}=\phi _{xp}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c753dd9260c256946628699673894bf853444f20)
(4) ![{\displaystyle \phi _{py}=\phi _{yp}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a03a0e4413292548956ec5e510feea9122f3d135)
(5) ![{\displaystyle \phi _{xy}=\phi _{yx}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/78f63eeb6fbbea88d7abc74577005b15fe4fb7b4)
Solving (1) for
and taking the partial with respect to
yields
![{\displaystyle \phi _{px}=f_{x}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a679fd1032e2e4fa690beeaa05fd7b345f896db3)
Then solving (2) for
and taking the partial with respect to
yields
![{\displaystyle \phi _{xp}=g_{p}-\phi _{yp}p-\phi _{y}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b337512bd81163377ad66c1aa1fb9f24e595dfe5)
Equating these two by (3) and solving for
yields
![{\displaystyle \phi _{y}=-f_{x}+g_{p}-\phi _{yp}p}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9ee9e40232b67ce8da641c70ed167d04496df625)
Performing a similar task, take (1), solve for
and take the partial with respect to ![{\displaystyle y}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b8a6208ec717213d4317e666f1ae872e00620a0d)
![{\displaystyle \phi _{py}=f_{y}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/614ae15d323b361c68f7d83242bd47cde0252aee)
Then solving (2) for
and taking the partial with respect to
yields
![{\displaystyle \phi _{yp}={\frac {g_{p}}{p}}-{\frac {g}{p^{2}}}-{\frac {\phi _{xp}}{p}}+{\frac {\phi _{x}}{p^{2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/49837535b2559626c92f2b8e6397c8c277bd4ab7)
Equating these two by (4) and solving for
yields
![{\displaystyle \phi _{x}=f_{y}p^{2}-g_{p}p+g+\phi _{xp}p}](https://wikimedia.org/api/rest_v1/media/math/render/svg/42fd42d518f3b54259de4d67e422dca19fcc8f52)
Finally,
and
as previously shown and so the two equations are
(a) ![{\displaystyle \phi _{y}=-f_{x}+g_{p}-f_{y}p}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a7d5bace269449f8353a65681b6c7e5a9fa375c9)
(b) ![{\displaystyle \phi _{x}=f_{y}p^{2}-g_{p}p+g+f_{x}p}](https://wikimedia.org/api/rest_v1/media/math/render/svg/14fda55ca48c1343398f6b24251adbaa1e619ccc)
Taking the partial with respect to
for (a) and the partial with respect to
for (b) and equating based on (5) yields
![{\displaystyle -f_{xx}+g_{px}-f_{yx}p=f_{yy}p^{2}-g_{py}p+g_{y}+f_{xy}p}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8fb9e22d0e500ed02584e219e4e308d4fd31174e)
Noting
,
,
, and moving the f terms to the left and the g terms to the right
You have performed much unnecessary rearranging of
terms to get here. You need to rearrange until the end. Equate
using the appropriate relations and your solution will be must more straight forward.--Egm6321.f09.TA 03:01, 28 September 2009 (UTC)
Problem Statement: Equations 4&5 on p.(10-2).
Given:
Show that the Nonlinear 2nd-Order ODE is exact
Second Condition of Exactness: 10-2 Eq (4),(5)
The Following partial derivatives are then identified:
Page 10-2 Eq(4) is then satisfied.
Applying the results in Page 10-2 Eq(5):
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,
is Exactness Second ODE is an exact, nonlinear second ODE
Problem Statement: Verify the exactness of the ODE
.
Given
Prove that the equation is not exact.
First Condition of Exactness Page10-1 Eq(1):
This satisfy satisfies the First Condition of Exactness.
The Second Exactness Condition for a second order ODE is as follows:
The following partial derivatives are found:
Using these values in Eq(4):
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:
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.
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)