[Débutant] Recherche d'erreur dans un code

Version imprimable

Bonjour,
je suis en train d'écrire un programme pour calculer C_{n}^{p}
J'ai eu 2 erreurs mais je n'arrive pas à m'en débarasser, pouvez-vous m'aider ?

Voilà mon 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 program expo implicit none integer n,p real sum, COMBI,comb print* ,'entrer le nombre n:' read*,n print* ,'entrer le nombre p' read*,p COMBI=comb(n,p) print *, COMBI end function fact(m) real fact integer m, p p = 1 do i = 1, m p = p * i end do fact = p end function comb(n,p) real comp integer n,p real q q=fact(n)/(fact(p)*fact(n-p)) comp=q end
Je rappelle C_{n}^{p}=n!/(p!.(n-p)!
Merci pour votre aide

17/07/2007, 10h12
genteur slayer

quel est le message d'erreur?

pourquoi tu n'utilise pas des double precision (ou même quadruple precision) les factorielles sont rapidement très grosses.... essaye juste d'écrire factorielle 20 déjà et tu va rire!!!

je crois qu'il y a d'autres algo pour touver le nombre de combinaisons; et il ne faut pas passé par les factorielles mais par le triangle de pascal, en gros on définit les combi par récurence.
17/07/2007, 13h03
Sylvain Bergeron

Question d'insister : Quel est le message d'erreur ?

En plus de ce qu'a dit genteur slayer, le "implicit none" doit être répété à chaque programme, fonction et subroutine. Quand je l'ajoute après function fact(m) et après function comb(n,p), j'obtiens 3 erreurs et 1 avertissement...

un petit rajout:

dans la fonction fact, le "i" n'est jamais défini... on sait pas ce que c de plus, fact n'est pas un real, mais un integer (bon c pas trop grave), question d'économie de mémoire, pas besoin de la variable p:
Code:

1 2 3 4 5 6 7 8 9 10 11 integer(kind=8) function fact(m) !pour le kind, je te conseil le max que ton compilo supporte implicit none integer,intent(in) :: m !pour etre sûr de pas le modifié dans la fonction integer i; fact=1 do i=2,m !pas besoin de commencer à 1, à causs de la ligne precedente fact=facf*i end do return !trèèèèèèèès important, faut retourné la valeur de fact!!!! end function fact
ensuite, dans la fonction combi, elle doit elle aussi retourné une valeur (donc fo bien collé sont "return"
de plus tu as confondu comp et comb.... pas malin
enfin, le nombre de combi est un entier, pas un reel!!!

donc si tu veux le faire avec la formule des factorielles (pas conseillé du tout) cela donne:
Code:

1 2 3 4 5 6 7 integer function comb(n,p) implicit none integer,intent(in) :: n,p comb=fact(n)/(fact(p)*fact(n-p)) return end function comb
enfin, dans ton prog principal, ne déclare pas de real comb alors que tu as une fonction qui a déjà ce nom, ce n'est pas une faute de language, c'est parfaitement autorisé (tout comme appeler un reel "if" par exemple) mais c fortement déconseillé car tu sais plus qui est quoi....
Code:

1 2 3 4 5 6 7 8 9 10 11 program expo implicit none integer n,p combi print * ,'entrer le nombre n:' read *,n print * ,'entrer le nombre p' read *,p COMBI=comb(n,p) print *, COMBI end program expo

17/07/2007, 19h31
feynman

Bonsoir,
Merci pour vos reponses,
mais je ne comprends pas le role de return, pourquoi c est tres important de le mettre?
17/07/2007, 19h35
Sylvain Bergeron

Le return est facultatif dans la mesure où il y a implicitement un return sur le end. Le return est requis quand on veut sortir avant la fin d'une fonction ou d'une routine. Tu peux donc les omettre dans ton programme.

Bonsoir,
J ai essaye le programme , mais dans la compilation j ai 2 erreurs :
voila mon 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
program expo
  implicit none
  integer n,p, combi
    print * ,'entrer le nombre n:'
    read *,n
    print * ,'entrer le nombre p'
    read *,p
    COMBI=comb(n,p)
    print *, COMBI
end program expo
 
integer function fact(m) !pour le kind, je te conseil le max que ton compilo supporte
  implicit none
  integer,intent(in) :: m !pour etre sûr de pas le modifié dans la fonction
  integer i;
    fact=1
    do i=2,m   !pas besoin de commencer à 1, à causs de la ligne precedente
      fact=fact*i
    end do
  return      !trèèèèèèèès important, faut retourné la valeur de fact!!!!
end function fact
 
integer function comb(n,p)
 implicit none
  integer,intent(in) :: n,p
    comb=fact(n)/(fact(p)*fact(n-p))
  return
end function comb

COMB must appear in a type declaration because implicit none has been used.
alors comment faire?
merci

Dans les units utilisant des fonctions externes, tu dois déclarer ces fonctions externes. Tu dois donc ajouter des lignes du genre :

« Type », external :: NomDeLaFonction

Ton programme devient donc :

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
30
program expo
implicit none
integer n,p, combi
integer, external :: comb
print * ,'entrer le nombre n:'
read *,n
print * ,'entrer le nombre p'
read *,p
COMBI=comb(n,p)
print *, COMBI
end program expo
 
integer function fact(m) !pour le kind, je te conseil le max que ton compilo supporte
implicit none
integer,intent(in) :: m !pour etre sûr de pas le modifié dans la fonction
integer i;
fact=1
do i=2,m !pas besoin de commencer à 1, à causs de la ligne precedente
fact=fact*i
end do
return !trèèèèèèèès important, faut retourné la valeur de fact!!!!
end function fact
 
integer function comb(n,p)
implicit none
integer,intent(in) :: n,p
integer, external :: fact
comb=fact(n)/(fact(p)*fact(n-p))
return
end function comb

La méthode alternative (que j'utilise toujours) est d'utiliser un module :

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
30
31
32
33
34
 
module Factoriels
  implicit none
  private
  public fact,comb
contains
 
  integer function fact(m) 
    integer,intent(in) :: m
    integer i;
    fact=1
    do i=2,m !pas besoin de commencer à 1, à causs de la ligne precedente
      fact=fact*i
    end do
    end function fact
 
  integer function comb(n,p)
    integer,intent(in) :: n,p
    comb=fact(n)/(fact(p)*fact(n-p))
    end function comb
 
end module
 
program expo
use Factoriels
implicit none
integer n,p, combi
print * ,'entrer le nombre n:'
read *,n
print * ,'entrer le nombre p'
read *,p
COMBI=comb(n,p)
print *, COMBI
end program expo

Avec le module, il y a validation automatique des arguments, des types, etc. Le retour est très rapide pour les 3 minutes investies...