Différences

Ci-dessous, les différences entre deux révisions de la page.

--- python:programmes_de_physique_chimie_en_python:parabole [2019/10/20 15:10] – physix
+++ python:programmes_de_physique_chimie_en_python:parabole [2020/07/24 03:31] (Version actuelle) – modification externe 127.0.0.1
@@ Ligne 1: / Ligne 1: @@
 ====== Parabole ======
-<code>
+<code python>
-especes=["CH4","O2","CO2","H2O"]
+#Conditions initiales:
-coefficients=[-1, -2, 1, 2]
-etat_init = [3., 20., 0., 0.] #quantites initiales en moles
-etat_inter = etat_init[:]#on initialise l'état intermediaire
-dksi = 0.0001 #pas
-x = 0 #avancement initial
+x = 0.
-while etat_inter[0]>0 and etat_inter[1]>0: #tant qu'il y a des réactifs
-   x = x + dksi
-    etat_inter[0] = etat_init[0] + coefficients[0]*x
-    etat_inter[1] = etat_init[1] + coefficients[1]*x
-    etat_inter[2] = etat_init[2] + coefficients[2]*x
-    etat_inter[3] = etat_init[3] + coefficients[3]*x
-print("====== Etat Final ======")
+y = 2.
-for i in range(4):
-    print(especes[i],": {:.3f}".format(etat_inter[i])) #on affiche le résultat avec 3 chiffres après la virgule
+vx = 3.5
-print("")
-print("Avancement final : xf = {:.3f} mol".format(x))
+vy = 0.
+ax = 0.
+ay = -9.81
+t = 0.
+#pas de temps
+dt = 0.001
+#Simulation
+while y> 0. :
+x = x + vx * dt
+y = y + vy * dt
+vx = vx + ax * dt
+vy = vy + ay * dt
+t = t + dt
+print("La balle touche le sol à t =", t, "s et à la position x =", x, "m")
 </code>