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