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