Résolution d'ODE avec scilab

[invite25e4f388]

Salut,
J'éprouve quelques difficultés avec un exercice que je dois résoudre. Il s'agit de modéliser un globe terrestre en sustentation magnétique avec Scilab. Nous avons un aimant qui attire vers le haut un globe terrestre qui est aussi soumis à la force de gravitation.
On nous donne les équations qui décrivent le système et nous devons utiliser ode pour les résoudre.



J'ai donc :

m=0.1;
lambda=0.0001;
e0=0.01;
R=0.1;
g=9.81;

i0= sqrt((2*g*2*m*e0^2)/lambda);
u0=R*i0;

function [xdot]=f(t,x)
  xdot(1)=x(2);
  xdot(2)=(lambda*x(3)^2)/(2*m*(e0-x(1))^2)-g;
  xdot(3)=-(x(3)*x(2))/(e0-x(1))+(u0-R*x(3))*(e0-x(1))/lambda;
endfunction

x0=[0.001; 0; i0];
  
t=0:0.001:10;
solution=ode(x0,0,t,f);
plot(t,solution(1));

Le problème est que la variable solution censée contenir toutes les valeurs de x[i] sur l'intervalle de temps t n'est qu'en fait qu'un vecteur ligne à 3 colonnes.
Je ne comprends pas le problème. Y aurait-t-il quelqu'un pour m'aider??
L'url du sujet : http://tice.utc.fr/moodle/mod/resour...w.php?id=13986. Je bloque dans la deuxième question

Image passée en pièce jointe : pas d'image sur un serveur extérieur.
JPL, modérateur
-----

[invite4b11d5c9]

Vive NF03 hein :P.