Solution of First Order DE using Integrating Factors

maverick280857
Messages
1,774
Reaction score
5
Hello everyone

I'm stuck with a problem from Simmons (not HW). The problem asks to solve the following equation by finding an integrating factor

(xy-1)dx + (x^2 - xy)dy = 0

I found by hit and trial that there can be no integrating factor that is a function of x alone or y alone. So how do I find an integrating factor for this equation? Any ideas?

Thanks
Vivek
 
Physics news on Phys.org
I'm not sure that integrating factors will be of any use here. But you can still solve it, i.e. by making the substitution

y = x + \sqrt{v}

this will give you

1 + \frac{1}{2}\frac{v^{\prime}}{\sqrt{v}} = \frac{1 - x^2 - x \sqrt{v}}{-x \sqrt{v}} = (x - \frac{1}{x}) \frac{1}{\sqrt{v}} + 1

which reduces to

v^{\prime} = 2(x - \frac{1}{x})

leaving you

v(x) = x^2 - 2 \ln{x} + C

where C is a constant. Feed this back into y to get

y(x) = x + \sqrt{x^2 - 2 \ln{x} + C}.
 
Last edited:
Thanks Matthew, interesting method.

I wonder why the questions asks to find an integrating factor though.
 
maverick280857 said:
Hello everyone

I'm stuck with a problem from Simmons (not HW). The problem asks to solve the following equation by finding an integrating factor

(xy-1)dx + (x^2 - xy)dy = 0

I found by hit and trial that there can be no integrating factor that is a function of x alone or y alone. So how do I find an integrating factor for this equation? Any ideas?

Thanks
Vivek

There is no obvious integrating factor, but there is one. If my understanding is correct, the problem you are given is worded wrong. You need to treat this problem as an Exact/Not Exact problem. Here's how you solve it:

M = (xy-1)
N = (x^2 - xy)

Now take the partial with respect to y of M, and the partial with respect to x of N:

M_y = x
N_x = 2x - y

Since

M_y is not equal to N_x, this equation is NOT exact, so you must find an integrating factor that makes it exact.
Check if the integrating factor can be only a function on x:

\frac{M_y - N_x}{N}\mu = \frac{d\mu}{dx}

Plugging in the corresponding values, you get:

\frac{x - 2x +y}{x^2-xy}\mu = \frac{d\mu}{dx}

\frac{-x +y}{-x ( -x + y )}\mu = \frac{d\mu}{dx}

-\frac{1}{x}\mu = \frac{d\mu}{dx}

-\frac{1}{x}dx = \frac{d\mu}{\mu}

-\ln(x) = \ln(mu)

e^{-\ln(x)} = e^{\ln(\mu)}

\mu = \frac{1}{x}

Now multiply both sides of the equation by the integrating factor:

(y-\frac{1}{x})dx + (x - y)dy = 0

Now check again if the equation is exact ( it must be, but just to make sure ):

M = y-\frac{1}{x}
N = x - y


M_y = 1
N_x = 1

Now it's exact:

(1)\psi(x,y) = y-\frac{1}{x}
(2) \psi(x,y) = x - y

\int(1) dx
(3)\int y-\frac{1}{x} = xy - \ln{x} + h(y)

Now differentiate (3):

(**)\frac{d}{dy}(y - \ln{x} + h(y)) = 1 + h\prime(y)

Set (2) = (**) :

x - y = 1 + h\prime(y)
h\prime(y) = x - y -1

Since h\prime(y) = x - y -1,

h(y) = xy - \frac{1}{2}y^2 - y

Now plug in h(y) into (3):

\psi(x,y) = xy - \ln{x} + xy - \frac{1}{2}y^2 - y
\psi(x,y) = 2xy - \ln{x} + \frac{1}{2}y^2 - y

So

xy - \ln{x} + xy - \frac{1}{2}y^2 - y = C

is a constant solution. I tried separating y from everything, but it looks nasty. Here, i'll post it anyway:

y = \sqrt(-2\ln{x} + x^2 - 2x - 2 (C - 2)) + x - 2

I hope I didn't mess up anywhere, it took sometime to type.
Hopefully this helps. Good Luck!
 
maverick280857 said:
Hello everyone

I'm stuck with a problem from Simmons (not HW). The problem asks to solve the following equation by finding an integrating factor

(xy-1)dx + (x^2 - xy)dy = 0


Thanks
Vivek

Answer given by the following online ODE solver at

http://www.fundp.ac.be/sciences/physique/administration/convode.html

i) integrating factor is 1/x
ii) possible solution set
{y=sqrt( - 2*log(x) - 2*cte + x ^2) + x, y=-sqrt( - 2*log(x) - 2*cte + x ^2) + x}
 
Last edited by a moderator:
Thread 'Direction Fields and Isoclines'
I sketched the isoclines for $$ m=-1,0,1,2 $$. Since both $$ \frac{dy}{dx} $$ and $$ D_{y} \frac{dy}{dx} $$ are continuous on the square region R defined by $$ -4\leq x \leq 4, -4 \leq y \leq 4 $$ the existence and uniqueness theorem guarantees that if we pick a point in the interior that lies on an isocline there will be a unique differentiable function (solution) passing through that point. I understand that a solution exists but I unsure how to actually sketch it. For example, consider a...

Similar threads

Back
Top