dsolve(diff(diff(y(x),x),x) = 2/x*(a*x+2*b)/(a*x+b)*diff(y(x),x)-(2*a*x+6*b)/(a*x+b)/x^2*y(x),y(x))