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