Bonjour, je souhaite écrire un programme qui fait un calcul sur trois sommes avec des coefficients binomiaux. Voici ce que j'ai fait, mais j'ai un doute sur çà. Je vous serai très reconnaissant si vous pouvez m'aider. Merci.
Pour vous donner une idée de la somme que je veux calculer: Somme = ∑i ∑j ∑k [ bino(i)*bino(j)*bino(k)*exp(x*y*(2*i+j+k)) ]
integer i,j,k,n
real bino1,bino2,bino3,Somme
a=fact(n)
Somme=0.
do 30 i=0,5
bino=fact(5)/(fact(5-i)*fact(i))
do 30 j=0,8
bino2=fact(8)/(fact(8-j)*fact(j))
do 30 k=0,6
Somme=Somme + bino1*bino2*bino3*exp(x*y*(2*i+j+k))
30 continue
stop
end
function fact(nn)
fact=1
do l=1,nn
if(nn.eq.0.)then
fact=1
else
fact=fact*l
endif
enddo
end
Partager