dsolve(diff(diff(y(x),x),x)-3*y(x)*diff(y(x),x)-3*a*y(x)^2-4*a^2*y(x)-b=0,y(x))