#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Mon Apr 28 15:59:10 2025

@author: caussin-jb
"""

import numpy as np
import matplotlib.pyplot as plt


g = 9.81  # Accélération de la pensanteur (m/s**2)
OmegaT = 2 * np.pi / 86400  # Vitesse angulaire de la Terre (rad/s)
H = 158   # Hauteur initiale
latitude = np.pi*51/180   # Latitude
dt = 0.001  # Pas de temps numérique

# Conditions initiales
vx,vy,vz = 0,0,0
x,y,z = 0,0,H
t = 0


# Calcul de la solution

T = [t]
Vx,Vy,Vz = [vx],[vy],[vz]
X,Y,Z = [x],[y],[z]

while      # A COMPLETER
    t = t + dt
    T.append(t)
    
    # A COMPLETER
    vx = 
    vy = 
    vz = 
    
    x = 
    y = 
    z = 
    
    Vx.append(vx)
    Vy.append(vy)
    Vz.append(vz)
    X.append(x)
    Y.append(y)
    Z.append(z)

# Tracé des courbes

plt.figure()
plt.plot(     )      # A COMPLETER
plt.xlabel('t (s)')
plt.ylabel('y (m)')
plt.title('Déviation horizontale')
plt.show()

plt.figure()
plt.plot(     )      # A COMPLETER
plt.xlabel('y (m)')
plt.ylabel('z (m)')
plt.title('Trajectoire')
plt.show()