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