* * $Id: d107s164.F,v 1.1.1.1 1996/04/01 15:03:08 mclareni Exp $ * * $Log: d107s164.F,v $ * Revision 1.1.1.1 1996/04/01 15:03:08 mclareni * Mathlib gen * * #include "sys/CERNLIB_machine.h" #include "_gen/pilot.h" SUBROUTINE D107S1(MODE,F,A,B,N,X,W) #include "imp64.inc" CHARACTER*6 NAME(2) DIMENSION X(*),W(*),KTBA(97),U(273),V(273) PARAMETER (Z1 = 1, HF = Z1/2) CHARACTER*80 ERRTXT #if defined(CERNLIB_DOUBLE) DATA NAME /'DGQUAD','DGSET'/ #endif #if !defined(CERNLIB_DOUBLE) DATA NAME /'RGQUAD','RGSET'/ #endif DATA KTBA 1/0,1,2,4,6,9,12,16,20,25,30,36,42,49,56,64,3*0,72,3*0,82,7*0,94, 2 7*0,110,7*0,130,15*0,154,15*0,186,15*0,226,0/ C N=2. DATA U(1) /5.7735026918962576D-1/, V(1) /1/ C N=3. DATA U(2) /7.7459666924148338D-1/, V(2) /5.5555555555555556D-1/ DATA U(3) /0/ , V(3) /8.8888888888888889D-1/ C N=4. DATA U(4) /8.6113631159405258D-1/, V(4) /3.4785484513745386D-1/ DATA U(5) /3.3998104358485626D-1/, V(5) /6.5214515486254614D-1/ C N=5. DATA U(6) /9.0617984593866399D-1/, V(6) /2.3692688505618909D-1/ DATA U(7) /5.3846931010568309D-1/, V(7) /4.7862867049936647D-1/ DATA U(8) /0/, V(8) /5.6888888888888889D-1/ C N=6. DATA U(9) /9.3246951420315203D-1/, V(9) /1.7132449237917035D-1/ DATA U(10) /6.6120938646626451D-1/, V(10) /3.6076157304813861D-1/ DATA U(11) /2.3861918608319691D-1/, V(11) /4.6791393457269105D-1/ C N=7. DATA U(12) /9.4910791234275852D-1/, V(12) /1.2948496616886969D-1/ DATA U(13) /7.4153118559939444D-1/, V(13) /2.7970539148927667D-1/ DATA U(14) /4.0584515137739717D-1/, V(14) /3.8183005050511894D-1/ DATA U(15) /0/ , V(15) /4.1795918367346939D-1/ C N=8. DATA U(16) /9.6028985649753623D-1/, V(16) /1.0122853629037626D-1/ DATA U(17) /7.9666647741362674D-1/, V(17) /2.2238103445337447D-1/ DATA U(18) /5.2553240991632899D-1/, V(18) /3.1370664587788729D-1/ DATA U(19) /1.8343464249564980D-1/, V(19) /3.6268378337836198D-1/ C N=9. DATA U(20) /9.6816023950762609D-1/, V(20) /8.1274388361574412D-2/ DATA U(21) /8.3603110732663579D-1/, V(21) /1.8064816069485740D-1/ DATA U(22) /6.1337143270059040D-1/, V(22) /2.6061069640293546D-1/ DATA U(23) /3.2425342340380893D-1/, V(23) /3.1234707704000284D-1/ DATA U(24) /0/ , V(24) /3.3023935500125976D-1/ C N=10. DATA U(25) /9.7390652851717172D-1/, V(25) /6.6671344308688138D-2/ DATA U(26) /8.6506336668898451D-1/, V(26) /1.4945134915058059D-1/ DATA U(27) /6.7940956829902441D-1/, V(27) /2.1908636251598204D-1/ DATA U(28) /4.3339539412924719D-1/, V(28) /2.6926671930999636D-1/ DATA U(29) /1.4887433898163121D-1/, V(29) /2.9552422471475287D-1/ C N=11. DATA U(30) /9.7822865814605699D-1/, V(30) /5.5668567116173666D-2/ DATA U(31) /8.8706259976809530D-1/, V(31) /1.2558036946490462D-1/ DATA U(32) /7.3015200557404932D-1/, V(32) /1.8629021092773425D-1/ DATA U(33) /5.1909612920681182D-1/, V(33) /2.3319376459199048D-1/ DATA U(34) /2.6954315595234497D-1/, V(34) /2.6280454451024666D-1/ DATA U(35) /0/ , V(35) /2.7292508677790063D-1/ C N=12. DATA U(36) /9.8156063424671925D-1/, V(36) /4.7175336386511827D-2/ DATA U(37) /9.0411725637047486D-1/, V(37) /1.0693932599531843D-1/ DATA U(38) /7.6990267419430469D-1/, V(38) /1.6007832854334623D-1/ DATA U(39) /5.8731795428661745D-1/, V(39) /2.0316742672306592D-1/ DATA U(40) /3.6783149899818019D-1/, V(40) /2.3349253653835481D-1/ DATA U(41) /1.2523340851146892D-1/, V(41) /2.4914704581340279D-1/ C N=13. DATA U(42) /9.8418305471858815D-1/, V(42) /4.0484004765315880D-2/ DATA U(43) /9.1759839922297797D-1/, V(43) /9.2121499837728448D-2/ DATA U(44) /8.0157809073330991D-1/, V(44) /1.3887351021978724D-1/ DATA U(45) /6.4234933944034022D-1/, V(45) /1.7814598076194574D-1/ DATA U(46) /4.4849275103644685D-1/, V(46) /2.0781604753688850D-1/ DATA U(47) /2.3045831595513479D-1/, V(47) /2.2628318026289724D-1/ DATA U(48) /0/ , V(48) /2.3255155323087391D-1/ C N=14. DATA U(49) /9.8628380869681234D-1/, V(49) /3.5119460331751863D-2/ DATA U(50) /9.2843488366357352D-1/, V(50) /8.0158087159760210D-2/ DATA U(51) /8.2720131506976499D-1/, V(51) /1.2151857068790318D-1/ DATA U(52) /6.8729290481168547D-1/, V(52) /1.5720316715819353D-1/ DATA U(53) /5.1524863635815409D-1/, V(53) /1.8553839747793781D-1/ DATA U(54) /3.1911236892788976D-1/, V(54) /2.0519846372129560D-1/ DATA U(55) /1.0805494870734366D-1/, V(55) /2.1526385346315779D-1/ C N=15. DATA U(56) /9.8799251802048543D-1/, V(56) /3.0753241996117268D-2/ DATA U(57) /9.3727339240070590D-1/, V(57) /7.0366047488108125D-2/ DATA U(58) /8.4820658341042722D-1/, V(58) /1.0715922046717194D-1/ DATA U(59) /7.2441773136017005D-1/, V(59) /1.3957067792615431D-1/ DATA U(60) /5.7097217260853885D-1/, V(60) /1.6626920581699393D-1/ DATA U(61) /3.9415134707756337D-1/, V(61) /1.8616100001556221D-1/ DATA U(62) /2.0119409399743452D-1/, V(62) /1.9843148532711158D-1/ DATA U(63) /0/ , V(63) /2.0257824192556127D-1/ C N=16. DATA U(64) /9.8940093499164993D-1/, V(64) /2.7152459411754095D-2/ DATA U(65) /9.4457502307323258D-1/, V(65) /6.2253523938647893D-2/ DATA U(66) /8.6563120238783174D-1/, V(66) /9.5158511682492785D-2/ DATA U(67) /7.5540440835500303D-1/, V(67) /1.2462897125553387D-1/ DATA U(68) /6.1787624440264375D-1/, V(68) /1.4959598881657673D-1/ DATA U(69) /4.5801677765722739D-1/, V(69) /1.6915651939500254D-1/ DATA U(70) /2.8160355077925891D-1/, V(70) /1.8260341504492359D-1/ DATA U(71) /9.5012509837637440D-2/, V(71) /1.8945061045506850D-1/ C N=20. DATA U(72) /9.9312859918509492D-1/, V(72) /1.7614007139152118D-2/ DATA U(73) /9.6397192727791379D-1/, V(73) /4.0601429800386941D-2/ DATA U(74) /9.1223442825132591D-1/, V(74) /6.2672048334109064D-2/ DATA U(75) /8.3911697182221882D-1/, V(75) /8.3276741576704749D-2/ DATA U(76) /7.4633190646015079D-1/, V(76) /1.0193011981724044D-1/ DATA U(77) /6.3605368072651503D-1/, V(77) /1.1819453196151842D-1/ DATA U(78) /5.1086700195082710D-1/, V(78) /1.3168863844917663D-1/ DATA U(79) /3.7370608871541956D-1/, V(79) /1.4209610931838205D-1/ DATA U(80) /2.2778585114164508D-1/, V(80) /1.4917298647260374D-1/ DATA U(81) /7.6526521133497334D-2/, V(81) /1.5275338713072585D-1/ C N=24. DATA U(82) /9.9518721999702136D-1/, V(82) /1.2341229799987200D-2/ DATA U(83) /9.7472855597130950D-1/, V(83) /2.8531388628933663D-2/ DATA U(84) /9.3827455200273276D-1/, V(84) /4.4277438817419806D-2/ DATA U(85) /8.8641552700440103D-1/, V(85) /5.9298584915436781D-2/ DATA U(86) /8.2000198597390292D-1/, V(86) /7.3346481411080306D-2/ DATA U(87) /7.4012419157855436D-1/, V(87) /8.6190161531953276D-2/ DATA U(88) /6.4809365193697557D-1/, V(88) /9.7618652104113888D-2/ DATA U(89) /5.4542147138883954D-1/, V(89) /1.0744427011596563D-1/ DATA U(90) /4.3379350762604514D-1/, V(90) /1.1550566805372560D-1/ DATA U(91) /3.1504267969616337D-1/, V(91) /1.2167047292780339D-1/ DATA U(92) /1.9111886747361631D-1/, V(92) /1.2583745634682830D-1/ DATA U(93) /6.4056892862605626D-2/, V(93) /1.2793819534675216D-1/ C N=32. DATA U(94) /9.9726386184948156D-1/, V(94) /7.0186100094700966D-3/ DATA U(95) /9.8561151154526834D-1/, V(95) /1.6274394730905671D-2/ DATA U(96) /9.6476225558750643D-1/, V(96) /2.5392065309262059D-2/ DATA U(97) /9.3490607593773969D-1/, V(97) /3.4273862913021433D-2/ DATA U(98) /8.9632115576605212D-1/, V(98) /4.2835898022226681D-2/ DATA U(99) /8.4936761373256997D-1/, V(99) /5.0998059262376176D-2/ DATA U(100)/7.9448379596794241D-1/, V(100)/5.8684093478535547D-2/ DATA U(101)/7.3218211874028968D-1/, V(101)/6.5822222776361847D-2/ DATA U(102)/6.6304426693021520D-1/, V(102)/7.2345794108848506D-2/ DATA U(103)/5.8771575724076233D-1/, V(103)/7.8193895787070306D-2/ DATA U(104)/5.0689990893222939D-1/, V(104)/8.3311924226946755D-2/ DATA U(105)/4.2135127613063535D-1/, V(105)/8.7652093004403811D-2/ DATA U(106)/3.3186860228212765D-1/, V(106)/9.1173878695763885D-2/ DATA U(107)/2.3928736225213707D-1/, V(107)/9.3844399080804566D-2/ DATA U(108)/1.4447196158279649D-1/, V(108)/9.5638720079274859D-2/ DATA U(109)/4.8307665687738316D-2/, V(109)/9.6540088514727801D-2/ C N=40. DATA U(110)/9.9823770971055920D-1/, V(110)/4.5212770985331913D-3/ DATA U(111)/9.9072623869945701D-1/, V(111)/1.0498284531152814D-2/ DATA U(112)/9.7725994998377426D-1/, V(112)/1.6421058381907889D-2/ DATA U(113)/9.5791681921379166D-1/, V(113)/2.2245849194166957D-2/ DATA U(114)/9.3281280827867653D-1/, V(114)/2.7937006980023401D-2/ DATA U(115)/9.0209880696887430D-1/, V(115)/3.3460195282547847D-2/ DATA U(116)/8.6595950321225950D-1/, V(116)/3.8782167974472018D-2/ DATA U(117)/8.2461223083331166D-1/, V(117)/4.3870908185673272D-2/ DATA U(118)/7.7830565142651939D-1/, V(118)/4.8695807635072232D-2/ DATA U(119)/7.2731825518992710D-1/, V(119)/5.3227846983936824D-2/ DATA U(120)/6.7195668461417955D-1/, V(120)/5.7439769099391551D-2/ DATA U(121)/6.1255388966798024D-1/, V(121)/6.1306242492928939D-2/ DATA U(122)/5.4946712509512820D-1/, V(122)/6.4804013456601038D-2/ DATA U(123)/4.8307580168617871D-1/, V(123)/6.7912045815233904D-2/ DATA U(124)/4.1377920437160500D-1/, V(124)/7.0611647391286780D-2/ DATA U(125)/3.4199409082575847D-1/, V(125)/7.2886582395804059D-2/ DATA U(126)/2.6815218500725368D-1/, V(126)/7.4723169057968264D-2/ DATA U(127)/1.9269758070137110D-1/, V(127)/7.6110361900626242D-2/ DATA U(128)/1.1608407067525521D-1/, V(128)/7.7039818164247966D-2/ DATA U(129)/3.8772417506050822D-2/, V(129)/7.7505947978424811D-2/ C N=48. DATA U(130)/9.9877100725242612D-1/, V(130)/3.1533460523058386D-3/ DATA U(131)/9.9353017226635076D-1/, V(131)/7.3275539012762621D-3/ DATA U(132)/9.8412458372282686D-1/, V(132)/1.1477234579234539D-2/ DATA U(133)/9.7059159254624725D-1/, V(133)/1.5579315722943849D-2/ DATA U(134)/9.5298770316043086D-1/, V(134)/1.9616160457355528D-2/ DATA U(135)/9.3138669070655433D-1/, V(135)/2.3570760839324379D-2/ DATA U(136)/9.0587913671556967D-1/, V(136)/2.7426509708356948D-2/ DATA U(137)/8.7657202027424789D-1/, V(137)/3.1167227832798089D-2/ DATA U(138)/8.4358826162439353D-1/, V(138)/3.4777222564770439D-2/ DATA U(139)/8.0706620402944263D-1/, V(139)/3.8241351065830706D-2/ DATA U(140)/7.6715903251574034D-1/, V(140)/4.1545082943464749D-2/ DATA U(141)/7.2403413092381465D-1/, V(141)/4.4674560856694280D-2/ DATA U(142)/6.7787237963266391D-1/, V(142)/4.7616658492490475D-2/ DATA U(143)/6.2886739677651362D-1/, V(143)/5.0359035553854475D-2/ DATA U(144)/5.7722472608397270D-1/, V(144)/5.2890189485193667D-2/ DATA U(145)/5.2316097472223303D-1/, V(145)/5.5199503699984163D-2/ DATA U(146)/4.6690290475095840D-1/, V(146)/5.7277292100403216D-2/ DATA U(147)/4.0868648199071673D-1/, V(147)/5.9114839698395636D-2/ DATA U(148)/3.4875588629216074D-1/, V(148)/6.0704439165893880D-2/ DATA U(149)/2.8736248735545558D-1/, V(149)/6.2039423159892664D-2/ DATA U(150)/2.2476379039468906D-1/, V(150)/6.3114192286254026D-2/ DATA U(151)/1.6122235606889172D-1/, V(151)/6.3924238584648187D-2/ DATA U(152)/9.7004699209462699D-2/, V(152)/6.4466164435950082D-2/ DATA U(153)/3.2380170962869362D-2/, V(153)/6.4737696812683923D-2/ C N=64. DATA U(154)/9.9930504173577214D-1/, V(154)/1.7832807216964329D-3/ DATA U(155)/9.9634011677195528D-1/, V(155)/4.1470332605624676D-3/ DATA U(156)/9.9101337147674432D-1/, V(156)/6.5044579689783629D-3/ DATA U(157)/9.8333625388462596D-1/, V(157)/8.8467598263639477D-3/ DATA U(158)/9.7332682778991096D-1/, V(158)/1.1168139460131129D-2/ DATA U(159)/9.6100879965205372D-1/, V(159)/1.3463047896718643D-2/ DATA U(160)/9.4641137485840282D-1/, V(160)/1.5726030476024719D-2/ DATA U(161)/9.2956917213193958D-1/, V(161)/1.7951715775697343D-2/ DATA U(162)/9.1052213707850281D-1/, V(162)/2.0134823153530209D-2/ DATA U(163)/8.8931544599511412D-1/, V(163)/2.2270173808383254D-2/ DATA U(164)/8.6599939815409282D-1/, V(164)/2.4352702568710873D-2/ DATA U(165)/8.4062929625258036D-1/, V(165)/2.6377469715054659D-2/ DATA U(166)/8.1326531512279756D-1/, V(166)/2.8339672614259483D-2/ DATA U(167)/7.8397235894334141D-1/, V(167)/3.0234657072402479D-2/ DATA U(168)/7.5281990726053190D-1/, V(168)/3.2057928354851554D-2/ DATA U(169)/7.1988185017161083D-1/, V(169)/3.3805161837141609D-2/ DATA U(170)/6.8523631305423324D-1/, V(170)/3.5472213256882384D-2/ DATA U(171)/6.4896547125465734D-1/, V(171)/3.7055128540240046D-2/ DATA U(172)/6.1115535517239325D-1/, V(172)/3.8550153178615629D-2/ DATA U(173)/5.7189564620263403D-1/, V(173)/3.9953741132720341D-2/ DATA U(174)/5.3127946401989455D-1/, V(174)/4.1262563242623529D-2/ DATA U(175)/4.8940314570705296D-1/, V(175)/4.2473515123653589D-2/ DATA U(176)/4.4636601725346409D-1/, V(176)/4.3583724529323453D-2/ DATA U(177)/4.0227015796399160D-1/, V(177)/4.4590558163756563D-2/ DATA U(178)/3.5722015833766812D-1/, V(178)/4.5491627927418144D-2/ DATA U(179)/3.1132287199021096D-1/, V(179)/4.6284796581314417D-2/ DATA U(180)/2.6468716220876742D-1/, V(180)/4.6968182816210017D-2/ DATA U(181)/2.1742364374000708D-1/, V(181)/4.7540165714830309D-2/ DATA U(182)/1.6964442042399282D-1/, V(182)/4.7999388596458308D-2/ DATA U(183)/1.2146281929612055D-1/, V(183)/4.8344762234802957D-2/ DATA U(184)/7.2993121787799039D-2/, V(184)/4.8575467441503427D-2/ DATA U(185)/2.4350292663424433D-2/, V(185)/4.8690957009139720D-2/ C N=80. DATA U(186)/9.9955382265163063D-1/, V(186)/1.1449500031869415D-3/ DATA U(187)/9.9764986439823769D-1/, V(187)/2.6635335895126817D-3/ DATA U(188)/9.9422754096568828D-1/, V(188)/4.1803131246948952D-3/ DATA U(189)/9.8929130249975553D-1/, V(189)/5.6909224514031986D-3/ DATA U(190)/9.8284857273862907D-1/, V(190)/7.1929047681173128D-3/ DATA U(191)/9.7490914058572779D-1/, V(191)/8.6839452692608584D-3/ DATA U(192)/9.6548508904379925D-1/, V(192)/1.0161766041103065D-2/ DATA U(193)/9.5459076634363491D-1/, V(193)/1.1624114120797827D-2/ DATA U(194)/9.4224276130987267D-1/, V(194)/1.3068761592401339D-2/ DATA U(195)/9.2845987717244580D-1/, V(195)/1.4493508040509076D-2/ DATA U(196)/9.1326310257175765D-1/, V(196)/1.5896183583725688D-2/ DATA U(197)/8.9667557943877068D-1/, V(197)/1.7274652056269306D-2/ DATA U(198)/8.7872256767821383D-1/, V(198)/1.8626814208299031D-2/ DATA U(199)/8.5943140666311110D-1/, V(199)/1.9950610878141999D-2/ DATA U(200)/8.3883147358025528D-1/, V(200)/2.1244026115782006D-2/ DATA U(201)/8.1695413868146347D-1/, V(201)/2.2505090246332462D-2/ DATA U(202)/7.9383271750460545D-1/, V(202)/2.3731882865930101D-2/ DATA U(203)/7.6950242013504137D-1/, V(203)/2.4922535764115491D-2/ DATA U(204)/7.4400029758359727D-1/, V(204)/2.6075235767565118D-2/ DATA U(205)/7.1736518536209988D-1/, V(205)/2.7188227500486381D-2/ DATA U(206)/6.8963764434202760D-1/, V(206)/2.8259816057276862D-2/ DATA U(207)/6.6085989898611980D-1/, V(207)/2.9288369583267848D-2/ DATA U(208)/6.3107577304687197D-1/, V(208)/3.0272321759557981D-2/ DATA U(209)/6.0033062282975174D-1/, V(209)/3.1210174188114702D-2/ DATA U(210)/5.6867126812270978D-1/, V(210)/3.2100498673487773D-2/ DATA U(211)/5.3614592089713193D-1/, V(211)/3.2941939397645401D-2/ DATA U(212)/5.0280411188878499D-1/, V(212)/3.3733214984611523D-2/ DATA U(213)/4.6869661517054448D-1/, V(213)/3.4473120451753929D-2/ DATA U(214)/4.3387537083175609D-1/, V(214)/3.5160529044747593D-2/ DATA U(215)/3.9839340588196923D-1/, V(215)/3.5794393953416055D-2/ DATA U(216)/3.6230475349948732D-1/, V(216)/3.6373749905835978D-2/ DATA U(217)/3.2566437074770191D-1/, V(217)/3.6897714638276009D-2/ DATA U(218)/2.8852805488451185D-1/, V(218)/3.7365490238730490D-2/ DATA U(219)/2.5095235839227212D-1/, V(219)/3.7776364362001397D-2/ DATA U(220)/2.1299450285766613D-1/, V(220)/3.8129711314477638D-2/ DATA U(221)/1.7471229183264681D-1/, V(221)/3.8424993006959423D-2/ DATA U(222)/1.3616402280914389D-1/, V(222)/3.8661759774076463D-2/ DATA U(223)/9.7408398441584599D-2/, V(223)/3.8839651059051969D-2/ DATA U(224)/5.8504437152420669D-2/, V(224)/3.8958395962769531D-2/ DATA U(225)/1.9511383256793998D-2/, V(225)/3.9017813656306655D-2/ C N=96. DATA U(226)/9.9968950388323077D-1/, V(226)/7.9679206555201243D-4/ DATA U(227)/9.9836437586318168D-1/, V(227)/1.8539607889469217D-3/ DATA U(228)/9.9598184298720929D-1/, V(228)/2.9107318179349464D-3/ DATA U(229)/9.9254390032376262D-1/, V(229)/3.9645543384446867D-3/ DATA U(230)/9.8805412632962380D-1/, V(230)/5.0142027429275177D-3/ DATA U(231)/9.8251726356301468D-1/, V(231)/6.0585455042359617D-3/ DATA U(232)/9.7593917458513647D-1/, V(232)/7.0964707911538653D-3/ DATA U(233)/9.6832682846326421D-1/, V(233)/8.1268769256987592D-3/ DATA U(234)/9.5968829144874254D-1/, V(234)/9.1486712307833866D-3/ DATA U(235)/9.5003271778443764D-1/, V(235)/1.0160770535008416D-2/ DATA U(236)/9.3937033975275522D-1/, V(236)/1.1162102099838499D-2/ DATA U(237)/9.2771245672230869D-1/, V(237)/1.2151604671088320D-2/ DATA U(238)/9.1507142312089807D-1/, V(238)/1.3128229566961573D-2/ DATA U(239)/9.0146063531585234D-1/, V(239)/1.4090941772314861D-2/ DATA U(240)/8.8689451740242042D-1/, V(240)/1.5038721026994938D-2/ DATA U(241)/8.7138850590929650D-1/, V(241)/1.5970562902562291D-2/ DATA U(242)/8.5495903343460146D-1/, V(242)/1.6885479864245172D-2/ DATA U(243)/8.3762351122818712D-1/, V(243)/1.7782502316045261D-2/ DATA U(244)/8.1940031073793168D-1/, V(244)/1.8660679627411467D-2/ DATA U(245)/8.0030874413914082D-1/, V(245)/1.9519081140145022D-2/ DATA U(246)/7.8036904386743322D-1/, V(246)/2.0356797154333325D-2/ DATA U(247)/7.5960234117664750D-1/, V(247)/2.1172939892191299D-2/ DATA U(248)/7.3803064374440013D-1/, V(248)/2.1966644438744349D-2/ DATA U(249)/7.1567681234896763D-1/, V(249)/2.2737069658329374D-2/ DATA U(250)/6.9256453664217156D-1/, V(250)/2.3483399085926220D-2/ DATA U(251)/6.6871831004391615D-1/, V(251)/2.4204841792364691D-2/ DATA U(252)/6.4416340378496712D-1/, V(252)/2.4900633222483610D-2/ DATA U(253)/6.1892584012546857D-1/, V(253)/2.5570036005349361D-2/ DATA U(254)/5.9303236477757208D-1/, V(254)/2.6212340735672414D-2/ DATA U(255)/5.6651041856139717D-1/, V(255)/2.6826866725591762D-2/ DATA U(256)/5.3938810832435744D-1/, V(256)/2.7412962726029243D-2/ DATA U(257)/5.1169417715466767D-1/, V(257)/2.7970007616848334D-2/ DATA U(258)/4.8345797392059636D-1/, V(258)/2.8497411065085386D-2/ DATA U(259)/4.5470942216774301D-1/, V(259)/2.8994614150555237D-2/ DATA U(260)/4.2547898840730055D-1/, V(260)/2.9461089958167906D-2/ DATA U(261)/3.9579764982890860D-1/, V(261)/2.9896344136328386D-2/ DATA U(262)/3.6569686147231364D-1/, V(262)/3.0299915420827594D-2/ DATA U(263)/3.3520852289262542D-1/, V(263)/3.0671376123669149D-2/ DATA U(264)/3.0436494435449635D-1/, V(264)/3.1010332586313837D-2/ DATA U(265)/2.7319881259104914D-1/, V(265)/3.1316425596861356D-2/ DATA U(266)/2.4174315616384001D-1/, V(266)/3.1589330770727167D-2/ DATA U(267)/2.1003131046056720D-1/, V(267)/3.1828758894411006D-2/ DATA U(268)/1.7809688236761860D-1/, V(268)/3.2034456231992663D-2/ DATA U(269)/1.4597371465489694D-1/, V(269)/3.2206204794030251D-2/ DATA U(270)/1.1369585011066592D-1/, V(270)/3.2343822568575928D-2/ DATA U(271)/8.1297495464425559D-2/, V(271)/3.2447163714064269D-2/ DATA U(272)/4.8812985136049731D-2/, V(272)/3.2516118713868836D-2/ DATA U(273)/1.6276744849602970D-2/, V(273)/3.2550614492363166D-2/ IF(KTBA(MIN(MAX(1,N),97)) .EQ. 0) THEN X(1)=0 WRITE(ERRTXT,101) N CALL MTLPRT(NAME(MODE),'D107.1',ERRTXT) RETURN ENDIF ALFA=HF*(B+A) BETA=HF*(B-A) IF(MODE .EQ. 1) THEN SUM=0 J1=MOD(N,2) J2=KTBA(N)+(N-1)/2 DO 1 J = KTBA(N),J2-J1 DELTA=BETA*U(J) SUM=SUM+V(J)*(F(ALFA+DELTA)+F(ALFA-DELTA)) 1 CONTINUE IF(J1 .EQ. 1) SUM=SUM+V(J2)*F(ALFA) X(1)=BETA*SUM ELSE J1=KTBA(N)-1 J2=N+1 DO 2 J=1,J2/2 WTEMP=BETA*V(J1+J) DELTA=BETA*U(J1+J) X(J)=ALFA-DELTA W(J)=WTEMP X(J2-J)=ALFA+DELTA W(J2-J)=WTEMP 2 CONTINUE ENDIF RETURN 101 FORMAT('N = ',I5,' IS NON-PERMISSIBLE') END