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