Exact differential equation explained

In mathematics, an exact differential equation or total differential equation is a certain kind of ordinary differential equation which is widely used in physics and engineering.

Definition

Given a simply connected and open subset D of

R2

and two functions I and J which are continuous on D, an implicit first-order ordinary differential equation of the form

I(x,y)dx+J(x,y)dy=0,

is called an exact differential equation if there exists a continuously differentiable function F, called the potential function,[1] [2] so that

\partialF
\partialx

=I

and
\partialF
\partialy

=J.

An exact equation may also be presented in the following form:

I(x,y)+J(x,y)y'(x)=0

where the same constraints on I and J apply for the differential equation to be exact.

The nomenclature of "exact differential equation" refers to the exact differential of a function. For a function

F(x0,x1,...,xn-1,xn)

, the exact or total derivative with respect to

x0

is given by
dF=
dx0
\partialF
\partialx0
n
+\sum
i=1
\partialF
\partialxi
dxi
dx0

.

Example

The function

F:R2\toR

given by

F(x,y)=

1
2

(x2+y2)+c

is a potential function for the differential equation

xdx+ydy=0.

First order exact differential equations

Identifying first order exact differential equations

Let the functions M, N, M_y, and N_x, where the subscripts denote the partial derivative with respect to the relative variable, be continuous in the region R: \alpha < x < \beta, \gamma < y < \delta. Then the differential equation

M(x,y)+N(x,y)

dy
dx

=0

is exact if and only if

My(x,y)=Nx(x,y)

That is, there exists a function

\psi(x,y)

, called a potential function, such that

\psix(x,y)=M(x,y)and\psiy(x,y)=N(x,y)

So, in general:

My(x,y)=Nx(x,y)\iff \begin{cases} \exists\psi(x,y)\\ \psix(x,y)=M(x,y)\\ \psiy(x,y)=N(x,y) \end{cases}

Proof

The proof has two parts.

First, suppose there is a function

\psi(x,y)

such that

\psix(x,y)=M(x,y)and\psiy(x,y)=N(x,y)

It then follows that

My(x,y)=\psixy(x,y)andNx(x,y)=\psiyx(x,y)

Since

My

and

Nx

are continuous, then

\psixy

and

\psiyx

are also continuous which guarantees their equality.

The second part of the proof involves the construction of

\psi(x,y)

and can also be used as a procedure for solving first-order exact differential equations. Suppose that

My(x,y)=Nx(x,y)

and let there be a function

\psi(x,y)

for which

\psix(x,y)=M(x,y)and\psiy(x,y)=N(x,y)

Begin by integrating the first equation with respect to

x

. In practice, it doesn't matter if you integrate the first or the second equation, so long as the integration is done with respect to the appropriate variable.

\frac(x, y) = M(x, y)\psi(x, y) = \int + h(y)\psi(x, y) = Q(x, y) + h(y)

where

Q(x,y)

is any differentiable function such that

Qx=M

. The function

h(y)

plays the role of a constant of integration, but instead of just a constant, it is function of

y

, since

M

is a function of both

x

and

y

and we are only integrating with respect to

x

.

Now to show that it is always possible to find an

h(y)

such that

\psiy=N

.\psi(x, y) = Q(x, y) + h(y)

Differentiate both sides with respect to

y

.\frac(x, y) = \frac(x, y) + h'(y)

Set the result equal to

N

and solve for

h'(y)

.h'(y) = N(x, y) - \frac(x, y)

In order to determine

h'(y)

from this equation, the right-hand side must depend only on

y

. This can be proven by showing that its derivative with respect to

x

is always zero, so differentiate the right-hand side with respect to

x

. \frac(x, y) - \frac\frac(x, y) \iff \frac(x, y) - \frac\frac(x, y)

Since

Qx=M

,\frac(x, y) - \frac(x, y)Now, this is zero based on our initial supposition that

My(x,y)=Nx(x,y)

Therefore,h'(y) = N(x, y) - \frac(x, y)h(y) = \int

\psi(x, y) = Q(x, y) + \int + C

And this completes the proof.

Solutions to first order exact differential equations

First order exact differential equations of the formM(x, y) + N(x, y)\frac = 0

can be written in terms of the potential function

\psi(x,y)

\frac + \frac\frac = 0

where\begin\psi _x(x, y) = M(x, y)\\\psi _y(x, y) = N(x, y)\end

This is equivalent to taking the exact differential of

\psi(x,y)

.\frac + \frac\frac = 0 \iff \frac\psi(x, y(x)) = 0

The solutions to an exact differential equation are then given by\psi(x, y(x)) = c

and the problem reduces to finding

\psi(x,y)

.

This can be done by integrating the two expressions

M(x,y)dx

and

N(x,y)dy

and then writing down each term in the resulting expressions only once and summing them up in order to get

\psi(x,y)

.

The reasoning behind this is the following. Since\begin\psi _x(x, y) = M(x, y)\\\psi _y(x, y) = N(x, y)\end

it follows, by integrating both sides, that\begin\psi(x, y) = \int + h(y) = Q(x, y) + h(y)\\\psi(x, y) = \int + g(x) = P(x, y) + g(x)\end

Therefore,Q(x, y) + h(y) = P(x, y) + g(x)

where

Q(x,y)

and

P(x,y)

are differentiable functions such that

Qx=M

and

Py=N

.

In order for this to be true and for both sides to result in the exact same expression, namely

\psi(x,y)

, then

h(y)

must be contained within the expression for

P(x,y)

because it cannot be contained within

g(x)

, since it is entirely a function of

y

and not

x

and is therefore not allowed to have anything to do with

x

. By analogy,

g(x)

must be contained within the expression

Q(x,y)

.

Ergo,Q(x, y) = g(x) + f(x, y) \text P(x, y) = h(y) + d(x, y)

for some expressions

f(x,y)

and

d(x,y)

.Plugging in into the above equation, we find thatg(x) + f(x, y) + h(y) = h(y) + d(x, y) + g(x) \Rightarrow f(x, y) = d(x, y)and so

f(x,y)

and

d(x,y)

turn out to be the same function. Therefore,Q(x, y) = g(x) + f(x, y) \text P(x, y) = h(y) + f(x, y)

Since we already showed that\begin\psi(x, y) = Q(x, y) + h(y)\\\psi(x, y) = P(x, y) + g(x)\end

it follows that\psi(x, y) = g(x) + f(x, y) + h(y)

So, we can construct

\psi(x,y)

by doing

\int{M(x,y)dx}

and

\int{N(x,y)dy}

and then taking the common terms we find within the two resulting expressions (that would be

f(x,y)

) and then adding the terms which are uniquely found in either one of them -

g(x)

and

h(y)

.

Second order exact differential equations

The concept of exact differential equations can be extended to second order equations.[3] Consider starting with the first-order exact equation:

I\left(x,y\right)+J\left(x,y\right){dy\overdx}=0

Since both functions

J\left(x,y\right)

are functions of two variables, implicitly differentiating the multivariate function yields

{dI\overdx}+\left({dJ\overdx}\right){dy\overdx}+{d2y\overdx2}\left(J\left(x,y\right)\right)=0

Expanding the total derivatives gives that

{dI\overdx}={\partialI\over\partialx}+{\partialI\over\partialy}{dy\overdx}

and that

{dJ\overdx}={\partialJ\over\partialx}+{\partialJ\over\partialy}{dy\overdx}

Combining the terms gives

{\partialI\over\partialx}+{dy\overdx}\left({\partialI\over\partialy}+{\partialJ\over\partialx}+{\partialJ\over\partialy}{dy\overdx}\right)+{d2y\overdx2}\left(J\left(x,y\right)\right)=0

If the equation is exact, then Additionally, the total derivative of

J\left(x,y\right)

is equal to its implicit ordinary derivative . This leads to the rewritten equation

{\partialI\over\partialx}+{dy\overdx}\left({\partialJ\over\partialx}+{dJ\overdx}\right)+{d2y\overdx2}\left(J\left(x,y\right)\right)=0

Now, let there be some second-order differential equation

f\left(x,y\right)+g\left(x,y,{dy\overdx}\right){dy\overdx}+{d2y\overdx2}\left(J\left(x,y\right)\right)=0

If

{\partialJ\over\partialx}={\partialI\over\partialy}

for exact differential equations, then

\int\left({\partialI\over\partialy}\right)dy=\int\left({\partialJ\over\partialx}\right)dy

and

\int\left({\partialI\over\partialy}\right)dy=\int\left({\partialJ\over\partialx}\right)dy=I\left(x,y\right)-h\left(x\right)

where

h\left(x\right)

is some arbitrary function only of

x

that was differentiated away to zero upon taking the partial derivative of

I\left(x,y\right)

with respect to

y

. Although the sign on

h\left(x\right)

could be positive, it is more intuitive to think of the integral's result as

I\left(x,y\right)

that is missing some original extra function

h\left(x\right)

that was partially differentiated to zero.

Next, if

{dI\overdx}={\partialI\over\partialx}+{\partialI\over\partialy}{dy\overdx}

then the term

{\partialI\over\partialx}

should be a function only of

x

and

y

, since partial differentiation with respect to

x

will hold

y

constant and not produce any derivatives of

y

. In the second order equation

f\left(x,y\right)+g\left(x,y,{dy\overdx}\right){dy\overdx}+{d2y\overdx2}\left(J\left(x,y\right)\right)=0

only the term

f\left(x,y\right)

is a term purely of

x

and

y

. Let

{\partialI\over\partialx}=f\left(x,y\right)

. If

{\partialI\over\partialx}=f\left(x,y\right)

, then

f\left(x,y\right)={dI\overdx}-{\partialI\over\partialy}{dy\overdx}

Since the total derivative of

I\left(x,y\right)

with respect to

x

is equivalent to the implicit ordinary derivative

{dI\overdx}

, then

f\left(x,y\right)+{\partialI\over\partialy}{dy\overdx}={dI\overdx}={d\overdx}\left(I\left(x,y\right)-h\left(x\right)\right)+{dh\left(x\right)\overdx}

So,

{dh\left(x\right)\overdx}=f\left(x,y\right)+{\partialI\over\partialy}{dy\overdx}-{d\overdx}\left(I\left(x,y\right)-h\left(x\right)\right)

and

h\left(x\right)=\int\left(f\left(x,y\right)+{\partialI\over\partialy}{dy\overdx}-{d\overdx}\left(I\left(x,y\right)-h\left(x\right)\right)\right)dx

Thus, the second order differential equation

f\left(x,y\right)+g\left(x,y,{dy\overdx}\right){dy\overdx}+{d2y\overdx2}\left(J\left(x,y\right)\right)=0

is exact only if

g\left(x,y,{dy\overdx}\right)={dJ\overdx}+{\partialJ\over\partialx}={dJ\overdx}+{\partialJ\over\partialx}

and only if the below expression

\int\left(f\left(x,y\right)+{\partialI\over\partialy}{dy\overdx}-{d\overdx}\left(I\left(x,y\right)-h\left(x\right)\right)\right)dx=\int\left(f\left(x,y\right)-{\partial\left(I\left(x,y\right)-h\left(x\right)\right)\over\partialx}\right)dx

is a function solely of

x

. Once

h\left(x\right)

is calculated with its arbitrary constant, it is added to

I\left(x,y\right)-h\left(x\right)

to make

I\left(x,y\right)

. If the equation is exact, then we can reduce to the first order exact form which is solvable by the usual method for first-order exact equations.

I\left(x,y\right)+J\left(x,y\right){dy\overdx}=0

Now, however, in the final implicit solution there will be a

C1x

term from integration of

h\left(x\right)

with respect to

x

twice as well as a

C2

, two arbitrary constants as expected from a second-order equation.

Example

Given the differential equation

\left(1-x2\right)y''-4xy'-2y=0

one can always easily check for exactness by examining the

y''

term. In this case, both the partial and total derivative of

1-x2

with respect to

x

are

-2x

, so their sum is

-4x

, which is exactly the term in front of

y'

. With one of the conditions for exactness met, one can calculate that

\int\left(-2x\right)dy=I\left(x,y\right)-h\left(x\right)=-2xy

Letting

f\left(x,y\right)=-2y

, then

\int\left(-2y-2xy'-{d\overdx}\left(-2xy\right)\right)dx=\int\left(-2y-2xy'+2xy'+2y\right)dx=\int\left(0\right)dx=h\left(x\right)

So,

h\left(x\right)

is indeed a function only of

x

and the second order differential equation is exact. Therefore,

h\left(x\right)=C1

and

I\left(x,y\right)=-2xy+C1

. Reduction to a first-order exact equation yields
2\right)y'=0
-2xy+C
1+\left(1-x

Integrating

I\left(x,y\right)

with respect to

x

yields
2y+C
-x
1x+i\left(y\right)=0

where

i\left(y\right)

is some arbitrary function of

y

. Differentiating with respect to

y

gives an equation correlating the derivative and the

y'

term.

-x2+i'\left(y\right)=1-x2

So,

i\left(y\right)=y+C2

and the full implicit solution becomes

C1x+C

2y=0
2+y-x

Solving explicitly for

y

yields

y=

C1x+C2
1-x2

Higher order exact differential equations

The concepts of exact differential equations can be extended to any order. Starting with the exact second order equation

{d2y\overdx2}\left(J\left(x,y\right)\right)+{dy\overdx}\left({dJ\overdx}+{\partialJ\over\partialx}\right)+f\left(x,y\right)=0

it was previously shown that equation is defined such that

f\left(x,y\right)={dh\left(x\right)\overdx}+{d\overdx}\left(I\left(x,y\right)-h\left(x\right)\right)-{\partialJ\over\partialx}{dy\overdx}

Implicit differentiation of the exact second-order equation

n

times will yield an

\left(n+2\right)

th order differential equation with new conditions for exactness that can be readily deduced from the form of the equation produced. For example, differentiating the above second-order differential equation once to yield a third-order exact equation gives the following form

{d3y\overdx3}\left(J\left(x,y\right)\right)+{d2y\overdx2}{dJ\overdx}+{d2y\overdx2}\left({dJ\overdx}+{\partialJ\over\partialx}\right)+{dy\overdx}\left({d2J\overdx2}+{d\overdx}\left({\partialJ\over\partialx}\right)\right)+{df\left(x,y\right)\overdx}=0

where

{df\left(x,y\right)\overdx}={d2h\left(x\right)\overdx2}+{d2\overdx2}\left(I\left(x,y\right)-h\left(x\right)\right)-{d2y\overdx2}{\partialJ\over\partialx}-{dy\overdx}{d\overdx}\left({\partialJ\over\partialx}\right)=F\left(x,y,{dy\overdx}\right)

and where

F\left(x,y,{dy\overdx}\right)

is a function only of

x,y

and

{dy\overdx}

. Combining all

{dy\overdx}

and

{d2y\overdx2}

terms not coming from

F\left(x,y,{dy\overdx}\right)

gives

{d3y\overdx3}\left(J\left(x,y\right)\right)+{d2y\overdx2}\left(2{dJ\overdx}+{\partialJ\over\partialx}\right)+{dy\overdx}\left({d2J\overdx2}+{d\overdx}\left({\partialJ\over\partialx}\right)\right)+F\left(x,y,{dy\overdx}\right)=0

Thus, the three conditions for exactness for a third-order differential equation are: the

{d2y\overdx2}

term must be

2{dJ\overdx}+{\partialJ\over\partialx}

, the

{dy\overdx}

term must be

{d2J\overdx2}+{d\overdx}\left({\partialJ\over\partialx}\right)

and

F\left(x,y,{dy\overdx}\right)-{d2\overdx2}\left(I\left(x,y\right)-h\left(x\right)\right)+{d2y\overdx2}{\partialJ\over\partialx}+{dy\overdx}{d\overdx}\left({\partialJ\over\partialx}\right)

must be a function solely of

x

.

Example

Consider the nonlinear third-order differential equation

yy'''+3y'y''+12x2=0

If

J\left(x,y\right)=y

, then

y''\left(2{dJ\overdx}+{\partialJ\over\partialx}\right)

is

2y'y''

and

y'\left({d2J\overdx2}+{d\overdx}\left({\partialJ\over\partialx}\right)\right)=y'y''

which together sum to

3y'y''

. Fortunately, this appears in our equation. For the last condition of exactness,

F\left(x,y,{dy\overdx}\right)-{d2\overdx2}\left(I\left(x,y\right)-h\left(x\right)\right)+{d2y\overdx2}{\partialJ\over\partialx}+{dy\overdx}{d\overdx}\left({\partialJ\over\partialx}\right)=12x2-0+0+0=12x2

which is indeed a function only of

x

. So, the differential equation is exact. Integrating twice yields that
4+C
h\left(x\right)=x
1x+C

2=I\left(x,y\right)

. Rewriting the equation as a first-order exact differential equation yields
4+C
x
1x+C

2+yy'=0

Integrating

I\left(x,y\right)

with respect to

x

gives that

{x5\over

2+C
5}+C
2x+i\left(y\right)=0
. Differentiating with respect to

y

and equating that to the term in front of

y'

in the first-order equation gives that

i'\left(y\right)=y

and that

i\left(y\right)={y2\over2}+C3

. The full implicit solution becomes

{x5\over

2+C
5}+C
2x+C
2\over
3+{y

2}=0

The explicit solution, then, is

2+C
y=\pm\sqrt{C
2x+C
3-2x5
5
}

See also

Further reading

Notes and References

  1. Book: Wolfgang Walter. Ordinary Differential Equations. 11 March 2013. Springer Science & Business Media. 978-1-4612-0601-9.
  2. Book: Vladimir A. Dobrushkin. Applied Differential Equations: The Primary Course. 16 December 2014. CRC Press. 978-1-4987-2835-5.
  3. Book: Ordinary Differential Equations: An Elementary Textbook for Students of Mathematics, Engineering and the Sciences. limited. Tenenbaum. Morris. Pollard. Harry. Dover. 1963. 0-486-64940-7. New York. 248. Solution of the Linear Differential Equation with Nonconstant Coefficients. Reduction of Order Method..