[[_1st_order, _with_linear_symmetries]]
Book solution method
Homogeneous equation, isobaric equation
Mathematica ✓
cpu = 0.0737663 (sec), leaf count = 25
Maple ✓
cpu = 0.072 (sec), leaf count = 47
DSolve[x^3 + y'[x] == x*Sqrt[x^4 + 4*y[x]],y[x],x]
Mathematica raw output
{{y[x] -> 2*E^(2*C[1])*(2*E^(2*C[1]) + x^2)}}
Maple raw input
dsolve(diff(y(x),x)+x^3 = x*(x^4+4*y(x))^(1/2), y(x),'implicit')
Maple raw output
((-_C1+y(x))*(x^4+4*y(x))^(1/2)-x^2*(_C1+y(x)))/(x^2+(x^4+4*y(x))^(1/2)) = 0