Implémentation méthode de pascual-meeker

Bonjour,

Je travail actuellement sur l'estimation d'une limite d'endurance sur des données de nickel. Pour cela j'étudie l'article de pascual et Meeker intitulé "Estimating Fatigue curves with the random Fatigue-Limit Model".
J'ai implémenté leur modèle mais celui-ci a soit des problèmes au niveau des intégrales soit au niveau des méthodes d'optimisation.
J'ai essayé beaucoup de choses mais rien ne marche. J'aimerais donc un regard nouveau sur mon code pour trouver ses erreurs.
Merci d'avance;

Florine