Scilab : Problème

Voir le sujet précédent Voir le sujet suivant Aller en bas

Scilab : Problème

Message par octo75 le Jeu 12 Avr - 15:58

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 Crying or Very sad

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.

octo75

Nombre de messages : 2
Age : 37
Date d'inscription : 11/04/2012

Revenir en haut Aller en bas

Re: Scilab : Problème

Message par domi le Lun 16 Avr - 16:01

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!

_________________
"Donne un poisson à un homme, tu le nourris pour un jour. Apprends-lui à pêcher, tu le nourris pour toujours." Lao-Tseu
avatar
domi
Admin

Nombre de messages : 1348
Age : 59
Localisation : Noyon
Date d'inscription : 22/10/2008

http://www.tangenteX.com

Revenir en haut Aller en bas

Voir le sujet précédent Voir le sujet suivant Revenir en haut

- Sujets similaires

 
Permission de ce forum:
Vous ne pouvez pas répondre aux sujets dans ce forum