Télécharger int3p2.eso

Retour à la liste

Numérotation des lignes :

  1. C INT3P2 SOURCE CHAT 05/01/13 00:39:29 5004
  2. SUBROUTINE INT3P2(P, RES)
  3. C
  4. C Procedure d'integration de polynomes de 2 variables sur le
  5. C domaine triangle (0., 1.)x(0., 1.)
  6. C Entree : P(9) : polynome a integrer : les coefficients sont
  7. C ranges comme suit :
  8. C indice : 1 2 3 4 5 6 7 8 9
  9. C terme : 1 T2 T1 T1*T2 T2^2 T1^2 T1*T2^2 T1^2*T2 T1^2*T2^2
  10. C Sortie : RES : REAL*8 : resultat de l'integration
  11. C
  12. IMPLICIT INTEGER(I-N)
  13. IMPLICIT REAL*8(A-H,O-Z)
  14. DIMENSION P(9), C1(9)
  15. DATA C1/.5D0, 0.166666666666667D0, 0.166666666666667D0,
  16. 1 0.208333333333333D0, 0.083333333333333D0,
  17. 2 0.083333333333333D0, 0.016666666666667D0,
  18. 3 0.016666666666667D0, 0.005555555555556D0/
  19. RES = 0.D0
  20. DO 30 I = 1, 9
  21. 30 RES = RES + C1(I)*P(I)
  22. RETURN
  23. END
  24.  
  25.  
  26.  
  27.  

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