Télécharger derip2.eso

Retour à la liste

Numérotation des lignes :

  1. C DERIP2 SOURCE CHAT 05/01/12 22:42:57 5004
  2. SUBROUTINE DERIP2(P, IVAR, DP)
  3. C
  4. C Procedure de multiplication de 2 polynomes de 2 variables
  5. C Entree : P(9) : polynome a deriver : les coefficients sont
  6. C ranges comme suit :
  7. C indice : 1 2 3 4 5 6 7 8 9
  8. C terme : 1 T2 T1 T1*T2 T2^2 T1^2 T1*T2^2 T1^2*T2 T1^2*T2^2
  9. C IVAR : numero de la variable par rapport a laquelle il
  10. C faut deriver
  11. C Sortie : DP(9) : polynome resultat de la derivation
  12. C
  13. IMPLICIT INTEGER(I-N)
  14. IMPLICIT REAL*8(A-H,O-Z)
  15. DIMENSION P(9), DP(9), INXD1(9), INXD2(9), C1(9), C2(9)
  16. DATA INXD1/ 1, 1, 1, 2, 1, 3, 5, 4, 7/
  17. DATA C1/0., 0., 1., 1., 0., 2., 1., 2., 2./
  18. DATA INXD2/ 1, 1, 1, 3, 2, 1, 4, 6, 8/
  19. DATA C2/0., 1., 0., 1., 2., 0., 2., 1., 2./
  20. C
  21. C Initialisation du polynome resultat
  22. C
  23. DO 70 I = 1, 9
  24. 70 DP(I) = 0.D0
  25. C
  26. C Derivation de P
  27. C
  28. IF (IVAR .EQ. 1) THEN
  29. DO 80 I = 2, 9
  30. 80 DP(INXD1(I)) = DP(INXD1(I)) + C1(I)*P(I)
  31. ELSE
  32. IF (IVAR .EQ. 2) THEN
  33. DO 90 I = 2, 9
  34. 90 DP(INXD2(I)) = DP(INXD2(I)) + C2(I)*P(I)
  35.  
  36. ENDIF
  37. ENDIF
  38. RETURN
  39. END
  40.  
  41.  
  42.  

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