Télécharger la plaquette Cast3M


fr en

Télécharger caluhg.eso

Retour à la liste

Numérotation des lignes : 

caluhg

  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