Scilab : Problème

Bonjour,

J'ai quelques soucis avec Scilab (j'suis débutant) et je ne vois pas où se situe le soucis ... J'ai pourtant calqué des exemples dans l'aide, en vain ... J'obtiens juste un graphique avec "Soleil" au centre, sans courbe

Code:

function [f]=LoiKepler1(t, u)
f(1)= du(1)
f(2) du(2)
df(1)=-u(1)/(u(1)^2+u(2)^2)^(3/2)
df(2)=-u(2)/(u(1)^2+u(2)^2)^(3/2)
endfunction

PositionInitiale1=[1;0];
VitesseInitiale1=[0;0.5];
init1=[PositionInitiale1;VitesseInitiale1];


T=24*3600*365;temps=0:1000:T;

sol1=ode(init1,0,temps,LoiKepler1);

// Plot de la solution
plot(0,0,'X')
xstring(0,0,'Soleil')
plot(sol1(1,:),sol1(2,:),'g')

Merci d'avance.

Bonjour,
Plusieurs choses m'étonnent dans ton code:
je n'ai pas connaissance que les lois de Kepler impliquent la résolution d'une équation différentielle.
je vois apparaître une vitesse initiale mais aucune grandeur qui ressemble à une accélération, ce qui est curieux lorsqu'il s'agit d'obtenir une vitesse par intégration.

Pour tout dire, ce code ressemble étrangement (sic!) à un code déjà vu sur le net mais mal recopié, et pas forcément compris...

Alors commençons par le commencement : quelles équations différentielles veux-tu utiliser pour simuler la trajectoire de la Terre autour du Soleil? Pour répondre à une de tes observations faites par ailleurs, l'excentricité de l'ellipse est telle qu'elle ressemble très fort à un cercle!