4.18.36 (3x+1)y(x)23(y(x)+2)y(x)+9=0

ODE
(3x+1)y(x)23(y(x)+2)y(x)+9=0 ODE Classification

[[_1st_order, _with_linear_symmetries], _Clairaut]

Book solution method
Clairaut’s equation and related types, f(yxy,y)=0

Mathematica
cpu = 0.386552 (sec), leaf count = 150

{{y(x)(sinh(2c1)+cosh(2c1))((3x34)cosh(c12)3(x12)sinh(c12))2+8sinh(c1)+8cosh(c1)18x294sinh(c1)+cosh(c1)36},{y(x)(sinh(2c1)+cosh(2c1))((3x34)cosh(c12)3(x12)sinh(c12))28sinh(c1)8cosh(c1)+18x+294sinh(c1)+cosh(c1)36}}

Maple
cpu = 0.03 (sec), leaf count = 37

{(y(x))2+4y(x)12x=0,y(x)=9+(1+3x)_C126_C13_C1} Mathematica raw input

DSolve[9 - 3*(2 + y[x])*y'[x] + (1 + 3*x)*y'[x]^2 == 0,y[x],x]

Mathematica raw output

{{y[x] -> -((-294 - 18*x + 8*Cosh[C[1]] + 8*Sinh[C[1]] + Sqrt[((-34 + 3*x)*Cosh[
C[1]/2] - 3*(-12 + x)*Sinh[C[1]/2])^2*(Cosh[2*C[1]] + Sinh[2*C[1]])])/(-36 + Cos
h[C[1]] + Sinh[C[1]]))}, {y[x] -> (294 + 18*x - 8*Cosh[C[1]] - 8*Sinh[C[1]] + Sq
rt[((-34 + 3*x)*Cosh[C[1]/2] - 3*(-12 + x)*Sinh[C[1]/2])^2*(Cosh[2*C[1]] + Sinh[
2*C[1]])])/(-36 + Cosh[C[1]] + Sinh[C[1]])}}

Maple raw input

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

Maple raw output

y(x)^2+4*y(x)-12*x = 0, y(x) = 1/3*(9+(1+3*x)*_C1^2-6*_C1)/_C1