dsolve(diff(diff(y(x),x),x) = 1/x*(7*a*x^2+5)/(a*x^2+1)*diff(y(x),x)-(15*a*x^2+5)/x^2/(a*x^2+1)*y(x),y(x))