Ça peut sembler évident, mais mieux vaut le répéter : les formules "doivent toutes avoir pour caractéristique fondamentale de retrouver DU(t) = c*(M/N)(t) pour « N stable »."
Cela me conforte dans l’idée que DUA(t+dt) = max[DUA(t) ; c*M(t+dt)/N(t+dt)]
Peut-être qu’il y a un problème d’homogénéité dans les conventions, un truc du genre :
- pour DUA, DU(t+dt) sert à calculer M(t+dt)
- pour DUĞ, DU(t) sert à calculer M(t+dt)
En tout cas, à ceux qui veulent se lancer dans des simulations : ces histoires d’indices sont très piégeux…