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