DSolve[Derivative[1][y][x] == (2 + Sqrt[1 + 3*x] + y[x])^(-1),y[x],x]