How can I simulate a discrete- time dynamic system in command window?

Hi everyone, I want to simulate a discrete system, perhaps there is something similar to the continuous case:

[t,y]=ode45(@(t,x)func(x),t0tf,x0);

plot(y(:,1),y(:,2))

How could I get the 'discrete version' of the above expressions?.

to avoid problems let's say

func(x)=[x(1)-x(2);x(1)+x(1)*x(2)+x(2)];

to begin with I have no idea how to discretize a nonlinear system in state-space.

If not possible, what would be the right way to do it?