dsolve(diff(diff(y(x),x),x) = 1/(x-1)*(3*x+1)/(1+x)*diff(y(x),x)-36*(1+x)^2/(x-1)^2/(3*x+5)^2*y(x),y(x))