Télécharger caluhg.eso

Retour à la liste

Numérotation des lignes :

  1. C CALUHG SOURCE CHAT 05/01/12 21:48:06 5004
  2. SUBROUTINE CALUHG(U,H,NG)
  3. C************************************************************************
  4. C CALCUL DES POIDS ET ABCISSES DES POINTS DE GAUSS SUR LE SEGMENT
  5. C (-1,1) DE 1 A 6 RACINES
  6. C
  7. C NG NOMBRE DE RACINES
  8. C U ABCISSES
  9. C H POIDS
  10. C************************************************************************
  11. IMPLICIT INTEGER(I-N)
  12. IMPLICIT REAL*8 (A-H,O-Z)
  13. REAL*8 U(NG),H(NG),USQ
  14. C***
  15.  
  16.  
  17.  
  18.  
  19. GO TO (10,20,30,40,50,60),NG
  20. CALL ARRET(0)
  21. C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  22. 10 CONTINUE
  23. U(1)=0.D0
  24. H(1)=2.D0
  25. RETURN
  26. C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  27. 20 CONTINUE
  28. U(1)=-SQRT(3.D0)/3.D0
  29. U(2)=-U(1)
  30. H(1)=1.D0
  31. H(2)=1.D0
  32. RETURN
  33. C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  34. 30 CONTINUE
  35. U(1)=-SQRT(0.6D0)
  36. U(2)=0.D0
  37. U(3)=-U(1)
  38. H(1)=5.D0/9.D0
  39. H(3)=H(1)
  40. H(2)=8.D0/9.D0
  41. RETURN
  42. C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  43. 40 CONTINUE
  44. USQ=SQRT(4.8D0)
  45. U(1)=(3.D0+USQ)/7.D0
  46. U(1)=-SQRT(U(1))
  47. U(4)=-U(1)
  48. U(2)=(3.D0-USQ)/7.D0
  49. U(2)=-SQRT(U(2))
  50. U(3)=-U(2)
  51. H(1)=(7.2D0-USQ)/14.4D0
  52. H(2)=(7.2D0+USQ)/14.4D0
  53. H(3)=H(2)
  54. H(4)=H(1)
  55. RETURN
  56. C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  57. 50 CONTINUE
  58. USQ=SQRT(1120.D0)
  59. U(1)=(70.D0+USQ)/126.D0
  60. U(1)=-SQRT(U(1))
  61. U(5)=-U(1)
  62. U(3)=0.D0
  63. U(2)=(70.D0-USQ)/126.D0
  64. U(2)=-SQRT(U(2))
  65. U(4)=-U(2)
  66. H(1)=21.D0/5.D0*(28.D0+5.D0*USQ)/(1120.D0+70.D0*USQ)
  67. H(2)=21.D0/5.D0*(28.D0-5.D0*USQ)/(1120.D0-70.D0*USQ)
  68. H(3)=128.D0/225.D0
  69. H(4)=H(2)
  70. H(5)=H(1)
  71. RETURN
  72. C>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  73. 60 CONTINUE
  74. U(1)=-0.932469514203152D0
  75. U(2)=-0.661209386466265D0
  76. U(3)=-0.238619186083197D0
  77. U(4)=-U(3)
  78. U(5)=-U(2)
  79. U(6)=-U(1)
  80. H(1)=0.171324492379170D0
  81. H(2)=0.360761573048139D0
  82. H(3)=0.467913934572691D0
  83. H(4)=H(3)
  84. H(5)=H(2)
  85. H(6)=H(1)
  86. RETURN
  87. END
  88.  
  89.  
  90.  

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