bonjour à tous j'ai un projet a faire sur le paradoxe des anniversaires sur scilab et la question est : A l’aide d’un programme Scilab déterminer le rang no à partir duquel pn > 1/2.
Je suis rendu la dans le programme :
i=0
n=0
x1=prod(1:n)/prod(1:365-n)
x2=1/365^n
x3=1-(x1*x2)
for i = 1:n
while x3 < 1/2
n=n+1
i=i+1
end
end
disp(i);

Et je dois avoir un résultat egal à 23 mais toute les valeurs que j'ai trouvée de correspondent pas.
Merci de votre réponse

Imaginons qu'on arrive a cette partie du code:

for i = 1:n
  while x3 < 1/2
    n=n+1
    i=i+1
  end
end

n=0, donc on ne rentre meme pas dans la boucle for. Ensuite, si on y rentrait, on a un while(condition sur x3) dans lequel on ne change jamais x3, donc si au depart on n'a pas la condition x3 < 1/2 on se retrouve avec une boucle infinie.

Ah et si prod(1:365-n) represente factoriel de 365-n, avec au depart n=0, tu demandes a Scilab de calculer 365!. C'est un nombre qui vaut a peu pres 10^800, ca va jamais marcher

La proba d'avoir au moins 2 personnes nées le même jour c'est 1 - prod(365-n+1:365)/365^n.

Ton i et ta structure for ne servent à rien.