dsolve(diff(y(x),x) = 1/(y(x)+2+(3*x+1)^(1/2)),y(x))