29
30#include "implicit_f.inc"
31
32
33
34
35
36
37
38 infinity = 1e20
39 zero = 0.
40 one = 1.
41 two = 2.
42 three = 3.
43 four = 4.
44 five = 5.
45 six = 6.
46 seven = 7.
47 eight = 8.
48 nine = 9.
49 ten = 10.
50 eleven = 11.
51 twelve = 12.
52 fourteen = 14.
53 fifteen = 15.
54 sixteen = 16.
55 twenty = 20.
56 twenty4 = 24.
57 twenty5 = 25.
58 twenty7 = 27.
59 thirty = 30.
60 thirty2 = 32.
61 thirty6 = 36.
62 fourty = 40.
63 fourty5 = 45.
64 fourty8 = 48.
65 fifty = 50.
66 sixty4 = 64.
67 eighty = 80.
68 ninty = 90.
69 eighty16 = 96.
70 eighty19 = 99.
71 hundred = 100.
72 hundred28 = 128.
73 hundred50 = 150.
74 hundred80 = 180.
75 three100 = 300.
76 twohundred25 = 225.
77 twohundred95 = 295.
78 twohundred98 = 298.
79 twohundred93 = 293.
80 threehundred24 = 324.
81 fivehundred12 = 512.
82 three1000 = 3000.
83
84
85
86 ep02 = hundred
87 ep03 = ep02 * ten
88 ep04 = ep03 * ten
89 ep05 = ep04 * ten
90 ep06 = ep05 * ten
91 ep07 = ep06 * ten
92 ep08 = ep07 * ten
93 ep09 = ep08 * ten
94 ep10 = ep09 * ten
95 ep11 = ep10 * ten
96 ep12 = ep11 * ten
97 ep13 = ep12 * ten
98 ep14 = ep13 * ten
99 ep15 = ep14 * ten
100 ep16 = ep15 * ten
101 ep17 = ep16 * ten
102 ep18 = ep17 * ten
103 ep19 = ep18 * ten
104 ep20 = ep19 * ten
105 ep29 = ep20 * ep09
106 ep30 = ep20 * ep10
107 ep31 = ep30 * ten
108 ep35 = ep30 * ep05
109 ep36 = ep30 * ep06
110 ep37 = ep30 * ep07
111
112
113
114 em01 = one / ten
115 em02 = one / ep02
116 em03 = one / ep03
117 em04 = one / ep04
118 em05 = one / ep05
119 em06 = one / ep06
120 em07 = one / ep07
121 em08 = one / ep08
122 em09 = one / ep09
123 em10 = one / ep10
124 em11 = one / ep11
125 em12 = one / ep12
126 em13 = one / ep13
127 em14 = one / ep14
128 em15 = one / ep15
129 em16 = one / ep16
130 em17 = one / ep17
131 em18 = one / ep18
132 em19 = one / ep19
133 em20 = one / ep20
134 em30 = one / ep30
135 em35 = one / ep35
136
137
138
139 zep1 = em01
140 zep2 = two / ten
141 zep3 = three / ten
142 zep4 = four / ten
143 zep5 = five / ten
144 zep6 = one - zep4
145 zep7 = one - zep3
146 zep8 = one - zep2
147 zep9 = one - zep1
148
149 zep01 = em02
150 zep02 = two / ep02
151 zep03 = three / ep02
152 zep04 = four / ep02
153 zep05 = five / ep02
154 zep06 = zep1 - zep04
155 zep07 = zep1 - zep03
156 zep08 = zep1 - zep02
157 zep09 = zep1 - zep01
158
159 twoem3 = two / ep03
160 threeem3 = three / ep03
161 fourem3 = four / ep03
162 fiveem3 = five / ep03
163 sixem3 = zep01 - fourem3
164 sevenem3 = zep01 - threeem3
165 eightem3 = zep01 - twoem3
166 nineem3 = zep01 - em03
167
168 twoem4 = two / ep04
169 threeem4 = three / ep04
170 fourem4 = four / ep04
171 fiveem4 = five / ep04
172 sixem4 = em03 - fourem4
173 sevenem4 = em03 - threeem4
174 eightem4 = em03 - twoem4
175 nineem4 = em03 - em04
176
177 twoem5 = two / ep05
178 threeem5 = three / ep05
179 fourem5 = four / ep05
180 fiveem5 = five / ep05
181 sixem5 = em04 - fourem5
182 sevenem5 = em04 - threeem5
183 eightem5 = em04 - twoem5
184 nineem5 = em04 - em05
185
186 twoem6 = two / ep06
187 threeem6 = three / ep06
188 fourem6 = four / ep06
189 fiveem6 = five / ep06
190 sixem6 = em05 - fourem6
191 sevenem6 = em05 - threeem6
192 eightem6 = em05 - twoem6
193 nineem6 = em05 - em06
194
195
196
197 twoep30 = two * ep30
198 threeep30 = three * ep30
199 fourep30 = four * ep30
200 zep9ep30 = nine *ep29
201
202
203
204 half = one/two
205 third = one/three
206 fourth = one/four
207 one_fifth = one/five
208 one_over_6 = one/six
209 one_over_8 = one/eight
210 one_over_9 = one/nine
211 one_over_12 = one/twelve
212 one_over_16 = one/sixteen
213 one_over_64 = one/sixty4
214 one_over_512 = one/fivehundred12
215 two_third = two/three
216 three_half = three/two
217 three_over_4 = three/four
218 trhee_over_14 = three/fourteen
219 four_over_3 = four/three
220 four_over_5 = four/five
221 five_over_3 = five/three
222 five_over_4 = five/four
223 five_over_6 = five/six
224 six_over_5 = six/five
225 seven_half = seven/two
226 seven_over_5 = seven/five
227 seven_over_9 = seven/nine
228 nine_over_2 = nine/two
229 nine_over_20 = nine/twenty
230
231
232
233 zep11 = zep1 + zep01
234 zep12 = zep1 + zep02
235 zep13 = zep1 + zep03
236 zep14 = zep1 + zep04
237 zep15 = zep1 + zep05
238 zep16 = zep1 + zep06
239 zep17 = zep1 + zep07
240 zep18 = zep1 + zep08
241 zep19 = zep1 + zep09
242
243 zep21 = zep2 + zep01
244 zep23 = zep2 + zep03
245 zep25 = zep2 + zep05
246 zep26 = zep2 + zep06
247 zep213 = zep21 + threeem3
248 zep29 = zep2 + zep09
249
250 zep33 = zep3 + zep03
251 zep34 = zep3 + zep04
252 zep35 = zep3 + zep05
253 zep37 = zep3 + zep07
254 zep38 = zep3 + zep08
255 zep333 = zep33 + threeem3
256 zep373 = zep37 + threeem3
257
258 zep41 = zep4 + zep01
259 zep4187 = zep41 + eightem3 + sevenem4
260 zep43 = zep4 + zep03
261 zep45 = zep4 + zep05
262 zep46 = zep4 + zep06
263 zep47 = zep4 + zep07
264 zep475 = zep4 + zep07 + fiveem3
265 zep499 = zep4 + zep09 + nineem3
266
267 zep62 = one - zep38
268 zep65 = one - zep35
269 zep66 = one - zep34
270 zep67 = one - zep33
271
272 zep71 = one - zep29
273 zep74 = one - zep26
274 zep707 = one - zep3 + sevenem3
275
276 zep81 = one - zep19
277 zep82 = one - zep18
278 zep84 = one - zep16
279 zep85 = one - zep15
280 zep87 = one - zep13
281 zep89 = one - zep11
282
283 zep92 = one - zep08
284 zep93 = one - zep07
285 zep94 = one - zep06
286 zep95 = one - zep05
287 zep96 = one - zep04
288 zep97 = one - zep03
289 zep98 = one - zep02
290 zep99 = one - zep01
291 zep995 = one - fiveem3
292 zep997 = one - threeem3
293 zep999 = one - em03
294 zep9999 = one - em04
295 zep99999 = one - em05
296 zep931 = zep93 + em03
297
298 zep015 = zep15 / ten
299 zep015625 = one / sixty4
300 zep025 = one / fourty
301 zep037 = zep03 + sevenem3
302 zep072169 = zep07 + twoem3 + em04 + sixem5 + nineem6
303 zep075 = three / fourty
304
305 onep01 = one + zep01
306 onep001 = one + em03
307 onep0001 = one + em04
308 onep00001 = one + em05
309 onep02 = one + zep02
310 onep03 = one + zep03
311 onep04 = one + zep04
312 onep05 = one + zep05
313 onep06 = one + zep06
314 onep08 = one + zep08
315 onep09 = one + zep09
316
317 onep1 = one + zep1
318 onep11 = onep1 + zep01
319 onep12 = onep1 + zep02
320 onep13 = onep1 + zep03
321 onep14 = onep1 + zep04
322 onep15 = onep1 + zep05
323 onep16 = onep1 + zep06
324 onep17 = onep1 + zep07
325 onep18 = onep1 + zep08
326 onep19 = onep1 + zep09
327 onep136 = onep13 + sixem3
328 onep138 = onep13 + eightem3
329
330 onep2 = one + zep2
331 onep21 = onep2 + zep01
332 onep23 = onep2 + zep03
333 onep24 = onep2 + zep04
334 onep25 = five / four
335 onep26 = onep2 + zep06
336 onep27 = onep2 + zep07
337 onep28 = onep2 + zep08
338 onep251 = onep25 + em03
339 onep252 = onep25 + twoem3
340
341 onep3 = one + zep3
342 onep31 = onep3 + zep01
343 onep32 = onep3 + zep02
344 onep33 = onep3 + zep03
345 onep34 = onep3 + zep04
346 onep35 = onep3 + zep05
347 onep36 = onep3 + zep06
348 onep37 = onep3 + zep07
349 onep38 = onep3 + zep08
350 onep39 = onep3 + zep09
351 onep333 = onep33 + threeem3
352 onep3348= onep33 + fourem3 + eightem4
353 onep369 = onep36 + nineem3
354
355 onep4 = one + zep4
356 onep41 = onep4 + zep01
357 onep42 = onep4 + zep02
358 onep44 = onep4 + zep04
359 onep45 = onep4 + zep05
360 onep46 = onep4 + zep06
361 onep48 = onep4 + zep08
362 onep49 = onep4 + zep09
363 onep414 = onep41 + fourem3
364 onep484 = onep48 + fourem3
365
366 onep5 = three_half
367 onep51 = onep5 + zep01
368 onep52 = onep5 + zep02
369 onep57 = onep5 + zep07
370 onep59 = onep5 + zep09
371
372 onep6 = one + zep6
373 onep62 = onep6 + zep02
374 onep65 = onep6 + zep05
375 onep66 = onep6 + zep06
376 onep67 = onep6 + zep07
377 onep68 = onep6 + zep08
378
379 onep7 = one + zep7
380 onep72 = onep7 + zep02
381 onep73 = onep7 + zep03
382 onep75 = onep7 + zep05
383 onep78 = onep7 + zep08
384 onep732 = onep73 + twoem3
385
386 onep8 = one + zep8
387 onep82 = onep8 + zep02
388 onep84 = onep8 + zep04
389 onep88 = onep8 + zep08
390 onep833 = onep8 + zep03 + threeem3
391 onep8688 = onep8 + zep06 + eightem3 + eightem4
392
393 onep9 = one + zep9
394 onep92 = onep9 + zep02
395 onep96 = onep9 + zep06
396 onep98 = onep9 + zep08
397
398 twop03 = two + zep03
399 twop05 = two + zep05
400 twop06 = two + zep06
401 twop07 = two + zep07
402 twop08 = two + zep08
403 twop1 = two + zep1
404 twop2 = two + zep2
405 twop30 = two + zep3
406 twop4 = two + zep4
407 twop5 = five / two
408 twop60 = two + zep6
409 twop8 = two + zep8
410
411 twop11 = twop1 + zep01
412 twop12 = twop1 + zep02
413 twop13 = twop1 + zep03
414 twop16 = twop1 + zep06
415 twop19 = twop1 + zep09
416 twop23 = twop2 + zep03
417 twop24 = twop2 + zep04
418 twop25 = twop2 + zep05
419 twop33 = twop30 + zep03
420 twop38 = twop30 + zep08
421 twop45 = twop4 + zep05
422 twop47 = twop4 + zep07
423 twop52 = twop5 + zep02
424 twop83 = twop8 + zep03
425 twop84 = twop8 + zep04
426 twop85 = twop8 + zep05
427
428 twop185 = twop1 + zep08 + fiveem3
429 twop261 = twop2 + zep06 + em03
430 twop645 = twop60 + zep04 + fiveem3
431
432 threep02 = three + zep02
433 threep06 = three + zep06
434 threep14 = three + zep14
435 threep15 = three + zep15
436 threep25 = three + zep2 + zep05
437 threep27 = three + zep2 + zep07
438 threep44 = three + zep4 + zep04
439 threep46 = three + zep4 + zep06
440 threep52 = three + zep5 + zep02
441 threep72 = three + zep7 + zep02
442 threep76 = three + zep7 + zep06
443 threep8 = three + zep8
444 threep94 = three + zep9 + zep04
445 threep266 = three + zep2 + zep06 + sixem3
446 threep464 = three + zep4 + zep06 + fourem3
447 threep4641 = threep464 + em04
448
449 fourp186 = four + zep1 + zep08 + sixem3
450 fourp40 = four + zep4
451 fourp76 = four + zep7 + zep06
452
453 fivep5 = five + zep5
454 fivep56 = five + zep5 + zep06
455 fivep76 = five + zep7 + zep06
456
457 sip2832 = six + zep2 + zep08 + threeem3 + twoem4
458
459 ninep24 = nine + zep2 + zep04
460 ninep38 = nine + zep3 + zep08
461 ninep793 = nine + zep7 + zep09 + threeem3
462
463 elevenp5 = eleven + zep5
464 twelvep6 = twelve + zep6
465 fourteenp2 = fourteen + zep2
466 sixteenp5 = sixteen + zep5
467 o88p9844 = 88. + zep9 + zep08 + fourem3 + fourem4
468
469 onep3em5 = onep3 * em05
470 onep001em32 = 1001. * em35
471 ninep637em5 = 9637. * em08
472 eightp314em5 = 8314. * em08
473
474 onep1ep30 = onep1 * ep30
475 onep2ep30 = onep2 * ep30
476
477 zep999ep31 = zep999 * ep31
478
479
480
481 zep426667 = hundred28 / three100
482 zep222 = ten / fourty5
483 zep00666666667 = two / three100
484 zep01666666667 = five / three100
485 zep444 = hundred / twohundred25
486 onep222 = one + zep222
487 onep8333 = one + twenty5 / thirty
488 twop444 = two + zep444
489 twop6666666667 = eight / three
490 twop6667 = twop6666666667
491 fivep333 = sixteen / three
492
493
494 zep54 = one - zep46
495 zep1381966 = zep13 + eightem3 + em04 + nineem5 + sixem6 + six*em07
496 zep5854102 = zep5 + zep08 + fiveem3 + fourem4 + em05 + two*em07
497
498
499
500 r_igc = 8.314472
501 stefboltz = 5.6704 / ep08
502
503
504
505 sqr2 = sqrt(two)
506 sqr3 = sqrt(three)
507 sqr5 = sqrt(five)
508 sqr6 = sqrt(six)
509 sqr7 = sqrt(seven)
510 sqr3_2 = sqrt(three_half)
511 sqr4_3 = two * sqrt(third)
512 sqr16_3 = two * sqr4_3
513 pi = atan2(zero,-one)
514 invpi = one / pi
515
516
517
518 ep9 = ep09
519 zep80 = zep8
520 zep90 = zep9
521 unp20 = onep2
522 unp30 = onep3
523 unp70 = onep7
524 unp90 = onep9
525 unem3 = em03
526 unem4 = em04
527 unem5 = em05
528 unem6 = em06
529 em3 = em03
530 em4 = em04
531 em5 = em05
532 em6 = em06
533 em7 = em07
534 em8 = em08
535 em9 = em09
536 twoem2 = zep02
537 fiveem2 = zep05
538 nineem2 = zep09
539
540 RETURN
1.15.0