Bonjour,
Après avoir fait le programme pour calculer une intégrale sur un domaine ]a,b[ par la méthode de Simpson, je me casse la tête pour savoir où est l'erreur.
Pouvez-vous alors m'aider ?
Merci
Voilà mon code :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
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
 
          program simpsone
          real inte
          CALL simpson(0.0,1.0,100,F,inte)
          print*,'integrale entre 0 et 1 de F =',inte
          end
 
          real function F(X)
          real X
          F = X
          return
          end
 
 
 
          subroutine simpson(A,B,m,F,REP)
          real A,B,REP,S,V
          integer m,j,i
          S=0.0
          V=0.0
          do i=1,m
             S=S+F(A+(2*i-1)*(B-A)/(2*m))
          enddo
          do j=1,m
             V=V+F(A+j*(B-A)/m)
          enddo
          REP=((B-A)/(3*m))*(F(A)+F(B)+2*S+V)
          end