dsolve((3*x+1)*diff(y(x),x)^2-3*(y(x)+2)*diff(y(x),x)+9 = 0,y(x))