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