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