I have a system of nonlinear differential equations with initial conditions: sys := diff(x(t), t) = -4*y(t)*x(t), diff(y(t), t) = -y(t)+2*y(t)^2-4*z(t)*y(t), diff(z(t), t) = -z(t)-2*u(t)*z(t), diff(u(t), t) = -u(t)+2*y(t)*u(t)-4*u(t)*z(t)+2*u(t)^2 ics := x(0) = 1, y(0) = 1, z(0) = 1, u(0) = 1 If I've got the equations right, I should get x(t) = 0.9019 as t -> infinity. How would I solve this system of equations on maple to get an answer for x(t)? I've tried this: dsolve({sys,ics}, numeric) but what I get is this: