Télécharger int6p2.eso

Retour à la liste

Numérotation des lignes :

  1. C INT6P2 SOURCE CHAT 05/01/13 00:39:32 5004
  2. SUBROUTINE INT6P2(P, RES)
  3. C
  4. C Procedure d'integration de polynomes de 2 variables sur le
  5. C domaine (-1., 1.)x(-1., 1.)
  6. C Entree : P(30) : premier polynome : les coefficients sont
  7. C ranges comme suit :
  8. C indice : 1 2 3 4 5 6 7 8 9
  9. C terme : 1 T1 T2 T1*T2 T1^2 T2^2 T1*T2^2 T1^2*T2 T1^3
  10. C indice : 10 11 12 13 14 15
  11. C terme : T2^3 T1*T2^3 T1^2*T2^2 T1^3*T2 T1^4 T2^4
  12. C indice : 16 17 18 19 20 21
  13. C terme : T1*T2^4 T1^2*T2^3 T1^3*T2^2 T1^4*T2 T1^5 T2^5
  14. C indice : 22 23 24 25 26
  15. C terme : T1*T2^5 T1^2*T2^4 T1^3*T2^3 T1^4*T2^2 T1^5*T2
  16. C indice : 27 28 29 30
  17. C terme : T1^2*T2^5 T1^3*T2^4 T1^4*T2^3 T1^5*T2^2
  18. C Sortie : RES : REAL*8 : resultat de l'integration
  19. C
  20. IMPLICIT INTEGER(I-N)
  21. IMPLICIT REAL*8(A-H,O-Z)
  22. DIMENSION P(30), C1(30)
  23. DATA C1/0.5, 0.166666666666667D0, 0.166666666666667D0,
  24. 1 0.041666666666667D0, 0.083333333333333D0, 0.083333333333333D0,
  25. 2 0.016666666666667D0, 0.016666666666667D0, 0.05D0, 0.05D0,
  26. 3 0.008333333333333D0, 0.005555555555556D0, 0.008333333333333D0,
  27. 4 0.033333333333333D0, 0.033333333333333D0, 0.004761904761904D0,
  28. 5 0.002380952380952D0, 0.002380952380952D0, 0.004761904761904D0,
  29. 6 0.023809523809524D0, 0.023809523809524D0, 0.002976190476D0,
  30. 7 0.001190476190D0, 0.0008928571429D0, 0.001190476190D0,
  31. 8 0.002976190476D0, 0.0006613756614D0, 0.0003968253968D0,
  32. 9 0.0003968253968D0, 0.0006613756614D0/
  33. RES = 0.D0
  34. DO 30 I = 1, 30
  35. 30 RES = RES + C1(I)*P(I)
  36. RETURN
  37. END
  38.  
  39.  
  40.  
  41.  

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