[MatLab] Equazioni Differenziali..aiuto

[rdiepo]

Ciao a tutti...non capisco una cosa!!!
Ho un sistema di due equazioni:
dx/dt = y
dy/dt= -x

Si tratta quindi di un sistema lineare della forma;
dz/dt = A*z

dove A=[0 1;-1 0]
dove gli autovalori sono pertanto +j e -j.

Il problema è questo:
Se utilizzo la funzione ode45 e faccio il grafico ottengo un andamento oscillante ma stabile delle due variabili (come dovrebbe essere)
MA se utilizzo il seguente codice(che è quello che devo usare):

ft=30;

dt=0.1;

T=0:dt:ft;
X0=[1 1]';
X=X0;
XArray=X
A=[0 1;-1 0]

for k=dt:dt:ft
X=X+(A*X)*dt;
XArray=[XArray X];
end

plot(T,XArray)

ottengo un andamento delle variabili che oscillano ma instabili (l'ampiezza aumenta all'infinito)
perchè?????? aiutooo