ODE
\[ (4 y(x)+x) y'(x)-y(x)+4 x=0 \] ODE Classification
[[_homogeneous, `class A`], _rational, [_Abel, `2nd type`, `class A`]]
Book solution method
Equation linear in the variables, \(y'(x)=f\left ( \frac {X_1}{X_2} \right ) \)
Mathematica ✓
cpu = 0.0375279 (sec), leaf count = 31
\[\text {Solve}\left [2 \log \left (\frac {y(x)^2}{x^2}+1\right )+\tan ^{-1}\left (\frac {y(x)}{x}\right )+4 \log (x)=c_1,y(x)\right ]\]
Maple ✓
cpu = 0.013 (sec), leaf count = 35
\[ \left \{ -{\frac {1}{2}\ln \left ( {\frac {{x}^{2}+ \left ( y \left ( x \right ) \right ) ^{2}}{{x}^{2}}} \right ) }-{\frac {1}{4}\arctan \left ( {\frac {y \left ( x \right ) }{x}} \right ) }-\ln \left ( x \right ) -{\it \_C1}=0 \right \} \] Mathematica raw input
DSolve[4*x - y[x] + (x + 4*y[x])*y'[x] == 0,y[x],x]
Mathematica raw output
Solve[ArcTan[y[x]/x] + 4*Log[x] + 2*Log[1 + y[x]^2/x^2] == C[1], y[x]]
Maple raw input
dsolve((x+4*y(x))*diff(y(x),x)+4*x-y(x) = 0, y(x),'implicit')
Maple raw output
-1/2*ln((x^2+y(x)^2)/x^2)-1/4*arctan(y(x)/x)-ln(x)-_C1 = 0