Télécharger d106xx.eso

Retour à la liste

Numérotation des lignes :

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

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