intégrale sur fortran [Résolu]

[suzanne1307]

Bonsoir,
j'ai programmé un programme sur fortran qui calcule une intégrale simple
il ne génère aucune erreur mais il donne toujours comme résultat 0.
svp aidez moi à détecter l'erreur je bloque
voici le programme:

Code :

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
external GAUS,TRAP
	double precision y
	y=TRAP(1,10,GAUS,1)
	print *, y
	end

	double precision FUNCTION GAUS (T)
	GAUS=EXP(T)
	return
	end

	double precision FUNCTION TRAP (A,B,F,DX)
	! A et B sont les bornes de l'intégrale
	!F est la fonction à intégrer
	!DX est le pas d'intégration
	!N est le nombre de pas nécessaires
	implicit none
	integer N,I
	double precision B,A,DX,F,X
	N=(B-A)/DX-.5
	TRAP=F(A)+F(B)/2.0
	X=A
	DO 10 I=1,N
	X=X+DX
	TRAP=TRAP+F(X)
10    CONTINUE 
      TRAP=TRAP*DX
	RETURN
	end

merci d'avance

[genteur slayer]

as-tu essayer de remplacer F par GAUS dans la seconde fonction (en mettant GAUS en external)

[Sylvain Bergeron]

On ne le répète jamais assez souvent :

implicit none, implicit none, implicit non...

La déclaration implicite des variables est une erreur de design de Fortran. On doit toujours débuter un programme ou sous-programme par un implicit none. Le temps de débogage sauvé en moyenne est nettement supérieur au temps requis pour taper cette instruction, même dans un petit programme.

Ton programme ne fonctionne pas car rien n'indique dans le programme principal que TRAP et GAUS sont double precision. Par défaut, ils sont donc considérés real simple précision.

[suzanne1307]

Bonsoir,
j'ai remplacé F dans la fonction TRAP, j'ai mis implicit none au début duprogramme et j'ai déclaré TRAP et GAUS comme double precision comme vous me l'aviez dit mais sans résultat, aucun changement
svp aidez moi

[mipicard]

Chez moi cela marche si tu appelles la fonctionde la façon suivante

Code :

1
2
3
4
5
6
       REAL*8 A,B,DX

         A= 1.0
         B= 10.0
         DX = 1.0
         y=TRAP(A,B,GAUS,DX)

l'appel simple y=TRAP(1,10,GAUS,1) ne respecte pas les déclarations de la fonction TRAP.

M.P.

[suzanne1307]

Je vous remercie beaucoup pour votre aide