dsolve(diff(y(x),x) = -y(x)*(tan(x)+ln(2*x)*x-ln(2*x)*x^2*y(x))/x/tan(x),y(x))