Télécharger la plaquette Cast3M


fr en

Télécharger d106xx.eso

Retour à la liste

Numérotation des lignes : 

d106xx

  1. C D106XX    SOURCE    FANDEUR   13/01/29    21:15:26     7683
  2.       SUBROUTINE D106XX(MODE,F,A,B,N,X,W)
  3.       IMPLICIT INTEGER(I-N)
  4.       IMPLICIT REAL*8(A-H,O-Z)
  5.       EXTERNAL F
  6.  
  7. -INC PPARAM
  8. -INC CCOPTIO
  9. *      SAVE NW
  10.       DIMENSION X(*),W(*)
  11. C
  12. C     ******************************************************************
  13. C
  14. C     MODE=1 ... SETS X(1) EQUAL TO THE INTEGRAL OF THE FUNCTION F OVER
  15. C     THE INTERVAL (A,B) AS APPROXIMATED BY THE GAUSS-LEGENDRE QUADRAT-
  16. C     URE FORMULA OF ORDER N.
  17. C
  18. C     MODE=2 ... STORES IN ARRAYS X AND W, RESPECTIVELY, THE ABSCISSAE
  19. C     AND WEIGHTS FOR N-POINT GAUSS-LEGENDRE QUADRATURE OVER THE INTER-
  20. C     VAL (A,B).
  21. C
  22. C     ******************************************************************
  23. C
  24.       DIMENSION KTAB(36),XT(273),WT(273)
  25. C
  26. C  TABLE OF INITIAL SUBSCRIPTS.
  27.       DATA KTAB(1) /0/
  28.       DATA KTAB(2) /  1/
  29.       DATA KTAB(3) /  2/
  30.       DATA KTAB(4) /  4/
  31.       DATA KTAB(5) /  6/
  32.       DATA KTAB(6) /  9/
  33.       DATA KTAB(7) / 12/
  34.       DATA KTAB(8) / 16/
  35.       DATA KTAB(9) / 20/
  36.       DATA KTAB(10)/ 25/
  37.       DATA KTAB(11)/ 30/
  38.       DATA KTAB(12)/ 36/
  39.       DATA KTAB(13)/ 42/
  40.       DATA KTAB(14)/ 49/
  41.       DATA KTAB(15)/ 56/
  42.       DATA KTAB(16)/ 64/
  43.       DATA KTAB(17)/ 72/
  44.       DATA KTAB(18)/ 82/
  45.       DATA KTAB(19)/ 0/
  46.       DATA KTAB(20)/ 94/
  47.       DATA KTAB(21) /0/
  48.       DATA KTAB(22)/110/
  49.       DATA KTAB(23) /0/
  50.       DATA KTAB(24)/130/
  51.       DATA KTAB(25) /0/
  52.       DATA KTAB(26) /0/
  53.       DATA KTAB(27) /0/
  54.       DATA KTAB(28)/154/
  55.       DATA KTAB(29) /0/
  56.       DATA KTAB(30) /0/
  57.       DATA KTAB(31) /0/
  58.       DATA KTAB(32)/186/
  59.       DATA KTAB(33) /0/
  60.       DATA KTAB(34) /0/
  61.       DATA KTAB(35) /0/
  62.       DATA KTAB(36)/226/
  63. C
  64. C  TABLE OF ABSCISSAE X AND WEIGHTS W FOR INTERVAL (-1,+1).
  65. C
  66. C     (EACH DECIMAL CONSTANT HAS MAXIMUM CDC NO-DIAGNOSTIC PRECISION.)
  67. C
  68. C  N=2.
  69.       DATA XT(1)  /5.7735026918963 D-01/, WT(1)  /1.00000000000000D+00/
  70. C  N=3.
  71.       DATA XT(2)  /7.7459666924148 D-01/, WT(2)  /5.5555555555556 D-01/
  72.       DATA XT(3)  /0.                  /, WT(3)  /8.8888888888889 D-01/
  73. C  N=4.
  74.       DATA XT(4)  /8.6113631159405 D-01/, WT(4)  /3.4785484513745 D-01/
  75.       DATA XT(5)  /3.3998104358486 D-01/, WT(5)  /6.5214515486255 D-01/
  76. C  N=5.
  77.       DATA XT(6)  /9.0617984593866 D-01/, WT(6)  /2.36926885056189D-01/
  78.       DATA XT(7)  /5.3846931010568 D-01/, WT(7)  /4.7862867049937 D-01/
  79.       DATA XT(8)  /0.                  /, WT(8)  /5.6888888888889 D-01/
  80. C  N=6.
  81.       DATA XT(9)  /9.3246951420315 D-01/, WT(9)  /1.71324492379170D-01/
  82.       DATA XT(10) /6.6120938646626 D-01/, WT(10) /3.6076157304814 D-01/
  83.       DATA XT(11) /2.38619186083197D-01/, WT(11) /4.6791393457269 D-01/
  84. C  N=7.
  85.       DATA XT(12) /9.4910791234276 D-01/, WT(12) /1.29484966168870D-01/
  86.       DATA XT(13) /7.4153118559939 D-01/, WT(13) /2.79705391489277D-01/
  87.       DATA XT(14) /4.0584515137740 D-01/, WT(14) /3.8183005050512 D-01/
  88.       DATA XT(15) /0.                  /, WT(15) /4.1795918367347 D-01/
  89. C  N=8.
  90.       DATA XT(16) /9.6028985649754 D-01/, WT(16) /1.01228536290376D-01/
  91.       DATA XT(17) /7.9666647741363 D-01/, WT(17) /2.22381034453374D-01/
  92.       DATA XT(18) /5.2553240991633 D-01/, WT(18) /3.1370664587789 D-01/
  93.       DATA XT(19) /1.83434642495650D-01/, WT(19) /3.6268378337836 D-01/
  94. C  N=9.
  95.       DATA XT(20) /9.6816023950763 D-01/, WT(20) /8.1274388361574 D-02/
  96.       DATA XT(21) /8.3603110732664 D-01/, WT(21) /1.80648160694857D-01/
  97.       DATA XT(22) /6.1337143270059 D-01/, WT(22) /2.60610696402935D-01/
  98.       DATA XT(23) /3.2425342340381 D-01/, WT(23) /3.1234707704000 D-01/
  99.       DATA XT(24) /0.                  /, WT(24) /3.3023935500126 D-01/
  100. C  N=10.
  101.       DATA XT(25) /9.7390652851717 D-01/, WT(25) /6.6671344308688 D-02/
  102.       DATA XT(26) /8.6506336668898 D-01/, WT(26) /1.49451349150581D-01/
  103.       DATA XT(27) /6.7940956829902 D-01/, WT(27) /2.19086362515982D-01/
  104.       DATA XT(28) /4.3339539412925 D-01/, WT(28) /2.69266719309996D-01/
  105.       DATA XT(29) /1.48874338981631D-01/, WT(29) /2.9552422471475 D-01/
  106. C  N=11.
  107.       DATA XT(30) /9.7822865814606 D-01/, WT(30) /5.5668567116174 D-02/
  108.       DATA XT(31) /8.8706259976810 D-01/, WT(31) /1.25580369464905D-01/
  109.       DATA XT(32) /7.3015200557405 D-01/, WT(32) /1.86290210927734D-01/
  110.       DATA XT(33) /5.1909612920681 D-01/, WT(33) /2.33193764591990D-01/
  111.       DATA XT(34) /2.69543155952345D-01/, WT(34) /2.62804544510247D-01/
  112.       DATA XT(35) /0.                  /, WT(35) /2.72925086777901D-01/
  113. C  N=12.
  114.       DATA XT(36) /9.8156063424672 D-01/, WT(36) /4.7175336386512 D-02/
  115.       DATA XT(37) /9.0411725637047 D-01/, WT(37) /1.06939325995318D-01/
  116.       DATA XT(38) /7.6990267419430 D-01/, WT(38) /1.60078328543346D-01/
  117.       DATA XT(39) /5.8731795428662 D-01/, WT(39) /2.03167426723066D-01/
  118.       DATA XT(40) /3.6783149899818 D-01/, WT(40) /2.33492536538355D-01/
  119.       DATA XT(41) /1.25233408511469D-01/, WT(41) /2.49147045813403D-01/
  120. C  N=13.
  121.       DATA XT(42) /9.8418305471859 D-01/, WT(42) /4.0484004765316 D-02/
  122.       DATA XT(43) /9.1759839922298 D-01/, WT(43) /9.2121499837728 D-02/
  123.       DATA XT(44) /8.0157809073331 D-01/, WT(44) /1.38873510219787D-01/
  124.       DATA XT(45) /6.4234933944034 D-01/, WT(45) /1.78145980761946D-01/
  125.       DATA XT(46) /4.4849275103645 D-01/, WT(46) /2.07816047536889D-01/
  126.       DATA XT(47) /2.30458315955135D-01/, WT(47) /2.26283180262897D-01/
  127.       DATA XT(48) /0.                  /, WT(48) /2.32551553230874D-01/
  128. C  N=14.
  129.       DATA XT(49) /9.8628380869681 D-01/, WT(49) /3.5119460331752 D-02/
  130.       DATA XT(50) /9.2843488366357 D-01/, WT(50) /8.0158087159760 D-02/
  131.       DATA XT(51) /8.2720131506976 D-01/, WT(51) /1.21518570687903D-01/
  132.       DATA XT(52) /6.8729290481169 D-01/, WT(52) /1.57203167158194D-01/
  133.       DATA XT(53) /5.1524863635815 D-01/, WT(53) /1.85538397477938D-01/
  134.       DATA XT(54) /3.1911236892789 D-01/, WT(54) /2.05198463721296D-01/
  135.       DATA XT(55) /1.08054948707344D-01/, WT(55) /2.15263853463158D-01/
  136. C  N=15.
  137.       DATA XT(56) /9.8799251802049 D-01/, WT(56) /3.0753241996117 D-02/
  138.       DATA XT(57) /9.3727339240071 D-01/, WT(57) /7.0366047488108 D-02/
  139.       DATA XT(58) /8.4820658341043 D-01/, WT(58) /1.07159220467172D-01/
  140.       DATA XT(59) /7.2441773136017 D-01/, WT(59) /1.39570677926154D-01/
  141.       DATA XT(60) /5.7097217260854 D-01/, WT(60) /1.66269205816994D-01/
  142.       DATA XT(61) /3.9415134707756 D-01/, WT(61) /1.86161000015562D-01/
  143.       DATA XT(62) /2.01194093997435D-01/, WT(62) /1.98431485327112D-01/
  144.       DATA XT(63) /0.                  /, WT(63) /2.02578241925561D-01/
  145. C  N=16.
  146.       DATA XT(64) /9.8940093499165 D-01/, WT(64) /2.71524594117541D-02/
  147.       DATA XT(65) /9.4457502307323 D-01/, WT(65) /6.2253523938648 D-02/
  148.       DATA XT(66) /8.6563120238783 D-01/, WT(66) /9.5158511682493 D-02/
  149.       DATA XT(67) /7.5540440835500 D-01/, WT(67) /1.24628971255534D-01/
  150.       DATA XT(68) /6.1787624440264 D-01/, WT(68) /1.49595988816577D-01/
  151.       DATA XT(69) /4.5801677765723 D-01/, WT(69) /1.69156519395003D-01/
  152.       DATA XT(70) /2.8160355077926 D-01/, WT(70) /1.82603415044924D-01/
  153.       DATA XT(71) /9.5012509837637 D-02/, WT(71) /1.89450610455068D-01/
  154. C  N=20.
  155.       DATA XT(72) /9.9312859918509 D-01/, WT(72) /1.76140071391521D-02/
  156.       DATA XT(73) /9.6397192727791 D-01/, WT(73) /4.0601429800387 D-02/
  157.       DATA XT(74) /9.1223442825133 D-01/, WT(74) /6.2672048334109 D-02/
  158.       DATA XT(75) /8.3911697182222 D-01/, WT(75) /8.3276741576705 D-02/
  159.       DATA XT(76) /7.4633190646015 D-01/, WT(76) /1.01930119817240D-01/
  160.       DATA XT(77) /6.3605368072652 D-01/, WT(77) /1.18194531961518D-01/
  161.       DATA XT(78) /5.1086700195083 D-01/, WT(78) /1.31688638449177D-01/
  162.       DATA XT(79) /3.7370608871542 D-01/, WT(79) /1.42096109318382D-01/
  163.       DATA XT(80) /2.27785851141645D-01/, WT(80) /1.49172986472604D-01/
  164.       DATA XT(81) /7.6526521133497 D-02/, WT(81) /1.52753387130726D-01/
  165. C  N=24.
  166.       DATA XT(82) /9.9518721999702 D-01/, WT(82) /1.23412297999872D-02/
  167.       DATA XT(83) /9.7472855597131 D-01/, WT(83) /2.8531388628934 D-02/
  168.       DATA XT(84) /9.3827455200273 D-01/, WT(84) /4.4277438817420 D-02/
  169.       DATA XT(85) /8.8641552700440 D-01/, WT(85) /5.9298584915437 D-02/
  170.       DATA XT(86) /8.2000198597390 D-01/, WT(86) /7.3346481411080 D-02/
  171.       DATA XT(87) /7.4012419157855 D-01/, WT(87) /8.6190161531953 D-02/
  172.       DATA XT(88) /6.4809365193698 D-01/, WT(88) /9.7618652104114 D-02/
  173.       DATA XT(89) /5.4542147138884 D-01/, WT(89) /1.07444270115966D-01/
  174.       DATA XT(90) /4.3379350762605 D-01/, WT(90) /1.15505668053726D-01/
  175.       DATA XT(91) /3.1504267969616 D-01/, WT(91) /1.21670472927803D-01/
  176.       DATA XT(92) /1.91118867473616D-01/, WT(92) /1.25837456346828D-01/
  177.       DATA XT(93) /6.4056892862606 D-02/, WT(93) /1.27938195346752D-01/
  178. C  N=32.
  179.       DATA XT(94) /9.9726386184948 D-01/, WT(94) /7.0186100094701 D-03/
  180.       DATA XT(95) /9.8561151154527 D-01/, WT(95) /1.62743947309057D-02/
  181.       DATA XT(96) /9.6476225558751 D-01/, WT(96) /2.53920653092621D-02/
  182.       DATA XT(97) /9.3490607593774 D-01/, WT(97) /3.4273862913021 D-02/
  183.       DATA XT(98) /8.9632115576605 D-01/, WT(98) /4.2835898022227 D-02/
  184.       DATA XT(99) /8.4936761373257 D-01/, WT(99) /5.0998059262376 D-02/
  185.       DATA XT(100)/7.9448379596794 D-01/, WT(100)/5.8684093478536 D-02/
  186.       DATA XT(101)/7.3218211874029 D-01/, WT(101)/6.5822222776362 D-02/
  187.       DATA XT(102)/6.6304426693022 D-01/, WT(102)/7.2345794108849 D-02/
  188.       DATA XT(103)/5.8771575724076 D-01/, WT(103)/7.8193895787070 D-02/
  189.       DATA XT(104)/5.0689990893223 D-01/, WT(104)/8.3311924226947 D-02/
  190.       DATA XT(105)/4.2135127613064 D-01/, WT(105)/8.7652093004404 D-02/
  191.       DATA XT(106)/3.3186860228213 D-01/, WT(106)/9.1173878695764 D-02/
  192.       DATA XT(107)/2.39287362252137D-01/, WT(107)/9.3844399080805 D-02/
  193.       DATA XT(108)/1.44471961582796D-01/, WT(108)/9.5638720079275 D-02/
  194.       DATA XT(109)/4.8307665687738 D-02/, WT(109)/9.6540088514728 D-02/
  195. C  N=40.
  196.       DATA XT(110)/9.9823770971056 D-01/, WT(110)/4.5212770985332 D-03/
  197.       DATA XT(111)/9.9072623869946 D-01/, WT(111)/1.04982845311528D-02/
  198.       DATA XT(112)/9.7725994998377 D-01/, WT(112)/1.64210583819079D-02/
  199.       DATA XT(113)/9.5791681921379 D-01/, WT(113)/2.22458491941670D-02/
  200.       DATA XT(114)/9.3281280827868 D-01/, WT(114)/2.79370069800234D-02/
  201.       DATA XT(115)/9.0209880696887 D-01/, WT(115)/3.3460195282548 D-02/
  202.       DATA XT(116)/8.6595950321226 D-01/, WT(116)/3.8782167974472 D-02/
  203.       DATA XT(117)/8.2461223083331 D-01/, WT(117)/4.3870908185673 D-02/
  204.       DATA XT(118)/7.7830565142652 D-01/, WT(118)/4.8695807635072 D-02/
  205.       DATA XT(119)/7.2731825518993 D-01/, WT(119)/5.3227846983937 D-02/
  206.       DATA XT(120)/6.7195668461418 D-01/, WT(120)/5.7439769099392 D-02/
  207.       DATA XT(121)/6.1255388966798 D-01/, WT(121)/6.1306242492929 D-02/
  208.       DATA XT(122)/5.4946712509513 D-01/, WT(122)/6.4804013456601 D-02/
  209.       DATA XT(123)/4.8307580168618 D-01/, WT(123)/6.7912045815234 D-02/
  210.       DATA XT(124)/4.1377920437161 D-01/, WT(124)/7.0611647391287 D-02/
  211.       DATA XT(125)/3.4199409082576 D-01/, WT(125)/7.2886582395804 D-02/
  212.       DATA XT(126)/2.68152185007254D-01/, WT(126)/7.4723169057968 D-02/
  213.       DATA XT(127)/1.92697580701371D-01/, WT(127)/7.6110361900626 D-02/
  214.       DATA XT(128)/1.16084070675255D-01/, WT(128)/7.7039818164248 D-02/
  215.       DATA XT(129)/3.8772417506051 D-02/, WT(129)/7.7505947978425 D-02/
  216. C  N=48.
  217.       DATA XT(130)/9.9877100725243 D-01/, WT(130)/3.1533460523058 D-03/
  218.       DATA XT(131)/9.9353017226635 D-01/, WT(131)/7.3275539012763 D-03/
  219.       DATA XT(132)/9.8412458372283 D-01/, WT(132)/1.14772345792345D-02/
  220.       DATA XT(133)/9.7059159254625 D-01/, WT(133)/1.55793157229438D-02/
  221.       DATA XT(134)/9.5298770316043 D-01/, WT(134)/1.96161604573555D-02/
  222.       DATA XT(135)/9.3138669070655 D-01/, WT(135)/2.35707608393244D-02/
  223.       DATA XT(136)/9.0587913671557 D-01/, WT(136)/2.74265097083569D-02/
  224.       DATA XT(137)/8.7657202027425 D-01/, WT(137)/3.1167227832798 D-02/
  225.       DATA XT(138)/8.4358826162439 D-01/, WT(138)/3.4777222564770 D-02/
  226.       DATA XT(139)/8.0706620402944 D-01/, WT(139)/3.8241351065831 D-02/
  227.       DATA XT(140)/7.6715903251574 D-01/, WT(140)/4.1545082943465 D-02/
  228.       DATA XT(141)/7.2403413092381 D-01/, WT(141)/4.4674560856694 D-02/
  229.       DATA XT(142)/6.7787237963266 D-01/, WT(142)/4.7616658492490 D-02/
  230.       DATA XT(143)/6.2886739677651 D-01/, WT(143)/5.0359035553854 D-02/
  231.       DATA XT(144)/5.7722472608397 D-01/, WT(144)/5.2890189485194 D-02/
  232.       DATA XT(145)/5.2316097472223 D-01/, WT(145)/5.5199503699984 D-02/
  233.       DATA XT(146)/4.6690290475096 D-01/, WT(146)/5.7277292100403 D-02/
  234.       DATA XT(147)/4.0868648199072 D-01/, WT(147)/5.9114839698396 D-02/
  235.       DATA XT(148)/3.4875588629216 D-01/, WT(148)/6.0704439165894 D-02/
  236.       DATA XT(149)/2.8736248735546 D-01/, WT(149)/6.2039423159893 D-02/
  237.       DATA XT(150)/2.24763790394689D-01/, WT(150)/6.3114192286254 D-02/
  238.       DATA XT(151)/1.61222356068892D-01/, WT(151)/6.3924238584648 D-02/
  239.       DATA XT(152)/9.7004699209463 D-02/, WT(152)/6.4466164435950 D-02/
  240.       DATA XT(153)/3.2380170962869 D-02/, WT(153)/6.4737696812684 D-02/
  241. C  N=64.
  242.       DATA XT(154)/9.9930504173577 D-01/, WT(154)/1.78328072169643D-03/
  243.       DATA XT(155)/9.9634011677196 D-01/, WT(155)/4.1470332605625 D-03/
  244.       DATA XT(156)/9.9101337147674 D-01/, WT(156)/6.5044579689784 D-03/
  245.       DATA XT(157)/9.8333625388463 D-01/, WT(157)/8.8467598263639 D-03/
  246.       DATA XT(158)/9.7332682778991 D-01/, WT(158)/1.11681394601311D-02/
  247.       DATA XT(159)/9.6100879965205 D-01/, WT(159)/1.34630478967186D-02/
  248.       DATA XT(160)/9.4641137485840 D-01/, WT(160)/1.57260304760247D-02/
  249.       DATA XT(161)/9.2956917213194 D-01/, WT(161)/1.79517157756973D-02/
  250.       DATA XT(162)/9.1052213707850 D-01/, WT(162)/2.01348231535302D-02/
  251.       DATA XT(163)/8.8931544599511 D-01/, WT(163)/2.22701738083833D-02/
  252.       DATA XT(164)/8.6599939815409 D-01/, WT(164)/2.43527025687109D-02/
  253.       DATA XT(165)/8.4062929625258 D-01/, WT(165)/2.63774697150547D-02/
  254.       DATA XT(166)/8.1326531512280 D-01/, WT(166)/2.8339672614259 D-02/
  255.       DATA XT(167)/7.8397235894334 D-01/, WT(167)/3.0234657072402 D-02/
  256.       DATA XT(168)/7.5281990726053 D-01/, WT(168)/3.2057928354852 D-02/
  257.       DATA XT(169)/7.1988185017161 D-01/, WT(169)/3.3805161837142 D-02/
  258.       DATA XT(170)/6.8523631305423 D-01/, WT(170)/3.5472213256882 D-02/
  259.       DATA XT(171)/6.4896547125466 D-01/, WT(171)/3.7055128540240 D-02/
  260.       DATA XT(172)/6.1115535517239 D-01/, WT(172)/3.8550153178616 D-02/
  261.       DATA XT(173)/5.7189564620263 D-01/, WT(173)/3.9953741132720 D-02/
  262.       DATA XT(174)/5.3127946401989 D-01/, WT(174)/4.1262563242624 D-02/
  263.       DATA XT(175)/4.8940314570705 D-01/, WT(175)/4.2473515123654 D-02/
  264.       DATA XT(176)/4.4636601725346 D-01/, WT(176)/4.3583724529323 D-02/
  265.       DATA XT(177)/4.0227015796399 D-01/, WT(177)/4.4590558163757 D-02/
  266.       DATA XT(178)/3.5722015833767 D-01/, WT(178)/4.5491627927418 D-02/
  267.       DATA XT(179)/3.1132287199021 D-01/, WT(179)/4.6284796581314 D-02/
  268.       DATA XT(180)/2.64687162208767D-01/, WT(180)/4.6968182816210 D-02/
  269.       DATA XT(181)/2.17423643740007D-01/, WT(181)/4.7540165714830 D-02/
  270.       DATA XT(182)/1.69644420423993D-01/, WT(182)/4.7999388596458 D-02/
  271.       DATA XT(183)/1.21462819296121D-01/, WT(183)/4.8344762234803 D-02/
  272.       DATA XT(184)/7.2993121787799 D-02/, WT(184)/4.8575467441503 D-02/
  273.       DATA XT(185)/2.43502926634244D-02/, WT(185)/4.8690957009140 D-02/
  274. C  N=80.
  275.       DATA XT(186)/9.9955382265163 D-01/, WT(186)/1.14495000318694D-03/
  276.       DATA XT(187)/9.9764986439824 D-01/, WT(187)/2.66353358951268D-03/
  277.       DATA XT(188)/9.9422754096569 D-01/, WT(188)/4.1803131246949 D-03/
  278.       DATA XT(189)/9.8929130249976 D-01/, WT(189)/5.6909224514032 D-03/
  279.       DATA XT(190)/9.8284857273863 D-01/, WT(190)/7.1929047681173 D-03/
  280.       DATA XT(191)/9.7490914058573 D-01/, WT(191)/8.6839452692609 D-03/
  281.       DATA XT(192)/9.6548508904380 D-01/, WT(192)/1.01617660411031D-02/
  282.       DATA XT(193)/9.5459076634363 D-01/, WT(193)/1.16241141207978D-02/
  283.       DATA XT(194)/9.4224276130987 D-01/, WT(194)/1.30687615924013D-02/
  284.       DATA XT(195)/9.2845987717245 D-01/, WT(195)/1.44935080405091D-02/
  285.       DATA XT(196)/9.1326310257176 D-01/, WT(196)/1.58961835837257D-02/
  286.       DATA XT(197)/8.9667557943877 D-01/, WT(197)/1.72746520562693D-02/
  287.       DATA XT(198)/8.7872256767821 D-01/, WT(198)/1.86268142082990D-02/
  288.       DATA XT(199)/8.5943140666311 D-01/, WT(199)/1.99506108781420D-02/
  289.       DATA XT(200)/8.3883147358026 D-01/, WT(200)/2.12440261157821D-02/
  290.       DATA XT(201)/8.1695413868146 D-01/, WT(201)/2.25050902463325D-02/
  291.       DATA XT(202)/7.9383271750461 D-01/, WT(202)/2.37318828659301D-02/
  292.       DATA XT(203)/7.6950242013504 D-01/, WT(203)/2.49225357641155D-02/
  293.       DATA XT(204)/7.4400029758360 D-01/, WT(204)/2.60752357675651D-02/
  294.       DATA XT(205)/7.1736518536210 D-01/, WT(205)/2.71882275004864D-02/
  295.       DATA XT(206)/6.8963764434203 D-01/, WT(206)/2.8259816057277 D-02/
  296.       DATA XT(207)/6.6085989898612 D-01/, WT(207)/2.9288369583268 D-02/
  297.       DATA XT(208)/6.3107577304687 D-01/, WT(208)/3.0272321759558 D-02/
  298.       DATA XT(209)/6.0033062282975 D-01/, WT(209)/3.1210174188115 D-02/
  299.       DATA XT(210)/5.6867126812271 D-01/, WT(210)/3.2100498673488 D-02/
  300.       DATA XT(211)/5.3614592089713 D-01/, WT(211)/3.2941939397645 D-02/
  301.       DATA XT(212)/5.0280411188878 D-01/, WT(212)/3.3733214984612 D-02/
  302.       DATA XT(213)/4.6869661517054 D-01/, WT(213)/3.4473120451754 D-02/
  303.       DATA XT(214)/4.3387537083176 D-01/, WT(214)/3.5160529044748 D-02/
  304.       DATA XT(215)/3.9839340588197 D-01/, WT(215)/3.5794393953416 D-02/
  305.       DATA XT(216)/3.6230475349949 D-01/, WT(216)/3.6373749905836 D-02/
  306.       DATA XT(217)/3.2566437074770 D-01/, WT(217)/3.6897714638276 D-02/
  307.       DATA XT(218)/2.8852805488451 D-01/, WT(218)/3.7365490238730 D-02/
  308.       DATA XT(219)/2.50952358392272D-01/, WT(219)/3.7776364362001 D-02/
  309.       DATA XT(220)/2.12994502857666D-01/, WT(220)/3.8129711314478 D-02/
  310.       DATA XT(221)/1.74712291832647D-01/, WT(221)/3.8424993006959 D-02/
  311.       DATA XT(222)/1.36164022809144D-01/, WT(222)/3.8661759774076 D-02/
  312.       DATA XT(223)/9.7408398441585 D-02/, WT(223)/3.8839651059052 D-02/
  313.       DATA XT(224)/5.8504437152421 D-02/, WT(224)/3.8958395962770 D-02/
  314.       DATA XT(225)/1.95113832567940D-02/, WT(225)/3.9017813656307 D-02/
  315. C  N=96.
  316.       DATA XT(226)/9.9968950388323 D-01/, WT(226)/7.9679206555201 D-04/
  317.       DATA XT(227)/9.9836437586318 D-01/, WT(227)/1.85396078894692D-03/
  318.       DATA XT(228)/9.9598184298721 D-01/, WT(228)/2.9107318179349 D-03/
  319.       DATA XT(229)/9.9254390032376 D-01/, WT(229)/3.9645543384447 D-03/
  320.       DATA XT(230)/9.8805412632962 D-01/, WT(230)/5.0142027429275 D-03/
  321.       DATA XT(231)/9.8251726356301 D-01/, WT(231)/6.0585455042360 D-03/
  322.       DATA XT(232)/9.7593917458514 D-01/, WT(232)/7.0964707911539 D-03/
  323.       DATA XT(233)/9.6832682846326 D-01/, WT(233)/8.1268769256988 D-03/
  324.       DATA XT(234)/9.5968829144874 D-01/, WT(234)/9.1486712307834 D-03/
  325.       DATA XT(235)/9.5003271778444 D-01/, WT(235)/1.01607705350084D-02/
  326.       DATA XT(236)/9.3937033975276 D-01/, WT(236)/1.11621020998385D-02/
  327.       DATA XT(237)/9.2771245672231 D-01/, WT(237)/1.21516046710883D-02/
  328.       DATA XT(238)/9.1507142312090 D-01/, WT(238)/1.31282295669616D-02/
  329.       DATA XT(239)/9.0146063531585 D-01/, WT(239)/1.40909417723149D-02/
  330.       DATA XT(240)/8.8689451740242 D-01/, WT(240)/1.50387210269949D-02/
  331.       DATA XT(241)/8.7138850590930 D-01/, WT(241)/1.59705629025623D-02/
  332.       DATA XT(242)/8.5495903343460 D-01/, WT(242)/1.68854798642452D-02/
  333.       DATA XT(243)/8.3762351122819 D-01/, WT(243)/1.77825023160453D-02/
  334.       DATA XT(244)/8.1940031073793 D-01/, WT(244)/1.86606796274115D-02/
  335.       DATA XT(245)/8.0030874413914 D-01/, WT(245)/1.95190811401450D-02/
  336.       DATA XT(246)/7.8036904386743 D-01/, WT(246)/2.03567971543333D-02/
  337.       DATA XT(247)/7.5960234117665 D-01/, WT(247)/2.11729398921913D-02/
  338.       DATA XT(248)/7.3803064374440 D-01/, WT(248)/2.19666444387443D-02/
  339.       DATA XT(249)/7.1567681234897 D-01/, WT(249)/2.27370696583294D-02/
  340.       DATA XT(250)/6.9256453664217 D-01/, WT(250)/2.34833990859262D-02/
  341.       DATA XT(251)/6.6871831004392 D-01/, WT(251)/2.42048417923647D-02/
  342.       DATA XT(252)/6.4416340378497 D-01/, WT(252)/2.49006332224836D-02/
  343.       DATA XT(253)/6.1892584012547 D-01/, WT(253)/2.55700360053494D-02/
  344.       DATA XT(254)/5.9303236477757 D-01/, WT(254)/2.62123407356724D-02/
  345.       DATA XT(255)/5.6651041856140 D-01/, WT(255)/2.68268667255918D-02/
  346.       DATA XT(256)/5.3938810832436 D-01/, WT(256)/2.74129627260292D-02/
  347.       DATA XT(257)/5.1169417715467 D-01/, WT(257)/2.79700076168483D-02/
  348.       DATA XT(258)/4.8345797392060 D-01/, WT(258)/2.8497411065085 D-02/
  349.       DATA XT(259)/4.5470942216774 D-01/, WT(259)/2.8994614150555 D-02/
  350.       DATA XT(260)/4.2547898840730 D-01/, WT(260)/2.9461089958168 D-02/
  351.       DATA XT(261)/3.9579764982891 D-01/, WT(261)/2.9896344136328 D-02/
  352.       DATA XT(262)/3.6569686147231 D-01/, WT(262)/3.0299915420828 D-02/
  353.       DATA XT(263)/3.3520852289263 D-01/, WT(263)/3.0671376123669 D-02/
  354.       DATA XT(264)/3.0436494435450 D-01/, WT(264)/3.1010332586314 D-02/
  355.       DATA XT(265)/2.73198812591049D-01/, WT(265)/3.1316425596861 D-02/
  356.       DATA XT(266)/2.41743156163840D-01/, WT(266)/3.1589330770727 D-02/
  357.       DATA XT(267)/2.10031310460567D-01/, WT(267)/3.1828758894411 D-02/
  358.       DATA XT(268)/1.78096882367619D-01/, WT(268)/3.2034456231993 D-02/
  359.       DATA XT(269)/1.45973714654897D-01/, WT(269)/3.2206204794030 D-02/
  360.       DATA XT(270)/1.13695850110666D-01/, WT(270)/3.2343822568576 D-02/
  361.       DATA XT(271)/8.1297495464426 D-02/, WT(271)/3.2447163714064 D-02/
  362.       DATA XT(272)/4.8812985136050 D-02/, WT(272)/3.2516118713869 D-02/
  363.       DATA XT(273)/1.62767448496030D-02/, WT(273)/3.2550614492363 D-02/
  364. C
  365. C     ******************************************************************
  366. C
  367. C  START.  IF N HAS CHANGED VALUE, SET K TO INITIAL SUBSCRIPT OF ARRAYS
  368. C  XT AND WT.
  369. C
  370. *      IF(N.EQ.NW) GO TO 2
  371. *      NWSAVE=NW
  372.       if(n.le.1.or.n.gt.96) go to 100
  373.       L=N
  374.       IF(N.LE.16) GO TO 1
  375.       L=(N/4)+12
  376. C
  377.     1 K=KTAB(L)
  378.       if(K.eq.0) go to 100
  379.       M=N/2
  380.       JMAX=K+M-1
  381.       JMID=0
  382.       IF(N-M-M.NE.0) JMID=K+M
  383. C
  384.     2 ALPHA=0.5D0*(B+A)
  385.       BETA=0.5D0*(B-A)
  386.       IF(MODE.EQ.2) GO TO 4
  387. C
  388. C  (MODE=1. EVALUATE QUDRATURE SUM.)
  389. C
  390.       SUM=0.D0
  391.       DO 3 J=K,JMAX
  392.          DELTA=BETA*XT(J)
  393.          SUM=SUM+WT(J)*(F(ALPHA+DELTA)+F(ALPHA-DELTA))
  394.     3 CONTINUE
  395.       IF(JMID.NE.0) SUM=SUM+WT(JMID)*F(ALPHA)
  396.       X(1)=BETA*SUM
  397.       RETURN
  398. C
  399. C  (MODE=2. SET ARRAYS X AND W.)
  400. C
  401. C    jamais appele en mode 2   TC
  402.     4 JLAST=(N+1)/2
  403.       DO 5 J=1,JLAST
  404.          JTAB=K-1+J
  405.          WTEMP=BETA*WT(JTAB)
  406.          DELTA=BETA*XT(JTAB)
  407.          X(J)=ALPHA-DELTA
  408.          W(J)=WTEMP
  409.          JP=NW+1-J
  410.          X(JP)=ALPHA+DELTA
  411.          W(JP)=WTEMP
  412.     5 CONTINUE
  413.       RETURN
  414. C
  415. C  ERROR MESSAGE FOR NON-PERMISSIBLE N.
  416. C
  417.  100  IF(MODE.EQ.1) X(1)=0.D0
  418.       INTERR(1)=N
  419.       CALL ERREUR(531)
  420. *      NW=NWSAVE
  421.       RETURN
  422.       END
  423.  
  424.  
  425.  

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