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