## Des cyclotrons pour la médecine##

import numpy as np
import matplotlib.pyplot as plt

plt.close('all')

## Définitions des grandeurs du problème ##

m =      # Masse de la particule H- (kg)
q =            # Charge de la partcicule H-
B =           # Champ magéntique (T)
V0 =        # Potentiel max entre les Dees (V)
d =        # Distance entre les Dees (m)
fc =      # Fréquence cyclotron (Hz)
Rmax =       # Rayon maximal (m)

## Définition de l'accélération##

def acceleration (x,y,vx,vy,t):
    if np.abs(x) < d/2:
        E=
    else :
        E=
    
    
    return ax,ay

## Définition de la fonction cyclotron##

def cyclotron(dt, tmax, Rmax):
    x = []
    y = []
    vx= []
    vy= []
    t = []
    k = 0
    while (np.sqrt(x[k]**2+y[k]**2) < Rmax) and (t[k]< tmax):
        ax,ay= acceleration (x[k],y[k],vx[k],vy[k],t[k])
        ...
    return np.array(x), np.array(y), np.array(vx), np.array(vy), np.array(t)

## Caractéristiques de la simulation ##

tmax =2e-6       # Durée de la simulation (s)
dt = 1e-10       # pas de temps (s)

## Tracé de la trajectoire ##

x,y,vx,vy,t = 

plt.axis ('equal')   # Même échelle sur les deux axes
plt.plot(....)
plt.show()