Bonjour à tous!

Nous sommes en train de travailler sur un projet scilab : attracteur de lorenz.
On doit résoudre le système d'equa diff avec la méthode RK4.
Voici le système:  dx/dt = 10(y-x)
                   dy/dt = -xz+px-y
                   dz/dt = xy-(8/3)z
On fixe p=28 ; et les conditions initiales sont:
x(0)=-3  ;  y(0)=-6   ;   z(0)=12 . On travaille sur l'intervalle [0;20].

On a trouvé ce programme qui correspond à ce qu'on veut:
http://www.umpa.ens-lyon.fr/~agregmath/Matlab/rklorenzae.m

Cependant, il est écrit pour Matlab, et on ne sait pas le transcrire en Scilab...
Quelqu'un peut nous aider?

Merci beaucoup!