dsolve({diff(x(t),t) = 1/2*x(t)^2-1/24*y(t), diff(y(t),t) = 2*x(t)*y(t)-3*z(t), diff(z(t),t) = 3*x(t)*z(t)-1/6*y(t)^2})