#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Sat Oct  7 15:58:46 2023

@author: fabrice
"""


"""particule chargée"""

###Entête###

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

ax = plt.figure().add_subplot(projection='3d')# pour pourvoir tracer en 3D



###Vecteur temps####

t=np.linspace(1,50,1000)


"""Premier cas"""

##Fonction de Cauchy  (à compléter)

def F(X,t):
     , , , , , = 
    B=  ##champ magnétique éventuellement variable dans l'espace
    return np.array([, , , ,, ])


## Conditions initiales (à compléter)
X0=[ , , , , , ]


##Récupération du calcul (à compléter)
Sol=
X=
Y=
Z=

##Tracé en 3D (C'est fonctionnel mais vous pouvez l'améliorer)

ax.plot(X, Y, Z, label='Courbe')  # Tracé de la courbe 3D
plt.title("Trajectoire d'un proton dans un champ B selon z")
plt.xlabel('x')
plt.axis('auto')
plt.ylabel('y')
plt.legend('trajectoire')
plt.show()

"""Second cas"""

##...

"""Troisième cas"""



##...
