| 21886bb6 |
1 | #include "pdf/pilot.h" |
| 2 | SUBROUTINE SFDSQA(X,SCALE,UPV,DNV,USEA,DSEA,STR,CHM,BOT,GLU) |
| 3 | |
| 4 | C THIS IS THE NEW "A" FIT -- May 1994 -- low Q^2 range |
| 5 | C DIS DIS DIS DIS DIS DIS DIS DIS DIS DIS DIS DIS DIS DIS |
| 6 | C |
| 7 | #include "pdf/impdp.inc" |
| 8 | C |
| 9 | PARAMETER(NX=47) |
| 10 | PARAMETER(NTENTH=21) |
| 11 | DIMENSION F(8,NX,8),G(8),XX(NX),N0(8) |
| 12 | DATA ((F(I,1,K),I=1,5),F(7,1,K),F(6,1,K),F(8,1,K),K=1,7)/ |
| 13 | + .00353, .00627, 4.23812, 1.87623, .00000, |
| 14 | + .00000, .51265, 1.87021, |
| 15 | + .00377, .00645, 6.42705, 1.94786, .00000, |
| 16 | + .00000, .61236, 1.94247, |
| 17 | + .00399, .00660, 9.43898, 2.04451, .00000, |
| 18 | + .00000, .73330, 2.03973, |
| 19 | + .00420, .00675, 13.16286, 2.16589, .00000, |
| 20 | + .00000, .87600, 2.16169, |
| 21 | + .00440, .00688, 17.54780, 2.30621, .00000, |
| 22 | + .00000, 1.03533, 2.30254, |
| 23 | + .00460, .00702, 22.55753, 2.46430, .06825, |
| 24 | + .00000, 1.21072, 2.46116, |
| 25 | + .00477, .00714, 28.12229, 2.63824, .27055, |
| 26 | + .00000, 1.40039, 2.63556 / |
| 27 | DATA ((F(I,2,K),I=1,5),F(7,2,K),F(6,2,K),F(8,2,K),K=1,7)/ |
| 28 | + .00538, .00802, 3.22632, 1.51244, .00000, |
| 29 | + .00000, .40761, 1.50993, |
| 30 | + .00571, .00826, 5.15995, 1.57471, .00000, |
| 31 | + .00000, .49238, 1.57262, |
| 32 | + .00602, .00846, 7.71045, 1.65616, .00000, |
| 33 | + .00000, .59325, 1.65446, |
| 34 | + .00632, .00866, 10.78281, 1.75593, .00000, |
| 35 | + .00000, .71002, 1.75459, |
| 36 | + .00659, .00884, 14.34431, 1.86953, .00000, |
| 37 | + .00000, .83878, 1.86852, |
| 38 | + .00686, .00902, 18.35565, 1.99595, .05595, |
| 39 | + .00000, .97896, 1.99526, |
| 40 | + .00712, .00919, 22.75467, 2.13362, .21642, |
| 41 | + .00000, 1.12915, 2.13323 / |
| 42 | DATA ((F(I,3,K),I=1,5),F(7,3,K),F(6,3,K),F(8,3,K),K=1,7)/ |
| 43 | + .00806, .01027, 2.45661, 1.21901, .00000, |
| 44 | + .00000, .32420, 1.21891, |
| 45 | + .00853, .01059, 4.14420, 1.27292, .00000, |
| 46 | + .00000, .39616, 1.27308, |
| 47 | + .00896, .01087, 6.30016, 1.34140, .00000, |
| 48 | + .00000, .48020, 1.34180, |
| 49 | + .00937, .01113, 8.83462, 1.42328, .00000, |
| 50 | + .00000, .57568, 1.42390, |
| 51 | + .00975, .01137, 11.72666, 1.51512, .00000, |
| 52 | + .00000, .67963, 1.51595, |
| 53 | + .01013, .01162, 14.93678, 1.61607, .04591, |
| 54 | + .00000, .79155, 1.61709, |
| 55 | + .01048, .01184, 18.41091, 1.72486, .17321, |
| 56 | + .00000, .91032, 1.72606 / |
| 57 | DATA ((F(I,4,K),I=1,5),F(7,4,K),F(6,4,K),F(8,4,K),K=1,7)/ |
| 58 | + .01014, .01188, 2.09495, 1.07443, .00000, |
| 59 | + .00000, .28369, 1.07538, |
| 60 | + .01071, .01225, 3.64657, 1.12390, .00000, |
| 61 | + .00000, .34904, 1.12504, |
| 62 | + .01125, .01259, 5.59929, 1.18568, .00000, |
| 63 | + .00000, .42450, 1.18701, |
| 64 | + .01174, .01289, 7.86366, 1.25856, .00000, |
| 65 | + .00000, .50933, 1.26006, |
| 66 | + .01220, .01319, 10.42364, 1.33960, .00000, |
| 67 | + .00000, .60099, 1.34125, |
| 68 | + .01266, .01348, 13.24045, 1.42800, .04089, |
| 69 | + .00000, .69902, 1.42979, |
| 70 | + .01308, .01374, 16.26478, 1.52271, .15196, |
| 71 | + .00000, .80246, 1.52463 / |
| 72 | DATA ((F(I,5,K),I=1,5),F(7,5,K),F(6,5,K),F(8,5,K),K=1,7)/ |
| 73 | + .01193, .01319, 1.87136, .98235, .00000, |
| 74 | + .00000, .25813, .98392, |
| 75 | + .01258, .01359, 3.33089, 1.02884, .00000, |
| 76 | + .00000, .31915, 1.03058, |
| 77 | + .01318, .01397, 5.15064, 1.08626, .00000, |
| 78 | + .00000, .38906, 1.08815, |
| 79 | + .01374, .01433, 7.24090, 1.15331, .00000, |
| 80 | + .00000, .46701, 1.15535, |
| 81 | + .01428, .01465, 9.58839, 1.22740, .00000, |
| 82 | + .00000, .55080, 1.22955, |
| 83 | + .01479, .01496, 12.15519, 1.30782, .03767, |
| 84 | + .00000, .63999, 1.31010, |
| 85 | + .01529, .01527, 14.89528, 1.39358, .13842, |
| 86 | + .00000, .73373, 1.39596 / |
| 87 | DATA ((F(I,6,K),I=1,5),F(7,6,K),F(6,6,K),F(8,6,K),K=1,7)/ |
| 88 | + .01351, .01430, 1.71468, .91640, .00000, |
| 89 | + .00000, .23995, .91839, |
| 90 | + .01423, .01475, 3.10549, .96070, .00000, |
| 91 | + .00000, .29779, .96283, |
| 92 | + .01490, .01516, 4.82815, 1.01490, .00000, |
| 93 | + .00000, .36364, 1.01717, |
| 94 | + .01552, .01555, 6.79254, 1.07773, .00000, |
| 95 | + .00000, .43666, 1.08013, |
| 96 | + .01611, .01591, 8.98726, 1.14683, .00000, |
| 97 | + .00000, .51480, 1.14934, |
| 98 | + .01668, .01626, 11.37514, 1.22153, .03529, |
| 99 | + .00000, .59766, 1.22414, |
| 100 | + .01724, .01658, 13.91271, 1.30091, .12850, |
| 101 | + .00000, .68447, 1.30361 / |
| 102 | DATA ((F(I,7,K),I=1,5),F(7,7,K),F(6,7,K),F(8,7,K),K=1,7)/ |
| 103 | + .01973, .01843, 1.30792, .73836, .00000, |
| 104 | + .00000, .19165, .74150, |
| 105 | + .02071, .01901, 2.50161, .77638, .00000, |
| 106 | + .00000, .24064, .77965, |
| 107 | + .02164, .01956, 3.95367, .82157, .00000, |
| 108 | + .00000, .29525, .82495, |
| 109 | + .02250, .02006, 5.57292, .87276, .00000, |
| 110 | + .00000, .35466, .87624, |
| 111 | + .02331, .02054, 7.35237, .92821, .00000, |
| 112 | + .00000, .41740, .93177, |
| 113 | + .02411, .02099, 9.25794, .98738, .02895, |
| 114 | + .00000, .48318, .99104, |
| 115 | + .02486, .02143, 11.25363, 1.04960, .10236, |
| 116 | + .00000, .55137, 1.05333 / |
| 117 | DATA ((F(I,8,K),I=1,5),F(7,8,K),F(6,8,K),F(8,8,K),K=1,7)/ |
| 118 | + .02860, .02387, .99986, .59496, .00000, |
| 119 | + .00000, .15390, .59922, |
| 120 | + .02994, .02463, 2.02268, .62745, .00000, |
| 121 | + .00000, .19534, .63186, |
| 122 | + .03120, .02534, 3.24642, .66488, .00000, |
| 123 | + .00000, .24044, .66941, |
| 124 | + .03238, .02600, 4.58011, .70627, .00000, |
| 125 | + .00000, .28853, .71093, |
| 126 | + .03349, .02660, 6.02017, .75042, .00000, |
| 127 | + .00000, .33861, .75519, |
| 128 | + .03458, .02721, 7.53631, .79695, .02377, |
| 129 | + .00000, .39051, .80182, |
| 130 | + .03562, .02777, 9.09952, .84533, .08133, |
| 131 | + .00000, .44376, .85029 / |
| 132 | DATA ((F(I,9,K),I=1,5),F(7,9,K),F(6,9,K),F(8,9,K),K=1,7)/ |
| 133 | + .03541, .02785, .85607, .52448, .00000, |
| 134 | + .00000, .13594, .52948, |
| 135 | + .03702, .02872, 1.79177, .55403, .00000, |
| 136 | + .00000, .17343, .55918, |
| 137 | + .03853, .02954, 2.89947, .58737, .00000, |
| 138 | + .00000, .21363, .59270, |
| 139 | + .03994, .03030, 4.08932, .62376, .00000, |
| 140 | + .00000, .25595, .62921, |
| 141 | + .04128, .03102, 5.35969, .66221, .00000, |
| 142 | + .00000, .29970, .66779, |
| 143 | + .04259, .03172, 6.68285, .70242, .02116, |
| 144 | + .00000, .34469, .70814, |
| 145 | + .04384, .03239, 8.03367, .74396, .07087, |
| 146 | + .00000, .39058, .74981 / |
| 147 | DATA ((F(I,10,K),I=1,5),F(7,10,K),F(6,10,K),F(8,10,K),K=1,7)/ |
| 148 | + .04114, .03109, .76776, .47972, .00000, |
| 149 | + .00000, .12483, .48529, |
| 150 | + .04297, .03205, 1.64767, .50719, .00000, |
| 151 | + .00000, .15966, .51298, |
| 152 | + .04468, .03295, 2.68028, .53784, .00000, |
| 153 | + .00000, .19662, .54379, |
| 154 | + .04629, .03380, 3.77709, .57092, .00000, |
| 155 | + .00000, .23521, .57706, |
| 156 | + .04781, .03460, 4.93803, .60569, .00000, |
| 157 | + .00000, .27487, .61197, |
| 158 | + .04931, .03538, 6.13732, .64187, .01945, |
| 159 | + .00000, .31546, .64830, |
| 160 | + .05072, .03614, 7.35252, .67905, .06418, |
| 161 | + .00000, .35667, .68564 / |
| 162 | DATA ((F(I,11,K),I=1,5),F(7,11,K),F(6,11,K),F(8,11,K),K=1,7)/ |
| 163 | + .04628, .03389, .70475, .44767, .00000, |
| 164 | + .00000, .11695, .45366, |
| 165 | + .04828, .03494, 1.54754, .47364, .00000, |
| 166 | + .00000, .14979, .47983, |
| 167 | + .05016, .03594, 2.52696, .50227, .00000, |
| 168 | + .00000, .18433, .50868, |
| 169 | + .05194, .03687, 3.55628, .53299, .00000, |
| 170 | + .00000, .22019, .53956, |
| 171 | + .05361, .03774, 4.63730, .56508, .00000, |
| 172 | + .00000, .25686, .57183, |
| 173 | + .05525, .03861, 5.74604, .59832, .01817, |
| 174 | + .00000, .29424, .60523, |
| 175 | + .05683, .03941, 6.86227, .63238, .05918, |
| 176 | + .00000, .33205, .63944 / |
| 177 | DATA ((F(I,12,K),I=1,5),F(7,12,K),F(6,12,K),F(8,12,K),K=1,7)/ |
| 178 | + .06634, .04460, .54955, .36164, .00000, |
| 179 | + .00000, .09697, .36963, |
| 180 | + .06902, .04597, 1.28832, .38302, .00000, |
| 181 | + .00000, .12399, .39128, |
| 182 | + .07156, .04728, 2.11904, .40586, .00000, |
| 183 | + .00000, .15169, .41438, |
| 184 | + .07395, .04850, 2.96136, .42982, .00000, |
| 185 | + .00000, .17986, .43859, |
| 186 | + .07622, .04966, 3.82264, .45449, .00000, |
| 187 | + .00000, .20829, .46348, |
| 188 | + .07843, .05079, 4.68465, .47971, .01465, |
| 189 | + .00000, .23691, .48893, |
| 190 | + .08055, .05186, 5.53362, .50526, .04585, |
| 191 | + .00000, .26555, .51468 / |
| 192 | DATA ((F(I,13,K),I=1,5),F(7,13,K),F(6,13,K),F(8,13,K),K=1,7)/ |
| 193 | + .09459, .05933, .45006, .29360, .00000, |
| 194 | + .00000, .08276, .30437, |
| 195 | + .09814, .06112, 1.10456, .31050, .00000, |
| 196 | + .00000, .10425, .32162, |
| 197 | + .10155, .06285, 1.80790, .32807, .00000, |
| 198 | + .00000, .12580, .33952, |
| 199 | + .10473, .06446, 2.49068, .34620, .00000, |
| 200 | + .00000, .14739, .35798, |
| 201 | + .10777, .06600, 3.16635, .36464, .00000, |
| 202 | + .00000, .16891, .37671, |
| 203 | + .11072, .06748, 3.82305, .38326, .01152, |
| 204 | + .00000, .19031, .39562, |
| 205 | + .11353, .06890, 4.45308, .40190, .03481, |
| 206 | + .00000, .21149, .41453 / |
| 207 | DATA ((F(I,14,K),I=1,5),F(7,14,K),F(6,14,K),F(8,14,K),K=1,7)/ |
| 208 | + .11598, .07030, .42083, .26052, .00000, |
| 209 | + .00000, .07654, .27338, |
| 210 | + .12018, .07245, 1.03073, .27483, .00000, |
| 211 | + .00000, .09494, .28811, |
| 212 | + .12420, .07450, 1.66542, .28956, .00000, |
| 213 | + .00000, .11322, .30324, |
| 214 | + .12798, .07641, 2.26414, .30469, .00000, |
| 215 | + .00000, .13143, .31874, |
| 216 | + .13156, .07822, 2.84358, .31997, .00000, |
| 217 | + .00000, .14946, .33439, |
| 218 | + .13503, .07996, 3.39582, .33531, .00979, |
| 219 | + .00000, .16730, .35006, |
| 220 | + .13832, .08161, 3.91637, .35057, .02918, |
| 221 | + -.00001, .18482, .36564 / |
| 222 | DATA ((F(I,15,K),I=1,5),F(7,15,K),F(6,15,K),F(8,15,K),K=1,7)/ |
| 223 | + .13377, .07939, .41558, .23959, .00000, |
| 224 | + .00000, .07283, .25410, |
| 225 | + .13849, .08184, .99416, .25213, .00000, |
| 226 | + .00000, .08912, .26710, |
| 227 | + .14304, .08416, 1.58125, .26497, .00000, |
| 228 | + .00000, .10522, .28042, |
| 229 | + .14730, .08633, 2.12318, .27815, .00000, |
| 230 | + .00000, .12124, .29401, |
| 231 | + .15132, .08837, 2.63867, .29140, .00000, |
| 232 | + .00000, .13703, .30768, |
| 233 | + .15519, .09033, 3.12247, .30466, .00857, |
| 234 | + .00000, .15259, .32132, |
| 235 | + .15886, .09215, 3.57229, .31776, .02549, |
| 236 | + -.00001, .16780, .33479 / |
| 237 | DATA ((F(I,16,K),I=1,5),F(7,16,K),F(6,16,K),F(8,16,K),K=1,7)/ |
| 238 | + .14903, .08698, .42079, .22390, .00000, |
| 239 | + .00000, .06996, .23937, |
| 240 | + .15428, .08974, .97466, .23519, .00000, |
| 241 | + .00000, .08469, .25125, |
| 242 | + .15929, .09235, 1.52498, .24674, .00000, |
| 243 | + .00000, .09927, .26336, |
| 244 | + .16398, .09477, 2.02410, .25857, .00000, |
| 245 | + .00000, .11375, .27570, |
| 246 | + .16840, .09703, 2.49209, .27045, .00000, |
| 247 | + .00000, .12799, .28805, |
| 248 | + .17264, .09917, 2.92573, .28225, .00757, |
| 249 | + .00000, .14197, .30029, |
| 250 | + .17663, .10116, 3.32426, .29388, .02270, |
| 251 | + .00000, .15557, .31233 / |
| 252 | DATA ((F(I,17,K),I=1,5),F(7,17,K),F(6,17,K),F(8,17,K),K=1,7)/ |
| 253 | + .20755, .11696, .49179, .18309, .00000, |
| 254 | + .00000, .06331, .20391, |
| 255 | + .21483, .12090, .96223, .19043, .00000, |
| 256 | + .00000, .07349, .21220, |
| 257 | + .22167, .12450, 1.39560, .19812, .00000, |
| 258 | + .00000, .08371, .22074, |
| 259 | + .22795, .12773, 1.76565, .20605, .00000, |
| 260 | + .00000, .09388, .22941, |
| 261 | + .23378, .13067, 2.09508, .21395, .00000, |
| 262 | + .00000, .10382, .23798, |
| 263 | + .23925, .13335, 2.38648, .22169, .00451, |
| 264 | + .00000, .11345, .24633, |
| 265 | + .24430, .13580, 2.64317, .22922, .01502, |
| 266 | + .00000, .12270, .25439 / |
| 267 | DATA ((F(I,18,K),I=1,5),F(7,18,K),F(6,18,K),F(8,18,K),K=1,7)/ |
| 268 | + .28677, .15844, .65201, .14832, .00000, |
| 269 | + .00000, .05678, .17806, |
| 270 | + .29713, .16379, 1.01515, .15204, .00000, |
| 271 | + .00000, .06325, .18312, |
| 272 | + .30643, .16834, 1.31959, .15631, .00000, |
| 273 | + .00000, .06998, .18849, |
| 274 | + .31466, .17220, 1.56052, .16082, .00000, |
| 275 | + .00000, .07671, .19391, |
| 276 | + .32205, .17554, 1.76050, .16531, .00000, |
| 277 | + .00000, .08320, .19917, |
| 278 | + .32867, .17842, 1.92619, .16961, .00128, |
| 279 | + .00000, .08935, .20410, |
| 280 | + .33465, .18092, 2.06321, .17369, .00833, |
| 281 | + -.00001, .09512, .20870 / |
| 282 | DATA ((F(I,19,K),I=1,5),F(7,19,K),F(6,19,K),F(8,19,K),K=1,7)/ |
| 283 | + .34542, .18867, .77868, .12710, .00000, |
| 284 | + .00000, .05139, .16412, |
| 285 | + .35819, .19454, 1.06122, .12899, .00000, |
| 286 | + .00000, .05624, .16744, |
| 287 | + .36909, .19926, 1.28184, .13152, .00000, |
| 288 | + .00000, .06139, .17101, |
| 289 | + .37837, .20302, 1.44554, .13426, .00000, |
| 290 | + .00000, .06644, .17453, |
| 291 | + .38641, .20610, 1.57317, .13698, .00000, |
| 292 | + .00000, .07123, .17784, |
| 293 | + .39332, .20854, 1.67235, .13951, .00001, |
| 294 | + .00000, .07566, .18081, |
| 295 | + .39934, .21054, 1.74894, .14183, .00482, |
| 296 | + .00000, .07971, .18346 / |
| 297 | DATA ((F(I,20,K),I=1,5),F(7,20,K),F(6,20,K),F(8,20,K),K=1,7)/ |
| 298 | + .39417, .21269, .87132, .11000, .00000, |
| 299 | + .00000, .04597, .15304, |
| 300 | + .40857, .21843, 1.08735, .11080, .00000, |
| 301 | + .00000, .04998, .15496, |
| 302 | + .42039, .22263, 1.24495, .11220, .00000, |
| 303 | + .00000, .05416, .15715, |
| 304 | + .42999, .22574, 1.35369, .11380, .00000, |
| 305 | + .00000, .05818, .15920, |
| 306 | + .43796, .22808, 1.43230, .11540, .00000, |
| 307 | + .00000, .06191, .16105, |
| 308 | + .44452, .22971, 1.48820, .11679, .00001, |
| 309 | + .00000, .06523, .16258, |
| 310 | + .44997, .23090, 1.52687, .11804, .00257, |
| 311 | + -.00001, .06820, .16384 / |
| 312 | DATA ((F(I,21,K),I=1,5),F(7,21,K),F(6,21,K),F(8,21,K),K=1,7)/ |
| 313 | + .43740, .23195, .93597, .09467, .00000, |
| 314 | + .00000, .04051, .14281, |
| 315 | + .45273, .23699, 1.09539, .09482, .00000, |
| 316 | + .00000, .04393, .14355, |
| 317 | + .46465, .24034, 1.20287, .09550, .00000, |
| 318 | + .00000, .04744, .14447, |
| 319 | + .47385, .24250, 1.26983, .09636, .00000, |
| 320 | + .00000, .05068, .14526, |
| 321 | + .48116, .24382, 1.31262, .09719, .00000, |
| 322 | + .00000, .05359, .14591, |
| 323 | + .48675, .24447, 1.33802, .09787, .00001, |
| 324 | + .00000, .05612, .14629, |
| 325 | + .49116, .24467, 1.35083, .09841, .00096, |
| 326 | + .00000, .05830, .14647 / |
| 327 | DATA ((F(I,22,K),I=1,5),F(7,22,K),F(6,22,K),F(8,22,K),K=1,7)/ |
| 328 | + .48495, .25086, .98361, .07726, .00000, |
| 329 | + .00000, .03390, .12954, |
| 330 | + .50004, .25430, 1.08372, .07710, .00000, |
| 331 | + .00000, .03696, .12905, |
| 332 | + .51094, .25600, 1.14164, .07738, .00000, |
| 333 | + .00000, .03990, .12871, |
| 334 | + .51874, .25651, 1.16898, .07772, .00000, |
| 335 | + .00000, .04250, .12835, |
| 336 | + .52440, .25628, 1.17918, .07799, .00000, |
| 337 | + .00000, .04472, .12789, |
| 338 | + .52821, .25543, 1.17779, .07813, .00001, |
| 339 | + .00000, .04656, .12725, |
| 340 | + .53076, .25425, 1.16862, .07816, .00001, |
| 341 | + -.00001, .04809, .12650 / |
| 342 | DATA ((F(I,23,K),I=1,5),F(7,23,K),F(6,23,K),F(8,23,K),K=1,7)/ |
| 343 | + .52685, .26533, 1.00111, .06198, .00000, |
| 344 | + .00000, .02725, .11577, |
| 345 | + .54044, .26638, 1.05264, .06174, .00000, |
| 346 | + .00000, .03024, .11435, |
| 347 | + .54931, .26589, 1.07218, .06180, .00000, |
| 348 | + .00000, .03289, .11316, |
| 349 | + .55485, .26445, 1.07015, .06185, .00000, |
| 350 | + .00000, .03506, .11194, |
| 351 | + .55816, .26246, 1.05720, .06185, .00000, |
| 352 | + .00000, .03682, .11070, |
| 353 | + .55957, .26006, 1.03756, .06171, .00001, |
| 354 | + .00000, .03820, .10935, |
| 355 | + .55978, .25746, 1.01406, .06150, .00001, |
| 356 | + -.00003, .03931, .10794 / |
| 357 | DATA ((F(I,24,K),I=1,5),F(7,24,K),F(6,24,K),F(8,24,K),K=1,7)/ |
| 358 | + .56454, .27445, .99514, .04789, .00000, |
| 359 | + .00000, .02161, .10191, |
| 360 | + .57507, .27304, 1.00701, .04803, .00000, |
| 361 | + .00000, .02440, .09988, |
| 362 | + .58075, .27034, .99689, .04819, .00000, |
| 363 | + .00000, .02674, .09801, |
| 364 | + .58312, .26698, .97320, .04827, .00000, |
| 365 | + .00000, .02854, .09624, |
| 366 | + .58342, .26328, .94399, .04821, .00000, |
| 367 | + .00000, .02995, .09448, |
| 368 | + .58201, .25935, .91213, .04801, .00001, |
| 369 | + .00000, .03101, .09271, |
| 370 | + .57958, .25541, .87953, .04770, .00001, |
| 371 | + -.00001, .03178, .09099 / |
| 372 | DATA ((F(I,25,K),I=1,5),F(7,25,K),F(6,25,K),F(8,25,K),K=1,7)/ |
| 373 | + .59691, .27992, .97091, .03603, .00000, |
| 374 | + .00000, .01614, .08820, |
| 375 | + .60334, .27563, .95082, .03651, .00000, |
| 376 | + .00000, .01898, .08577, |
| 377 | + .60506, .27058, .91809, .03687, .00000, |
| 378 | + .00000, .02114, .08357, |
| 379 | + .60374, .26525, .87878, .03699, .00000, |
| 380 | + .00000, .02277, .08146, |
| 381 | + .60062, .25988, .83847, .03696, .00000, |
| 382 | + .00000, .02394, .07943, |
| 383 | + .59609, .25451, .79881, .03678, .00001, |
| 384 | + .00000, .02476, .07747, |
| 385 | + .59086, .24929, .76083, .03650, .00001, |
| 386 | + -.00002, .02533, .07557 / |
| 387 | DATA ((F(I,26,K),I=1,5),F(7,26,K),F(6,26,K),F(8,26,K),K=1,7)/ |
| 388 | + .62371, .28203, .93107, .02649, .00000, |
| 389 | + .00000, .01124, .07481, |
| 390 | + .62519, .27467, .88680, .02720, .00000, |
| 391 | + .00000, .01413, .07232, |
| 392 | + .62239, .26718, .83766, .02763, .00000, |
| 393 | + .00000, .01625, .07001, |
| 394 | + .61707, .25983, .78764, .02781, .00000, |
| 395 | + .00000, .01772, .06790, |
| 396 | + .61031, .25286, .74029, .02785, .00000, |
| 397 | + .00000, .01879, .06577, |
| 398 | + .60253, .24619, .69616, .02770, .00001, |
| 399 | + .00000, .01949, .06376, |
| 400 | + .59439, .23988, .65555, .02744, .00001, |
| 401 | + -.00002, .01993, .06188 / |
| 402 | DATA ((F(I,27,K),I=1,5),F(7,27,K),F(6,27,K),F(8,27,K),K=1,7)/ |
| 403 | + .64685, .27834, .87914, .01677, .00000, |
| 404 | + .00000, .00808, .06246, |
| 405 | + .64193, .26874, .81691, .01836, .00000, |
| 406 | + .00000, .01072, .06006, |
| 407 | + .63371, .25949, .75668, .01935, .00000, |
| 408 | + .00000, .01258, .05778, |
| 409 | + .62377, .25070, .70017, .01991, .00000, |
| 410 | + .00000, .01381, .05566, |
| 411 | + .61311, .24249, .64913, .02015, .00000, |
| 412 | + .00000, .01465, .05364, |
| 413 | + .60199, .23483, .60318, .02021, .00001, |
| 414 | + .00000, .01518, .05174, |
| 415 | + .59088, .22773, .56203, .02013, .00001, |
| 416 | + .00002, .01554, .04995 / |
| 417 | DATA ((F(I,28,K),I=1,5),F(7,28,K),F(6,28,K),F(8,28,K),K=1,7)/ |
| 418 | + .66107, .27515, .82378, .01211, .00000, |
| 419 | + .00000, .00374, .05101, |
| 420 | + .65020, .26245, .74720, .01349, .00000, |
| 421 | + .00000, .00665, .04889, |
| 422 | + .63701, .25088, .67896, .01434, .00000, |
| 423 | + .00000, .00861, .04683, |
| 424 | + .62284, .24040, .61850, .01480, .00000, |
| 425 | + .00000, .00991, .04486, |
| 426 | + .60857, .23090, .56585, .01502, .00000, |
| 427 | + .00000, .01076, .04303, |
| 428 | + .59439, .22225, .51975, .01504, .00001, |
| 429 | + .00000, .01129, .04132, |
| 430 | + .58066, .21436, .47938, .01495, .00001, |
| 431 | + -.00014, .01159, .03972 / |
| 432 | DATA ((F(I,29,K),I=1,5),F(7,29,K),F(6,29,K),F(8,29,K),K=1,7)/ |
| 433 | + .67179, .26585, .75762, .00606, .00000, |
| 434 | + .00000, .00178, .04082, |
| 435 | + .65387, .25121, .67378, .00798, .00000, |
| 436 | + .00000, .00448, .03902, |
| 437 | + .63509, .23825, .60220, .00918, .00000, |
| 438 | + .00000, .00622, .03725, |
| 439 | + .61643, .22673, .54096, .00988, .00000, |
| 440 | + .00000, .00733, .03555, |
| 441 | + .59846, .21645, .48891, .01024, .00000, |
| 442 | + .00000, .00805, .03397, |
| 443 | + .58120, .20725, .44429, .01048, .00001, |
| 444 | + .00000, .00851, .03240, |
| 445 | + .56495, .19892, .40590, .01053, .00001, |
| 446 | + -.00002, .00877, .03101 / |
| 447 | DATA ((F(I,30,K),I=1,5),F(7,30,K),F(6,30,K),F(8,30,K),K=1,7)/ |
| 448 | + .67527, .25529, .69187, .00245, .00000, |
| 449 | + .00000, .00001, .03191, |
| 450 | + .65063, .23870, .60286, .00452, .00000, |
| 451 | + .00000, .00236, .03056, |
| 452 | + .62661, .22441, .52998, .00578, .00000, |
| 453 | + .00000, .00405, .02908, |
| 454 | + .60372, .21204, .46960, .00656, .00000, |
| 455 | + .00000, .00510, .02766, |
| 456 | + .58244, .20113, .41941, .00701, .00000, |
| 457 | + .00000, .00575, .02631, |
| 458 | + .56258, .19146, .37721, .00720, .00001, |
| 459 | + .00000, .00615, .02511, |
| 460 | + .54411, .18290, .34143, .00732, .00001, |
| 461 | + -.00004, .00641, .02390 / |
| 462 | DATA ((F(I,31,K),I=1,5),F(7,31,K),F(6,31,K),F(8,31,K),K=1,7)/ |
| 463 | + .67251, .24277, .62540, .00010, .00000, |
| 464 | + .00000, .00001, .02445, |
| 465 | + .64126, .22455, .53420, .00214, .00000, |
| 466 | + .00000, .00078, .02338, |
| 467 | + .61229, .20924, .46215, .00332, .00000, |
| 468 | + .00000, .00236, .02228, |
| 469 | + .58571, .19616, .40404, .00407, .00000, |
| 470 | + .00000, .00338, .02110, |
| 471 | + .56145, .18486, .35676, .00452, .00000, |
| 472 | + .00000, .00399, .02008, |
| 473 | + .53927, .17498, .31763, .00477, .00001, |
| 474 | + .00000, .00440, .01905, |
| 475 | + .51898, .16628, .28494, .00492, .00001, |
| 476 | + .00000, .00462, .01810 / |
| 477 | DATA ((F(I,32,K),I=1,5),F(7,32,K),F(6,32,K),F(8,32,K),K=1,7)/ |
| 478 | + .66391, .22706, .55633, .00001, .00000, |
| 479 | + .00000, .00001, .01820, |
| 480 | + .62628, .20811, .46726, .00001, .00000, |
| 481 | + .00000, .00003, .01752, |
| 482 | + .59273, .19238, .39852, .00116, .00000, |
| 483 | + .00000, .00139, .01675, |
| 484 | + .56270, .17922, .34417, .00202, .00000, |
| 485 | + .00000, .00225, .01585, |
| 486 | + .53591, .16787, .30062, .00253, .00000, |
| 487 | + .00000, .00280, .01501, |
| 488 | + .51187, .15809, .26506, .00284, .00001, |
| 489 | + .00000, .00311, .01421, |
| 490 | + .49013, .14952, .23573, .00306, .00001, |
| 491 | + .00004, .00331, .01340 / |
| 492 | DATA ((F(I,33,K),I=1,5),F(7,33,K),F(6,33,K),F(8,33,K),K=1,7)/ |
| 493 | + .64879, .21198, .49256, .00001, .00000, |
| 494 | + .00000, .00001, .01307, |
| 495 | + .60557, .19218, .40616, .00001, .00000, |
| 496 | + .00000, .00001, .01281, |
| 497 | + .56813, .17607, .34118, .00036, .00000, |
| 498 | + .00000, .00038, .01218, |
| 499 | + .53536, .16266, .29090, .00107, .00000, |
| 500 | + .00000, .00120, .01162, |
| 501 | + .50661, .15134, .25128, .00154, .00000, |
| 502 | + .00000, .00172, .01096, |
| 503 | + .48111, .14168, .21942, .00184, .00001, |
| 504 | + .00000, .00204, .01034, |
| 505 | + .45837, .13328, .19346, .00198, .00001, |
| 506 | + .00001, .00221, .00988 / |
| 507 | DATA ((F(I,34,K),I=1,5),F(7,34,K),F(6,34,K),F(8,34,K),K=1,7)/ |
| 508 | + .62786, .19433, .42790, .00001, .00000, |
| 509 | + .00000, .00001, .00916, |
| 510 | + .57984, .17452, .34752, .00001, .00000, |
| 511 | + .00000, .00001, .00903, |
| 512 | + .53927, .15865, .28809, .00001, .00000, |
| 513 | + .00000, .00001, .00871, |
| 514 | + .50443, .14561, .24287, .00016, .00000, |
| 515 | + .00000, .00065, .00832, |
| 516 | + .47415, .13472, .20769, .00067, .00000, |
| 517 | + .00000, .00105, .00783, |
| 518 | + .44781, .12546, .17973, .00096, .00001, |
| 519 | + .00000, .00131, .00743, |
| 520 | + .42452, .11746, .15717, .00116, .00001, |
| 521 | + .00003, .00148, .00700 / |
| 522 | DATA ((F(I,35,K),I=1,5),F(7,35,K),F(6,35,K),F(8,35,K),K=1,7)/ |
| 523 | + .60128, .17617, .36867, .00001, .00000, |
| 524 | + .00000, .00001, .00600, |
| 525 | + .54953, .15670, .29445, .00001, .00000, |
| 526 | + .00000, .00001, .00617, |
| 527 | + .50651, .14136, .24076, .00001, .00000, |
| 528 | + .00000, .00001, .00596, |
| 529 | + .47024, .12886, .20061, .00001, .00000, |
| 530 | + .00000, .00024, .00577, |
| 531 | + .43920, .11844, .16982, .00005, .00000, |
| 532 | + .00000, .00059, .00551, |
| 533 | + .41244, .10970, .14560, .00035, .00001, |
| 534 | + .00000, .00083, .00521, |
| 535 | + .38903, .10226, .12625, .00054, .00001, |
| 536 | + .00004, .00095, .00490 / |
| 537 | DATA ((F(I,36,K),I=1,5),F(7,36,K),F(6,36,K),F(8,36,K),K=1,7)/ |
| 538 | + .57046, .15882, .31457, .00001, .00000, |
| 539 | + .00000, .00001, .00370, |
| 540 | + .51571, .13983, .24747, .00001, .00000, |
| 541 | + .00000, .00001, .00397, |
| 542 | + .47116, .12491, .19955, .00001, .00000, |
| 543 | + .00000, .00001, .00406, |
| 544 | + .43400, .11296, .16422, .00001, .00000, |
| 545 | + .00000, .00001, .00391, |
| 546 | + .40271, .10316, .13748, .00001, .00000, |
| 547 | + .00000, .00027, .00374, |
| 548 | + .37598, .09498, .11669, .00001, .00001, |
| 549 | + .00000, .00048, .00356, |
| 550 | + .35274, .08802, .10031, .00018, .00001, |
| 551 | + .00003, .00061, .00334 / |
| 552 | DATA ((F(I,37,K),I=1,5),F(7,37,K),F(6,37,K),F(8,37,K),K=1,7)/ |
| 553 | + .53556, .14173, .26561, .00001, .00000, |
| 554 | + .00000, .00001, .00219, |
| 555 | + .47905, .12344, .20521, .00001, .00000, |
| 556 | + .00000, .00001, .00262, |
| 557 | + .43373, .10934, .16317, .00001, .00000, |
| 558 | + .00000, .00001, .00266, |
| 559 | + .39671, .09815, .13264, .00001, .00000, |
| 560 | + .00000, .00001, .00257, |
| 561 | + .36558, .08902, .10991, .00001, .00000, |
| 562 | + .00000, .00001, .00250, |
| 563 | + .33924, .08151, .09246, .00007, .00001, |
| 564 | + .00000, .00015, .00237, |
| 565 | + .31668, .07519, .07885, .00019, .00001, |
| 566 | + -.00002, .00027, .00222 / |
| 567 | DATA ((F(I,38,K),I=1,5),F(7,38,K),F(6,38,K),F(8,38,K),K=1,7)/ |
| 568 | + .49683, .12381, .22084, .00001, .00000, |
| 569 | + .00000, .00001, .00113, |
| 570 | + .43944, .10672, .16749, .00001, .00000, |
| 571 | + .00000, .00001, .00150, |
| 572 | + .39428, .09371, .13118, .00001, .00000, |
| 573 | + .00000, .00001, .00168, |
| 574 | + .35761, .08355, .10540, .00001, .00000, |
| 575 | + .00000, .00001, .00159, |
| 576 | + .32742, .07528, .08641, .00001, .00000, |
| 577 | + .00000, .00001, .00161, |
| 578 | + .30203, .06857, .07204, .00001, .00001, |
| 579 | + .00000, .00009, .00150, |
| 580 | + .28045, .06294, .06092, .00001, .00001, |
| 581 | + .00002, .00018, .00141 / |
| 582 | DATA ((F(I,39,K),I=1,5),F(7,39,K),F(6,39,K),F(8,39,K),K=1,7)/ |
| 583 | + .45431, .10635, .17814, .00001, .00000, |
| 584 | + .00000, .00001, .00036, |
| 585 | + .39783, .09091, .13388, .00001, .00000, |
| 586 | + .00000, .00001, .00072, |
| 587 | + .35381, .07923, .10374, .00001, .00000, |
| 588 | + .00000, .00001, .00084, |
| 589 | + .31856, .07014, .08246, .00001, .00000, |
| 590 | + .00000, .00001, .00096, |
| 591 | + .28960, .06287, .06697, .00001, .00000, |
| 592 | + .00000, .00001, .00093, |
| 593 | + .26556, .05692, .05530, .00001, .00001, |
| 594 | + .00000, .00004, .00092, |
| 595 | + .24521, .05198, .04637, .00001, .00001, |
| 596 | + .00004, .00011, .00086 / |
| 597 | DATA ((F(I,40,K),I=1,5),F(7,40,K),F(6,40,K),F(8,40,K),K=1,7)/ |
| 598 | + .41254, .09148, .14592, .00001, .00000, |
| 599 | + .00000, .00001, .00003, |
| 600 | + .35692, .07724, .10695, .00001, .00000, |
| 601 | + .00000, .00001, .00033, |
| 602 | + .31428, .06668, .08140, .00001, .00000, |
| 603 | + .00000, .00001, .00047, |
| 604 | + .28057, .05850, .06381, .00001, .00000, |
| 605 | + .00000, .00001, .00053, |
| 606 | + .25332, .05206, .05122, .00001, .00000, |
| 607 | + .00000, .00001, .00056, |
| 608 | + .23091, .04683, .04190, .00001, .00001, |
| 609 | + .00000, .00001, .00054, |
| 610 | + .21210, .04253, .03482, .00001, .00001, |
| 611 | + .00000, .00001, .00053 / |
| 612 | DATA ((F(I,41,K),I=1,5),F(7,41,K),F(6,41,K),F(8,41,K),K=1,7)/ |
| 613 | + .36734, .07563, .11355, .00001, .00000, |
| 614 | + .00000, .00001, .00001, |
| 615 | + .31446, .06338, .08237, .00001, .00000, |
| 616 | + .00000, .00001, .00001, |
| 617 | + .27437, .05433, .06203, .00001, .00000, |
| 618 | + .00000, .00001, .00017, |
| 619 | + .24304, .04740, .04812, .00001, .00000, |
| 620 | + .00000, .00001, .00024, |
| 621 | + .21785, .04189, .03825, .00001, .00000, |
| 622 | + .00000, .00001, .00027, |
| 623 | + .19734, .03751, .03100, .00001, .00001, |
| 624 | + .00000, .00001, .00030, |
| 625 | + .18021, .03390, .02554, .00001, .00001, |
| 626 | + .00003, .00001, .00029 / |
| 627 | DATA ((F(I,42,K),I=1,5),F(7,42,K),F(6,42,K),F(8,42,K),K=1,7)/ |
| 628 | + .27911, .05064, .06747, .00001, .00000, |
| 629 | + .00000, .00001, .00001, |
| 630 | + .23354, .04145, .04718, .00001, .00000, |
| 631 | + .00000, .00001, .00001, |
| 632 | + .19979, .03487, .03442, .00001, .00000, |
| 633 | + .00000, .00001, .00001, |
| 634 | + .17401, .02990, .02599, .00001, .00000, |
| 635 | + .00000, .00001, .00005, |
| 636 | + .15364, .02606, .02017, .00001, .00000, |
| 637 | + .00000, .00001, .00006, |
| 638 | + .13734, .02302, .01600, .00001, .00001, |
| 639 | + .00000, .00001, .00008, |
| 640 | + .12395, .02056, .01294, .00001, .00001, |
| 641 | + .00000, .00001, .00008 / |
| 642 | DATA ((F(I,43,K),I=1,5),F(7,43,K),F(6,43,K),F(8,43,K),K=1,7)/ |
| 643 | + .19610, .03075, .03606, .00001, .00000, |
| 644 | + .00000, .00001, .00001, |
| 645 | + .16013, .02462, .02425, .00001, .00000, |
| 646 | + .00000, .00001, .00001, |
| 647 | + .13421, .02028, .01712, .00001, .00000, |
| 648 | + .00000, .00001, .00001, |
| 649 | + .11479, .01712, .01257, .00001, .00000, |
| 650 | + .00000, .00001, .00001, |
| 651 | + .09976, .01469, .00952, .00001, .00000, |
| 652 | + .00000, .00001, .00001, |
| 653 | + .08792, .01280, .00739, .00000, .00001, |
| 654 | + .00000, .00001, .00001, |
| 655 | + .07838, .01129, .00586, .00001, .00001, |
| 656 | + .00000, .00001, .00000 / |
| 657 | DATA ((F(I,44,K),I=1,5),F(7,44,K),F(6,44,K),F(8,44,K),K=1,7)/ |
| 658 | + .12471, .01656, .01668, .00001, .00000, |
| 659 | + .00000, .00001, .00001, |
| 660 | + .09913, .01295, .01080, .00001, .00000, |
| 661 | + .00000, .00001, .00001, |
| 662 | + .08124, .01044, .00737, .00001, .00000, |
| 663 | + .00000, .00001, .00001, |
| 664 | + .06814, .00865, .00525, .00001, .00000, |
| 665 | + .00000, .00001, .00001, |
| 666 | + .05819, .00728, .00387, .00000, .00000, |
| 667 | + .00000, .00001, .00000, |
| 668 | + .05052, .00627, .00294, .00001, .00001, |
| 669 | + .00000, .00001, .00001, |
| 670 | + .04438, .00545, .00227, .00001, .00001, |
| 671 | + .00000, .00001, .00001 / |
| 672 | DATA ((F(I,45,K),I=1,5),F(7,45,K),F(6,45,K),F(8,45,K),K=1,7)/ |
| 673 | + .06789, .00717, .00626, .00001, .00000, |
| 674 | + .00000, .00001, .00001, |
| 675 | + .05238, .00549, .00388, .00001, .00000, |
| 676 | + .00000, .00001, .00001, |
| 677 | + .04184, .00436, .00255, .00001, .00000, |
| 678 | + .00000, .00001, .00001, |
| 679 | + .03429, .00357, .00176, .00001, .00000, |
| 680 | + .00000, .00001, .00001, |
| 681 | + .02872, .00295, .00126, .00001, .00000, |
| 682 | + .00000, .00001, .00001, |
| 683 | + .02448, .00252, .00093, .00001, .00001, |
| 684 | + .00000, .00001, .00001, |
| 685 | + .02116, .00218, .00070, .00001, .00001, |
| 686 | + .00000, .00001, .00001 / |
| 687 | DATA ((F(I,46,K),I=1,5),F(7,46,K),F(6,46,K),F(8,46,K),K=1,7)/ |
| 688 | + .00856, .00063, .00037, .00001, .00000, |
| 689 | + .00000, .00001, .00001, |
| 690 | + .00610, .00043, .00017, .00001, .00000, |
| 691 | + .00000, .00001, .00001, |
| 692 | + .00453, .00032, .00010, .00001, .00000, |
| 693 | + .00000, .00001, .00001, |
| 694 | + .00349, .00023, .00006, .00001, .00000, |
| 695 | + .00000, .00001, .00001, |
| 696 | + .00277, .00019, .00004, .00001, .00000, |
| 697 | + .00000, .00001, .00001, |
| 698 | + .00226, .00014, .00003, .00001, .00001, |
| 699 | + .00000, .00001, .00001, |
| 700 | + .00187, .00011, .00002, .00001, .00001, |
| 701 | + .00000, .00001, .00001 / |
| 702 | DATA XX/1.D-5,2.D-5,4.D-5,6.D-5,8.D-5, |
| 703 | . 1.D-4,2.D-4,4.D-4,6.D-4,8.D-4, |
| 704 | . 1.D-3,2.D-3,4.D-3,6.D-3,8.D-3, |
| 705 | . 1.D-2,2.D-2,4.D-2,6.D-2,8.D-2, |
| 706 | . .1D0,.125D0,.15D0,.175D0,.2D0,.225D0,.25D0,.275D0, |
| 707 | . .3D0,.325D0,.35D0,.375D0,.4D0,.425D0,.45D0,.475D0, |
| 708 | . .5D0,.525D0,.55D0,.575D0,.6D0,.65D0,.7D0,.75D0, |
| 709 | . .8D0,.9D0,1.D0/ |
| 710 | DATA XMIN,XMAX,QSQMIN,QSQMAX/1.D-5,1.D0,0.625D0,5.D0/ |
| 711 | DATA N0/2,5,5,9,0,0,9,9/ |
| 712 | DATA ZEROD/0.D0/,PONED/0.1D0/,ONED/1.D0/,ONEDO/1.1D0/,TWOD/2.D0/ |
| 713 | DATA INIT/0/ |
| 714 | C don't let x be altered if it's out of range!! |
| 715 | xsave=x |
| 716 | C |
| 717 | IF(INIT.NE.0) GOTO 10 |
| 718 | INIT=1 |
| 719 | DO 20 N=1,NX-1 |
| 720 | DO 20 M=1,7 |
| 721 | C READ(36,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), |
| 722 | C . F(6,N,M),F(8,N,M) |
| 723 | C 1=UV 2=DV 3=GLUE 4=UBAR 5=CBAR 7=BBAR 6=SBAR 8=DBAR |
| 724 | DO 25 I=1,8 |
| 725 | 25 F(I,N,M)=F(I,N,M)/(ONED-XX(N))**N0(I) |
| 726 | 20 CONTINUE |
| 727 | DO 31 J=1,NTENTH-1 |
| 728 | XX(J)= LOG10(XX(J))+ONEDO |
| 729 | DO 31 I=1,8 |
| 730 | IF(I.EQ.7.OR.I.EQ.5) GO TO 31 |
| 731 | DO 30 K=1,7 |
| 732 | 30 F(I,J,K)= LOG(F(I,J,K))*F(I,NTENTH,K)/ LOG(F(I,NTENTH,K)) |
| 733 | 31 CONTINUE |
| 734 | 50 FORMAT(8F10.5) |
| 735 | DO 40 I=1,8 |
| 736 | DO 40 M=1,7 |
| 737 | 40 F(I,nx,M)=ZEROD |
| 738 | 10 CONTINUE |
| 739 | IF(X.LT.XMIN) X=XMIN |
| 740 | IF(X.GT.XMAX) X=XMAX |
| 741 | QSQ=SCALE**2 |
| 742 | IF(QSQ.LT.QSQMIN) QSQ=QSQMIN |
| 743 | IF(QSQ.GT.QSQMAX) QSQ=QSQMAX |
| 744 | XXX=X |
| 745 | IF(X.LT.PONED) XXX= LOG10(X)+ONEDO |
| 746 | N=0 |
| 747 | 70 N=N+1 |
| 748 | IF(XXX.GT.XX(N+1)) GOTO 70 |
| 749 | A=(XXX-XX(N))/(XX(N+1)-XX(N)) |
| 750 | RM= LOG(QSQ/QSQMIN)/ LOG(TWOD)*TWOD |
| 751 | B=RM-AINT(RM) |
| 752 | M=1+ INT(RM) |
| 753 | DO 60 I=1,8 |
| 754 | G(I)= (ONED-A)*(ONED-B)*F(I,N,M)+(ONED-A)*B*F(I,N,M+1) |
| 755 | . + A*(ONED-B)*F(I,N+1,M) + A*B*F(I,N+1,M+1) |
| 756 | IF(N.GE.ntenth) GOTO 65 |
| 757 | IF(I.EQ.7.or.i.eq.5) GOTO 65 |
| 758 | FAC=(ONED-B)*F(I,ntenth,M)+B*F(I,ntenth,M+1) |
| 759 | G(I)=FAC**(G(I)/FAC) |
| 760 | 65 CONTINUE |
| 761 | G(I)=G(I)*(ONED-X)**N0(I) |
| 762 | 60 CONTINUE |
| 763 | UPV=G(1) |
| 764 | DNV=G(2) |
| 765 | USEA=G(4) |
| 766 | DSEA=G(8) |
| 767 | STR=G(6) |
| 768 | CHM=G(5) |
| 769 | GLU=G(3) |
| 770 | BOT=G(7) |
| 771 | C restore x |
| 772 | x=xsave |
| 773 | C |
| 774 | RETURN |
| 775 | END |
git://git.uio.no / u / mrichter / AliRoot.git / blame