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