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