1 Introduction

This report gives the result of running the computer algebra independent integration problems.

The listing of the problems used by this report are

The above zip ﬁles are maintained by and were downloaded from Albert Rich Rubi web site.

The current number of problems in this test suite is [69522].

1.1 Listing of CAS systems tested

The following systems were tested at this time.

1. Mathematica 11.3 (64 bit).
2. Rubi 4.15.2 in Mathematica 11.3.
3. Maple 2018.1 (64 bit).
4. Maxima 5.41 Using Lisp ECL 16.1.2.
5. Fricas 1.3.4.
6. Sympy 1.1.1 under Python 3.6.5 using Anaconda distribution.
7. Giac/Xcas 1.4.9.

Maxima, Fricas and Giac/Xcas were called from inside SageMath version 8.3. This was done using SageMath integrate command by changing the name of the algorithm to use the diﬀerent CAS systems.

Sympy was called directly using Python.

1.2 Design of the test system

The following diagram gives a high level view of the current test build system.

1.3 Timing

The command AboluteTiming[] was used in Mathematica to obtain the elapsed time for each integrate call. In Maple, the command Usage was used as in the following example

cpu_time := Usage(assign ('result_of _int',int(expr,x)),output='realtime'

For all other CAS systems, the elapsed time to complete each integral was found by taking the diﬀerence between the time after the call has completed from the time before the call was made. This was done using Python’s time.time() call.

All elapsed times shown are in seconds. A time limit of 3 minutes was used for each integral. If the integrate command did not complete within this time limit, the integral was aborted and considered to have failed and assigned an F grade. The time used by failed integrals due to time out is not counted in the ﬁnal statistics. (May be I should change this and make it count?)

1.4 Veriﬁcation

A veriﬁcation phase was applied on the result of integration for Rubi and Mathematica. Future version of this report will implement veriﬁcation for the other CAS systems. For the integrals whose result was not run through a veriﬁcation phase, it is assumed that the antiderivative produced was correct.

Veriﬁcation phase has 3 minutes time out. An integral whose result was not veriﬁed could still be correct. Further investigation is needed on those integrals which failed veriﬁcations. Such integrals are marked in the summary table below and also in each integral separate section so they are easy to identify and locate.

1.5 Important notes about some of the results

Important note about Maxima results Since these integrals are run in a batch mode, using an automated script, and by using sagemath (SageMath uses Maxima), then any integral where Maxima needs an interactive response from the user to answer a question during evaluation of the integral in order to complete the integration, will fail and is counted as failed.

The exception raised is ValueError. Therefore Maxima result below is lower than what could result if Maxima was run directly and each question Maxima asks was answered correctly.

The percentage of such failures were not counted for each test ﬁle, but for an example, for the Timofeev test ﬁle, there were about 30 such integrals out of total 705, or about 4 percent. This pecrentage can be higher or lower depending on the speciﬁc input test ﬁle.

Such integrals can be indentiﬁed by looking at the output of the integration in each section for Maxima. If the output was an exception ValueError then this is most likely due to this reason.

Maxima integrate was run using SageMath with the following settings set by default

'besselexpand : true'
'display2d : false'
'domain : complex'
'keepfloat : true'

SageMath loading of Maxima abs_integrate was found to cause some problem. So the following code was added to disable this eﬀect.

from sage.interfaces.maxima_lib import maxima_lib
maxima_lib.set('extra_definite_integration_methods', '[]')
maxima_lib.set('extra_integration_methods', '[]')

Important note about FriCAS and Giac/XCAS results There are Few integrals which failed due to SageMath not able to translate the result back to SageMath syntax and not because these CAS system were not able to do the integrations.

These will fail With error Exception raised: NotImplementedError

The number of such cases seems to be very small. About 1 or 2 percent of all integrals.

Hopefully the next version of SageMath will have complete translation of FriCAS and XCAS syntax and I will re-run all the tests again when this happens.

Important note about ﬁnding leaf size of antiderivative For Mathematica, Rubi and Maple, the buildin system function LeafSize is used to ﬁnd the leaf size of each antiderivative.

The other CAS systems (SageMath and Sympy) do not have special buildin function for this purpose at this time. Therefore the leaf size is determined as follows.

For Fricas, Giac and Maxima (all called via sagemath) the following code is used

#see https://stackoverflow.com/questions/25202346/how-to-obtain-leaf-count-expression-size-in-sage

def tree(expr):
if expr.operator() is None:
return expr
else:
return [expr.operator()]+map(tree, expr.operands())

try:
# 1.35 is a fudge factor since this estimate of leaf count is bit lower than
#what it should be compared to Mathematica's
leafCount = round(1.35*len(flatten(tree(anti))))
except Exception as ee:
leafCount =1

For Sympy, called directly from Python, the following code is used

try:
# 1.7 is a fudge factor since it is low side from actual leaf count
leafCount = round(1.7*count_ops(anti))

except Exception as ee:
leafCount =1

When these cas systems have a buildin function to ﬁnd the leaf size of expressions, it will be used instead, and these tests run again.

The table below summarizes the grading of each CAS system.

Important note: A number of problems in this test suite have no antiderivative in closed form. This means the antiderivative of these integrals can not be expressed in terms of elementary, special functions or Hypergeometric2F1 functions. RootSum and RootOf are not allowed.

If a CAS returns the above integral unevaluated within the time limit, then the result is counted as passed and assigned an A grade.

However, if CAS times out, then it is assigned an F grade even if the integral is not integrable, as this implies CAS could not determine that the integral is not integrable in the time limit.

If a CAS returns an antiderivative to such an integral, it is assigned an A grade automatically and this special result is listed in the introduction section of each individual test report to make it easy to identify as this can be important result to investigate.

The results given in in the table below reﬂects the above.

 System solved Failed Rubi % 99.9 ( 69452 ) % 0.1 ( 70 ) Mathematica % 97.33 ( 67666 ) % 2.67 ( 1856 ) Maple % 83.8 ( 58257 ) % 16.2 ( 11265 ) Maxima % 45.54 ( 31662 ) % 54.46 ( 37860 ) Fricas % 66.27 ( 46074 ) % 33.73 ( 23448 ) Sympy % 30.29 ( 21057 ) % 69.71 ( 48465 ) Giac % 53.12 ( 36930 ) % 46.88 ( 32592 )

Table 1: Time and leaf size performance for each CAS

The table below gives additional break down of the grading of quality of the antiderivatives generated by each CAS. The grading is given using the letters A,B,C and F with A being the best quality. The grading is accomplished by comparing the antiderivative generated with the optimal antiderivatives included in the test suite. The following table describes the meaning of these grades.

 grade description A Integral was solved and antiderivative is optimal in quality and leaf size. B Integral was solved and antiderivative is optimal in quality but leaf size is larger than twice the optimal antiderivatives leaf size. C Integral was solved and antiderivative is non-optimal in quality. This can be due to one or more of the following reasons antiderivative contains a hypergeometric function and the optimal antiderivative does not. antiderivative contains a special function and the optimal antiderivative does not. antiderivative contains the imaginary unit and the optimal antiderivative does not. F Integral was not solved. Either the integral was returned unevaluated within the time limit, or it timed out, or CAS hanged or crashed or an exception was raised.

Table 2: Description of grading applied to integration result

Currently grading is implemented only for Mathematica, Rubi and Maple. For other CAS systems (Maxima, Fricas, Sympy, Giac), the grading function is not yet implemented. For these systems, a grade of A is assigned if the integrate command completes successfully and a grade of F otherwise.

Based on the above, the following table summarizes the grading for this test suite.

 System % A grade % B grade % C grade % F grade Rubi 99.76 0.08 0.06 0.1 Mathematica 75.37 8.46 15.81 2.67 Maple 52.49 22.13 9.18 16.2 Maxima 45.54 0. 0. 54.46 Fricas 66.27 0. 0. 33.73 Sympy 30.29 0. 0. 69.71 Giac 53.12 0. 0. 46.88

Table 3: Antiderivative Grade distribution for each CAS

The following is a Bar chart illustration of the data in the above table.

The Mathematica code used to generate the above ﬁgure is contained in this small notebook make_3d_barchart.nb

The ﬁgure below compares the CAS systems for each grade level.

1.7 Time and leaf size Performance

The table below summarizes the performance of each CAS system in terms of time used and leaf size of results.

 System Mean time (sec) Mean size Normalized mean Median size Normalized median Rubi 0.53 154.16 1. 107. 1. Mathematica 3.02 1871.7 6.4 104. 1. Maple 0.21 40868.6 461.8 130. 1.26 Maxima 1.81 270.68 2.53 96. 1.36 Fricas 1.02 335.95 2.32 77. 1.27 Sympy 9.82 247.65 2.65 68. 1.14 Giac 0.9 175.9 1.8 104.5 1.4

1.8 Performance per integrand type

The following are the diﬀerent integrand types the test suite contains.

1. Algebraic Binomial problems (products involving powers of binomials and monomials).
2. Algebraic Trinomial problems (products involving powers of trinomials, binomials and monomials).
3. Miscellaneous Algebraic functions.
4. Exponentials.
5. Logarithms.
6. Trigonometric.
7. Inverse Trigonometric.
8. Hyperbolic functions.
9. Inverse Hyperbolic functions.
10. Special functions.
11. Independent tests.

The following table gives percentage solved of each CAS per integrand type.

 Integrand type problems Rubi Mathematica Maple Maxima Fricas Sympy Giac Independent tests 1867 98.8 98.6 93.1 82.2 92.6 69.4 83.3 Algebraic Binomial 13952 100. 99. 83.8 43.8 71.6 58.2 64. Algebraic Trinomial 10279 100. 98.5 90.4 39.5 77.5 39.9 64.5 Algebraic Miscellaneous 1379 99.3 94.8 88. 49.6 72.3 47.2 58.8 Exponentials 948 99.8 93.8 80.8 65.8 90.3 40.8 46.1 Logarithms 2505 99.8 96.2 53.5 52.8 51.3 25.9 46. Trigonometric 22422 99.9 96.6 84.4 46.1 57.3 10. 43.8 Inverse Trigonometric 4505 100. 96.5 83.7 31.8 46.4 23.3 45.7 Hyperbolic 5043 99.9 97.3 83.4 57.2 85.4 19.5 59.9 Inverse Hyperbolic 6421 100. 96.3 80.5 40.2 61.4 24.4 39.1 Special functions 201 100. 96.5 71.6 44.3 31.3 8. 0.5

Table 4: Percentage solved per integrand type

In addition to the above table, for each type of integrand listed above, 3D chart is made which shows how each CAS performed on that speciﬁc integrand type.

These charts and the table above can be used to show where each CAS relative strength or weakness in the area of integration.

1.9 Maximum leaf size ratio for each CAS against the optimal result

The following table gives the largest ratio found in each test ﬁle, between each CAS antiderivative and the optimal antiderivative.

For each test input ﬁle, the problem with the largest ratio $$\frac{\text{CAS leaf size}}{\text{Optimal leaf size}}$$ is recorded with the corresponding problem number.

In each column in the table below, the ﬁrst number is the maximum leaf size ratio, and the number that follows inside the parentheses is the problem number in that speciﬁc ﬁle where this maximum ratio was found. This ratio is determined only when CAS solved the the problem and also when an optimal antiderivative is known.

If it happens that a CAS was not able to solve all the integrals in the input test ﬁle, or if it was not possible to obtain leaf size for the CAS result for all the problems in the ﬁle, then a zero is used for the ratio and -1 is used for the problem number.

This makes it easy to locate the problem. In the future, a direct link will be added as well.

 Table 5: Maximum leaf size ratio for each CAS against the optimal result ﬁle # Rubi Mathematica Maple Maxima FriCAS Sympy Giac 1 1. (1) 3.9 (50) 16.6 (114) 5.1 (169) 7.2 (47) 4.5 (48) 15.5 (55) 2 7.3 (21) 7.7 (14) 3.6 (17) 2.6 (4) 13.6 (24) 16.7 (5) 6.2 (2) 3 1. (1) 15.7 (6) 16.5 (6) 14.9 (7) 2.4 (5) 1.9 (5) 2.6 (5) 4 6.8 (5) 149.4 (44) 40.8 (46) 5.1 (40) 11.6 (8) 4.4 (40) 7.2 (1) 5 1. (38) 54.7 (278) 12737.8 (278) 11. (280) 10.4 (280) 73.4 (175) 22.9 (161) 6 1. (1) 1. (2) 2.2 (4) 2.6 (1) 1.9 (7) 0.8 (6) 3.2 (5) 7 2.2 (3) 5.6 (7) 1.8 (3) 3.8 (3) 3.4 (2) 3.6 (8) 2.6 (3) 8 1.6 (50) 5.3 (31) 7.9 (70) 9. (11) 7.7 (42) 26.4 (71) 7. (70) 9 1. (1) 7.2 (80) 4.3 (341) 16.3 (328) 7.2 (142) 191.1 (251) 48.1 (368) 10 3.2 (335) 242.6 (327) 3343.5 (327) 125.1 (408) 106.2 (593) 76.7 (217) 20. (494) 11 579. (80) 100. (80) 317. (80) 3.7 (2) 95. (80) 41.3 (17) 3.7 (2) 12 1.8 (6) 2.3 (4) 1.2 (8) 2. (2) 2.4 (2) 3.4 (3) 2.1 (2) 13 7.1 (369) 23.8 (1323) 30.9 (1323) 44.4 (1323) 44.4 (1323) 166. (484) 29.8 (1488) 14 2. (858) 198.6 (3104) 22.6 (1084) 30. (1699) 29.4 (1084) 84.5 (366) 13.7 (141) 15 3.3 (97) 56. (147) 28.5 (100) 2.7 (154) 9.2 (131) 49.2 (119) 2.6 (27) 16 2.6 (69) 28.1 (136) 68. (142) 2.7 (69) 11.5 (25) 60.1 (27) 8.7 (28) 17 8.2 (664) 8.7 (663) 7.9 (196) 13.5 (196) 13.5 (196) 55.3 (528) 14.8 (416) 18 1.6 (251) 11.2 (346) 147.4 (69) 5.7 (73) 11.9 (221) 18.8 (32) 6.4 (73) 19 1. (656) 3.3 (972) 46.7 (754) 4.2 (313) 12. (491) 28.8 (325) 11.6 (553) 20 1.3 (64) 3.4 (88) 15.2 (57) 1.5 (18) 0. (-1) 3. (21) 3.2 (63) 21 1. (1) 2.8 (47) 10.4 (15) 0. (-1) 9.5 (15) 43. (16) 8.4 (18) 22 1.2 (173) 1.9 (45) 2. (162) 3.4 (163) 4.6 (43) 18.3 (93) 5.3 (157) 23 8.4 (2679) 13.4 (2906) 141.8 (2906) 17.8 (2278) 15.2 (2584) 84.8 (1709) 17. (2807) 24 1.4 (243) 42.1 (30) 17.9 (166) 2.6 (253) 13.8 (106) 28.1 (139) 8.1 (185) 25 1.4 (629) 12.3 (882) 77.4 (546) 39.3 (894) 25.4 (317) 48.4 (810) 14.9 (882) 26 1.2 (46) 6. (14) 51.1 (15) 0. (-1) 38.3 (15) 0. (-1) 7.8 (4) 27 1.2 (316) 3.8 (45) 10.4 (43) 4.3 (161) 31.8 (156) 43.6 (2) 12.6 (579) 28 1.3 (278) 10. (328) 51.5 (297) 11.3 (328) 13.8 (348) 10. (328) 13.8 (348) 29 1. (1) 4.5 (277) 4.9 (269) 3.1 (279) 5.4 (269) 21.6 (269) 9.6 (269) 30 2.8 (83) 4.8 (25) 5.8 (74) 3.5 (98) 4.7 (103) 16.4 (63) 4.5 (74) 31 2. (2408) 14.9 (2447) 70.8 (2340) 19.2 (1486) 48.8 (2282) 115.7 (1974) 49.7 (971) 32 1.3 (1471) 16.3 (1088) 91. (1626) 29. (2015) 61.5 (1452) 92.2 (2054) 17.7 (2163) 33 2.1 (833) 46.1 (890) 116.1 (801) 8.1 (579) 45.9 (616) 70.9 (361) 9. (357) 34 1. (1) 11.1 (14) 424.6 (78) 3.7 (95) 40.7 (112) 1.2 (19) 5.4 (121) 35 1. (129) 62.2 (37) 14197.2 (12) 8.9 (27) 40.7 (117) 6.9 (13) 7.8 (37) 36 1.8 (76) 18.1 (265) 421. (278) 120.1 (278) 11. (367) 114.1 (278) 17.2 (41) 37 1.7 (636) 11.8 (1093) 13.9 (882) 7.3 (515) 33.2 (1077) 28.5 (1105) 8.5 (785) 38 1.7 (216) 18.1 (4) 50.7 (224) 9. (96) 72.9 (117) 19.6 (222) 10. (96) 39 1.9 (319) 15.1 (208) 25.9 (136) 7.5 (70) 57.4 (388) 47.5 (220) 7.5 (70) 40 1. (59) 3.6 (103) 31.6 (57) 1.8 (111) 3.5 (46) 58.1 (87) 5. (108) 41 1.6 (135) 1.8 (51) 18.6 (128) 2.1 (131) 65. (60) 27.3 (39) 4.3 (131) 42 1. (1) 3. (5) 2.4 (6) 0. (-1) 5.6 (5) 0. (-1) 1.6 (13) 43 2.1 (154) 32.7 (602) 54.7 (609) 8.5 (609) 28.6 (632) 26.3 (438) 13.2 (598) 44 1. (1) 36. (86) 2.7 (37) 1.8 (26) 41.1 (41) 42.2 (68) 9.1 (67) 45 1. (67) 25.1 (143) 2909.3 (93) 93. (94) 76.1 (96) 82.9 (93) 13.9 (101) 46 1. (1) 11. (17) 1.7 (11) 2.8 (16) 3. (16) 3.2 (11) 5.7 (11) 47 1. (1) 2.7 (104) 8. (92) 1.2 (71) 12.8 (102) 20.7 (24) 8.2 (72) 48 6.2 (424) 11.5 (160) 1223.1 (192) 57.1 (63) 22.2 (215) 84.3 (192) 41. (195) 49 4.1 (859) 172.1 (872) 3059.3 (872) 6.8 (451) 54.5 (715) 57.7 (157) 14.1 (754) 50 1. (1) 1.2 (41) 9.5 (87) 3. (2) 2.7 (81) 2.6 (2) 4.5 (40) 51 1. (1) 2.5 (49) 13. (46) 2.2 (49) 6.1 (58) 2.4 (25) 13.6 (35) 52 1.2 (638) 5.8 (39) 34.4 (267) 283.8 (255) 38.5 (292) 42.2 (352) 23.4 (565) 53 1.5 (769) 12.4 (657) 504.9 (785) 14. (654) 28.5 (388) 21.2 (329) 13.3 (550) 54 1.3 (368) 12.6 (472) 161.9 (62) 19.4 (335) 11.9 (398) 14.9 (398) 9.2 (399) 55 1.5 (390) 46.7 (131) 54.3 (175) 18.6 (380) 6.3 (204) 13.6 (321) 8.2 (329) 56 1. (74) 49.4 (142) 35. (31) 17.2 (118) 7.6 (83) 5.2 (36) 6. (83) 57 1. (1) 14.3 (5) 2190.9 (52) 19.2 (75) 9. (169) 15.9 (158) 7.6 (128) 58 1. (1) 114.1 (497) 33.3 (493) 121.9 (284) 9.5 (301) 33.1 (355) 7.6 (105) 59 1. (1) 13.4 (250) 7.6 (83) 26.2 (132) 18.6 (209) 10.7 (114) 15.8 (5) 60 1. (1) 9.2 (12) 4.3 (51) 42.1 (5) 4.3 (6) 17.1 (49) 5.8 (7) 61 1. (1) 8.6 (38) 7.7 (65) 28.8 (45) 2.1 (103) 2.2 (12) 16. (6) 62 1. (1) 4.2 (36) 7.8 (201) 227.2 (37) 10.8 (36) 8.4 (115) 7.2 (155) 63 2. (615) 577.9 (348) 447.6 (605) 33.6 (209) 15.8 (415) 10.1 (415) 11.5 (186) 64 1. (1) 1.1 (10) 1.4 (29) 20.7 (33) 1.5 (10) 0. (-1) 3.2 (30) 65 1.6 (103) 536.8 (138) 3.6 (200) 8.5 (92) 5.5 (105) 2.6 (40) 8.1 (105) 66 1.9 (621) 363.7 (767) 8120. (795) 47.1 (596) 8.8 (256) 15.9 (4) 24.3 (336) 67 1.6 (1108) 14860. (1237) 101. (174) 42. (41) 7. (913) 31.3 (921) 18.1 (174) 68 1.3 (12) 3375. (37) 687.4 (48) 9.7 (16) 4. (13) 3.4 (1) 8.3 (13) 69 1.2 (205) 234.6 (201) 8076. (351) 47.4 (47) 5.8 (47) 8.7 (39) 12. (118) 70 1. (1) 6.7 (10) 3.9 (2) 16.8 (1) 3.2 (2) 13.1 (3) 7.3 (12) 71 1.4 (32) 177.5 (23) 4.4 (33) 4.4 (20) 2.9 (18) 2.3 (32) 1.3 (32) 72 1.8 (236) 1209.5 (555) 24. (357) 23. (475) 55.8 (257) 108.5 (89) 20.7 (29) 73 1. (1) 2.2 (2) 2.1 (4) 1.7 (2) 2.7 (2) 6.7 (2) 3.1 (2) 74 1. (1) 2.3 (14) 6. (13) 1.4 (19) 34.1 (1) 4.4 (19) 2.5 (14) 75 1. (1) 38. (284) 8.3 (12) 39. (30) 4.2 (272) 2.7 (64) 16.9 (64) 76 1. (1) 4.5 (129) 8.3 (76) 75. (91) 13.2 (33) 4.1 (9) 15. (5) 77 1. (1) 2.4 (61) 3.4 (50) 30.7 (12) 2.8 (6) 6.2 (32) 1.4 (19) 78 1. (1) 2.1 (3) 4.2 (26) 6.8 (40) 2.2 (86) 6. (61) 4.2 (92) 79 4.3 (11) 4.6 (60) 13.2 (78) 4.4 (3) 5. (11) 9.5 (1) 5. (11) 80 1. (1) 1. (10) 1.4 (29) 20.7 (32) 1.5 (10) 0. (-1) 3.1 (30) 81 1. (1) 3.2 (1) 3.4 (3) 5.5 (3) 3. (3) 0. (-1) 4.1 (3) 82 1.4 (370) 35.3 (907) 9.3 (642) 50.1 (213) 3.7 (781) 12.8 (781) 10.4 (782) 83 1. (1) 21.5 (4) 2.9 (2) 0. (-1) 0. (-1) 0. (-1) 0. (-1) 84 1. (1) 3. (1) 1.8 (1) 4.7 (1) 2.3 (1) 0. (-1) 2.2 (1) 85 1.1 (40) 77.4 (639) 14.3 (436) 61.1 (445) 3.5 (285) 7.4 (40) 5.4 (80) 86 1. (1) 159.1 (370) 8. (29) 40. (29) 2.8 (12) 3. (9) 7.1 (35) 87 1.4 (940) 60.6 (1445) 18. (1154) 52.8 (190) 3.6 (296) 6.8 (41) 5.6 (81) 88 1.2 (81) 1735.1 (96) 6.9 (70) 13.2 (64) 94.2 (82) 3.6 (44) 5. (91) 89 1. (1) 2.1 (9) 7.6 (21) 1.4 (2) 2. (2) 1.4 (3) 5.7 (13) 90 1. (20) 1.9 (5) 10.5 (13) 1.1 (12) 33.8 (13) 2.8 (3) 2.5 (5) 91 1. (1) 341.4 (356) 326.6 (179) 16.3 (124) 10.7 (251) 3.1 (376) 28.2 (4) 92 1. (1) 4.5 (44) 6. (54) 15. (49) 4.3 (54) 1.4 (8) 7.3 (8) 93 1. (1) 3.8 (44) 1.5 (21) 10.6 (52) 3.9 (44) 8.3 (15) 2.6 (21) 94 1.5 (562) 101.2 (604) 173.9 (617) 110.2 (143) 8.7 (569) 1.5 (512) 14.1 (32) 95 1. (1) 15.2 (67) 4.1 (61) 3.9 (67) 10.2 (75) 5.3 (2) 10.3 (14) 96 1.4 (891) 2373.9 (1245) 9943.8 (636) 48.9 (1043) 122.7 (551) 14. (1209) 13.6 (426) 97 1. (1) 1046.6 (336) 17625.3 (454) 53.1 (689) 101.6 (340) 12. (271) 14.5 (232) 98 1. (86) 6480.5 (100) 172.7 (79) 4. (39) 15.4 (39) 20.1 (27) 15.2 (9) 99 1. (1) 1638.1 (398) 2905.5 (351) 7.6 (277) 27. (288) 23.8 (232) 24.5 (186) 100 1. (1) 1561.2 (27) 31765.8 (14) 1.1 (32) 0. (-1) 0. (-1) 0. (-1) 101 1. (1) 166.2 (47) 288.2 (43) 1.7 (4) 33.6 (30) 2.6 (1) 5.7 (3) 102 1. (1) 4.5 (42) 10.1 (27) 25.1 (47) 7.5 (37) 2.1 (8) 4.4 (25) 103 1. (1) 2.5 (11) 3.4 (16) 4.4 (11) 5.1 (14) 4.2 (2) 3.5 (7) 104 1. (1) 2.4 (5) 4.3 (9) 5.8 (7) 4.4 (15) 5.3 (2) 3.7 (6) 105 1. (1) 4. (94) 69.3 (103) 2.6 (94) 69.7 (104) 9. (92) 3.2 (94) 106 1. (1) 1775.8 (62) 35. (29) 18.4 (8) 14.4 (9) 21.5 (6) 6.4 (38) 107 1. (1) 1997.4 (22) 31765.8 (3) 0. (-1) 0. (-1) 0. (-1) 0. (-1) 108 1. (1) 35.1 (280) 9.6 (259) 35. (47) 7.7 (42) 3.3 (1) 15.7 (42) 109 1. (1) 5.9 (39) 4.1 (29) 20. (16) 9.8 (36) 0. (-1) 8.1 (18) 110 1. (1) 3.5 (5) 5.9 (73) 162.5 (20) 10.3 (18) 2. (53) 2.7 (12) 111 1.4 (423) 116.7 (711) 14.7 (578) 70.7 (255) 5.7 (267) 2.6 (5) 6.2 (95) 112 1. (1) 74.2 (140) 12.6 (284) 3.9 (65) 7.8 (227) 0. (-1) 7.7 (197) 113 1.7 (340) 172.5 (353) 53.6 (339) 5. (67) 6.7 (306) 17.1 (90) 6.3 (301) 114 1.3 (115) 33.2 (241) 1120.6 (153) 50.6 (109) 2.9 (22) 0. (-1) 5.6 (63) 115 2.2 (197) 105.6 (175) 7.1 (238) 58.3 (130) 4.1 (130) 3.1 (170) 5.9 (256) 116 1.3 (265) 88.8 (383) 16.3 (633) 71. (259) 3.3 (47) 2.2 (47) 7.3 (124) 117 1. (1) 4.6 (41) 19.4 (25) 18.2 (25) 3.6 (58) 1.3 (32) 3.6 (41) 118 1.3 (775) 145. (857) 19.8 (970) 62. (1289) 3.3 (309) 3.1 (930) 7.6 (488) 119 1. (1) 66.8 (138) 544.5 (434) 37.4 (297) 9.2 (256) 9.8 (460) 15.3 (366) 120 1. (1) 5.6 (42) 12.4 (21) 45.1 (39) 5.2 (42) 3.1 (1) 4.2 (41) 121 1. (1) 4.7 (10) 5. (74) 53.1 (15) 9.7 (18) 2. (53) 4. (61) 122 1. (1) 5.3 (36) 19.6 (18) 8.7 (13) 4.6 (3) 0. (-1) 18.4 (15) 123 1. (1) 2.5 (8) 4. (9) 6.6 (8) 2.1 (7) 0. (-1) 3. (8) 124 1.3 (20) 3.3 (10) 2.3 (22) 4.8 (1) 3.6 (18) 0. (-1) 3. (10) 125 1. (1) 2.7 (3) 2.2 (8) 3.4 (8) 3.1 (9) 4.9 (18) 4.5 (12) 126 1. (1) 1.2 (1) 1.8 (1) 0. (-1) 0. (-1) 0. (-1) 0. (-1) 127 1. (12) 3.6 (2) 26.8 (15) 32.8 (18) 11.1 (25) 0. (-1) 4.2 (1) 128 1. (1) 123.5 (187) 2946325.1 (77) 114.8 (57) 13.7 (186) 14.9 (225) 33.2 (29) 129 3.3 (23) 45. (29) 5.9 (146) 12.8 (209) 11.5 (143) 13.2 (275) 6.7 (220) 130 1.1 (281) 16.6 (196) 14.6 (80) 78.9 (391) 24.5 (273) 9.2 (80) 31.4 (42) 131 1. (1) 2.7 (1) 3.4 (5) 0.6 (5) 26.4 (4) 1.1 (5) 0.9 (5) 132 4.3 (179) 14. (130) 13.1 (179) 122.7 (150) 17.4 (151) 6.9 (12) 23.2 (88) 133 19.2 (34) 9.1 (133) 43.8 (34) 109.8 (34) 5.6 (63) 9.8 (39) 12.9 (140) 134 10.8 (759) 718.9 (434) 708.4 (860) 241.3 (64) 37.4 (503) 223.4 (478) 82.5 (796) 135 1.4 (107) 2.5 (209) 4.8 (156) 2.3 (155) 2.4 (7) 2.3 (11) 13.3 (145) 136 1.7 (100) 4. (196) 19.9 (90) 5.6 (195) 1.9 (7) 2.5 (100) 16.8 (7) 137 1. (1) 2.3 (56) 7.7 (48) 2.1 (67) 2.1 (7) 2. (9) 6.8 (9) 138 1.9 (147) 15.5 (85) 13.9 (55) 3.8 (186) 4.3 (206) 8.1 (206) 12.3 (346) 139 1. (1) 4.9 (41) 2.8 (156) 2.4 (155) 3.6 (7) 2.3 (11) 17. (146) 140 1. (1) 4.5 (28) 5.6 (4) 4. (11) 3.9 (147) 2. (46) 17.7 (44) 141 1.3 (231) 6.4 (508) 85.8 (225) 7. (297) 4.2 (67) 10.8 (1417) 5.1 (67) 142 1. (1) 6.5 (55) 48.6 (37) 6.9 (1) 4.1 (30) 8.7 (51) 5. (7) 143 1.1 (117) 3.3 (75) 25.3 (49) 7.3 (67) 8.2 (83) 4.1 (10) 3.5 (10) 144 2. (101) 4.5 (108) 13.6 (86) 10.4 (132) 13.5 (132) 5.3 (133) 11.3 (132) 145 1. (1) 3.4 (203) 13. (215) 5.7 (207) 9.1 (188) 64.6 (177) 6.1 (211) 146 1.3 (109) 22.3 (110) 72.1 (110) 7.2 (164) 8.2 (190) 7.1 (105) 4.6 (129) 147 1. (1) 1.2 (7) 1. (2) 1.3 (2) 1.5 (5) 1.9 (3) 1.6 (2) 148 1.2 (68) 3.3 (25) 11.9 (105) 4.6 (31) 9. (56) 1.2 (9) 7.2 (42) 149 1. (1) 3.3 (42) 4.7 (26) 2.3 (14) 6.7 (18) 0.9 (10) 3. (25) 150 1.4 (51) 3. (114) 11.9 (112) 2.6 (22) 9. (44) 0. (-1) 3. (29) 151 1. (1) 3.8 (31) 5.7 (26) 2.2 (13) 6.7 (17) 0. (-1) 3. (24) 152 1. (1) 14.8 (45) 7.5 (379) 4.8 (1) 25.1 (329) 17.6 (262) 7.9 (260) 153 1. (1) 5.9 (53) 3.4 (98) 18.8 (90) 9.3 (20) 1.9 (10) 5.2 (95) 154 1. (1) 1.5 (24) 1.9 (28) 8.1 (7) 7.2 (21) 0. (-1) 2.2 (29) 155 1. (1) 8.7 (358) 8.4 (198) 18. (208) 35.4 (157) 43.5 (162) 11. (324) 156 1.3 (16) 6.8 (26) 33.4 (55) 29.6 (315) 85.6 (214) 28.1 (183) 11. (392) 157 1. (1) 14.8 (48) 7.4 (28) 4.8 (1) 21.6 (36) 4.1 (8) 4.8 (1) 158 1. (1) 7.9 (106) 3.6 (79) 4.7 (5) 6.6 (107) 2.3 (12) 4.6 (38) 159 1. (1) 2.1 (3) 3.4 (64) 18.9 (56) 9.3 (20) 1.9 (10) 5.2 (61) 160 1. (1) 1.5 (12) 1.9 (28) 8.1 (7) 7.2 (21) 0. (-1) 2.2 (29) 161 1. (1) 8.7 (321) 7.4 (165) 15.4 (196) 35.4 (139) 22. (152) 11. (287) 162 1.3 (60) 2.5 (11) 27.9 (34) 7.9 (40) 49. (72) 7.7 (1) 7. (38) 163 1. (1) 3.8 (61) 12.5 (46) 4.3 (8) 31.7 (11) 1.7 (8) 4.5 (8) 164 1.2 (109) 1353.9 (162) 6.6 (102) 17. (150) 50.9 (102) 8.4 (49) 5.7 (102) 165 1.3 (257) 500.2 (259) 34.4 (127) 24.4 (104) 94.7 (157) 25.2 (183) 7.4 (159) 166 1. (1) 4.7 (48) 12.7 (27) 5. (8) 27. (47) 11.7 (27) 4.3 (8) 167 1. (1) 615.2 (167) 6.2 (35) 17.5 (155) 44.1 (161) 11.6 (148) 7.9 (113) 168 1. (1) 426.7 (53) 9.2 (24) 8.7 (10) 31.9 (48) 6.4 (5) 12. (10) 169 1. (1) 3.8 (6) 7.4 (3) 3.3 (7) 21.6 (9) 0. (-1) 3.6 (7) 170 1. (1) 4.3 (18) 5.7 (77) 2.7 (15) 28.8 (15) 0. (-1) 2.7 (31) 171 3.5 (179) 10.9 (178) 12.7 (179) 11.6 (59) 45.5 (113) 2.2 (119) 4.4 (86) 172 1.4 (54) 18.2 (43) 28.3 (169) 17.8 (124) 91.1 (167) 5.4 (142) 6.7 (167) 173 1. (1) 6.9 (26) 5. (29) 4.2 (7) 21.1 (9) 0. (-1) 3.7 (7) 174 1. (1) 5.9 (18) 3.9 (78) 3.1 (15) 38.9 (15) 0. (-1) 2.5 (5) 175 3.3 (153) 10. (24) 22.3 (24) 12. (91) 45.7 (124) 0. (-1) 12. (24) 176 1.1 (12) 3.6 (2) 42.5 (8) 8.4 (1) 107.3 (19) 0. (-1) 4. (1) 177 1.9 (192) 682.6 (777) 141.6 (860) 35.1 (100) 232.6 (1047) 187.7 (808) 13.7 (1050) 178 1. (1) 1.8 (30) 2.7 (38) 2.3 (15) 7.6 (11) 1. (22) 3.5 (19) 179 1. (1) 4.3 (33) 10. (242) 5.4 (154) 8.8 (25) 2.2 (154) 4.2 (118) 180 1. (1) 2.4 (44) 2.5 (7) 2.1 (11) 2.8 (11) 1.9 (1) 3.3 (11) 181 1. (1) 7.2 (46) 6.2 (151) 15.6 (276) 29. (125) 8.2 (147) 11. (115) 182 1. (1) 3.6 (118) 2.3 (18) 1.8 (135) 5.3 (140) 1.1 (135) 3.1 (19) 183 1.5 (138) 11.9 (170) 16.1 (93) 3.7 (22) 7. (25) 1.6 (22) 3.5 (22) 184 1.2 (16) 6.6 (54) 13.8 (48) 1.9 (4) 8. (27) 1.7 (63) 2.9 (63) 185 1.3 (73) 5.7 (49) 26. (288) 3.4 (99) 20.5 (100) 7.6 (119) 8.2 (90) 186 1.6 (70) 3.6 (70) 885.7 (659) 21.7 (315) 6.7 (329) 15.2 (622) 6.1 (56) 187 1. (43) 27. (61) 63.3 (46) 7. (15) 6.7 (37) 5.2 (13) 6. (31) 188 1.7 (81) 24. (319) 24.6 (312) 5.8 (72) 12.1 (332) 2.6 (270) 7.3 (133) 189 1.5 (439) 8.9 (108) 13.7 (486) 37.8 (270) 16.5 (486) 30.7 (115) 13. (486) 190 2. (172) 2.9 (440) 16.4 (866) 13.5 (870) 12. (368) 16. (518) 8.2 (235) 191 1.2 (78) 27.2 (82) 9792.1 (185) 5.2 (95) 12.1 (251) 5.3 (92) 0.9 (198) 192 1.9 (172) 3.3 (736) 11.6 (22) 4.8 (37) 5.3 (130) 12.1 (767) 8.2 (235) 193 1. (1) 6.6 (161) 19.3 (124) 1.9 (47) 4.3 (35) 1.5 (35) 0. (-1) 194 2.8 (38) 5.5 (18) 40.2 (80) 1.3 (34) 7.6 (100) 0. (-1) 4.1 (47) 195 1.2 (75) 3.2 (114) 11.5 (112) 2.8 (10) 6.9 (50) 1.2 (9) 0. (-1) 196 1.6 (55) 8.2 (13) 5.2 (65) 3.6 (31) 9.8 (71) 2.3 (65) 3.4 (31) 197 1. (1) 9.8 (22) 2.4 (144) 5.5 (155) 2.3 (155) 0. (-1) 0. (-1)