DSolve[-((-1 + b)*b) + a*x^k + x^2*(y[x]^2 + Derivative[1][y][x]) == 0,y[x],x]