Blame view
fvn_fnlib/d9sifg.f
22.5 KB
38581db0c
git-svn-id: https...
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 |
subroutine d9sifg (x, f, g)
c december 1980 edition. w. fullerton, bell labs.
double precision x, f, g, f1cs(43), f2cs(99), g1cs(44),
1 g2cs(44), g3cs(56), xbnd, xbndg, xbig, xmaxf, xmaxg,
2 dcsevl, d1mach, dexp, dlog, dsqrt
external d1mach, dcsevl, initds
c
c series for f1 on the interval 2.00000e-02 to 6.25000e-02
c with weighted error 2.45e-32
c log weighted error 31.61
c significant figures required 30.42
c decimal places required 32.43
c
data f1 cs( 1) / -0.1191081969 0513636103 4820196582 8918d0/
data f1 cs( 2) / -0.0247823144 9962362475 9007415082 3133d0/
data f1 cs( 3) / 0.0011910281 4533578212 6812036305 4457d0/
data f1 cs( 4) / -0.0000927027 7143885617 4830860036 0706d0/
data f1 cs( 5) / 0.0000093373 1415682709 9686820458 2766d0/
data f1 cs( 6) / -0.0000011058 2878205571 4393897942 6306d0/
data f1 cs( 7) / 0.0000001464 7720714601 6216933655 0799d0/
data f1 cs( 8) / -0.0000000210 6944962876 8953260122 7548d0/
data f1 cs( 9) / 0.0000000032 2934923668 4823638285 7374d0/
data f1 cs( 10) / -0.0000000005 2065296175 2937582801 4986d0/
data f1 cs( 11) / 0.0000000000 8748788845 7027875026 8316d0/
data f1 cs( 12) / -0.0000000000 1521761870 5612366829 4574d0/
data f1 cs( 13) / 0.0000000000 0272571924 0541957390 0583d0/
data f1 cs( 14) / -0.0000000000 0050070530 7596855629 0255d0/
data f1 cs( 15) / 0.0000000000 0009402409 0272606851 1779d0/
data f1 cs( 16) / -0.0000000000 0001800144 4479180367 8336d0/
data f1 cs( 17) / 0.0000000000 0000350626 2143274178 5826d0/
data f1 cs( 18) / -0.0000000000 0000069352 8292676914 9709d0/
data f1 cs( 19) / 0.0000000000 0000013909 2513645421 6568d0/
data f1 cs( 20) / -0.0000000000 0000002824 8688507417 0585d0/
data f1 cs( 21) / 0.0000000000 0000000580 3130569357 9081d0/
data f1 cs( 22) / -0.0000000000 0000000120 4690157337 5820d0/
data f1 cs( 23) / 0.0000000000 0000000025 2505244365 5940d0/
data f1 cs( 24) / -0.0000000000 0000000005 3398026880 5594d0/
data f1 cs( 25) / 0.0000000000 0000000001 1385578627 4122d0/
data f1 cs( 26) / -0.0000000000 0000000000 2446286150 5259d0/
data f1 cs( 27) / 0.0000000000 0000000000 0529365932 0439d0/
data f1 cs( 28) / -0.0000000000 0000000000 0115318494 0277d0/
data f1 cs( 29) / 0.0000000000 0000000000 0025278656 8318d0/
data f1 cs( 30) / -0.0000000000 0000000000 0005573864 5378d0/
data f1 cs( 31) / 0.0000000000 0000000000 0001235824 5621d0/
data f1 cs( 32) / -0.0000000000 0000000000 0000275435 0842d0/
data f1 cs( 33) / 0.0000000000 0000000000 0000061690 6808d0/
data f1 cs( 34) / -0.0000000000 0000000000 0000013881 7443d0/
data f1 cs( 35) / 0.0000000000 0000000000 0000003137 5329d0/
data f1 cs( 36) / -0.0000000000 0000000000 0000000712 1249d0/
data f1 cs( 37) / 0.0000000000 0000000000 0000000162 2778d0/
data f1 cs( 38) / -0.0000000000 0000000000 0000000037 1206d0/
data f1 cs( 39) / 0.0000000000 0000000000 0000000008 5221d0/
data f1 cs( 40) / -0.0000000000 0000000000 0000000001 9633d0/
data f1 cs( 41) / 0.0000000000 0000000000 0000000000 4538d0/
data f1 cs( 42) / -0.0000000000 0000000000 0000000000 1052d0/
data f1 cs( 43) / 0.0000000000 0000000000 0000000000 0245d0/
c
c series for f2 on the interval 0.00000e+00 to 2.00000e-02
c with weighted error 2.38e-32
c log weighted error 31.62
c significant figures required 30.01
c decimal places required 32.62
c
data f2 cs( 1) / -0.0348409253 8970132330 8360497337 45577d0/
data f2 cs( 2) / -0.0166842205 6779596873 2467863122 78676d0/
data f2 cs( 3) / 0.0006752901 2412377385 0452078592 39727d0/
data f2 cs( 4) / -0.0000535066 6225447013 6287855775 57429d0/
data f2 cs( 5) / 0.0000062693 4217790075 2670507594 31626d0/
data f2 cs( 6) / -0.0000009526 6388019916 6806777904 14293d0/
data f2 cs( 7) / 0.0000001745 6292242509 8804255044 27666d0/
data f2 cs( 8) / -0.0000000368 7954030653 0933070976 46628d0/
data f2 cs( 9) / 0.0000000087 2026777051 3952640758 16938d0/
data f2 cs( 10) / -0.0000000022 6019703919 7387485304 23167d0/
data f2 cs( 11) / 0.0000000006 3246249765 2506125204 44877d0/
data f2 cs( 12) / -0.0000000001 8889118884 7178692409 11480d0/
data f2 cs( 13) / 0.0000000000 5967746729 9978133726 20472d0/
data f2 cs( 14) / -0.0000000000 1980443117 3722390111 96007d0/
data f2 cs( 15) / 0.0000000000 0686413954 7721033837 13264d0/
data f2 cs( 16) / -0.0000000000 0247310193 0701991060 74890d0/
data f2 cs( 17) / 0.0000000000 0092263594 5499414041 96042d0/
data f2 cs( 18) / -0.0000000000 0035523634 9992617844 97297d0/
data f2 cs( 19) / 0.0000000000 0014076049 6253515914 61820d0/
data f2 cs( 20) / -0.0000000000 0005726228 4997476527 94311d0/
data f2 cs( 21) / 0.0000000000 0002386537 5454131718 10106d0/
data f2 cs( 22) / -0.0000000000 0001017141 8907645971 42232d0/
data f2 cs( 23) / 0.0000000000 0000442594 5310783644 24968d0/
data f2 cs( 24) / -0.0000000000 0000196344 9330491897 61979d0/
data f2 cs( 25) / 0.0000000000 0000088688 7483148104 61024d0/
data f2 cs( 26) / -0.0000000000 0000040743 3450273115 46948d0/
data f2 cs( 27) / 0.0000000000 0000019016 8372156753 39859d0/
data f2 cs( 28) / -0.0000000000 0000009009 7072974780 42442d0/
data f2 cs( 29) / 0.0000000000 0000004329 2112740956 68667d0/
data f2 cs( 30) / -0.0000000000 0000002108 1444653224 79526d0/
data f2 cs( 31) / 0.0000000000 0000001039 6379070264 52274d0/
data f2 cs( 32) / -0.0000000000 0000000518 8910079489 31936d0/
data f2 cs( 33) / 0.0000000000 0000000261 9553248698 99371d0/
data f2 cs( 34) / -0.0000000000 0000000133 6903999513 01570d0/
data f2 cs( 35) / 0.0000000000 0000000068 9410577029 31664d0/
data f2 cs( 36) / -0.0000000000 0000000035 9053626104 37250d0/
data f2 cs( 37) / 0.0000000000 0000000018 8780772557 91706d0/
data f2 cs( 38) / -0.0000000000 0000000010 0161252655 94380d0/
data f2 cs( 39) / 0.0000000000 0000000005 3607256915 78228d0/
data f2 cs( 40) / -0.0000000000 0000000002 8931989749 44827d0/
data f2 cs( 41) / 0.0000000000 0000000001 5740651002 02625d0/
data f2 cs( 42) / -0.0000000000 0000000000 8630271064 31206d0/
data f2 cs( 43) / 0.0000000000 0000000000 4767156028 62288d0/
data f2 cs( 44) / -0.0000000000 0000000000 2652227399 98504d0/
data f2 cs( 45) / 0.0000000000 0000000000 1485828650 63866d0/
data f2 cs( 46) / -0.0000000000 0000000000 0837972359 23135d0/
data f2 cs( 47) / 0.0000000000 0000000000 0475659164 22711d0/
data f2 cs( 48) / -0.0000000000 0000000000 0271690733 53112d0/
data f2 cs( 49) / 0.0000000000 0000000000 0156127388 81686d0/
data f2 cs( 50) / -0.0000000000 0000000000 0090245550 78347d0/
data f2 cs( 51) / 0.0000000000 0000000000 0052460970 49119d0/
data f2 cs( 52) / -0.0000000000 0000000000 0030664508 18697d0/
data f2 cs( 53) / 0.0000000000 0000000000 0018019962 50957d0/
data f2 cs( 54) / -0.0000000000 0000000000 0010644430 50752d0/
data f2 cs( 55) / 0.0000000000 0000000000 0006319421 58881d0/
data f2 cs( 56) / -0.0000000000 0000000000 0003770138 12246d0/
data f2 cs( 57) / 0.0000000000 0000000000 0002259975 42918d0/
data f2 cs( 58) / -0.0000000000 0000000000 0001361008 44814d0/
data f2 cs( 59) / 0.0000000000 0000000000 0000823332 32003d0/
data f2 cs( 60) / -0.0000000000 0000000000 0000500259 86091d0/
data f2 cs( 61) / 0.0000000000 0000000000 0000305262 45684d0/
data f2 cs( 62) / -0.0000000000 0000000000 0000187051 64021d0/
data f2 cs( 63) / 0.0000000000 0000000000 0000115084 04393d0/
data f2 cs( 64) / -0.0000000000 0000000000 0000071087 14611d0/
data f2 cs( 65) / 0.0000000000 0000000000 0000044080 65533d0/
data f2 cs( 66) / -0.0000000000 0000000000 0000027437 60867d0/
data f2 cs( 67) / 0.0000000000 0000000000 0000017141 44851d0/
data f2 cs( 68) / -0.0000000000 0000000000 0000010747 68860d0/
data f2 cs( 69) / 0.0000000000 0000000000 0000006762 59777d0/
data f2 cs( 70) / -0.0000000000 0000000000 0000004269 81348d0/
data f2 cs( 71) / 0.0000000000 0000000000 0000002705 00637d0/
data f2 cs( 72) / -0.0000000000 0000000000 0000001719 33331d0/
data f2 cs( 73) / 0.0000000000 0000000000 0000001096 36138d0/
data f2 cs( 74) / -0.0000000000 0000000000 0000000701 32573d0/
data f2 cs( 75) / 0.0000000000 0000000000 0000000450 01784d0/
data f2 cs( 76) / -0.0000000000 0000000000 0000000289 63835d0/
data f2 cs( 77) / 0.0000000000 0000000000 0000000186 97009d0/
data f2 cs( 78) / -0.0000000000 0000000000 0000000121 04646d0/
data f2 cs( 79) / 0.0000000000 0000000000 0000000078 59065d0/
data f2 cs( 80) / -0.0000000000 0000000000 0000000051 16867d0/
data f2 cs( 81) / 0.0000000000 0000000000 0000000033 40627d0/
data f2 cs( 82) / -0.0000000000 0000000000 0000000021 86851d0/
data f2 cs( 83) / 0.0000000000 0000000000 0000000014 35340d0/
data f2 cs( 84) / -0.0000000000 0000000000 0000000009 44523d0/
data f2 cs( 85) / 0.0000000000 0000000000 0000000006 23117d0/
data f2 cs( 86) / -0.0000000000 0000000000 0000000004 12101d0/
data f2 cs( 87) / 0.0000000000 0000000000 0000000002 73208d0/
data f2 cs( 88) / -0.0000000000 0000000000 0000000001 81558d0/
data f2 cs( 89) / 0.0000000000 0000000000 0000000001 20934d0/
data f2 cs( 90) / -0.0000000000 0000000000 0000000000 80737d0/
data f2 cs( 91) / 0.0000000000 0000000000 0000000000 54022d0/
data f2 cs( 92) / -0.0000000000 0000000000 0000000000 36227d0/
data f2 cs( 93) / 0.0000000000 0000000000 0000000000 24348d0/
data f2 cs( 94) / -0.0000000000 0000000000 0000000000 16401d0/
data f2 cs( 95) / 0.0000000000 0000000000 0000000000 11074d0/
data f2 cs( 96) / -0.0000000000 0000000000 0000000000 07497d0/
data f2 cs( 97) / 0.0000000000 0000000000 0000000000 05091d0/
data f2 cs( 98) / -0.0000000000 0000000000 0000000000 03470d0/
data f2 cs( 99) / 0.0000000000 0000000000 0000000000 02377d0/
c
c series for g1 on the interval 2.00000e-02 to 6.25000e-02
c with weighted error 7.23e-32
c log weighted error 31.14
c significant figures required 30.35
c decimal places required 31.96
c
data g1 cs( 1) / -0.3040578798 2534959544 9972668209 1083d0/
data g1 cs( 2) / -0.0566890984 5971205877 3133915611 8269d0/
data g1 cs( 3) / 0.0039046158 1732756439 1998407155 4082d0/
data g1 cs( 4) / -0.0003746075 9592022606 1861933986 7489d0/
data g1 cs( 5) / 0.0000435431 5565598436 7955222084 0065d0/
data g1 cs( 6) / -0.0000057417 2944530250 4656197072 3475d0/
data g1 cs( 7) / 0.0000008282 5521045026 2974193761 6492d0/
data g1 cs( 8) / -0.0000001278 2458925946 4272788391 3223d0/
data g1 cs( 9) / 0.0000000207 9783529486 8788443925 7529d0/
data g1 cs( 10) / -0.0000000035 3132059219 9079804203 2682d0/
data g1 cs( 11) / 0.0000000006 2108242363 0895106863 1449d0/
data g1 cs( 12) / -0.0000000001 1252154744 4629264933 6987d0/
data g1 cs( 13) / 0.0000000000 2090889176 8442160526 7019d0/
data g1 cs( 14) / -0.0000000000 0397158317 3768172768 9158d0/
data g1 cs( 15) / 0.0000000000 0076904313 1427208993 9005d0/
data g1 cs( 16) / -0.0000000000 0015146967 4273161351 9826d0/
data g1 cs( 17) / 0.0000000000 0003028921 4655235968 4119d0/
data g1 cs( 18) / -0.0000000000 0000613997 0383470882 5400d0/
data g1 cs( 19) / 0.0000000000 0000126006 0582951093 3553d0/
data g1 cs( 20) / -0.0000000000 0000026150 2925093948 3683d0/
data g1 cs( 21) / 0.0000000000 0000005482 7884489179 6821d0/
data g1 cs( 22) / -0.0000000000 0000001160 3818212952 6571d0/
data g1 cs( 23) / 0.0000000000 0000000247 7165410712 9795d0/
data g1 cs( 24) / -0.0000000000 0000000053 3067275322 3389d0/
data g1 cs( 25) / 0.0000000000 0000000011 5566607559 8465d0/
data g1 cs( 26) / -0.0000000000 0000000002 5228054774 4957d0/
data g1 cs( 27) / 0.0000000000 0000000000 5542903855 0786d0/
data g1 cs( 28) / -0.0000000000 0000000000 1225220842 1297d0/
data g1 cs( 29) / 0.0000000000 0000000000 0272366431 8684d0/
data g1 cs( 30) / -0.0000000000 0000000000 0060870783 1422d0/
data g1 cs( 31) / 0.0000000000 0000000000 0013672487 4476d0/
data g1 cs( 32) / -0.0000000000 0000000000 0003085662 6806d0/
data g1 cs( 33) / 0.0000000000 0000000000 0000699521 2319d0/
data g1 cs( 34) / -0.0000000000 0000000000 0000159258 7569d0/
data g1 cs( 35) / 0.0000000000 0000000000 0000036405 1056d0/
data g1 cs( 36) / -0.0000000000 0000000000 0000008353 9465d0/
data g1 cs( 37) / 0.0000000000 0000000000 0000001924 0303d0/
data g1 cs( 38) / -0.0000000000 0000000000 0000000444 6816d0/
data g1 cs( 39) / 0.0000000000 0000000000 0000000103 1182d0/
data g1 cs( 40) / -0.0000000000 0000000000 0000000023 9887d0/
data g1 cs( 41) / 0.0000000000 0000000000 0000000005 5976d0/
data g1 cs( 42) / -0.0000000000 0000000000 0000000001 3100d0/
data g1 cs( 43) / 0.0000000000 0000000000 0000000000 3074d0/
data g1 cs( 44) / -0.0000000000 0000000000 0000000000 0723d0/
c
c series for g2 on the interval 5.00000e-03 to 2.00000e-02
c with weighted error 3.25e-32
c log weighted error 31.49
c significant figures required 30.32
c decimal places required 32.31
c
data g2 cs( 1) / -0.1211802894 7316462635 4183404685 8267d0/
data g2 cs( 2) / -0.0316761386 3949502867 0140792350 5610d0/
data g2 cs( 3) / 0.0013383199 7788626801 6381942949 2182d0/
data g2 cs( 4) / -0.0000895511 0113922524 2553190506 9518d0/
data g2 cs( 5) / 0.0000079155 5629617182 1311524946 7924d0/
data g2 cs( 6) / -0.0000008438 7933222415 2018141898 2080d0/
data g2 cs( 7) / 0.0000001029 9804256775 3014664722 7274d0/
data g2 cs( 8) / -0.0000000139 2957506051 8383579583 4444d0/
data g2 cs( 9) / 0.0000000020 4227039598 7598040067 7594d0/
data g2 cs( 10) / -0.0000000003 1965346942 0642703543 4752d0/
data g2 cs( 11) / 0.0000000000 5281478326 5726769861 5312d0/
data g2 cs( 12) / -0.0000000000 0913395546 7267103373 5289d0/
data g2 cs( 13) / 0.0000000000 0164262512 3896776044 4819d0/
data g2 cs( 14) / -0.0000000000 0030558970 3932266000 2410d0/
data g2 cs( 15) / 0.0000000000 0005856558 2578577971 7892d0/
data g2 cs( 16) / -0.0000000000 0001152291 9773094012 0563d0/
data g2 cs( 17) / 0.0000000000 0000232094 6911998853 7310d0/
data g2 cs( 18) / -0.0000000000 0000047743 5583417753 5025d0/
data g2 cs( 19) / 0.0000000000 0000010009 9676580018 0573d0/
data g2 cs( 20) / -0.0000000000 0000002135 3377808225 6704d0/
data g2 cs( 21) / 0.0000000000 0000000462 7719077736 7671d0/
data g2 cs( 22) / -0.0000000000 0000000101 7580741022 7657d0/
data g2 cs( 23) / 0.0000000000 0000000022 6765739988 4672d0/
data g2 cs( 24) / -0.0000000000 0000000005 1163077607 6426d0/
data g2 cs( 25) / 0.0000000000 0000000001 1676701491 3108d0/
data g2 cs( 26) / -0.0000000000 0000000000 2693542767 2470d0/
data g2 cs( 27) / 0.0000000000 0000000000 0627566584 1146d0/
data g2 cs( 28) / -0.0000000000 0000000000 0147588055 7531d0/
data g2 cs( 29) / 0.0000000000 0000000000 0035014531 4739d0/
data g2 cs( 30) / -0.0000000000 0000000000 0008375773 2152d0/
data g2 cs( 31) / 0.0000000000 0000000000 0002019181 5152d0/
data g2 cs( 32) / -0.0000000000 0000000000 0000490356 7705d0/
data g2 cs( 33) / 0.0000000000 0000000000 0000119912 3348d0/
data g2 cs( 34) / -0.0000000000 0000000000 0000029517 0610d0/
data g2 cs( 35) / 0.0000000000 0000000000 0000007311 3112d0/
data g2 cs( 36) / -0.0000000000 0000000000 0000001821 7843d0/
data g2 cs( 37) / 0.0000000000 0000000000 0000000456 5148d0/
data g2 cs( 38) / -0.0000000000 0000000000 0000000115 0151d0/
data g2 cs( 39) / 0.0000000000 0000000000 0000000029 1267d0/
data g2 cs( 40) / -0.0000000000 0000000000 0000000007 4125d0/
data g2 cs( 41) / 0.0000000000 0000000000 0000000001 8953d0/
data g2 cs( 42) / -0.0000000000 0000000000 0000000000 4868d0/
data g2 cs( 43) / 0.0000000000 0000000000 0000000000 1256d0/
data g2 cs( 44) / -0.0000000000 0000000000 0000000000 0325d0/
c
c series for g3 on the interval 0.00000e+00 to 5.00000e-03
c with weighted error 3.83e-32
c log weighted error 31.42
c significant figures required 29.71
c decimal places required 32.29
c
data g3 cs( 1) / -0.0280574367 8094729284 0281526433 5299d0/
data g3 cs( 2) / -0.0137271597 1622369754 0910050808 9556d0/
data g3 cs( 3) / 0.0002894032 6387602960 2744894127 3751d0/
data g3 cs( 4) / -0.0000114129 2393911971 4590874362 2517d0/
data g3 cs( 5) / 0.0000006813 9655907262 4299772020 7302d0/
data g3 cs( 6) / -0.0000000547 9522896046 5236366905 8052d0/
data g3 cs( 7) / 0.0000000055 2074299182 1252910940 6521d0/
data g3 cs( 8) / -0.0000000006 6414641993 2292002249 1428d0/
data g3 cs( 9) / 0.0000000000 9223736634 8704110856 4960d0/
data g3 cs( 10) / -0.0000000000 1442990888 8668286261 1718d0/
data g3 cs( 11) / 0.0000000000 0249639048 9203071024 8705d0/
data g3 cs( 12) / -0.0000000000 0047082406 7587524472 2971d0/
data g3 cs( 13) / 0.0000000000 0009572176 5921675998 8140d0/
data g3 cs( 14) / -0.0000000000 0002078899 6609580903 0537d0/
data g3 cs( 15) / 0.0000000000 0000478750 9997087743 1627d0/
data g3 cs( 16) / -0.0000000000 0000116190 7058337717 3759d0/
data g3 cs( 17) / 0.0000000000 0000029565 0896926783 6974d0/
data g3 cs( 18) / -0.0000000000 0000007852 9498825649 2025d0/
data g3 cs( 19) / 0.0000000000 0000002169 2226436825 6612d0/
data g3 cs( 20) / -0.0000000000 0000000621 1351583167 6342d0/
data g3 cs( 21) / 0.0000000000 0000000183 8456883845 0977d0/
data g3 cs( 22) / -0.0000000000 0000000056 1088748213 7276d0/
data g3 cs( 23) / 0.0000000000 0000000017 6186280528 0062d0/
data g3 cs( 24) / -0.0000000000 0000000005 6811105054 1451d0/
data g3 cs( 25) / 0.0000000000 0000000001 8778627958 2313d0/
data g3 cs( 26) / -0.0000000000 0000000000 6353169415 1124d0/
data g3 cs( 27) / 0.0000000000 0000000000 2196880236 8238d0/
data g3 cs( 28) / -0.0000000000 0000000000 0775466655 0395d0/
data g3 cs( 29) / 0.0000000000 0000000000 0279101835 6581d0/
data g3 cs( 30) / -0.0000000000 0000000000 0102317852 5247d0/
data g3 cs( 31) / 0.0000000000 0000000000 0038169340 3919d0/
data g3 cs( 32) / -0.0000000000 0000000000 0014476789 5606d0/
data g3 cs( 33) / 0.0000000000 0000000000 0005577951 2634d0/
data g3 cs( 34) / -0.0000000000 0000000000 0002181723 9071d0/
data g3 cs( 35) / 0.0000000000 0000000000 0000865664 6309d0/
data g3 cs( 36) / -0.0000000000 0000000000 0000348215 7895d0/
data g3 cs( 37) / 0.0000000000 0000000000 0000141918 8130d0/
data g3 cs( 38) / -0.0000000000 0000000000 0000058571 4314d0/
data g3 cs( 39) / 0.0000000000 0000000000 0000024466 0482d0/
data g3 cs( 40) / -0.0000000000 0000000000 0000010338 7099d0/
data g3 cs( 41) / 0.0000000000 0000000000 0000004417 7299d0/
data g3 cs( 42) / -0.0000000000 0000000000 0000001908 0079d0/
data g3 cs( 43) / 0.0000000000 0000000000 0000000832 6038d0/
data g3 cs( 44) / -0.0000000000 0000000000 0000000366 9553d0/
data g3 cs( 45) / 0.0000000000 0000000000 0000000163 2875d0/
data g3 cs( 46) / -0.0000000000 0000000000 0000000073 3357d0/
data g3 cs( 47) / 0.0000000000 0000000000 0000000033 2327d0/
data g3 cs( 48) / -0.0000000000 0000000000 0000000015 1906d0/
data g3 cs( 49) / 0.0000000000 0000000000 0000000007 0020d0/
data g3 cs( 50) / -0.0000000000 0000000000 0000000003 2539d0/
data g3 cs( 51) / 0.0000000000 0000000000 0000000001 5240d0/
data g3 cs( 52) / -0.0000000000 0000000000 0000000000 7193d0/
data g3 cs( 53) / 0.0000000000 0000000000 0000000000 3420d0/
data g3 cs( 54) / -0.0000000000 0000000000 0000000000 1638d0/
data g3 cs( 55) / 0.0000000000 0000000000 0000000000 0790d0/
data g3 cs( 56) / -0.0000000000 0000000000 0000000000 0383d0/
c
data nf1, nf2, ng1, ng2, ng3 / 5*0 /
data xbnd, xbndg, xbig, xmaxf, xmaxg / 5*0.0d0 /
c
if (nf1.ne.0) go to 10
tol = 0.1d0*d1mach(3)
nf1 = initds (f1cs, 43, tol)
nf2 = initds (f2cs, 99, tol)
ng1 = initds (g1cs, 44, tol)
ng2 = initds (g2cs, 44, tol)
ng3 = initds (g3cs, 56, tol)
c
xbig = dsqrt(1.0d0/d1mach(3))
xmaxf = dexp (dmin1 (-dlog(d1mach(1)), dlog(d1mach(2))) - .01d0)
xmaxg = 1.0d0/dsqrt(d1mach(1))
xbnd = dsqrt (50.0d0)
xbndg = dsqrt (200.d0)
c
10 if (x.lt.4.0d0) call seteru (
1 34hd9sifg approxs invalid for x lt 4, 34, 1, 2)
c
if (x.gt.xbnd) go to 20
f = (1.0d0 + dcsevl ((1.d0/x**2-0.04125d0)/.02125d0, f1cs, nf1))/x
g = (1.0d0 + dcsevl((1.d0/x**2-.04125d0)/.02125d0, g1cs,ng1))/x**2
return
c
20 if (x.gt.xbig) go to 30
f = (1.0d0 + dcsevl (100.d0/x**2-1.d0, f2cs, nf2))/x
if (x.le.xbndg) g = (1.0d0 + dcsevl ((10000.d0/x**2-125.d0)/75.d0,
1 g2cs, ng2))/x**2
if (x.gt.xbndg) g = (1.0d0 + dcsevl (400.d0/x**2-1.d0, g3cs,
1 ng3))/x**2
return
c
30 f = 0.d0
if (x.lt.xmaxf) f = 1.0d0/x
g = 0.d0
if (x.lt.xmaxg) g = 1.0d0/x**2
return
c
end
|