\[ \left \{x'(t)=x(t) \left (y(t)^2-z(t)^2\right ),y'(t)=-y(t) \left (x(t)^2+z(t)^2\right ),z'(t)=z(t) \left (x(t)^2+y(t)^2\right )\right \} \] ✗ Mathematica : cpu = 0.0439705 (sec), leaf count = 0 , could not solve
DSolve[{Derivative[1][x][t] == x[t]*(y[t]^2 - z[t]^2), Derivative[1][y][t] == -(y[t]*(x[t]^2 + z[t]^2)), Derivative[1][z][t] == (x[t]^2 + y[t]^2)*z[t]}, {x[t], y[t], z[t]}, t]
✓ Maple : cpu = 1.076 (sec), leaf count = 704
\[\left \{[\{x \relax (t ) = 0\}, \{y \relax (t ) = 0\}, \{z \relax (t ) = c_{1}\}], \left [\{x \relax (t ) = 0\}, \left \{y \relax (t ) = \frac {\sqrt {c_{1} {\mathrm e}^{2 c_{1} c_{2}} {\mathrm e}^{2 c_{1} t}-c_{1} {\mathrm e}^{4 c_{1} c_{2}} {\mathrm e}^{4 c_{1} t}}}{{\mathrm e}^{2 c_{1} c_{2}} {\mathrm e}^{2 c_{1} t}-1}, y \relax (t ) = -\frac {\sqrt {c_{1} {\mathrm e}^{2 c_{1} c_{2}} {\mathrm e}^{2 c_{1} t}-c_{1} {\mathrm e}^{4 c_{1} c_{2}} {\mathrm e}^{4 c_{1} t}}}{{\mathrm e}^{2 c_{1} c_{2}} {\mathrm e}^{2 c_{1} t}-1}\right \}, \left \{z \relax (t ) = \frac {\sqrt {-\left (\frac {d}{d t}y \relax (t )\right ) y \relax (t )}}{y \relax (t )}, z \relax (t ) = -\frac {\sqrt {-\left (\frac {d}{d t}y \relax (t )\right ) y \relax (t )}}{y \relax (t )}\right \}\right ], [\{x \relax (t ) = c_{1}\}, \{y \relax (t ) = i x \relax (t ), y \relax (t ) = -i x \relax (t )\}, \{z \relax (t ) = i x \relax (t ), z \relax (t ) = -i x \relax (t )\}], \left [\left \{x \relax (t ) = \frac {\sqrt {c_{1} {\mathrm e}^{2 c_{1} c_{2}} {\mathrm e}^{2 c_{1} t}-c_{1} {\mathrm e}^{4 c_{1} c_{2}} {\mathrm e}^{4 c_{1} t}}}{{\mathrm e}^{2 c_{1} c_{2}} {\mathrm e}^{2 c_{1} t}-1}, x \relax (t ) = -\frac {\sqrt {c_{1} {\mathrm e}^{2 c_{1} c_{2}} {\mathrm e}^{2 c_{1} t}-c_{1} {\mathrm e}^{4 c_{1} c_{2}} {\mathrm e}^{4 c_{1} t}}}{{\mathrm e}^{2 c_{1} c_{2}} {\mathrm e}^{2 c_{1} t}-1}\right \}, \{y \relax (t ) = 0\}, \left \{z \relax (t ) = \frac {\sqrt {-\left (\frac {d}{d t}x \relax (t )\right ) x \relax (t )}}{x \relax (t )}, z \relax (t ) = -\frac {\sqrt {-\left (\frac {d}{d t}x \relax (t )\right ) x \relax (t )}}{x \relax (t )}\right \}\right ], \left [\left \{x \relax (t ) = \RootOf \left (c_{3}+t -\left (\int _{}^{\textit {\_Z}}-\frac {2}{\sqrt {4 \textit {\_a}^{4}-16 c_{2} \textit {\_a}^{2}+16 c_{2}^{2}+c_{1}}\, \textit {\_a}}d \textit {\_a} \right )\right ), x \relax (t ) = \RootOf \left (c_{3}+t -\left (\int _{}^{\textit {\_Z}}\frac {2}{\sqrt {4 \textit {\_a}^{4}-16 c_{2} \textit {\_a}^{2}+16 c_{2}^{2}+c_{1}}\, \textit {\_a}}d \textit {\_a} \right )\right )\right \}, \left \{y \relax (t ) = -\frac {\sqrt {-2 \left (x \relax (t )^{3}-\frac {d}{d t}x \relax (t )-\sqrt {x \relax (t )^{6}-\left (\frac {d^{2}}{d t^{2}}x \relax (t )\right ) x \relax (t )+2 \left (\frac {d}{d t}x \relax (t )\right )^{2}}\right ) x \relax (t )}}{2 x \relax (t )}, y \relax (t ) = \frac {\sqrt {-2 \left (x \relax (t )^{3}-\frac {d}{d t}x \relax (t )-\sqrt {x \relax (t )^{6}-\left (\frac {d^{2}}{d t^{2}}x \relax (t )\right ) x \relax (t )+2 \left (\frac {d}{d t}x \relax (t )\right )^{2}}\right ) x \relax (t )}}{2 x \relax (t )}, y \relax (t ) = -\frac {\sqrt {-2 \left (x \relax (t )^{3}-\frac {d}{d t}x \relax (t )+\sqrt {x \relax (t )^{6}-\left (\frac {d^{2}}{d t^{2}}x \relax (t )\right ) x \relax (t )+2 \left (\frac {d}{d t}x \relax (t )\right )^{2}}\right ) x \relax (t )}}{2 x \relax (t )}, y \relax (t ) = \frac {\sqrt {-2 \left (x \relax (t )^{3}-\frac {d}{d t}x \relax (t )+\sqrt {x \relax (t )^{6}-\left (\frac {d^{2}}{d t^{2}}x \relax (t )\right ) x \relax (t )+2 \left (\frac {d}{d t}x \relax (t )\right )^{2}}\right ) x \relax (t )}}{2 x \relax (t )}\right \}, \left \{z \relax (t ) = \frac {\sqrt {\left (x \relax (t ) y \relax (t )^{2}-\frac {d}{d t}x \relax (t )\right ) x \relax (t )}}{x \relax (t )}, z \relax (t ) = -\frac {\sqrt {\left (x \relax (t ) y \relax (t )^{2}-\frac {d}{d t}x \relax (t )\right ) x \relax (t )}}{x \relax (t )}\right \}\right ]\right \}\]