Télécharger invaj3.eso

Retour à la liste

Numérotation des lignes :

  1. C INVAJ3 SOURCE CHAT 05/01/13 00:42:12 5004
  2.  
  3. SUBROUTINE INVARI_J3 (SIG,N,J3)
  4. IMPLICIT INTEGER(I-N)
  5. integer n,i
  6. real*8 sig(n),J3,desv(6),desvdesv(6)
  7. call Desviador(sig,desv,n)
  8. desvdesv(1)=desv(1)*desv(1)
  9. desvdesv(2)=desv(2)*desv(2)
  10. desvdesv(3)=desv(3)*desv(3)
  11. if (n.ge.4) then
  12. desvdesv(1)=desvdesv(1)+desv(4)*desv(4)
  13. desvdesv(2)=desvdesv(2)+desv(4)*desv(4)
  14. desvdesv(4)=(desv(1)+desv(2))*desv(4)
  15. if (n.eq.6) then
  16. desvdesv(1)=desvdesv(1)+desv(5)*desv(5)
  17. desvdesv(2)=desvdesv(2)+desv(6)*desv(6)
  18. desvdesv(3)=desvdesv(3)+desv(5)*desv(5)+desv(6)*desv(6)
  19. desvdesv(4)=desvdesv(4)+desv(5)*desv(6)
  20. desvdesv(5)=(desv(1)+desv(3))*desv(5)+desv(4)*desv(6)
  21. desvdesv(6)=(desv(2)+desv(3))*desv(6)+desv(4)*desv(5)
  22. endif
  23. endif
  24. J3=0.D0
  25. do i=1,3
  26. J3=J3+desvdesv(i)*desv(i)
  27. enddo
  28. J3=J3/3.D0
  29. if (n.lt.4) return
  30. do i=4,n
  31. J3=J3+2.D0*desvdesv(i)*desv(i)/3.D0
  32. enddo
  33. return
  34. end
  35.  
  36.  
  37.  
  38.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales