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