dsolve(diff(y(x),x) = (a*y(x)^2+b*x^2)^3/a^(5/2)*x/(a*y(x)^2+b*x^2+a)/y(x),y(x))