✓ Mathematica : cpu = 0.0062074 (sec), leaf count = 24
DSolve[y[x] - (x + 3*y[x])*Derivative[1][y][x] + (5 + 3*x)*Derivative[1][y][x]^2 == 0,y[x],x]
✓ Maple : cpu = 0.035 (sec), leaf count = 60
dsolve((3*x+5)*diff(y(x),x)^2-(3*y(x)+x)*diff(y(x),x)+y(x) = 0,y(x))