# C'est ici qu'on calcule les valeurs du diagramme de bifurcation
x = []
y = []
while r < rmax :
    for n in range (nmin ,nmax) :
        x.append(float(r))
        y.append(float(u(n,r,u0)))
    r=r+dr