Discussion :

[feynman]

Bonjour,
J ai devant moi a calculer la fonction de Bessel J0 pour un X donner et d arreter le programme quand le nouveau terme de la somme soit inferieur a 0.00001*(le terme totale),
Mais quand je la donne une valeur X , la reponse est NAN, alors ou est le truc? merci
voila le code

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
  program bessel
             real X,J
             print*,'donner X'
             read*,X
             CALL BESSELO(X,J)
             print*,'J0(X)=',J
             end
 
 
 
             SUBROUTINE BESSELO(Y,K)
             real Y,K
             K=0
             do i=1,10000
               K=K+((-1)**i)*((Y/2)**(2*i))/(fact(i))**2
               if (abs(K).LT.K) exit
             enddo
             end 
 
           function fact(m)
         real fact
       integer  m, p
       p = 1
              do i = 1, m
                   p = p * i
              end do
       fact = p
       end

Je rappel que J0(X)=somme(i=1,infini) (-1)^i/(i!)^2*(x/2)^(2*i)

[FR119492]

Salut.

Premièrement, si tu veux qu'on réponde à tes questions, il faudrait qu'on les comprenne. Alors, essaie de les écrire en français.

Ensuite, j'ai un peu mieux compris ton Fortran. Et je pense que ta formule est inutilisable pour la raison suivante:

[FR119492]

Désolé, mais je n'ai pas pressé la bonne touche!

Suite:
Je n'ai pas calculé combien vaut 10000!, mais je suis convaincu que c'est plus que le plus grand nombre représentable dans le processeur avec lequel tu travailles.

Jean-Marc Blanc

[feynman]

Je clarifie
tu sais que la fonction de Bessel est donnee par une serie .
le probleme dit : etant donne' un reel X, calculer J_0(X),
tu as J_0(X)= \sum_{m=0}^{\infty} U_m, ou U_m est une suite dont j ai donne' son terme general. il faut arreter le programme quand la difference entre \sum_{m=0}^{N} U_{m}-\sum_{m=0}^{N+1} U_{m} soit plus petite que 0.00001
Voila le probleme detailler
moi j ai fait le programme, mais il ne donne pas les reponses voulu

[afrancisco]

salut,

K=K+((-1)**i)*((Y/2)**(2*i))/(fact(i))**2

FR119492 a raison, ton proc verra rouge lorsqu'il devra évaluer fact(i) avec i se rapprochant de ... 100, alors imagine 10 000...

Juste par curiosité, que fais-tu comme études (travail) ?

[FR119492]

Tu utilises le développement en série de Taylor ou de MacLaurin qui est probablement une des plus mauvaises méthodes pour évaluer une fonction transcendante ailleurs que tout près de l'origine. As-tu choisi toi-même cette méthode ou t'a-t-elle été imposée par ton prof. ? Si c'est le cas, change de prof. au plus vite !

Je n'ai pas ma documentation sous la main, mais je vais essayer ces prochains jours de te trouver quelque-chose de mieux.

Jean-Marc Blanc

[feynman]

Pour afrancisco: mathématiques appliquées
pour Jean Marc: C est un exercice proposé dans un livre

[feynman]

J ai modifié le programme, mais quand meme je ne suis pas convaincu des résultats
voila mon code

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 bessel
             real X,J
             print*,'donner X'
             read*,X
             CALL BESSELO(X,J)
             print*,'J0(X)=',J
             end
 
 
             SUBROUTINE BESSELO(Y,K)
             real Y,K,Kold
             Kold=0
             do i=1,10000
               K=Kold+((-1)**i)*((Y/2)**(2*i))/(fact(i))**2
               Kold=K
               if (abs(K-Kold).LT.0.0001) exit
             enddo
             end 
 
             function fact(m)
              real fact
              integer  m, p
              p = 1
              do i = 1, m
                   p = p * i
              end do
               fact = p
             end

[Sylvain Bergeron]

Es-tu convaincu que les lignes suivantes sont dans le bon ordre ?

Code :

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

1
2
3
4
5
6
           .
...
               K=Kold+((-1)**i)*((Y/2)**(2*i))/(fact(i))**2
               Kold=K
               if (abs(K-Kold).LT.0.0001) exit
...

[feynman]

Bah oui, que penses tu?

[Sylvain Bergeron]

Comme Kold = K, comment veux-tu que abs(K - KOld) donne autre chose que 0 ?

[feynman]

t as raison j ai pas fait attention a ca,
alors quel est l equivalent? qu est ce qu il faut modifier?

[FR119492]

Salut !

Je t'avais suggéré de changer de prof. Je rectifie: change de livre.

D'autre part, à la venir, si tu veux qu'on t'aide efficacement, ce qu'on fait volontiers, copie textuellement dans tes messages l'énoncé complet de tes problèmes. Ce sera plus rapide que si on doit te demander trois fois des précisions.

Jean-Marc Blanc

[afrancisco]

salut,

D'une manière générale, lorque l'on calcul des quotients type coefficient du binôme, il faut arranger l'expression à implémenter :
si q=fact(i+1)/fact(i), il faut le remplacer lors de la programmation par q=i+1

Je sais, l'exemple est trivial, mais le principe s'applique à beaucoup de cas ... même si, parfois, il faut regrouper en plus des termes de la série ( notamment en présence d' un (-1)**i )

Il est utopique de vouloir systématiquement coder une formule telle quelle.

[DaZumba]

D'autre part, à la venir, si tu veux qu'on t'aide efficacement, ce qu'on fait volontiers, copie textuellement dans tes messages l'énoncé complet de tes problèmes. Ce sera plus rapide que si on doit te demander trois fois des précisions.

Certes, mais on n'est pas la pour faire ses devoirs non plus. Il faut qu'il apprenne par lui-meme. Pour l'instant, feynman lance un petit programme ecrit un peu au hasard sur le forum deux fois par jours, et attend la reponse. Ce n'est pas une bonne methode d'apprentissage.

[feynman]

Salut !

Je t'avais suggéré de changer de prof. Je rectifie: change de livre.

D'autre part, à la venir, si tu veux qu'on t'aide efficacement, ce qu'on fait volontiers, copie textuellement dans tes messages l'énoncé complet de tes problèmes. Ce sera plus rapide que si on doit te demander trois fois des précisions.

Jean-Marc Blanc

cher Jean Marc, j apprend le fortran seul, y a pas de prof.
voila l ennonce' en anglais:
write a program to approximate the value of the infinite series J_0(X)=1-X^2/2^2+X^4/(2^2*4^2)-X^6/(2^2*4^2*6^2)+..........
when computing the Bessel function , terminate the series when the absolute value of any new term changes the approximation by no more than 0.1%


J espere que je suis claire maintenant

[afrancisco]

Salut,

Malgré les difficultés que tu peux rencontrer (apprentissage du fortran et de l'anglais), tu dois réfléchir avant d'écrire ton programme : ce n'est pas parcequ'il t'est demandé de programmer une série qu'il faut recopier aveuglément la formule en fortran (ou autre ...)
Tu dois d'abord écrire l'algorithme sur papier et ensuite l'implémenter. On a l'impression que tu ponds un programme un peu vite, et que s'il ne fonctionne pas du premier coup, tu postes pour que l'on te le debugg.
Bref, "aide-toi et le ciel t'aideras" : la communauté est prête à t'aider, mais le premier niveau de debugg doit être assuré par toi-même.

Bon courage !