dsolve((y(x)^2/b+x^2/a)*(y(x)*diff(y(x),x)+x)+(a-b)/(a+b)*(y(x)*diff(y(x),x)-x) = 0,y(x))