Discussion :

[lucastof]

bonjour amis développeurs !!
voila je veut implémenter l'équation suivante en matlab mais le résultat et tjr NAN quand je calcule P(1,449) quand le N est grand quoi !!
alors que le calcule a la main me donne un résultats
je crois que je me suis trompé au niveau de l'algo.

Code :

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

1
2
3
function [p] = P(n,N)
p=factorial(N+n)/(((2*n)+1)*factorial((N-n)-1));
return;

[Jerome Briot]

Code :

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

1
2
3
4
5
6
7
8
9
10
11
>> factorial(449)
 
ans =
 
   Inf
 
>> Inf/Inf
 
ans =
 
   NaN

[lucastof]

svp comment faire pour que factorial(449) me retourne une valeur et non pas "INF" ???

[prgasp77]

Tu peux écrire à l'adresse suivante

MathWorks participates in the Better Business Bureau’s Safe Harbor dispute resolution process. European users may submit privacy complaints to MathWorks, or may file a complaint with the Better Business Bureau.

The Better Business Bureau may be contacted at:
Better Business Bureau EU Safe Harbor Program
4200 Wilson Blvd.
Arlington, VA 22203

[pseudocode]

Tu peux aussi écrire P(n,N+1) en fonction de P(n,N) et en déduire une relation de récurrence... Ou simplement factoriser le numérateur par (N-n-1)!.

P(1,449) = 30172800

[lucastof]

j'ai utilisé la foction VPI qui se trouve a ce lien
http://www.mathworks.com/matlabcentr...ger-arithmetic

donc problème résolu !!

[pseudocode]

j'ai utilisé la foction VPI qui se trouve a ce lien
http://www.mathworks.com/matlabcentr...ger-arithmetic

donc problème résolu !!

Ca me rappelle une histoire de marteau pilon et de mouche...

[lucastof]

pourquoi ??

[pseudocode]

pourquoi ??

C'est juste que ca me parait disproportionné de calculer des factorielles à 500 chiffres, alors qu'il suffit de simplifier les termes dans la fonction.

Code :

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

1
2
3
4
5
6
(N+n)! = 1 * 2 * 3 * ................................... * N+n
 
       = 1 * 2 * 3 * ... * (N-n-1) * (N-n) * (N-n+1) *...* (N+n)
         \_______________________/
 
       =          (N-n-1)!         * (N-n) * (N-n+1) *...* (N+n)

D'où

Code :

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

1
2
3
4
 
          (N-n)*...*(N+n)
 P(n,N) = ---------------
              (2n+1)

En particulier

P(1,449) = (448*449*450)/3 = 30172800

[lucastof]

oui c'est vrai !!
que va être la forme générale de ma fonction alors ??
et une petite question :
comment gérer 360! /450! en matlab ??

[pseudocode]

oui c'est vrai !!
que va être la forme générale de ma fonction alors ??

Après il faut demander a des spécialistes de Matlab comment cela se code.

[Nebulix]

comment gérer 360! /450! en matlab ??

Si pour ton problème cette quantité est significativement différente de ZERO, explique nous le, il doit être passionnant !

[nihao]

Bonjour,

Avec SciLab, il y a la fonction prod(x) et tu peux écrire :

Code :

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

1
2
 
prod(N-n:N+n)/(2*n+1)

Je ne connais pas en Mathlab
Remarque :

Tu va retomber sur le même problème si n devient grand.

170! est représentable en double précision IEEE mais 171! provoque un débordement (Numerical Recipes 3th edition p.257)

Si le contexte dans lequel tu utilise P(n,N) le permet, passe en log

Cordialement