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