// *********************************
// Etude de la chute d'un corps
// Tracé de la courbe z = f(t)
// Programme de mécanique de seconde
// Dominique LEFEBVRE
// 26 février 2007
// *********************************

clear;             // Initialisation de toutes les variables

// Initialisation des paramètres de simulation
deltaT = 0.01;     // base de temps = 0,1 seconde

// Initialsation des paramètres du mouvement
g  =  9.81;        // accélération de la pesanteur 9,81 N.m-2 (ou N.kg)
v0 =  2 ;          // vitesse initiale du corps (en m.s-1)
z0 =  10;          // ordonnée initiale du corps (en m)

// Initialisation des variables de calcul
z(1) = z0;         // ordonnée initiale = y0, paramètre de la simulation
vZ(1) = 0;         // vitesse verticale initiale du corps
t(1) = 0;
// Initialisation graphique
xbasc();

// Calcul du mouvement
i = 2;                          // indice de boucle 
while (z(i-1) > 0)              // la boucle est exécutée tant que le corps n'est pas sur le sol  
  vZ(i) = vZ(i-1) - g*deltaT;
  z(i) = z(i-1) + vZ(i)*deltaT; 
  t(i) = t(i-1) + deltaT;
  i = i + 1;                    // on passe à l'indice de tableau suivant
end

// Affichage du mouvement dans le repère choisi
plot2d(t,z,rect=[0,0,3,z0],style = 5);
xtitle('Chute d''un corps - Pas de temp: 0.01 ','Temps','Altitude');