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

@author: caussin-jb
"""

import numpy as np
from scipy.integrate import odeint
import matplotlib.pyplot as plt


# Paramètres
g = 9.81  # Accélération de la pensanteur (m/s**2)
OmegaT = 2 * np.pi / 86400  # Vitesse angulaire de la Terre (rad/s)
L = 67   # Hauteur initiale
latitude = np.pi*45.85/180   # Latitude

w0 =      # A COMPLETER
f =       # A COMPLETER

# Système à résoudre
def F(U,t):
    x,y,vx,vy = U
    ax =      # A COMPLETER
    ay =      # A COMPLETER
    return    # A COMPLETER


U0 = [1,0,0,0]    # Condition initiale
    
T = np.linspace(0,6*3600,10000)
R = odeint(F,U0,T)

plt.figure()
plt.plot(      )     # A COMPLETER
plt.xlabel('x (m)')
plt.ylabel('y (m)')
plt.title('Trajectoire du pendule')
plt.tight_layout()
plt.show()