ODE
\[ a x y(x) y'(x)=x^2+y(x)^2 \] ODE Classification
[[_homogeneous, `class A`], _rational, _Bernoulli]
Book solution method
Homogeneous equation
Mathematica ✓
cpu = 0.257285 (sec), leaf count = 68
\[\left \{\left \{y(x)\to -\frac {\sqrt {x^2+(a-1) c_1 x^{2/a}}}{\sqrt {a-1}}\right \},\left \{y(x)\to \frac {\sqrt {x^2+(a-1) c_1 x^{2/a}}}{\sqrt {a-1}}\right \}\right \}\]
Maple ✓
cpu = 0.034 (sec), leaf count = 80
\[\left [y \left (x \right ) = \frac {\sqrt {\left (a -1\right ) \left (x^{\frac {2}{a}} \textit {\_C1} a -x^{\frac {2}{a}} \textit {\_C1} +x^{2}\right )}}{a -1}, y \left (x \right ) = -\frac {\sqrt {\left (a -1\right ) \left (x^{\frac {2}{a}} \textit {\_C1} a -x^{\frac {2}{a}} \textit {\_C1} +x^{2}\right )}}{a -1}\right ]\] Mathematica raw input
DSolve[a*x*y[x]*y'[x] == x^2 + y[x]^2,y[x],x]
Mathematica raw output
{{y[x] -> -(Sqrt[x^2 + (-1 + a)*x^(2/a)*C[1]]/Sqrt[-1 + a])}, {y[x] -> Sqrt[x^2 + (-1 + a)*x^(2/a)*C[1]]/Sqrt[-1 + a]}}
Maple raw input
dsolve(a*x*y(x)*diff(y(x),x) = x^2+y(x)^2, y(x))
Maple raw output
[y(x) = 1/(a-1)*((a-1)*(x^(2/a)*_C1*a-x^(2/a)*_C1+x^2))^(1/2), y(x) = -1/(a-1)*((a-1)*(x^(2/a)*_C1*a-x^(2/a)*_C1+x^2))^(1/2)]