Bonjour, voila je suis en train de faire un programme sous fortran 90 et je dois calculer la somme suivante :

Som=((-1)**n/n**2)*exp((-(4*n**2*(pi**2)*alphac*deltato))/t0**2)*cos(n*pi)

avec :
  • alphac= 2.4185849e-05
  • deltato=4.2531456e-02
  • t0=1e-4
Voici le code que j’ai inclus dans mon programme :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
 
        som = 0.d0
	   do 30 n =1,1000
 
 som = som + ((-1.d0)**n/n**2.d0)*exp((-4.d0*n**2.d0*(pi**2.d0)* alphac*deltato)/t0**2.d0)*cos(n*pi)
 
30	continue   
write(*,*)”la somme est”,som
Le problème que je trouve c’est que la somme est égale toujours à 0, j’ai trouver la source de ce problème, j’ai changé la valeur de t0, en prenant par exemple t0=0.1
et j’ai vu que ça marche la somme égale à 1.53.
Alors je vous demande si il y a une solution à ce problème, je pense que c’est à cause de la très faible valeur pour la EXP.

Merci d’avance