# -*- coding: utf-8 -*-
"""
Created on Tue Apr  7 17:01:28 2026

@author: delph
"""

import numpy as np
import matplotlib.pyplot as plt

# Paramètres
m = 1.0  # masse
R = 1.0  # rayon
g = 9.81  # accélération gravitationnelle
a = 1.0   # accélération constante

# Fonction pour d(alpha)
def d(alpha):
    return # A compléter

# Fonction pour J
def J(m, R):
    return 0.5 * m * R**2

# Fonction pour E_c_max(alpha)
def E_c_max(alpha):
    d_alpha = d(alpha)
    J_alpha = J(m, R)
    return # A compléter

# Plage de valeurs pour alpha
alpha_values = np.linspace(0.1, np.pi, 1000)
E_c_values = [E_c_max(alpha) for alpha in alpha_values]

# Tracé
plt.plot(alpha_values, E_c_values)
plt.xlabel('Alpha (rad)')
plt.ylabel('E_c_max')
plt.title('Énergie cinétique maximale en fonction de alpha')
plt.grid()
plt.show()