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