Forum de mathématiques - Bibm@th.net

Bonjour,

Suite à la réponse de Mathrack, sans être un expert en Fortran, je donnerai les indications suivantes, valables quelque soit le langage de programmation.

1) D'une manière générale dans un programme de "calcul numérique", je pense par expérience qu' il vaut mieux pré calculer tout ce qu'on peut pré calculer. Par exemple les constantes :

INTEGER N
DOUBLE PRECISION PI,TWOPI,RTWOPI,SIGMA, DOUBLESIGMACARRE, DX
PARAMETER(PI=3.141592653589793D0,TWOPI=2D0*PI,RTWOPI=SQRT(TWOPI),SIGMA=.15D0, DOUBLESIGMACARRE=2D0*SIGMA*SIGMA, N=100, DX=1D0/(2D0*N+1)

2) Il vaut mieux utiliser les expressions algébriques simples et quand on peut éviter les puissances. Simplifier les expressions et éviter les divisions en virgule flottante à répétition c'est moins coûteux et plus précis suivant les machines:

abcs = (i-1) * DX - l/2.D0
U(i) = 1/(exp( (abcs*abcs) /DOUBLESIGMACARRE ) * (RTWOPI * SIGMA ))

Je ne garantis pas que ma syntaxe soit la bonne mais en gros  l'idée c'est de laisser le moins possible de travail  à l'optimiseur.

3) Si l'édition du graphe proposé par Mathrack ne fonctionne pas, je testerai avec quelques valeurs de N telle que 1, 5, 10.

Cordialement
Denis