Discussion :

[feynman]

Bonjour,
J'ai écrit un programme pour calculer l'angle entre deux vecteurs, par la formule Cos(alpha)=A.B/(|A|.|B|), mais la réponse c est toujours NAN,
alors où est l'erreur dans mon code ci-dessous ?
Merci

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
30
31
32
 
             program alphae
             real A(3),B(3)
             real alpha
             read*,(A(i),i=1,3)
             read*,(B(i),i=1,3)
             alpha=ACOS(DOTT(A,B,3)/(MODU(A,3)*MODU(B,3)))
             print*,'alpha=',alpha
             end
 
 
             function DOTT(H, P, L)
             integer L
             real H(L)
             real P(L),DOT 
             DOT = 0.0
             do 10 i  = 1,L
                DOT = DOT + H(i)*P(i)
10           continue
             return
             end  
 
             function MODU(V,N)
             integer N
             real V(N),Z
             Z=0.0
             do 20 i=1,N
                Z = Z+V(i)**2
20           continue
             Z=sqrt(Z) 
             return
             end

[phy4me]

real A(3),B(3)

Produit scalaire dans l'espace ??? les vecteurs A et B ne devront avoir que deux coordonnées je vois pas la raison de mettre A(3)

[feynman]

Oui, le produit scalaire peut etre dans dimension N>3
A.B= sum(i=1,N)A(i)*B(i)

[FR119492]

Comment veux-tu que ta fonction DOTT sache qu'elle doit prendre la valeur que tu as mise dans la variable locale DOT? Idem pour MODU et Z.

Jean-Marc Blanc

[feynman]

Tu as raison,
alors qu est ce qu il faut modifier?

[FR119492]

Ta fonction DOTT doit se terminer par une instruction
DOTT=DOT
et ta fonction MODU par
MODU=Z

JMB

[feynman]

Le alpha est en radian ou degre?

[afrancisco]

salut,

Les angles renvoyés par les fonctions réciproques trigo sont en radian. L'indéterminations liée à l'usage de acos ([0,pi]) doit être levée avec un asin et l'usage d'un produit vectoriel, si toutefois c'est un angle orienté que tu cherches ...

[feynman]

D accord
Merci