Problème de subroutine dans une boucle

[invite17d39b9f]

Bonjour à tous,
J'ai un problème de dimensionnement de réacteur en série à réaliser sur fortran (avec utilisation de la méthode de newton-raphson obligatoire). Lorsque j'exécute le programme il me sort un premier résultat qui est bon (j'ai vérifié) puis il rentre dans une boucle ou il appelle une subroutine. Lors de la première boucle pour i=1 il me sort le résultat attendu cependant pour i=2 mon programme plante et je ne comprend pas pourquoi. (mon critère d'arrêt sur la méthode newton raphson déconne).
La subroutine qui plante se trouve à la ligne 26 et le critère d'arrêt de newton raphson est ligne 166.

si quelqu'un souhaite m'aider je lui ai mis mon programme et le fichier que mon programme doit lire en pièce jointe normalement.
-----

[umfred]

j'ai du mal à voir l'intérêt du DO i=1,1 en ligne 165 (mais bon pourquoi pas) au lieu d'affecter simplement i à 1 ou remplacer i par 1 dans la formule.
j'afficherai la valeur de crit avant le test ligne 168 pour voir comment il évolue (peut-être une erreur de signe quelque part?)

[invite17d39b9f]

Merci de votre réponse
les boucles allant de i=1 à 1 c'est juste que j'ai repris un ancien code qui utilisait aussi newton-raphson mais ça ne sert à rien ici. Je viens de regarder les valeurs du crit via le débogueur de fortran 90. Donc j'ai vu que c'était xvar qui était mauvais car je vois que xvar=-nan alors que à la ligne 149 je lui fixe une valeur initiale donc je ne comprends toujours pas pourquoi j'ai ce problème.

[umfred]

tu lui rajoutes deltax ensuite ligne 162; donc vérifier sa valeur aussi

[A voir en vidéo sur Futura]

[invite17d39b9f]

Bon j'ai cherché l'erreur pendant un certains temps et j'ai enfin trouvé. Ligne 204 j'approxime une dérivé par cette formule (F(x+h)-F(x-h))/2*h. or lorsque je calcule F(x+h) et F(x-h), je trouve les mêmes valeurs car fortran ne possède pas assez de chiffre après la virgule pour pouvoir les différencier ce qui donne une dérivé nulle et lorsque je l'inverse pour procéder à la méthode de newton raphson cela donne une valeur infinie.

Merci de m'avoir aidé à trouver le problème