Discussion :

[feynman]

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 :

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
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

[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.

[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...

[genteur slayer]

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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

[feynman]

Bonsoir,
Merci pour vos reponses,
mais je ne comprends pas le role de return, pourquoi c est tres important de le mettre?

[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.