#include "pdf/pilot.h" SUBROUTINE SFDISA(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,GLU) C THIS IS THE NEW "A" FIT -- May 1994 -- standard Q^2 range C DIS DIS DIS DIS DIS DIS DIS DIS DIS DIS DIS DIS DIS DIS DIS #include "pdf/impdp.inc" PARAMETER(NX=47) PARAMETER(NTENTH=21) DIMENSION F(8,NX,20),G(8),XX(NX),N0(8) DATA ((F(I,1,K),I=1,5),F(7,1,K),F(6,1,K),F(8,1,K),K=1,7)/ + .00477, .00714, 28.15842, 2.63928, .27227, + .00000, 1.40153, 2.63661, + .00513, .00738, 40.95661, 3.03369, .71801, + .00000, 1.82413, 3.03197, + .00546, .00760, 55.69226, 3.48003, 1.20884, + .00000, 2.29496, 3.47918, + .00577, .00781, 71.98843, 3.97328, 1.74137, + .34007, 2.80989, 3.97324, + .00608, .00800, 89.79082, 4.51017, 2.31361, + .97946, 3.36634, 4.51093, + .00637, .00819, 108.92063, 5.08768, 2.92334, + 1.62786, 3.96169, 5.08919, + .00665, .00839, 129.17682, 5.70179, 3.56701, + 2.31175, 4.59222, 5.70402 / DATA ((F(I,1,K),I=1,5),F(7,1,K),F(6,1,K),F(8,1,K),K=8,14)/ + .00693, .00856, 150.39406, 6.34899, 4.24147, + 3.02502, 5.25458, 6.35191, + .00719, .00873, 172.43767, 7.02603, 4.94380, + 3.76499, 5.94570, 7.02963, + .00745, .00890, 195.20667, 7.73017, 5.67150, + 4.52908, 6.66297, 7.73443, + .00770, .00906, 218.62808, 8.45928, 6.42264, + 5.31542, 7.40439, 8.46416, + .00794, .00922, 242.57538, 9.21087, 7.19493, + 6.12166, 8.16754, 9.21635, + .00819, .00936, 266.79108, 9.98022, 7.98383, + 6.94318, 8.94784, 9.98631, + .00841, .00951, 291.48782, 10.76914, 8.79120, + 7.78204, 9.74709, 10.77580 / DATA ((F(I,1,K),I=1,5),F(7,1,K),F(6,1,K),F(8,1,K),K=15,19)/ + .00865, .00965, 316.54001, 11.57541, 9.61495, + 8.63616, 10.56316, 11.58263, + .00887, .00980, 341.62692, 12.39302, 10.44920, + 9.49945, 11.39009, 12.40078, + .00908, .00993, 367.17233, 13.22829, 11.30029, + 10.37867, 12.23424, 13.23658, + .00930, .01007, 392.56421, 14.07036, 12.15750, + 11.26270, 13.08482, 14.07917, + .00951, .01019, 418.32579, 14.92685, 13.02841, + 12.15962, 13.94940, 14.93615 / DATA ((F(I,2,K),I=1,5),F(7,2,K),F(6,2,K),F(8,2,K),K=1,7)/ + .00712, .00919, 22.78315, 2.13445, .21779, + .00000, 1.13005, 2.13406, + .00761, .00951, 32.73079, 2.44280, .56738, + .00000, 1.46077, 2.44298, + .00807, .00980, 43.96273, 2.78701, .94726, + .00000, 1.82436, 2.78773, + .00850, .01008, 56.16534, 3.16268, 1.35441, + .25995, 2.21716, 3.16390, + .00892, .01035, 69.29437, 3.56730, 1.78730, + .74521, 2.63720, 3.56898, + .00932, .01061, 83.20822, 3.99827, 2.24403, + 1.23206, 3.08224, 4.00042, + .00971, .01086, 97.75517, 4.45242, 2.72177, + 1.74079, 3.54929, 4.45499 / DATA ((F(I,2,K),I=1,5),F(7,2,K),F(6,2,K),F(8,2,K),K=8,14)/ + .01008, .01108, 112.81502, 4.92700, 3.21811, + 2.26672, 4.03578, 4.93000, + .01045, .01131, 128.29343, 5.41959, 3.73084, + 2.80787, 4.53940, 5.42298, + .01080, .01153, 144.12270, 5.92812, 4.25813, + 3.36238, 5.05822, 5.93190, + .01115, .01174, 160.25618, 6.45108, 4.79862, + 3.92896, 5.59081, 6.45523, + .01148, .01194, 176.60817, 6.98664, 5.35064, + 4.50594, 6.13542, 6.99115, + .01180, .01215, 193.00063, 7.53130, 5.91082, + 5.08989, 6.68860, 7.53615, + .01212, .01233, 209.59521, 8.08663, 6.48079, + 5.68264, 7.25199, 8.09182 / DATA ((F(I,2,K),I=1,5),F(7,2,K),F(6,2,K),F(8,2,K),K=15,19)/ + .01242, .01252, 226.30714, 8.65105, 7.05909, + 6.28271, 7.82405, 8.65657, + .01273, .01271, 242.91707, 9.22010, 7.64133, + 6.88566, 8.40037, 9.22595, + .01302, .01288, 259.72992, 9.79875, 8.23252, + 7.49675, 8.98593, 9.80490, + .01331, .01306, 276.32343, 10.37895, 8.82469, + 8.10778, 9.57272, 10.38540, + .01358, .01323, 293.06784, 10.96654, 9.42371, + 8.72496, 10.16662, 10.97328 / DATA ((F(I,3,K),I=1,5),F(7,3,K),F(6,3,K),F(8,3,K),K=1,7)/ + .01048, .01184, 18.43334, 1.72552, .17411, + .00000, .91103, 1.72672, + .01116, .01227, 26.15385, 1.96612, .44786, + .00000, 1.16941, 1.96765, + .01180, .01267, 34.69698, 2.23092, .74132, + .00000, 1.44959, 2.23278, + .01240, .01303, 43.80949, 2.51627, 1.05190, + .19835, 1.74850, 2.51842, + .01297, .01339, 53.46115, 2.82028, 1.37855, + .56595, 2.06469, 2.82270, + .01352, .01373, 63.54490, 3.14086, 1.71974, + .93070, 2.39636, 3.14354, + .01406, .01404, 73.95010, 3.47554, 2.07328, + 1.30822, 2.74120, 3.47846 / DATA ((F(I,3,K),I=1,5),F(7,3,K),F(6,3,K),F(8,3,K),K=8,14)/ + .01457, .01435, 84.59291, 3.82226, 2.43737, + 1.69498, 3.09728, 3.82543, + .01506, .01466, 95.41053, 4.17924, 2.81042, + 2.08959, 3.46294, 4.18263, + .01555, .01494, 106.36025, 4.54504, 3.19117, + 2.49080, 3.83679, 4.54864, + .01601, .01522, 117.41489, 4.91856, 3.57867, + 2.89773, 4.21788, 4.92238, + .01647, .01549, 128.51851, 5.29855, 3.97174, + 3.30925, 4.60493, 5.30256, + .01690, .01575, 139.54984, 5.68238, 4.36792, + 3.72285, 4.99544, 5.68659, + .01733, .01600, 150.63226, 6.07148, 4.76868, + 4.14017, 5.39084, 6.07588 / DATA ((F(I,3,K),I=1,5),F(7,3,K),F(6,3,K),F(8,3,K),K=15,19)/ + .01775, .01624, 161.71013, 6.46473, 5.17293, + 4.56016, 5.79004, 6.46930, + .01815, .01648, 172.63414, 6.85887, 5.57754, + 4.97961, 6.18984, 6.86362, + .01855, .01671, 183.62479, 7.25774, 5.98637, + 5.40262, 6.59411, 7.26267, + .01894, .01694, 194.39084, 7.65545, 6.39358, + 5.82317, 6.99696, 7.66054, + .01932, .01715, 205.19538, 8.05647, 6.80367, + 6.24603, 7.40290, 8.06173 / DATA ((F(I,4,K),I=1,5),F(7,4,K),F(6,4,K),F(8,4,K),K=1,7)/ + .01309, .01374, 16.28428, 1.52327, .15268, + .00000, .80308, 1.52520, + .01391, .01425, 22.93523, 1.73109, .38971, + .00000, 1.02647, 1.73327, + .01467, .01472, 30.20590, 1.95792, .64177, + .00000, 1.26675, 1.96034, + .01540, .01515, 37.87578, 2.20047, .90651, + .16916, 1.52116, 2.20310, + .01608, .01556, 45.92312, 2.45723, 1.18317, + .48132, 1.78853, 2.46004, + .01675, .01595, 54.25890, 2.72634, 1.47039, + .78902, 2.06732, 2.72934, + .01739, .01635, 62.79275, 3.00573, 1.76635, + 1.10565, 2.35557, 3.00891 / DATA ((F(I,4,K),I=1,5),F(7,4,K),F(6,4,K),F(8,4,K),K=8,14)/ + .01801, .01671, 71.45830, 3.29369, 2.06956, + 1.42830, 2.65168, 3.29703, + .01860, .01706, 80.20737, 3.58874, 2.37871, + 1.75584, 2.95428, 3.59225, + .01918, .01738, 89.00863, 3.88971, 2.69280, + 2.08729, 3.26227, 3.89338, + .01974, .01772, 97.84350, 4.19579, 3.01111, + 2.42202, 3.57489, 4.19959, + .02029, .01802, 106.66942, 4.50589, 3.33270, + 2.75909, 3.89116, 4.50985, + .02081, .01833, 115.38995, 4.81788, 3.65550, + 3.09646, 4.20893, 4.82198, + .02132, .01863, 124.11079, 5.13308, 3.98090, + 3.43564, 4.52959, 5.13730 / DATA ((F(I,4,K),I=1,5),F(7,4,K),F(6,4,K),F(8,4,K),K=15,19)/ + .02181, .01891, 132.78889, 5.45054, 4.30802, + 3.77580, 4.85223, 5.45490, + .02230, .01919, 141.30566, 5.76762, 4.63425, + 4.11428, 5.17421, 5.77208, + .02277, .01945, 149.84323, 6.08760, 4.96294, + 4.45463, 5.49885, 6.09220, + .02323, .01971, 158.16785, 6.40555, 5.28920, + 4.79183, 5.82127, 6.41026, + .02368, .01997, 166.49508, 6.72534, 5.61691, + 5.12997, 6.14531, 6.73016 / DATA ((F(I,5,K),I=1,5),F(7,5,K),F(6,5,K),F(8,5,K),K=1,7)/ + .01529, .01527, 14.91292, 1.39409, .13906, + .00000, .73429, 1.39648, + .01623, .01585, 20.89327, 1.58127, .35294, + .00000, .93565, 1.58388, + .01711, .01636, 27.37356, 1.78435, .57908, + .00000, 1.15099, 1.78713, + .01793, .01686, 34.15500, 2.00029, .81531, + .15102, 1.37771, 2.00325, + .01871, .01732, 41.22172, 2.22782, 1.06101, + .42883, 1.61490, 2.23094, + .01948, .01777, 48.49636, 2.46528, 1.31501, + .70137, 1.86115, 2.46855, + .02021, .01819, 55.90141, 2.71082, 1.57567, + .98067, 2.11473, 2.71422 / DATA ((F(I,5,K),I=1,5),F(7,5,K),F(6,5,K),F(8,5,K),K=8,14)/ + .02091, .01860, 63.38121, 2.96294, 1.84171, + 1.26417, 2.37427, 2.96649, + .02159, .01899, 70.89649, 3.22038, 2.11201, + 1.55092, 2.63855, 3.22406, + .02224, .01937, 78.42281, 3.48215, 2.38574, + 1.84011, 2.90669, 3.48595, + .02288, .01973, 85.94665, 3.74756, 2.66232, + 2.13124, 3.17803, 3.75147, + .02349, .02008, 93.43315, 4.01569, 2.94090, + 2.42353, 3.45172, 4.01971, + .02409, .02042, 100.80080, 4.28466, 3.21973, + 2.71520, 3.72596, 4.28880, + .02466, .02075, 108.14440, 4.55573, 3.50009, + 3.00767, 4.00197, 4.55997 / DATA ((F(I,5,K),I=1,5),F(7,5,K),F(6,5,K),F(8,5,K),K=15,19)/ + .02523, .02106, 115.42803, 4.82808, 3.78124, + 3.30025, 4.27900, 4.83242, + .02577, .02137, 122.55126, 5.09940, 4.06090, + 3.59061, 4.55476, 5.10383, + .02631, .02167, 129.67320, 5.37266, 4.34209, + 3.88197, 4.83224, 5.37719, + .02683, .02196, 136.59415, 5.64352, 4.62051, + 4.16990, 5.10714, 5.64814, + .02733, .02224, 143.50098, 5.91545, 4.89964, + 4.45808, 5.38290, 5.92015 / DATA ((F(I,6,K),I=1,5),F(7,6,K),F(6,6,K),F(8,6,K),K=1,7)/ + .01724, .01659, 13.92903, 1.30138, .12931, + .00000, .68499, 1.30409, + .01827, .01720, 19.43438, 1.47391, .32677, + .00000, .87072, 1.47679, + .01925, .01778, 25.35865, 1.66018, .53457, + .00000, 1.06840, 1.66324, + .02016, .01832, 31.51905, 1.85742, .75073, + .13825, 1.27568, 1.86062, + .02103, .01882, 37.90410, 2.06447, .97476, + .39196, 1.49172, 2.06781, + .02187, .01930, 44.44489, 2.27983, 1.20553, + .63992, 1.71525, 2.28330, + .02268, .01977, 51.07295, 2.50181, 1.44163, + .89320, 1.94471, 2.50540 / DATA ((F(I,6,K),I=1,5),F(7,6,K),F(6,6,K),F(8,6,K),K=8,14)/ + .02346, .02022, 57.74014, 2.72908, 1.68188, + 1.14949, 2.17887, 2.73280, + .02421, .02064, 64.41331, 2.96053, 1.92532, + 1.40800, 2.41667, 2.96434, + .02493, .02106, 71.07264, 3.19525, 2.17120, + 1.66801, 2.65731, 3.19918, + .02563, .02145, 77.70800, 3.43268, 2.41903, + 1.92912, 2.90027, 3.43671, + .02631, .02183, 84.28979, 3.67201, 2.66813, + 2.19064, 3.14477, 3.67613, + .02696, .02220, 90.74677, 3.91155, 2.91684, + 2.45100, 3.38917, 3.91577, + .02760, .02255, 97.16567, 4.15248, 3.16643, + 2.71154, 3.63468, 4.15678 / DATA ((F(I,6,K),I=1,5),F(7,6,K),F(6,6,K),F(8,6,K),K=15,19)/ + .02821, .02290, 103.51585, 4.39409, 3.41623, + 2.97165, 3.88065, 4.39848, + .02882, .02323, 109.70876, 4.63430, 3.66421, + 3.22926, 4.12497, 4.63876, + .02941, .02356, 115.88792, 4.87584, 3.91314, + 3.48733, 4.37044, 4.88039, + .02998, .02388, 121.87642, 5.11480, 4.15914, + 3.74185, 4.61314, 5.11943, + .03054, .02419, 127.84157, 5.35435, 4.40541, + 3.99622, 4.85625, 5.35905 / DATA ((F(I,7,K),I=1,5),F(7,7,K),F(6,7,K),F(8,7,K),K=1,7)/ + .02487, .02143, 11.26643, 1.04997, .10301, + .00000, .55178, 1.05371, + .02630, .02226, 15.51377, 1.18338, .25659, + .00000, .69583, 1.18726, + .02764, .02301, 19.98321, 1.32529, .41591, + .00000, .84691, 1.32931, + .02889, .02371, 24.53687, 1.47351, .57943, + .10483, 1.00319, 1.47765, + .03008, .02438, 29.17542, 1.62731, .74694, + .29572, 1.16421, 1.63157, + .03123, .02502, 33.85272, 1.78557, .91765, + .48005, 1.32902, 1.78994, + .03234, .02563, 38.52400, 1.94707, 1.09054, + .66639, 1.49653, 1.95154 / DATA ((F(I,7,K),I=1,5),F(7,7,K),F(6,7,K),F(8,7,K),K=8,14)/ + .03340, .02621, 43.16029, 2.11090, 1.26483, + .85312, 1.66586, 2.11546, + .03442, .02677, 47.74363, 2.27629, 1.43990, + 1.03974, 1.83635, 2.28094, + .03540, .02731, 52.26545, 2.44270, 1.61529, + 1.22587, 2.00748, 2.44743, + .03635, .02782, 56.72363, 2.60977, 1.79075, + 1.41131, 2.17896, 2.61458, + .03726, .02832, 61.10136, 2.77695, 1.96579, + 1.59564, 2.35028, 2.78185, + .03815, .02880, 65.35197, 2.94306, 2.13930, + 1.77777, 2.52029, 2.94804, + .03900, .02926, 69.54204, 3.10911, 2.31231, + 1.95884, 2.69001, 3.11416 / DATA ((F(I,7,K),I=1,5),F(7,7,K),F(6,7,K),F(8,7,K),K=15,19)/ + .03984, .02970, 73.65256, 3.27464, 2.48443, + 2.13848, 2.85899, 3.27975, + .04064, .03014, 77.62486, 3.43812, 2.65417, + 2.31520, 3.02577, 3.44331, + .04144, .03056, 81.56219, 3.60171, 2.82369, + 2.49131, 3.19248, 3.60696, + .04220, .03097, 85.34420, 3.76254, 2.99018, + 2.66389, 3.35628, 3.76784, + .04295, .03138, 89.08908, 3.92303, 3.15606, + 2.83554, 3.51960, 3.92838 / DATA ((F(I,8,K),I=1,5),F(7,8,K),F(6,8,K),F(8,8,K),K=1,7)/ + .03563, .02778, 9.10952, .84562, .08174, + .00000, .44408, .85058, + .03760, .02886, 12.36868, .94796, .20049, + .00000, .55507, .95310, + .03943, .02984, 15.71702, 1.05517, .32184, + .00000, .66975, 1.06049, + .04114, .03077, 19.05450, 1.16557, .44468, + .07908, .78671, 1.17105, + .04277, .03164, 22.39167, 1.27879, .56903, + .22192, .90580, 1.28439, + .04433, .03248, 25.70026, 1.39398, .69434, + .35811, 1.02632, 1.39973, + .04584, .03328, 28.95338, 1.51031, .81991, + .49428, 1.14752, 1.51618 / DATA ((F(I,8,K),I=1,5),F(7,8,K),F(6,8,K),F(8,8,K),K=8,14)/ + .04728, .03404, 32.13588, 1.62719, .94529, + .62934, 1.26888, 1.63317, + .04865, .03477, 35.24037, 1.74411, 1.07007, + .76302, 1.38994, 1.75022, + .04998, .03546, 38.26561, 1.86078, 1.19404, + .89517, 1.51045, 1.86699, + .05125, .03612, 41.21448, 1.97699, 1.31705, + 1.02573, 1.63025, 1.98330, + .05248, .03677, 44.07856, 2.09241, 1.43886, + 1.15449, 1.74903, 2.09884, + .05366, .03739, 46.82820, 2.20620, 1.55864, + 1.28068, 1.86598, 2.21272, + .05481, .03799, 49.51406, 2.31923, 1.67733, + 1.40529, 1.98198, 2.32584 / DATA ((F(I,8,K),I=1,5),F(7,8,K),F(6,8,K),F(8,8,K),K=15,19)/ + .05593, .03857, 52.12506, 2.43118, 1.79463, + 1.52809, 2.09674, 2.43787, + .05701, .03913, 54.62273, 2.54098, 1.90950, + 1.64804, 2.20922, 2.54776, + .05806, .03966, 57.08110, 2.65028, 2.02363, + 1.76692, 2.32106, 2.65715, + .05908, .04020, 59.41876, 2.75702, 2.13496, + 1.88262, 2.43022, 2.76396, + .06008, .04071, 61.71886, 2.86303, 2.24536, + 1.99713, 2.53854, 2.87005 / DATA ((F(I,9,K),I=1,5),F(7,9,K),F(6,9,K),F(8,9,K),K=1,7)/ + .04385, .03239, 8.04229, .74421, .07120, + .00000, .39086, .75006, + .04621, .03364, 10.82227, .83137, .17294, + .00000, .48574, .83745, + .04841, .03480, 13.63543, .92186, .27599, + .00000, .58289, .92815, + .05045, .03588, 16.40099, 1.01427, .37944, + .06681, .68112, 1.02073, + .05240, .03690, 19.13413, 1.10833, .48340, + .18688, .78043, 1.11498, + .05427, .03788, 21.81510, 1.20337, .58744, + .30047, .88025, 1.21020, + .05606, .03880, 24.42526, 1.29875, .69103, + .41327, .97998, 1.30574 / DATA ((F(I,9,K),I=1,5),F(7,9,K),F(6,9,K),F(8,9,K),K=8,14)/ + .05778, .03969, 26.95555, 1.39401, .79384, + .52443, 1.07922, 1.40114, + .05941, .04054, 29.40300, 1.48879, .89560, + .63382, 1.17771, 1.49607, + .06099, .04135, 31.76937, 1.58285, .99616, + .74134, 1.27521, 1.59028, + .06250, .04213, 34.05931, 1.67610, 1.09547, + .84703, 1.37167, 1.68366, + .06396, .04288, 36.26794, 1.76830, 1.19333, + .95075, 1.46686, 1.77597, + .06536, .04359, 38.37289, 1.85874, 1.28912, + 1.05189, 1.56014, 1.86655, + .06672, .04429, 40.41710, 1.94823, 1.38363, + 1.15135, 1.65229, 1.95615 / DATA ((F(I,9,K),I=1,5),F(7,9,K),F(6,9,K),F(8,9,K),K=15,19)/ + .06802, .04495, 42.39277, 2.03650, 1.47668, + 1.24895, 1.74309, 2.04455, + .06930, .04560, 44.27015, 2.12273, 1.56741, + 1.34387, 1.83169, 2.13087, + .07053, .04623, 46.10979, 2.20827, 1.65725, + 1.43762, 1.91951, 2.21651, + .07174, .04684, 47.84745, 2.29146, 1.74453, + 1.52848, 2.00487, 2.29981, + .07290, .04742, 49.55035, 2.37383, 1.83080, + 1.61813, 2.08931, 2.38228 / DATA ((F(I,10,K),I=1,5),F(7,10,K),F(6,10,K),F(8,10,K),K=1,7)/ + .05073, .03614, 7.36026, .67927, .06447, + .00000, .35692, .68586, + .05343, .03754, 9.83669, .75687, .15548, + .00000, .44164, .76373, + .05592, .03884, 12.31420, .83689, .24708, + .00000, .52782, .84399, + .05825, .04005, 14.72451, .91809, .33845, + .05913, .61444, .92541, + .06046, .04119, 17.08565, 1.00032, .42978, + .16500, .70150, 1.00785, + .06258, .04227, 19.38313, 1.08297, .52072, + .26457, .78856, 1.09070, + .06460, .04331, 21.60332, 1.16552, .61084, + .36296, .87513, 1.17344 / DATA ((F(I,10,K),I=1,5),F(7,10,K),F(6,10,K),F(8,10,K),K=8,14)/ + .06653, .04430, 23.74072, 1.24759, .69987, + .45946, .96092, 1.25570, + .06838, .04524, 25.79492, 1.32891, .78764, + .55399, 1.04566, 1.33720, + .07015, .04614, 27.76928, 1.40931, .87402, + .64654, 1.12925, 1.41776, + .07185, .04701, 29.66939, 1.48872, .95902, + .73716, 1.21164, 1.49733, + .07349, .04784, 31.49225, 1.56695, 1.04248, + .82576, 1.29265, 1.57572, + .07507, .04863, 33.21981, 1.64343, 1.12387, + .91183, 1.37176, 1.65234, + .07659, .04939, 34.89013, 1.71887, 1.20396, + .99621, 1.44966, 1.72791 / DATA ((F(I,10,K),I=1,5),F(7,10,K),F(6,10,K),F(8,10,K),K=15,19)/ + .07807, .05013, 36.49722, 1.79307, 1.28255, + 1.07877, 1.52620, 1.80225, + .07948, .05085, 38.01646, 1.86529, 1.35894, + 1.15878, 1.60064, 1.87460, + .08086, .05154, 39.50019, 1.93678, 1.43440, + 1.23761, 1.67425, 1.94622, + .08220, .05221, 40.89431, 2.00608, 1.50747, + 1.31377, 1.74556, 2.01564, + .08351, .05287, 42.25639, 2.07455, 1.57955, + 1.38873, 1.81595, 2.08423 / DATA ((F(I,11,K),I=1,5),F(7,11,K),F(6,11,K),F(8,11,K),K=1,7)/ + .05684, .03942, 6.86937, .63259, .05956, + .00000, .33227, .63965, + .05980, .04096, 9.12498, .70329, .14280, + .00000, .40965, .71065, + .06255, .04237, 11.35974, .77580, .22615, + .00000, .48793, .78343, + .06511, .04370, 13.51488, .84900, .30885, + .05360, .56618, .85686, + .06753, .04493, 15.61026, .92279, .39116, + .14928, .64450, .93087, + .06986, .04612, 17.63542, .99663, .47276, + .23882, .72247, 1.00494, + .07208, .04725, 19.58012, 1.07009, .55330, + .32693, .79970, 1.07860 / DATA ((F(I,11,K),I=1,5),F(7,11,K),F(6,11,K),F(8,11,K),K=8,14)/ + .07420, .04832, 21.44131, 1.14286, .63258, + .41301, .87593, 1.15156, + .07622, .04935, 23.22055, 1.21471, .71045, + .49702, .95099, 1.22359, + .07815, .05032, 24.92199, 1.28551, .78685, + .57899, 1.02477, 1.29458, + .08001, .05126, 26.55169, 1.35523, .86180, + .65901, 1.09729, 1.36446, + .08179, .05216, 28.10815, 1.42371, .93518, + .73699, 1.16837, 1.43311, + .08351, .05302, 29.57612, 1.49045, 1.00652, + .81253, 1.23757, 1.50000, + .08516, .05384, 30.99007, 1.55611, 1.07652, + .88638, 1.30554, 1.56581 / DATA ((F(I,11,K),I=1,5),F(7,11,K),F(6,11,K),F(8,11,K),K=15,19)/ + .08676, .05465, 32.34525, 1.62054, 1.14505, + .95844, 1.37216, 1.63037, + .08830, .05542, 33.62056, 1.68308, 1.21150, + 1.02810, 1.43677, 1.69306, + .08981, .05616, 34.86252, 1.74486, 1.27699, + 1.09659, 1.50054, 1.75497, + .09125, .05688, 36.02439, 1.80458, 1.34024, + 1.16258, 1.56215, 1.81482, + .09267, .05758, 37.15666, 1.86348, 1.40251, + 1.22741, 1.62284, 1.87384 / DATA ((F(I,12,K),I=1,5),F(7,12,K),F(6,12,K),F(8,12,K),K=1,7)/ + .08056, .05187, 5.53899, .50541, .04615, + .00000, .26572, .51484, + .08454, .05389, 7.20629, .55765, .10862, + .00000, .32344, .56746, + .08822, .05575, 8.80460, .61029, .17010, + .00000, .38083, .62046, + .09163, .05747, 10.30043, .66253, .23011, + .03913, .43725, .67304, + .09486, .05909, 11.71833, .71443, .28897, + .10827, .49288, .72524, + .09792, .06063, 13.05725, .76563, .34651, + .17194, .54750, .77674, + .10084, .06208, 14.31548, .81589, .40256, + .23372, .60086, .82726 / DATA ((F(I,12,K),I=1,5),F(7,12,K),F(6,12,K),F(8,12,K),K=8,14)/ + .10361, .06345, 15.49562, .86506, .45706, + .29330, .65291, .87669, + .10624, .06475, 16.60275, .91304, .50997, + .35074, .70354, .92492, + .10874, .06600, 17.64306, .95982, .56133, + .40615, .75279, .97193, + .11115, .06718, 18.62338, 1.00541, .61120, + .45967, .80069, 1.01775, + .11345, .06831, 19.54482, 1.04974, .65956, + .51132, .84719, 1.06230, + .11564, .06938, 20.39902, 1.09251, .70610, + .56081, .89199, 1.10525, + .11776, .07042, 21.21094, 1.13422, .75138, + .60879, .93563, 1.14717 / DATA ((F(I,12,K),I=1,5),F(7,12,K),F(6,12,K),F(8,12,K),K=15,19)/ + .11981, .07142, 21.97861, 1.17480, .79535, + .65521, .97803, 1.18794, + .12177, .07238, 22.68923, 1.21384, .83757, + .69967, 1.01878, 1.22714, + .12367, .07330, 23.37444, 1.25213, .87893, + .74308, 1.05873, 1.26562, + .12551, .07419, 24.00451, 1.28880, .91849, + .78453, 1.09697, 1.30245, + .12728, .07505, 24.61300, 1.32474, .95722, + .82500, 1.13441, 1.33856 / DATA ((F(I,13,K),I=1,5),F(7,13,K),F(6,13,K),F(8,13,K),K=1,7)/ + .11354, .06891, 4.45704, .40201, .03499, + .00000, .21161, .41464, + .11880, .07154, 5.65313, .43955, .08085, + .00000, .25367, .45270, + .12364, .07395, 6.75557, .47666, .12516, + .00000, .29469, .49028, + .12810, .07617, 7.75081, .51280, .16762, + .02790, .33427, .52684, + .13227, .07823, 8.66580, .54811, .20861, + .07664, .37266, .56254, + .13622, .08017, 9.50576, .58235, .24803, + .12069, .40972, .59715, + .13995, .08200, 10.27447, .61544, .28588, + .16278, .44541, .63060 / DATA ((F(I,13,K),I=1,5),F(7,13,K),F(6,13,K),F(8,13,K),K=8,14)/ + .14346, .08370, 10.97758, .64735, .32216, + .20275, .47969, .66283, + .14679, .08531, 11.62169, .67807, .35691, + .24076, .51262, .69384, + .14995, .08684, 12.21352, .70762, .39023, + .27696, .54422, .72369, + .15295, .08828, 12.75969, .73608, .42221, + .31150, .57460, .75243, + .15581, .08965, 13.26261, .76341, .45286, + .34443, .60375, .78002, + .15853, .09095, 13.71802, .78946, .48202, + .37563, .63149, .80632, + .16115, .09220, 14.14334, .81462, .51013, + .40558, .65826, .83171 / DATA ((F(I,13,K),I=1,5),F(7,13,K),F(6,13,K),F(8,13,K),K=15,19)/ + .16366, .09339, 14.53804, .83885, .53714, + .43426, .68400, .85616, + .16605, .09452, 14.89501, .86186, .56280, + .46143, .70845, .87938, + .16837, .09560, 15.23480, .88426, .58774, + .48776, .73225, .90199, + .17058, .09665, 15.53910, .90547, .61134, + .51262, .75476, .92339, + .17273, .09765, 15.82952, .92610, .63428, + .53673, .77665, .94421 / DATA ((F(I,14,K),I=1,5),F(7,14,K),F(6,14,K),F(8,14,K),K=1,7)/ + .13834, .08162, 3.91963, .35066, .02931, + -.00001, .18493, .36574, + .14448, .08469, 4.88547, .38108, .06713, + -.00001, .21940, .39675, + .15010, .08747, 5.75208, .41078, .10326, + -.00001, .25262, .42700, + .15524, .09000, 6.51529, .43936, .13748, + .02258, .28429, .45606, + .16004, .09233, 7.20213, .46695, .17014, + .06173, .31467, .48411, + .16455, .09452, 7.82023, .49342, .20125, + .09669, .34369, .51100, + .16879, .09656, 8.37516, .51871, .23080, + .12975, .37132, .53669 / DATA ((F(I,14,K),I=1,5),F(7,14,K),F(6,14,K),F(8,14,K),K=8,14)/ + .17276, .09846, 8.87344, .54287, .25888, + .16084, .39764, .56121, + .17651, .10025, 9.32211, .56590, .28555, + .19013, .42266, .58458, + .18004, .10191, 9.72756, .58785, .31089, + .21778, .44647, .60685, + .18339, .10349, 10.09574, .60880, .33503, + .24396, .46917, .62811, + .18658, .10499, 10.42926, .62879, .35800, + .26874, .49080, .64838, + .18959, .10639, 10.72539, .64764, .37966, + .29201, .51120, .66749, + .19248, .10773, 10.99823, .66571, .40040, + .31421, .53074, .68583 / DATA ((F(I,14,K),I=1,5),F(7,14,K),F(6,14,K),F(8,14,K),K=15,19)/ + .19524, .10901, 11.24775, .68300, .42021, + .33533, .54941, .70335, + .19785, .11021, 11.46847, .69928, .43888, + .35519, .56701, .71986, + .20039, .11137, 11.67636, .71504, .45694, + .37434, .58403, .73585, + .20279, .11246, 11.85815, .72983, .47390, + .39229, .60001, .75084, + .20511, .11353, 12.02992, .74414, .49031, + .40962, .61548, .76535 / DATA ((F(I,15,K),I=1,5),F(7,15,K),F(6,15,K),F(8,15,K),K=1,7)/ + .15889, .09217, 3.57509, .31783, .02560, + -.00001, .16789, .33487, + .16571, .09555, 4.39473, .34380, .05837, + -.00001, .19761, .36150, + .17193, .09861, 5.11463, .36890, .08940, + -.00001, .22600, .38721, + .17758, .10135, 5.73604, .39281, .11854, + .01928, .25279, .41166, + .18283, .10389, 6.28565, .41568, .14612, + .05251, .27827, .43504, + .18774, .10623, 6.77214, .43743, .17218, + .08191, .30240, .45723, + .19233, .10841, 7.20198, .45802, .19674, + .10948, .32517, .47824 / DATA ((F(I,15,K),I=1,5),F(7,15,K),F(6,15,K),F(8,15,K),K=8,14)/ + .19662, .11043, 7.58193, .47754, .21990, + .13522, .34670, .49814, + .20064, .11231, 7.91886, .49598, .24175, + .15931, .36703, .51696, + .20442, .11407, 8.21884, .51343, .26237, + .18188, .38624, .53476, + .20800, .11572, 8.48738, .52998, .28189, + .20313, .40444, .55162, + .21139, .11727, 8.72705, .54566, .30036, + .22312, .42165, .56759, + .21458, .11872, 8.93594, .56034, .31767, + .24177, .43777, .58254, + .21763, .12010, 9.12589, .57432, .33415, + .25947, .45314, .59678 / DATA ((F(I,15,K),I=1,5),F(7,15,K),F(6,15,K),F(8,15,K),K=15,19)/ + .22052, .12141, 9.29704, .58761, .34982, + .27624, .46774, .61033, + .22327, .12263, 9.44518, .60004, .36449, + .29190, .48141, .62298, + .22592, .12381, 9.58329, .61202, .37862, + .30695, .49457, .63518, + .22842, .12492, 9.70083, .62316, .39182, + .32097, .50684, .64654, + .23084, .12598, 9.81077, .63390, .40452, + .33445, .51867, .65748 / DATA ((F(I,16,K),I=1,5),F(7,16,K),F(6,16,K),F(8,16,K),K=1,7)/ + .17665, .10117, 3.32673, .29395, .02286, + .00000, .15565, .31240, + .18406, .10485, 4.04194, .31681, .05213, + .00000, .18208, .33603, + .19076, .10813, 4.65887, .33873, .07965, + -.00001, .20709, .35862, + .19684, .11106, 5.18219, .35943, .10529, + .01696, .23053, .37991, + .20244, .11374, 5.63806, .37907, .12939, + .04607, .25264, .40008, + .20766, .11620, 6.03570, .39757, .15199, + .07162, .27342, .41908, + .21251, .11848, 6.38188, .41497, .17315, + .09542, .29291, .43693 / DATA ((F(I,16,K),I=1,5),F(7,16,K),F(6,16,K),F(8,16,K),K=8,14)/ + .21703, .12057, 6.68365, .43132, .19299, + .11750, .31119, .45369, + .22125, .12251, 6.94738, .44668, .21158, + .13803, .32834, .46943, + .22521, .12431, 7.17876, .46112, .22904, + .15717, .34445, .48422, + .22894, .12600, 7.38304, .47472, .24548, + .17510, .35963, .49816, + .23246, .12757, 7.56283, .48752, .26095, + .19187, .37390, .51126, + .23577, .12904, 7.71659, .49942, .27535, + .20745, .38718, .52344, + .23890, .13044, 7.85442, .51069, .28903, + .22217, .39979, .53499 / DATA ((F(I,16,K),I=1,5),F(7,16,K),F(6,16,K),F(8,16,K),K=15,19)/ + .24189, .13175, 7.97654, .52135, .30196, + .23605, .41170, .54590, + .24470, .13297, 8.07988, .53125, .31401, + .24895, .42278, .55602, + .24741, .13414, 8.17520, .54075, .32557, + .26130, .43342, .56575, + .24996, .13524, 8.25361, .54951, .33630, + .27275, .44327, .57472, + .25243, .13629, 8.32609, .55794, .34660, + .28373, .45273, .58335 / DATA ((F(I,17,K),I=1,5),F(7,17,K),F(6,17,K),F(8,17,K),K=1,7)/ + .24433, .13581, 2.64474, .22926, .01515, + .00000, .12276, .25444, + .25356, .14018, 3.07957, .24369, .03515, + .00000, .14033, .26982, + .26171, .14396, 3.42895, .25709, .05345, + .00000, .15648, .28402, + .26893, .14722, 3.70512, .26930, .07002, + .01097, .17114, .29691, + .27545, .15010, 3.93073, .28052, .08522, + .02948, .18458, .30873, + .28138, .15267, 4.11511, .29074, .09911, + .04528, .19683, .31947, + .28678, .15496, 4.26501, .30004, .11179, + .05965, .20802, .32923 / DATA ((F(I,17,K),I=1,5),F(7,17,K),F(6,17,K),F(8,17,K),K=8,14)/ + .29172, .15701, 4.38661, .30853, .12341, + .07268, .21824, .33812, + .29627, .15887, 4.48505, .31626, .13408, + .08455, .22759, .34620, + .30045, .16055, 4.56467, .32332, .14388, + .09540, .23617, .35359, + .30434, .16210, 4.62914, .32979, .15292, + .10536, .24406, .36036, + .30794, .16350, 4.68024, .33571, .16127, + .11453, .25131, .36654, + .31127, .16477, 4.71772, .34103, .16889, + .12286, .25790, .37209, + .31440, .16594, 4.74756, .34597, .17600, + .13062, .26401, .37724 / DATA ((F(I,17,K),I=1,5),F(7,17,K),F(6,17,K),F(8,17,K),K=15,19)/ + .31732, .16703, 4.77000, .35049, .18259, + .13782, .26967, .38197, + .32002, .16801, 4.78330, .35455, .18862, + .14439, .27480, .38620, + .32261, .16894, 4.79354, .35837, .19432, + .15060, .27965, .39018, + .32499, .16978, 4.79590, .36175, .19950, + .15625, .28402, .39370, + .32726, .17056, 4.79639, .36494, .20441, + .16161, .28814, .39702 / DATA ((F(I,18,K),I=1,5),F(7,18,K),F(6,18,K),F(8,18,K),K=1,7)/ + .33468, .18093, 2.06404, .17371, .00840, + -.00001, .09516, .20873, + .34513, .18513, 2.27574, .18123, .02138, + -.00002, .10573, .21711, + .35397, .18849, 2.42600, .18783, .03280, + -.00002, .11501, .22437, + .36143, .19119, 2.52849, .19350, .04277, + .00649, .12304, .23050, + .36791, .19339, 2.60026, .19839, .05157, + .01722, .13007, .23577, + .37351, .19518, 2.64872, .20254, .05933, + .02607, .13618, .24022, + .37838, .19663, 2.67907, .20608, .06619, + .03387, .14151, .24396 / DATA ((F(I,18,K),I=1,5),F(7,18,K),F(6,18,K),F(8,18,K),K=8,14)/ + .38269, .19784, 2.69554, .20908, .07226, + .04075, .14616, .24712, + .38646, .19883, 2.70144, .21163, .07767, + .04683, .15023, .24977, + .38981, .19963, 2.69946, .21379, .08249, + .05224, .15380, .25200, + .39281, .20030, 2.69167, .21561, .08682, + .05710, .15695, .25387, + .39548, .20084, 2.67899, .21714, .09070, + .06145, .15971, .25542, + .39782, .20124, 2.66139, .21836, .09413, + .06532, .16208, .25663, + .39994, .20157, 2.64186, .21937, .09723, + .06883, .16419, .25763 / DATA ((F(I,18,K),I=1,5),F(7,18,K),F(6,18,K),F(8,18,K),K=15,19)/ + .40185, .20181, 2.62045, .22019, .10006, + .07202, .16605, .25843, + .40352, .20196, 2.59630, .22080, .10254, + .07485, .16763, .25897, + .40505, .20208, 2.57226, .22131, .10485, + .07748, .16906, .25942, + .40637, .20213, 2.54596, .22160, .10685, + .07980, .17025, .25965, + .40760, .20213, 2.52020, .22184, .10871, + .08196, .17133, .25982 / DATA ((F(I,19,K),I=1,5),F(7,19,K),F(6,19,K),F(8,19,K),K=1,7)/ + .39938, .21056, 1.74940, .14184, .00488, + .00000, .07974, .18348, + .40948, .21363, 1.85509, .14594, .01467, + .00000, .08691, .18801, + .41753, .21579, 1.91718, .14930, .02300, + .00001, .09293, .19159, + .42394, .21725, 1.94774, .15193, .03004, + .00449, .09787, .19429, + .42917, .21821, 1.95939, .15400, .03607, + .01176, .10201, .19633, + .43342, .21876, 1.95773, .15553, .04122, + .01761, .10542, .19775, + .43685, .21901, 1.94648, .15665, .04565, + .02265, .10824, .19873 / DATA ((F(I,19,K),I=1,5),F(7,19,K),F(6,19,K),F(8,19,K),K=8,14)/ + .43965, .21903, 1.92844, .15742, .04948, + .02698, .11057, .19933, + .44194, .21889, 1.90574, .15793, .05280, + .03075, .11250, .19960, + .44380, .21861, 1.88002, .15820, .05567, + .03403, .11409, .19963, + .44531, .21822, 1.85251, .15827, .05819, + .03692, .11539, .19949, + .44650, .21776, 1.82368, .15820, .06039, + .03946, .11645, .19915, + .44739, .21718, 1.79338, .15795, .06228, + .04166, .11725, .19863, + .44807, .21658, 1.76337, .15761, .06395, + .04363, .11790, .19803 / DATA ((F(I,19,K),I=1,5),F(7,19,K),F(6,19,K),F(8,19,K),K=15,19)/ + .44856, .21593, 1.73355, .15720, .06542, + .04538, .11841, .19732, + .44883, .21522, 1.70326, .15667, .06666, + .04690, .11875, .19652, + .44902, .21451, 1.67414, .15611, .06780, + .04830, .11901, .19569, + .44901, .21374, 1.64475, .15544, .06874, + .04948, .11913, .19474, + .44894, .21298, 1.61661, .15477, .06961, + .05058, .11920, .19379 / DATA ((F(I,20,K),I=1,5),F(7,20,K),F(6,20,K),F(8,20,K),K=1,7)/ + .45000, .23091, 1.52710, .11805, .00262, + -.00001, .06822, .16385, + .45865, .23234, 1.56946, .12007, .01049, + .00000, .07326, .16573, + .46503, .23290, 1.58186, .12152, .01698, + .00000, .07729, .16685, + .46961, .23282, 1.57417, .12242, .02229, + .00330, .08042, .16736, + .47296, .23235, 1.55566, .12293, .02674, + .00856, .08287, .16739, + .47526, .23153, 1.53020, .12307, .03043, + .01271, .08477, .16704, + .47679, .23047, 1.50033, .12294, .03352, + .01621, .08621, .16638 / DATA ((F(I,20,K),I=1,5),F(7,20,K),F(6,20,K),F(8,20,K),K=8,14)/ + .47770, .22929, 1.46792, .12261, .03613, + .01917, .08729, .16553, + .47815, .22798, 1.43427, .12211, .03833, + .02169, .08807, .16450, + .47823, .22662, 1.40030, .12150, .04019, + .02385, .08863, .16336, + .47799, .22521, 1.36666, .12079, .04177, + .02572, .08901, .16214, + .47751, .22378, 1.33354, .12000, .04312, + .02733, .08923, .16086, + .47678, .22228, 1.30069, .11912, .04423, + .02870, .08928, .15949, + .47590, .22081, 1.26915, .11822, .04520, + .02990, .08924, .15811 / DATA ((F(I,20,K),I=1,5),F(7,20,K),F(6,20,K),F(8,20,K),K=15,19)/ + .47487, .21934, 1.23873, .11730, .04602, + .03094, .08911, .15671, + .47368, .21785, 1.20894, .11633, .04669, + .03184, .08890, .15527, + .47247, .21639, 1.18072, .11536, .04728, + .03265, .08865, .15387, + .47110, .21490, 1.15308, .11435, .04775, + .03332, .08831, .15243, + .46971, .21346, 1.12691, .11337, .04815, + .03393, .08797, .15104 / DATA ((F(I,21,K),I=1,5),F(7,21,K),F(6,21,K),F(8,21,K),K=1,7)/ + .49119, .24467, 1.35090, .09841, .00099, + .00000, .05831, .14647, + .49755, .24418, 1.35189, .09909, .00750, + -.00001, .06187, .14636, + .50153, .24301, 1.33372, .09936, .01269, + -.00001, .06453, .14577, + .50373, .24131, 1.30385, .09924, .01684, + .00250, .06646, .14474, + .50473, .23936, 1.26890, .09886, .02022, + .00644, .06785, .14344, + .50474, .23717, 1.23142, .09824, .02296, + .00948, .06880, .14195, + .50405, .23484, 1.19300, .09746, .02520, + .01201, .06940, .14031 / DATA ((F(I,21,K),I=1,5),F(7,21,K),F(6,21,K),F(8,21,K),K=8,14)/ + .50287, .23248, 1.15472, .09655, .02704, + .01412, .06975, .13860, + .50131, .23010, 1.11726, .09558, .02855, + .01588, .06990, .13683, + .49946, .22775, 1.08108, .09455, .02980, + .01736, .06990, .13503, + .49741, .22542, 1.04646, .09348, .03085, + .01862, .06978, .13322, + .49519, .22311, 1.01335, .09240, .03171, + .01968, .06955, .13143, + .49278, .22080, .98146, .09128, .03239, + .02058, .06923, .12965, + .49031, .21857, .95137, .09016, .03296, + .02136, .06887, .12792 / DATA ((F(I,21,K),I=1,5),F(7,21,K),F(6,21,K),F(8,21,K),K=15,19)/ + .48779, .21639, .92285, .08907, .03343, + .02203, .06846, .12619, + .48515, .21423, .89548, .08797, .03379, + .02256, .06799, .12449, + .48256, .21214, .86981, .08689, .03409, + .02304, .06751, .12284, + .47985, .21007, .84513, .08581, .03431, + .02342, .06698, .12121, + .47721, .20808, .82195, .08476, .03449, + .02377, .06645, .11964 / DATA ((F(I,22,K),I=1,5),F(7,22,K),F(6,22,K),F(8,22,K),K=1,7)/ + .53077, .25425, 1.16857, .07816, .00001, + -.00001, .04809, .12649, + .53340, .25128, 1.13617, .07794, .00499, + -.00001, .05042, .12468, + .53377, .24785, 1.09503, .07741, .00909, + -.00001, .05197, .12267, + .53252, .24422, 1.04985, .07663, .01224, + .00180, .05295, .12047, + .53030, .24048, 1.00464, .07568, .01476, + .00462, .05350, .11823, + .52726, .23672, .96062, .07462, .01673, + .00676, .05374, .11588, + .52369, .23300, .91844, .07347, .01828, + .00851, .05375, .11357 / DATA ((F(I,22,K),I=1,5),F(7,22,K),F(6,22,K),F(8,22,K),K=8,14)/ + .51982, .22934, .87847, .07229, .01951, + .00992, .05360, .11131, + .51572, .22580, .84086, .07109, .02051, + .01107, .05331, .10909, + .51150, .22238, .80565, .06989, .02130, + .01204, .05293, .10691, + .50725, .21905, .77279, .06869, .02193, + .01285, .05249, .10483, + .50295, .21584, .74210, .06751, .02242, + .01352, .05200, .10281, + .49858, .21272, .71321, .06633, .02281, + .01405, .05146, .10084, + .49425, .20972, .68636, .06521, .02310, + .01450, .05090, .09895 / DATA ((F(I,22,K),I=1,5),F(7,22,K),F(6,22,K),F(8,22,K),K=15,19)/ + .48998, .20685, .66128, .06413, .02331, + .01489, .05034, .09711, + .48571, .20403, .63763, .06303, .02347, + .01520, .04975, .09534, + .48157, .20133, .61566, .06198, .02359, + .01546, .04917, .09365, + .47737, .19872, .59484, .06096, .02365, + .01567, .04857, .09198, + .47332, .19622, .57546, .05998, .02368, + .01584, .04798, .09040 / DATA ((F(I,23,K),I=1,5),F(7,23,K),F(6,23,K),F(8,23,K),K=1,7)/ + .55978, .25744, 1.01392, .06150, .00001, + -.00003, .03931, .10793, + .55800, .25193, .96111, .06081, .00327, + -.00003, .04084, .10509, + .55431, .24637, .90733, .05989, .00659, + -.00002, .04170, .10225, + .54938, .24087, .85494, .05880, .00907, + .00132, .04209, .09943, + .54376, .23556, .80589, .05763, .01095, + .00340, .04218, .09671, + .53757, .23041, .76041, .05643, .01239, + .00494, .04204, .09403, + .53116, .22546, .71846, .05520, .01351, + .00617, .04174, .09146 / DATA ((F(I,23,K),I=1,5),F(7,23,K),F(6,23,K),F(8,23,K),K=8,14)/ + .52464, .22077, .67991, .05398, .01437, + .00716, .04132, .08899, + .51813, .21628, .64453, .05278, .01503, + .00794, .04083, .08663, + .51169, .21198, .61211, .05161, .01554, + .00858, .04030, .08439, + .50536, .20790, .58240, .05046, .01592, + .00911, .03973, .08227, + .49911, .20402, .55506, .04936, .01621, + .00954, .03916, .08025, + .49297, .20028, .52976, .04828, .01643, + .00987, .03856, .07830, + .48702, .19672, .50653, .04725, .01657, + .01015, .03796, .07645 / DATA ((F(I,23,K),I=1,5),F(7,23,K),F(6,23,K),F(8,23,K),K=15,19)/ + .48121, .19333, .48508, .04625, .01666, + .01038, .03735, .07471, + .47547, .19005, .46511, .04528, .01669, + .01054, .03675, .07303, + .46996, .18695, .44670, .04437, .01672, + .01068, .03617, .07145, + .46453, .18394, .42945, .04347, .01670, + .01079, .03559, .06990, + .45929, .18108, .41349, .04262, .01666, + .01088, .03502, .06845 / DATA ((F(I,24,K),I=1,5),F(7,24,K),F(6,24,K),F(8,24,K),K=1,7)/ + .57956, .25538, .87934, .04769, .00001, + .00000, .03179, .09098, + .57288, .24758, .81474, .04690, .00204, + -.00001, .03274, .08754, + .56485, .24017, .75488, .04592, .00472, + -.00001, .03316, .08432, + .55611, .23310, .70015, .04481, .00667, + .00099, .03322, .08123, + .54706, .22650, .65100, .04367, .00812, + .00251, .03302, .07831, + .53788, .22021, .60687, .04250, .00921, + .00361, .03268, .07556, + .52874, .21433, .56723, .04134, .01000, + .00449, .03222, .07298 / DATA ((F(I,24,K),I=1,5),F(7,24,K),F(6,24,K),F(8,24,K),K=8,14)/ + .51977, .20878, .53158, .04022, .01061, + .00518, .03169, .07056, + .51106, .20354, .49950, .03911, .01104, + .00571, .03114, .06830, + .50259, .19865, .47056, .03805, .01137, + .00614, .03057, .06617, + .49439, .19405, .44440, .03706, .01162, + .00649, .02998, .06414, + .48650, .18967, .42069, .03611, .01180, + .00677, .02938, .06223, + .47880, .18550, .39898, .03516, .01189, + .00697, .02880, .06047, + .47142, .18158, .37923, .03427, .01195, + .00714, .02824, .05878 / DATA ((F(I,24,K),I=1,5),F(7,24,K),F(6,24,K),F(8,24,K),K=15,19)/ + .46428, .17787, .36118, .03343, .01198, + .00728, .02768, .05718, + .45738, .17430, .34454, .03261, .01196, + .00737, .02713, .05569, + .45080, .17090, .32930, .03183, .01192, + .00744, .02660, .05427, + .44433, .16767, .31518, .03110, .01188, + .00748, .02607, .05292, + .43815, .16460, .30217, .03039, .01182, + .00752, .02556, .05163 / DATA ((F(I,25,K),I=1,5),F(7,25,K),F(6,25,K),F(8,25,K),K=1,7)/ + .59083, .24926, .76061, .03650, .00001, + -.00002, .02534, .07556, + .57907, .23940, .69015, .03572, .00119, + -.00002, .02596, .07196, + .56672, .23035, .62849, .03480, .00337, + -.00003, .02613, .06862, + .55425, .22200, .57445, .03378, .00493, + .00073, .02599, .06554, + .54201, .21428, .52734, .03272, .00604, + .00187, .02567, .06274, + .52999, .20716, .48607, .03168, .00687, + .00268, .02524, .06012, + .51836, .20055, .44979, .03065, .00746, + .00332, .02474, .05769 / DATA ((F(I,25,K),I=1,5),F(7,25,K),F(6,25,K),F(8,25,K),K=8,14)/ + .50721, .19442, .41774, .02969, .00788, + .00380, .02422, .05542, + .49653, .18870, .38931, .02876, .00819, + .00417, .02367, .05332, + .48632, .18336, .36400, .02785, .00840, + .00446, .02311, .05140, + .47658, .17838, .34139, .02700, .00853, + .00469, .02256, .04961, + .46729, .17369, .32109, .02618, .00863, + .00487, .02202, .04795, + .45833, .16929, .30275, .02541, .00867, + .00500, .02150, .04634, + .44979, .16516, .28620, .02468, .00868, + .00510, .02099, .04486 / DATA ((F(I,25,K),I=1,5),F(7,25,K),F(6,25,K),F(8,25,K),K=15,19)/ + .44164, .16125, .27120, .02399, .00868, + .00517, .02048, .04349, + .43376, .15755, .25747, .02334, .00865, + .00523, .02001, .04217, + .42628, .15406, .24498, .02273, .00861, + .00527, .01956, .04093, + .41906, .15071, .23347, .02210, .00854, + .00528, .01911, .03980, + .41219, .14753, .22294, .02153, .00846, + .00528, .01868, .03874 / DATA ((F(I,26,K),I=1,5),F(7,26,K),F(6,26,K),F(8,26,K),K=1,7)/ + .59434, .23984, .65531, .02744, .00001, + -.00002, .01993, .06187, + .57754, .22821, .58311, .02679, .00060, + -.00001, .02036, .05835, + .56109, .21779, .52246, .02595, .00241, + -.00002, .02038, .05517, + .54509, .20846, .47102, .02511, .00367, + .00056, .02014, .05227, + .52986, .19998, .42728, .02422, .00453, + .00142, .01980, .04964, + .51529, .19225, .38971, .02335, .00514, + .00201, .01937, .04724, + .50151, .18514, .35721, .02251, .00557, + .00246, .01889, .04506 / DATA ((F(I,26,K),I=1,5),F(7,26,K),F(6,26,K),F(8,26,K),K=8,14)/ + .48847, .17863, .32891, .02168, .00589, + .00280, .01837, .04308, + .47616, .17260, .30416, .02091, .00609, + .00307, .01788, .04124, + .46455, .16700, .28238, .02016, .00623, + .00326, .01739, .03954, + .45357, .16183, .26312, .01948, .00632, + .00341, .01689, .03797, + .44309, .15706, .24601, .01884, .00637, + .00353, .01644, .03649, + .43316, .15255, .23069, .01820, .00636, + .00361, .01597, .03518, + .42377, .14832, .21697, .01762, .00635, + .00367, .01554, .03394 / DATA ((F(I,26,K),I=1,5),F(7,26,K),F(6,26,K),F(8,26,K),K=15,19)/ + .41487, .14436, .20461, .01707, .00632, + .00372, .01511, .03278, + .40629, .14063, .19339, .01655, .00628, + .00374, .01471, .03169, + .39819, .13714, .18323, .01606, .00623, + .00376, .01432, .03067, + .39039, .13380, .17395, .01561, .00617, + .00377, .01396, .02972, + .38300, .13066, .16549, .01517, .00610, + .00377, .01361, .02883 / DATA ((F(I,27,K),I=1,5),F(7,27,K),F(6,27,K),F(8,27,K),K=1,7)/ + .59081, .22768, .56179, .02013, .00001, + .00002, .01554, .04994, + .56917, .21473, .49090, .01969, .00027, + .00002, .01577, .04665, + .54874, .20343, .43319, .01913, .00172, + .00003, .01569, .04375, + .52961, .19340, .38555, .01843, .00267, + .00046, .01547, .04113, + .51178, .18442, .34581, .01775, .00336, + .00106, .01513, .03878, + .49505, .17634, .31224, .01706, .00382, + .00150, .01468, .03670, + .47948, .16899, .28363, .01636, .00414, + .00181, .01424, .03480 / DATA ((F(I,27,K),I=1,5),F(7,27,K),F(6,27,K),F(8,27,K),K=8,14)/ + .46493, .16226, .25905, .01571, .00434, + .00206, .01380, .03309, + .45131, .15614, .23781, .01509, .00448, + .00224, .01338, .03151, + .43854, .15051, .21931, .01453, .00457, + .00237, .01295, .03007, + .42656, .14532, .20310, .01399, .00463, + .00247, .01253, .02876, + .41532, .14052, .18885, .01349, .00465, + .00255, .01214, .02755, + .40468, .13602, .17617, .01299, .00464, + .00259, .01175, .02645, + .39466, .13185, .16489, .01255, .00461, + .00262, .01139, .02542 / DATA ((F(I,27,K),I=1,5),F(7,27,K),F(6,27,K),F(8,27,K),K=15,19)/ + .38518, .12796, .15481, .01215, .00458, + .00265, .01106, .02444, + .37621, .12429, .14571, .01172, .00454, + .00265, .01072, .02356, + .36776, .12085, .13751, .01134, .00448, + .00266, .01041, .02275, + .35964, .11764, .13007, .01098, .00443, + .00266, .01010, .02195, + .35197, .11461, .12332, .01065, .00437, + .00265, .00983, .02123 / DATA ((F(I,28,K),I=1,5),F(7,28,K),F(6,28,K),F(8,28,K),K=1,7)/ + .58058, .21431, .47915, .01495, .00001, + -.00014, .01159, .03971, + .55451, .20028, .41142, .01462, .00001, + -.00013, .01187, .03674, + .53072, .18822, .35785, .01409, .00100, + -.00012, .01182, .03423, + .50889, .17769, .31456, .01356, .00183, + .00021, .01160, .03193, + .48888, .16841, .27910, .01299, .00239, + .00069, .01128, .02994, + .47044, .16014, .24960, .01245, .00275, + .00102, .01095, .02811, + .45345, .15270, .22480, .01190, .00299, + .00127, .01060, .02652 / DATA ((F(I,28,K),I=1,5),F(7,28,K),F(6,28,K),F(8,28,K),K=8,14)/ + .43774, .14598, .20377, .01136, .00314, + .00143, .01022, .02512, + .42317, .13986, .18575, .01088, .00326, + .00158, .00986, .02384, + .40964, .13429, .17022, .01043, .00332, + .00168, .00951, .02263, + .39703, .12919, .15674, .01000, .00335, + .00174, .00918, .02154, + .38523, .12446, .14494, .00963, .00337, + .00179, .00885, .02056, + .37417, .12009, .13454, .00925, .00334, + .00183, .00855, .01964, + .36383, .11606, .12534, .00890, .00332, + .00185, .00826, .01880 / DATA ((F(I,28,K),I=1,5),F(7,28,K),F(6,28,K),F(8,28,K),K=15,19)/ + .35412, .11230, .11716, .00856, .00328, + .00186, .00799, .01804, + .34490, .10878, .10984, .00826, .00325, + .00186, .00773, .01733, + .33626, .10550, .10327, .00795, .00321, + .00186, .00749, .01668, + .32802, .10243, .09735, .00768, .00316, + .00185, .00725, .01606, + .32024, .09954, .09200, .00743, .00311, + .00185, .00703, .01549 / DATA ((F(I,29,K),I=1,5),F(7,29,K),F(6,29,K),F(8,29,K),K=1,7)/ + .56485, .19887, .40568, .01053, .00001, + -.00002, .00877, .03101, + .53479, .18424, .34266, .01035, .00001, + -.00003, .00897, .02852, + .50803, .17183, .29389, .01003, .00072, + -.00003, .00890, .02631, + .48390, .16119, .25527, .00964, .00134, + .00021, .00870, .02445, + .46220, .15188, .22417, .00921, .00175, + .00055, .00845, .02277, + .44244, .14363, .19863, .00880, .00201, + .00078, .00815, .02131, + .42438, .13633, .17743, .00841, .00220, + .00095, .00784, .01998 / DATA ((F(I,29,K),I=1,5),F(7,29,K),F(6,29,K),F(8,29,K),K=8,14)/ + .40784, .12976, .15963, .00805, .00233, + .00108, .00754, .01877, + .39262, .12384, .14454, .00769, .00241, + .00117, .00724, .01771, + .37861, .11845, .13166, .00735, .00245, + .00123, .00696, .01677, + .36565, .11352, .12056, .00702, .00246, + .00127, .00668, .01592, + .35363, .10901, .11090, .00671, .00244, + .00130, .00643, .01515, + .34238, .10485, .10245, .00644, .00242, + .00132, .00619, .01442, + .33190, .10104, .09503, .00618, .00240, + .00133, .00596, .01375 / DATA ((F(I,29,K),I=1,5),F(7,29,K),F(6,29,K),F(8,29,K),K=15,19)/ + .32209, .09751, .08848, .00594, .00236, + .00133, .00575, .01314, + .31282, .09418, .08264, .00574, .00233, + .00133, .00554, .01259, + .30416, .09110, .07742, .00556, .00230, + .00133, .00535, .01209, + .29599, .08822, .07275, .00535, .00225, + .00132, .00516, .01161, + .28832, .08554, .06854, .00515, .00221, + .00131, .00500, .01116 / DATA ((F(I,30,K),I=1,5),F(7,30,K),F(6,30,K),F(8,30,K),K=1,7)/ + .54401, .18285, .34122, .00732, .00001, + -.00004, .00641, .02389, + .51054, .16789, .28359, .00732, .00001, + -.00003, .00660, .02177, + .48139, .15539, .24004, .00709, .00042, + -.00004, .00656, .02002, + .45559, .14480, .20611, .00680, .00093, + .00014, .00639, .01846, + .43257, .13565, .17915, .00650, .00126, + .00040, .00619, .01707, + .41191, .12763, .15733, .00620, .00148, + .00058, .00595, .01589, + .39322, .12055, .13943, .00592, .00160, + .00070, .00569, .01485 / DATA ((F(I,30,K),I=1,5),F(7,30,K),F(6,30,K),F(8,30,K),K=8,14)/ + .37627, .11423, .12454, .00560, .00168, + .00078, .00547, .01392, + .36077, .10856, .11204, .00534, .00172, + .00083, .00524, .01307, + .34656, .10345, .10145, .00511, .00175, + .00087, .00502, .01230, + .33352, .09881, .09239, .00490, .00176, + .00091, .00481, .01160, + .32145, .09457, .08456, .00468, .00175, + .00093, .00462, .01100, + .31025, .09069, .07777, .00448, .00174, + .00094, .00442, .01045, + .29986, .08710, .07184, .00429, .00171, + .00095, .00425, .00994 / DATA ((F(I,30,K),I=1,5),F(7,30,K),F(6,30,K),F(8,30,K),K=15,19)/ + .29018, .08380, .06662, .00412, .00169, + .00095, .00409, .00950, + .28112, .08074, .06200, .00393, .00166, + .00094, .00392, .00907, + .27268, .07790, .05789, .00377, .00163, + .00093, .00378, .00867, + .26465, .07526, .05424, .00364, .00160, + .00092, .00365, .00830, + .25716, .07280, .05095, .00353, .00157, + .00091, .00353, .00796 / DATA ((F(I,31,K),I=1,5),F(7,31,K),F(6,31,K),F(8,31,K),K=1,7)/ + .51886, .16623, .28475, .00492, .00001, + .00000, .00462, .01809, + .48273, .15131, .23299, .00495, .00001, + .00000, .00479, .01643, + .45168, .13904, .19461, .00489, .00026, + .00000, .00478, .01493, + .42463, .12874, .16529, .00473, .00064, + .00013, .00466, .01369, + .40087, .11986, .14226, .00450, .00091, + .00032, .00447, .01263, + .37965, .11223, .12386, .00432, .00105, + .00043, .00430, .01166, + .36075, .10549, .10893, .00409, .00116, + .00051, .00411, .01084 / DATA ((F(I,31,K),I=1,5),F(7,31,K),F(6,31,K),F(8,31,K),K=8,14)/ + .34365, .09953, .09660, .00392, .00122, + .00058, .00393, .01010, + .32820, .09421, .08636, .00372, .00124, + .00062, .00375, .00945, + .31411, .08943, .07775, .00353, .00125, + .00064, .00358, .00887, + .30120, .08513, .07042, .00336, .00127, + .00066, .00342, .00836, + .28935, .08120, .06417, .00320, .00125, + .00067, .00327, .00790, + .27837, .07762, .05875, .00306, .00125, + .00067, .00313, .00748, + .26824, .07434, .05404, .00294, .00123, + .00067, .00299, .00709 / DATA ((F(I,31,K),I=1,5),F(7,31,K),F(6,31,K),F(8,31,K),K=15,19)/ + .25887, .07133, .04993, .00280, .00120, + .00068, .00286, .00675, + .25010, .06852, .04631, .00269, .00118, + .00067, .00276, .00643, + .24196, .06595, .04311, .00258, .00116, + .00067, .00266, .00613, + .23429, .06356, .04025, .00248, .00112, + .00065, .00254, .00585, + .22714, .06134, .03770, .00239, .00110, + .00064, .00245, .00560 / DATA ((F(I,32,K),I=1,5),F(7,32,K),F(6,32,K),F(8,32,K),K=1,7)/ + .49000, .14947, .23556, .00306, .00001, + .00004, .00331, .01339, + .45183, .13488, .18990, .00323, .00001, + .00005, .00346, .01213, + .41960, .12303, .15658, .00322, .00017, + .00004, .00345, .01099, + .39197, .11315, .13151, .00310, .00049, + .00013, .00336, .01001, + .36781, .10477, .11214, .00300, .00068, + .00026, .00324, .00919, + .34661, .09759, .09679, .00288, .00080, + .00034, .00310, .00842, + .32774, .09131, .08445, .00273, .00084, + .00040, .00295, .00781 / DATA ((F(I,32,K),I=1,5),F(7,32,K),F(6,32,K),F(8,32,K),K=8,14)/ + .31088, .08578, .07441, .00257, .00088, + .00043, .00280, .00728, + .29565, .08087, .06611, .00246, .00092, + .00046, .00266, .00678, + .28184, .07647, .05918, .00236, .00091, + .00048, .00253, .00633, + .26929, .07254, .05334, .00226, .00090, + .00049, .00242, .00594, + .25787, .06896, .04835, .00211, .00091, + .00049, .00229, .00561, + .24732, .06571, .04409, .00202, .00088, + .00049, .00220, .00528, + .23763, .06276, .04041, .00195, .00087, + .00048, .00210, .00499 / DATA ((F(I,32,K),I=1,5),F(7,32,K),F(6,32,K),F(8,32,K),K=15,19)/ + .22867, .06006, .03720, .00187, .00086, + .00048, .00201, .00472, + .22033, .05757, .03438, .00178, .00084, + .00047, .00192, .00447, + .21259, .05528, .03190, .00172, .00082, + .00046, .00184, .00425, + .20537, .05312, .02970, .00165, .00079, + .00046, .00177, .00405, + .19866, .05115, .02775, .00157, .00078, + .00046, .00169, .00386 / DATA ((F(I,33,K),I=1,5),F(7,33,K),F(6,33,K),F(8,33,K),K=1,7)/ + .45824, .13324, .19331, .00198, .00001, + .00001, .00221, .00988, + .41888, .11919, .15350, .00211, .00001, + .00000, .00238, .00877, + .38607, .10791, .12501, .00212, .00006, + .00001, .00239, .00795, + .35814, .09860, .10386, .00212, .00030, + .00008, .00236, .00721, + .33417, .09079, .08771, .00204, .00045, + .00017, .00227, .00653, + .31325, .08408, .07507, .00196, .00054, + .00023, .00216, .00599, + .29486, .07828, .06503, .00183, .00060, + .00027, .00206, .00552 / DATA ((F(I,33,K),I=1,5),F(7,33,K),F(6,33,K),F(8,33,K),K=8,14)/ + .27839, .07322, .05691, .00176, .00063, + .00030, .00194, .00512, + .26372, .06875, .05027, .00167, .00063, + .00032, .00185, .00476, + .25048, .06481, .04476, .00160, .00064, + .00033, .00176, .00441, + .23847, .06125, .04014, .00153, .00064, + .00033, .00166, .00412, + .22755, .05803, .03623, .00146, .00061, + .00032, .00159, .00387, + .21758, .05514, .03290, .00139, .00062, + .00033, .00151, .00363, + .20844, .05250, .03003, .00131, .00061, + .00033, .00144, .00343 / DATA ((F(I,33,K),I=1,5),F(7,33,K),F(6,33,K),F(8,33,K),K=15,19)/ + .20003, .05009, .02754, .00124, .00059, + .00033, .00135, .00326, + .19224, .04789, .02537, .00117, .00057, + .00032, .00130, .00308, + .18505, .04586, .02347, .00112, .00055, + .00032, .00125, .00291, + .17834, .04396, .02179, .00108, .00055, + .00032, .00120, .00277, + .17210, .04222, .02030, .00104, .00054, + .00032, .00116, .00264 / DATA ((F(I,34,K),I=1,5),F(7,34,K),F(6,34,K),F(8,34,K),K=1,7)/ + .42438, .11742, .15704, .00116, .00001, + .00003, .00148, .00700, + .38440, .10416, .12291, .00136, .00001, + .00003, .00163, .00627, + .35168, .09361, .09888, .00146, .00001, + .00002, .00165, .00558, + .32418, .08501, .08131, .00142, .00018, + .00007, .00160, .00507, + .30066, .07781, .06805, .00137, .00027, + .00012, .00153, .00461, + .28036, .07173, .05779, .00133, .00034, + .00016, .00146, .00421, + .26256, .06649, .04969, .00129, .00038, + .00018, .00139, .00386 / DATA ((F(I,34,K),I=1,5),F(7,34,K),F(6,34,K),F(8,34,K),K=8,14)/ + .24687, .06192, .04320, .00122, .00040, + .00020, .00131, .00356, + .23290, .05790, .03793, .00115, .00041, + .00021, .00124, .00328, + .22041, .05436, .03359, .00108, .00042, + .00022, .00117, .00305, + .20915, .05121, .02997, .00103, .00042, + .00022, .00112, .00283, + .19889, .04837, .02693, .00098, .00041, + .00022, .00106, .00266, + .18956, .04582, .02435, .00096, .00039, + .00021, .00099, .00249, + .18107, .04350, .02215, .00090, .00038, + .00021, .00094, .00234 / DATA ((F(I,34,K),I=1,5),F(7,34,K),F(6,34,K),F(8,34,K),K=15,19)/ + .17329, .04139, .02025, .00084, .00038, + .00022, .00090, .00223, + .16607, .03947, .01860, .00083, .00038, + .00021, .00086, .00209, + .15941, .03770, .01715, .00080, .00038, + .00021, .00082, .00196, + .15326, .03606, .01588, .00075, .00036, + .00020, .00078, .00187, + .14758, .03456, .01476, .00070, .00035, + .00020, .00075, .00178 / DATA ((F(I,35,K),I=1,5),F(7,35,K),F(6,35,K),F(8,35,K),K=1,7)/ + .38889, .10222, .12614, .00054, .00001, + .00004, .00095, .00490, + .34910, .08989, .09729, .00076, .00001, + .00004, .00108, .00438, + .31692, .08021, .07732, .00087, .00001, + .00003, .00112, .00389, + .29023, .07233, .06288, .00085, .00013, + .00006, .00109, .00352, + .26753, .06583, .05212, .00087, .00020, + .00010, .00105, .00316, + .24816, .06036, .04389, .00082, .00025, + .00013, .00100, .00288, + .23136, .05569, .03747, .00077, .00028, + .00014, .00095, .00262 / DATA ((F(I,35,K),I=1,5),F(7,35,K),F(6,35,K),F(8,35,K),K=8,14)/ + .21652, .05163, .03236, .00075, .00031, + .00015, .00091, .00240, + .20342, .04809, .02824, .00072, .00032, + .00016, .00085, .00223, + .19174, .04497, .02487, .00068, .00032, + .00016, .00080, .00206, + .18128, .04221, .02209, .00064, .00031, + .00016, .00075, .00191, + .17188, .03974, .01976, .00060, .00029, + .00016, .00072, .00179, + .16329, .03751, .01780, .00058, .00029, + .00016, .00069, .00167, + .15548, .03551, .01613, .00056, .00028, + .00016, .00064, .00156 / DATA ((F(I,35,K),I=1,5),F(7,35,K),F(6,35,K),F(8,35,K),K=15,19)/ + .14834, .03370, .01469, .00053, .00027, + .00016, .00060, .00147, + .14178, .03203, .01345, .00051, .00026, + .00015, .00058, .00139, + .13577, .03051, .01237, .00050, .00026, + .00014, .00055, .00130, + .13019, .02912, .01143, .00046, .00025, + .00014, .00053, .00124, + .12506, .02784, .01059, .00044, .00024, + .00014, .00049, .00117 / DATA ((F(I,36,K),I=1,5),F(7,36,K),F(6,36,K),F(8,36,K),K=1,7)/ + .35261, .08798, .10021, .00018, .00001, + .00003, .00061, .00334, + .31369, .07665, .07613, .00034, .00001, + .00003, .00073, .00296, + .28250, .06781, .05972, .00043, .00000, + .00001, .00077, .00266, + .25690, .06075, .04805, .00045, .00012, + .00005, .00077, .00236, + .23540, .05494, .03946, .00046, .00017, + .00008, .00074, .00213, + .21714, .05011, .03296, .00046, .00021, + .00010, .00071, .00190, + .20138, .04601, .02792, .00046, .00022, + .00011, .00065, .00174 / DATA ((F(I,36,K),I=1,5),F(7,36,K),F(6,36,K),F(8,36,K),K=8,14)/ + .18767, .04246, .02396, .00041, .00023, + .00012, .00064, .00159, + .17558, .03937, .02078, .00038, .00023, + .00013, .00060, .00148, + .16485, .03668, .01821, .00035, .00023, + .00013, .00057, .00137, + .15529, .03431, .01609, .00034, .00022, + .00012, .00054, .00127, + .14672, .03219, .01433, .00033, .00022, + .00013, .00049, .00115, + .13892, .03030, .01286, .00033, .00022, + .00012, .00048, .00107, + .13189, .02859, .01161, .00031, .00022, + .00012, .00045, .00102 / DATA ((F(I,36,K),I=1,5),F(7,36,K),F(6,36,K),F(8,36,K),K=15,19)/ + .12552, .02703, .01054, .00027, .00020, + .00012, .00044, .00095, + .11964, .02563, .00962, .00026, .00021, + .00012, .00040, .00089, + .11427, .02434, .00882, .00025, .00020, + .00011, .00038, .00083, + .10927, .02319, .00813, .00025, .00019, + .00011, .00037, .00080, + .10468, .02213, .00752, .00024, .00018, + .00011, .00036, .00076 / DATA ((F(I,37,K),I=1,5),F(7,37,K),F(6,37,K),F(8,37,K),K=1,7)/ + .31655, .07516, .07877, .00019, .00001, + -.00002, .00027, .00221, + .27901, .06489, .05899, .00032, .00001, + -.00001, .00037, .00201, + .24943, .05701, .04574, .00036, .00001, + -.00002, .00044, .00174, + .22521, .05074, .03644, .00040, .00001, + .00001, .00043, .00156, + .20515, .04561, .02966, .00039, .00007, + .00003, .00044, .00140, + .18824, .04138, .02458, .00036, .00009, + .00004, .00042, .00127, + .17375, .03780, .02069, .00035, .00010, + .00005, .00039, .00115 / DATA ((F(I,37,K),I=1,5),F(7,37,K),F(6,37,K),F(8,37,K),K=8,14)/ + .16119, .03473, .01764, .00034, .00011, + .00006, .00037, .00105, + .15022, .03208, .01522, .00032, .00011, + .00006, .00035, .00095, + .14049, .02977, .01327, .00028, .00012, + .00006, .00034, .00088, + .13185, .02773, .01167, .00028, .00013, + .00006, .00032, .00081, + .12412, .02595, .01036, .00027, .00012, + .00007, .00030, .00076, + .11718, .02433, .00926, .00025, .00011, + .00006, .00027, .00070, + .11090, .02290, .00833, .00024, .00010, + .00005, .00027, .00065 / DATA ((F(I,37,K),I=1,5),F(7,37,K),F(6,37,K),F(8,37,K),K=15,19)/ + .10520, .02161, .00754, .00023, .00010, + .00005, .00025, .00061, + .10001, .02043, .00686, .00023, .00010, + .00005, .00024, .00056, + .09528, .01936, .00628, .00022, .00010, + .00005, .00023, .00052, + .09092, .01839, .00577, .00019, .00010, + .00005, .00021, .00050, + .08691, .01750, .00532, .00018, .00010, + .00005, .00020, .00048 / DATA ((F(I,38,K),I=1,5),F(7,38,K),F(6,38,K),F(8,38,K),K=1,7)/ + .28032, .06291, .06086, .00001, .00001, + .00002, .00018, .00141, + .24487, .05384, .04492, .00006, .00001, + .00002, .00027, .00127, + .21708, .04694, .03440, .00018, .00001, + .00002, .00031, .00108, + .19475, .04150, .02713, .00018, .00002, + .00003, .00031, .00100, + .17628, .03710, .02188, .00019, .00006, + .00003, .00029, .00088, + .16086, .03347, .01798, .00019, .00006, + .00004, .00026, .00082, + .14774, .03043, .01502, .00019, .00008, + .00005, .00025, .00075 / DATA ((F(I,38,K),I=1,5),F(7,38,K),F(6,38,K),F(8,38,K),K=8,14)/ + .13640, .02784, .01273, .00017, .00009, + .00005, .00024, .00066, + .12652, .02561, .01092, .00018, .00008, + .00005, .00023, .00060, + .11786, .02369, .00947, .00018, .00007, + .00004, .00021, .00055, + .11022, .02200, .00830, .00016, .00008, + .00005, .00020, .00050, + .10343, .02049, .00733, .00013, .00007, + .00005, .00019, .00048, + .09732, .01915, .00653, .00013, .00007, + .00005, .00018, .00045, + .09180, .01797, .00585, .00014, .00008, + .00005, .00017, .00041 / DATA ((F(I,38,K),I=1,5),F(7,38,K),F(6,38,K),F(8,38,K),K=15,19)/ + .08684, .01691, .00528, .00014, .00008, + .00004, .00015, .00037, + .08230, .01594, .00479, .00014, .00007, + .00004, .00014, .00035, + .07819, .01507, .00437, .00013, .00007, + .00004, .00014, .00034, + .07442, .01428, .00401, .00014, .00007, + .00004, .00013, .00030, + .07098, .01356, .00369, .00012, .00006, + .00004, .00013, .00028 / DATA ((F(I,39,K),I=1,5),F(7,39,K),F(6,39,K),F(8,39,K),K=1,7)/ + .24509, .05196, .04632, .00001, .00001, + .00004, .00011, .00085, + .21210, .04406, .03370, .00001, .00001, + .00004, .00019, .00078, + .18655, .03810, .02549, .00003, .00001, + .00002, .00018, .00070, + .16611, .03346, .01987, .00009, .00003, + .00003, .00019, .00064, + .14945, .02973, .01588, .00006, .00004, + .00004, .00018, .00056, + .13556, .02667, .01294, .00007, .00006, + .00004, .00018, .00050, + .12385, .02414, .01074, .00010, .00007, + .00004, .00018, .00044 / DATA ((F(I,39,K),I=1,5),F(7,39,K),F(6,39,K),F(8,39,K),K=8,14)/ + .11382, .02197, .00903, .00008, .00007, + .00004, .00016, .00041, + .10511, .02012, .00770, .00008, .00008, + .00004, .00015, .00037, + .09751, .01852, .00665, .00009, .00007, + .00004, .00014, .00034, + .09083, .01713, .00580, .00009, .00006, + .00004, .00013, .00033, + .08491, .01591, .00510, .00008, .00006, + .00004, .00012, .00029, + .07962, .01482, .00452, .00007, .00006, + .00003, .00012, .00026, + .07490, .01386, .00404, .00007, .00006, + .00003, .00011, .00025 / DATA ((F(I,39,K),I=1,5),F(7,39,K),F(6,39,K),F(8,39,K),K=15,19)/ + .07064, .01300, .00363, .00007, .00005, + .00003, .00010, .00023, + .06675, .01222, .00329, .00007, .00005, + .00003, .00010, .00022, + .06323, .01153, .00299, .00006, .00005, + .00003, .00009, .00020, + .06003, .01089, .00274, .00005, .00004, + .00002, .00009, .00019, + .05711, .01031, .00251, .00005, .00004, + .00002, .00008, .00018 / DATA ((F(I,40,K),I=1,5),F(7,40,K),F(6,40,K),F(8,40,K),K=1,7)/ + .21200, .04251, .03478, .00001, .00001, + .00000, .00001, .00053, + .18163, .03571, .02492, .00002, .00001, + .00000, .00005, .00047, + .15838, .03066, .01862, .00006, .00001, + .00001, .00008, .00045, + .14006, .02670, .01436, .00006, .00001, + .00000, .00009, .00037, + .12519, .02356, .01137, .00007, .00002, + .00001, .00010, .00034, + .11293, .02103, .00920, .00008, .00002, + .00002, .00009, .00031, + .10265, .01893, .00757, .00007, .00002, + .00002, .00009, .00030 / DATA ((F(I,40,K),I=1,5),F(7,40,K),F(6,40,K),F(8,40,K),K=8,14)/ + .09386, .01715, .00633, .00007, .00003, + .00001, .00007, .00025, + .08629, .01564, .00537, .00008, .00003, + .00002, .00008, .00022, + .07974, .01434, .00461, .00008, .00003, + .00002, .00007, .00020, + .07397, .01322, .00400, .00007, .00002, + .00001, .00006, .00018, + .06889, .01223, .00351, .00008, .00002, + .00001, .00006, .00018, + .06439, .01137, .00310, .00005, .00003, + .00001, .00006, .00017, + .06038, .01059, .00276, .00004, .00003, + .00001, .00005, .00015 / DATA ((F(I,40,K),I=1,5),F(7,40,K),F(6,40,K),F(8,40,K),K=15,19)/ + .05675, .00989, .00248, .00006, .00003, + .00001, .00005, .00014, + .05347, .00929, .00224, .00006, .00003, + .00001, .00005, .00013, + .05052, .00872, .00203, .00006, .00002, + .00001, .00004, .00013, + .04784, .00821, .00185, .00004, .00002, + .00001, .00004, .00012, + .04539, .00776, .00169, .00005, .00002, + .00001, .00004, .00010 / DATA ((F(I,41,K),I=1,5),F(7,41,K),F(6,41,K),F(8,41,K),K=1,7)/ + .18011, .03388, .02551, .00001, .00001, + .00003, .00001, .00029, + .15273, .02821, .01800, .00001, .00001, + .00002, .00005, .00028, + .13211, .02402, .01327, .00001, .00001, + .00001, .00006, .00024, + .11593, .02078, .01013, .00001, .00001, + .00002, .00007, .00023, + .10291, .01823, .00794, .00001, .00001, + .00002, .00006, .00020, + .09228, .01618, .00637, .00002, .00002, + .00002, .00007, .00018, + .08339, .01448, .00520, .00002, .00002, + .00002, .00006, .00015 / DATA ((F(I,41,K),I=1,5),F(7,41,K),F(6,41,K),F(8,41,K),K=8,14)/ + .07588, .01305, .00432, .00003, .00002, + .00002, .00006, .00013, + .06943, .01185, .00364, .00003, .00003, + .00002, .00005, .00013, + .06386, .01083, .00311, .00003, .00003, + .00002, .00005, .00012, + .05902, .00993, .00269, .00003, .00002, + .00001, .00004, .00009, + .05474, .00916, .00235, .00003, .00003, + .00001, .00004, .00009, + .05098, .00847, .00207, .00003, .00002, + .00001, .00004, .00008, + .04764, .00787, .00183, .00002, .00002, + .00001, .00004, .00008 / DATA ((F(I,41,K),I=1,5),F(7,41,K),F(6,41,K),F(8,41,K),K=15,19)/ + .04464, .00733, .00164, .00002, .00002, + .00001, .00003, .00008, + .04194, .00685, .00148, .00002, .00002, + .00001, .00003, .00007, + .03950, .00643, .00134, .00002, .00002, + .00001, .00003, .00007, + .03730, .00604, .00122, .00002, .00001, + .00001, .00003, .00006, + .03530, .00569, .00111, .00001, .00001, + .00001, .00002, .00006 / DATA ((F(I,42,K),I=1,5),F(7,42,K),F(6,42,K),F(8,42,K),K=1,7)/ + .12388, .02055, .01292, .00001, .00001, + .00000, .00001, .00008, + .10291, .01676, .00883, .00001, .00001, + .00000, .00000, .00007, + .08742, .01400, .00634, .00001, .00001, + .00000, .00001, .00006, + .07549, .01195, .00472, .00001, .00001, + .00000, .00001, .00006, + .06607, .01034, .00363, .00001, .00000, + .00000, .00001, .00006, + .05851, .00905, .00286, .00000, .00001, + .00000, .00002, .00005, + .05228, .00801, .00230, .00001, .00001, + .00000, .00002, .00005 / DATA ((F(I,42,K),I=1,5),F(7,42,K),F(6,42,K),F(8,42,K),K=8,14)/ + .04705, .00715, .00189, .00001, .00001, + .00000, .00002, .00005, + .04265, .00643, .00158, .00001, .00001, + .00000, .00001, .00003, + .03889, .00581, .00133, .00001, .00001, + .00000, .00001, .00003, + .03563, .00529, .00114, .00001, .00001, + .00000, .00001, .00004, + .03279, .00484, .00099, .00001, .00001, + .00000, .00001, .00002, + .03029, .00445, .00086, .00001, .00001, + .00000, .00001, .00002, + .02811, .00410, .00076, .00001, .00001, + .00000, .00001, .00002 / DATA ((F(I,42,K),I=1,5),F(7,42,K),F(6,42,K),F(8,42,K),K=15,19)/ + .02616, .00379, .00068, .00001, .00001, + .00000, .00001, .00002, + .02442, .00353, .00061, .00001, .00001, + .00000, .00001, .00002, + .02286, .00328, .00055, .00001, .00001, + .00000, .00001, .00001, + .02145, .00307, .00049, .00001, .00001, + .00000, .00001, .00001, + .02018, .00287, .00045, .00001, .00001, + .00000, .00001, .00001 / DATA ((F(I,43,K),I=1,5),F(7,43,K),F(6,43,K),F(8,43,K),K=1,7)/ + .07832, .01128, .00585, .00001, .00001, + .00000, .00001, .00000, + .06364, .00901, .00386, .00001, .00001, + .00000, .00000, .00002, + .05299, .00740, .00269, .00001, .00001, + .00000, .00001, .00002, + .04501, .00620, .00196, .00001, .00001, + .00000, .00001, .00001, + .03881, .00529, .00148, .00001, .00001, + .00000, .00001, .00001, + .03389, .00457, .00114, .00001, .00001, + .00000, .00001, .00001, + .02993, .00400, .00091, .00001, .00001, + .00000, .00001, .00001 / DATA ((F(I,43,K),I=1,5),F(7,43,K),F(6,43,K),F(8,43,K),K=8,14)/ + .02664, .00351, .00073, .00001, .00001, + .00000, .00001, .00001, + .02388, .00314, .00060, .00001, .00001, + .00000, .00001, .00001, + .02155, .00281, .00051, .00001, .00001, + .00000, .00001, .00001, + .01956, .00253, .00043, .00001, .00001, + .00000, .00001, .00001, + .01785, .00230, .00037, .00001, .00001, + .00000, .00001, .00001, + .01636, .00209, .00032, .00001, .00001, + .00000, .00001, .00001, + .01505, .00192, .00028, .00001, .00001, + .00000, .00001, .00001 / DATA ((F(I,43,K),I=1,5),F(7,43,K),F(6,43,K),F(8,43,K),K=15,19)/ + .01390, .00176, .00025, .00001, .00001, + .00000, .00001, .00001, + .01289, .00162, .00022, .00001, .00001, + .00000, .00001, .00001, + .01199, .00150, .00020, .00001, .00001, + .00000, .00001, .00001, + .01119, .00139, .00018, .00001, .00001, + .00000, .00001, .00001, + .01047, .00130, .00016, .00001, .00001, + .00000, .00001, .00001 / DATA ((F(I,44,K),I=1,5),F(7,44,K),F(6,44,K),F(8,44,K),K=1,7)/ + .04435, .00544, .00227, .00001, .00001, + .00000, .00001, .00001, + .03516, .00424, .00145, .00001, .00001, + .00000, .00001, .00001, + .02865, .00341, .00098, .00001, .00001, + .00000, .00001, .00001, + .02390, .00281, .00069, .00001, .00001, + .00000, .00001, .00001, + .02024, .00235, .00051, .00001, .00001, + .00000, .00001, .00001, + .01741, .00200, .00039, .00001, .00001, + .00000, .00001, .00001, + .01516, .00173, .00030, .00001, .00001, + .00000, .00001, .00001 / DATA ((F(I,44,K),I=1,5),F(7,44,K),F(6,44,K),F(8,44,K),K=8,14)/ + .01332, .00151, .00024, .00001, .00001, + .00000, .00001, .00001, + .01180, .00132, .00020, .00001, .00001, + .00000, .00001, .00001, + .01053, .00118, .00016, .00001, .00001, + .00000, .00001, .00001, + .00946, .00105, .00014, .00001, .00001, + .00000, .00001, .00001, + .00856, .00094, .00012, .00001, .00001, + .00000, .00001, .00001, + .00777, .00085, .00010, .00001, .00001, + .00000, .00001, .00001, + .00708, .00077, .00009, .00001, .00001, + .00000, .00001, .00001 / DATA ((F(I,44,K),I=1,5),F(7,44,K),F(6,44,K),F(8,44,K),K=15,19)/ + .00649, .00070, .00008, .00001, .00001, + .00000, .00001, .00001, + .00598, .00064, .00007, .00001, .00001, + .00000, .00001, .00001, + .00551, .00058, .00006, .00001, .00001, + .00000, .00001, .00001, + .00510, .00054, .00006, .00001, .00001, + .00000, .00001, .00001, + .00474, .00050, .00005, .00001, .00001, + .00000, .00001, .00001 / DATA ((F(I,45,K),I=1,5),F(7,45,K),F(6,45,K),F(8,45,K),K=1,7)/ + .02114, .00218, .00070, .00001, .00001, + .00000, .00001, .00001, + .01629, .00164, .00043, .00001, .00001, + .00000, .00001, .00001, + .01295, .00129, .00028, .00001, .00001, + .00000, .00001, .00001, + .01056, .00106, .00019, .00001, .00001, + .00000, .00001, .00001, + .00877, .00085, .00014, .00001, .00001, + .00000, .00001, .00001, + .00742, .00071, .00010, .00001, .00001, + .00000, .00001, .00001, + .00635, .00062, .00008, .00001, .00001, + .00000, .00001, .00001 / DATA ((F(I,45,K),I=1,5),F(7,45,K),F(6,45,K),F(8,45,K),K=8,14)/ + .00549, .00053, .00006, .00001, .00001, + .00000, .00001, .00001, + .00480, .00045, .00005, .00001, .00001, + .00000, .00001, .00001, + .00424, .00039, .00004, .00001, .00001, + .00000, .00001, .00001, + .00376, .00034, .00004, .00001, .00001, + .00000, .00001, .00001, + .00335, .00030, .00003, .00001, .00001, + .00000, .00001, .00001, + .00302, .00027, .00003, .00001, .00001, + .00000, .00001, .00001, + .00272, .00024, .00002, .00001, .00001, + .00000, .00001, .00001 / DATA ((F(I,45,K),I=1,5),F(7,45,K),F(6,45,K),F(8,45,K),K=15,19)/ + .00248, .00022, .00002, .00001, .00001, + .00000, .00001, .00001, + .00225, .00019, .00002, .00001, .00001, + .00000, .00001, .00001, + .00207, .00018, .00002, .00001, .00001, + .00000, .00001, .00001, + .00190, .00016, .00001, .00001, .00001, + .00000, .00001, .00001, + .00174, .00015, .00001, .00001, .00001, + .00000, .00001, .00001 / DATA ((F(I,46,K),I=1,5),F(7,46,K),F(6,46,K),F(8,46,K),K=1,7)/ + .00187, .00011, .00002, .00001, .00001, + .00000, .00001, .00001, + .00132, .00008, .00001, .00001, .00001, + .00000, .00001, .00001, + .00099, .00006, .00001, .00001, .00001, + .00000, .00001, .00001, + .00075, .00004, .00000, .00000, .00000, + .00000, .00000, .00000, + .00060, .00004, .00000, .00000, .00000, + .00000, .00000, .00000, + .00048, .00003, .00000, .00000, .00000, + .00000, .00000, .00000, + .00039, .00002, .00000, .00000, .00000, + .00000, .00000, .00000 / DATA ((F(I,46,K),I=1,5),F(7,46,K),F(6,46,K),F(8,46,K),K=8,14)/ + .00032, .00002, .00000, .00000, .00000, + .00000, .00000, .00000, + .00028, .00002, .00000, .00000, .00000, + .00000, .00000, .00000, + .00023, .00002, .00000, .00000, .00000, + .00000, .00000, .00000, + .00020, .00002, .00000, .00000, .00000, + .00000, .00000, .00000, + .00017, .00002, .00000, .00000, .00000, + .00000, .00000, .00000, + .00015, .00002, .00000, .00000, .00000, + .00000, .00000, .00000, + .00013, .00001, .00000, .00000, .00000, + .00000, .00000, .00000 / DATA ((F(I,46,K),I=1,5),F(7,46,K),F(6,46,K),F(8,46,K),K=15,19)/ + .00012, .00001, .00000, .00000, .00000, + .00000, .00000, .00000, + .00011, .00001, .00000, .00000, .00000, + .00000, .00000, .00000, + .00009, .00001, .00000, .00000, .00000, + .00000, .00000, .00000, + .00008, .00001, .00000, .00000, .00000, + .00000, .00000, .00000, + .00008, .00001, .00000, .00000, .00000, + .00000, .00000, .00000 / DATA XX/1.D-5,2.D-5,4.D-5,6.D-5,8.D-5, . 1.D-4,2.D-4,4.D-4,6.D-4,8.D-4, . 1.D-3,2.D-3,4.D-3,6.D-3,8.D-3, . 1.D-2,2.D-2,4.D-2,6.D-2,8.D-2, . .1D0,.125D0,.15D0,.175D0,.2D0,.225D0,.25D0,.275D0, . .3D0,.325D0,.35D0,.375D0,.4D0,.425D0,.45D0,.475D0, . .5D0,.525D0,.55D0,.575D0,.6D0,.65D0,.7D0,.75D0, . .8D0,.9D0,1.D0/ DATA XMIN,XMAX,QSQMIN,QSQMAX/1.D-5,1.D0,5.D0,1310720.D0/ DATA N0/2,5,5,9,0,0,9,9/ DATA ZEROD/0.D0/,PONED/0.1D0/,ONED/1.D0/,ONEDO/1.1D0/,TWOD/2.D0/ DATA INIT/0/ C xsave=x C IF(INIT.NE.0) GOTO 10 INIT=1 DO 20 N=1,NX-1 DO 20 M=1,19 C READ(35,50)F(1,N,M),F(2,N,M),F(3,N,M),F(4,N,M),F(5,N,M),F(7,N,M), C . F(6,N,M),F(8,N,M) C 1=UV 2=DV 3=GLUE 4=UBAR 5=CBAR 7=BBAR 6=SBAR 8=DBAR DO 25 I=1,8 25 F(I,N,M)=F(I,N,M)/(ONED-XX(N))**N0(I) 20 CONTINUE DO 31 J=1,NTENTH-1 XX(J)= LOG10(XX(J))+ONEDO DO 31 I=1,8 IF(I.EQ.7) GO TO 31 DO 30 K=1,19 30 F(I,J,K)= LOG(F(I,J,K))*F(I,NTENTH,K)/ LOG(F(I,NTENTH,K)) 31 CONTINUE 50 FORMAT(8F10.5) DO 40 I=1,8 DO 40 M=1,19 40 F(I,NX,M)=ZEROD 10 CONTINUE IF(X.LT.XMIN) X=XMIN IF(X.GT.XMAX) X=XMAX QSQ=SCALE**2 IF(QSQ.LT.QSQMIN) QSQ=QSQMIN IF(QSQ.GT.QSQMAX) QSQ=QSQMAX XXX=X IF(X.LT.PONED) XXX= LOG10(X)+ONEDO N=0 70 N=N+1 IF(XXX.GT.XX(N+1)) GOTO 70 A=(XXX-XX(N))/(XX(N+1)-XX(N)) RM= LOG(QSQ/QSQMIN)/ LOG(TWOD) B=RM-AINT(RM) M=1+ INT(RM) DO 60 I=1,8 G(I)= (ONED-A)*(ONED-B)*F(I,N,M)+(ONED-A)*B*F(I,N,M+1) . + A*(ONED-B)*F(I,N+1,M) + A*B*F(I,N+1,M+1) IF(N.GE.NTENTH) GOTO 65 IF(I.EQ.7) GOTO 65 FAC=(ONED-B)*F(I,NTENTH,M)+B*F(I,NTENTH,M+1) G(I)=FAC**(G(I)/FAC) 65 CONTINUE G(I)=G(I)*(ONED-X)**N0(I) 60 CONTINUE UPV=G(1) DNV=G(2) USEA=G(4) DSEA=G(8) STR=G(6) CHM=G(5) GLU=G(3) BOT=G(7) x=xsave RETURN END