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