dsolve({diff(x(t),t) = x(t)*(y(t)^2-z(t)^2), diff(y(t),t) = -y(t)*(z(t)^2+x(t)^2), diff(z(t),t) = z(t)*(x(t)^2+y(t)^2)})