function X = DynSys(x0, T, h)
f = @(x) [x(1)*(1-x(2)); x(2)*(x(1)-1)];
N = T/h;
X = zeros(length(x0),1+N);
X(:,1) = x0;
for n=1:N
    X(:,n+1) =  X(:,n) + h*f(X(:,n));
end