Hello itsthewoo,
I finally found the solution. It kept running in my head and I couldn't put it aside, so here it is. First of all I used an advanced technique called Lie point transformations to find the answer. It comes down for this equation to find an integrating factor so the equation becomes exact. Writing the equation as:
\left(2x\right)\cdot dy-\left(-y+y^2\sqrt{x-x^2y^2}\right)\cdot dx=0
The integrating factor is the following:
M=\frac{1}{xy^2\sqrt{x-x^2y^2}}
The new (exact) equation is now:
\left(\frac{2x}{xy^2\sqrt{x-x^2y^2}}\right)\cdot dy+<br />
\left(\frac{1}{xy\sqrt{x-x^2y^2}}-\frac{1}{x}\right)\cdot dx=0
Which has as solution:
x=Ke^{-\frac{2\sqrt{x-x^2y^2}}{xy}}
This is the same solution as maxima (like maple or mathematica) gives. If any more info is necessary, I would recommend the following book:
"Ordinary differential equations, an elementary text-book with an introduction
to Lie's theory of the group of one parameter." written by James Morris Page
Hope this helps,
coomast