89:
32:
14884:
874:
14594:
14215:
14585:
13978:
14879:{\displaystyle {\mathcal {J}}={\begin{bmatrix}a_{0}&{\sqrt {b_{1}}}&0&\cdots &0\\{\sqrt {b_{1}}}&a_{1}&{\sqrt {b_{2}}}&\ddots &\vdots \\0&{\sqrt {b_{2}}}&\ddots &\ddots &0\\\vdots &\ddots &\ddots &a_{n-2}&{\sqrt {b_{n-1}}}\\0&\cdots &0&{\sqrt {b_{n-1}}}&a_{n-1}\end{bmatrix}}.}
7494:
5409:
2715:
14288:
6995:
746:
3247:
13423:
15691:
derivative, and furthermore the actual error may be much less than a bound established by the derivative. Another approach is to use two
Gaussian quadrature rules of different orders, and to estimate the error as the difference between the two results. For this purpose, Gauss–Kronrod quadrature rules
7272:
5095:
2486:
10300:
4843:
8092:
9661:
8664:
5824:
7157:
7758:
9202:
6729:
14210:{\displaystyle \mathbf {J} ={\begin{bmatrix}a_{0}&1&0&\cdots &0\\b_{1}&a_{1}&1&\ddots &\vdots \\0&b_{2}&\ddots &\ddots &0\\\vdots &\ddots &\ddots &a_{n-2}&1\\0&\cdots &0&b_{n-1}&a_{n-1}\end{bmatrix}}.}
2219:
17762:. Applied Mathematics Series. Vol. 55 (Ninth reprint with additional corrections of tenth original printing with corrections (December 1972); first ed.). Washington D.C.; New York: United States Department of Commerce, National Bureau of Standards; Dover Publications.
8291:
6772:
6410:
1937:
16483:
15353:
11959:
15680:
13841:
12299:
10720:
3597:
2866:
14580:{\displaystyle {\begin{aligned}{\mathcal {J}}_{k,i}=J_{k,i}&=a_{k-1}&k&=1,2,\ldots ,n\\{\mathcal {J}}_{k-1,i}={\mathcal {J}}_{k,k-1}={\sqrt {J_{k,k-1}J_{k-1,k}}}&={\sqrt {b_{k-1}}}&k&={\hphantom {1,\,}}2,\ldots ,n.\end{aligned}}}
13144:
3478:
3363:
7909:
541:
9366:
501:
13965:
4348:
15977:
15730:. This allows for computing higher-order estimates while re-using the function values of a lower-order estimate. The difference between a Gauss quadrature rule and its Kronrod extension is often used as an estimate of the approximation error.
12537:
18512:
Riener, Cordian; Schweighofer, Markus (2018). "Optimization approaches to quadrature: New characterizations of
Gaussian quadrature on the line and quadrature with few nodes on plane algebraic curves, on the plane and in higher dimensions".
10153:
6098:
4631:
3871:
7930:
11289:
9441:
8457:
5604:
7009:
98:
approximates the function with a linear function that coincides with the integrand at the endpoints of the interval and is represented by an orange dashed line. The approximation is apparently not good, so the error is large (the
8989:
6572:
17478:
17373:
15706:
If the interval is subdivided, the Gauss evaluation points of the new subintervals never coincide with the previous evaluation points (except at zero for odd numbers), and thus the integrand must be evaluated at every point.
11566:
7622:
1474:
1366:
2056:
8099:
334:
16237:
8398:
4991:
1788:
17145:
17046:
16244:
15156:
1724:
1617:
2396:
1108:
13712:
7489:{\displaystyle \prod _{\begin{smallmatrix}1\leq j\leq n\\j\neq i\end{smallmatrix}}\left(x-x_{j}\right)={\frac {\prod _{1\leq j\leq n}\left(x-x_{j}\right)}{x-x_{i}}}={\frac {p_{n}(x)}{a_{n}\left(x-x_{i}\right)}}}
5404:{\displaystyle \int _{a}^{b}\omega (x)\,r(x)\,dx=\int _{a}^{b}\omega (x)\,\sum _{i=1}^{n}l_{i}(x)\,r(x_{i})\,dx=\sum _{i=1}^{n}\,r(x_{i})\,\int _{a}^{b}\omega (x)\,l_{i}(x)\,dx=\sum _{i=1}^{n}\,r(x_{i})\,w_{i},}
2710:{\displaystyle {\begin{aligned}\int _{8}^{30}{f(t)dt}&={\frac {30-8}{2}}\int _{-1}^{1}{f\left({\frac {30-8}{2}}x+{\frac {30+8}{2}}\right){dx}}\\&=11\int _{-1}^{1}{f\left(11x+19\right){dx}}\end{aligned}}}
18121:
12966:
11699:
5595:
15498:
9959:
13553:
12058:
15035:
10544:
14293:
11401:
8894:
13092:
3520:
15480:
17972:
5088:
18287:
Kabir, Hossein; Matikolaei, Sayed Amir
Hossein Hassanpour (2017). "Implementing an Accurate Generalized Gaussian Quadrature Solution to Find the Elastic Field in a Homogeneous Anisotropic Media".
4127:
10112:
1997:
15130:
6221:
3372:
3257:
2871:
2491:
132:
The
Gaussian quadrature chooses more suitable points instead, so even a linear function approximates the function better (the black dashed line). As the integrand is the polynomial of degree 3 (
18041:
3678:
3367:
3252:
1770:
1663:
17521:
17416:
13717:
6990:{\displaystyle \int _{a}^{b}\omega (x)r(x)dx=\sum _{i=1}^{n}r(x_{i})\int _{a}^{b}\omega (x)\prod _{\begin{smallmatrix}1\leq j\leq n\\j\neq i\end{smallmatrix}}{\frac {x-x_{j}}{x_{i}-x_{j}}}dx}
12774:
9209:
385:
5965:
5861:. Any complex conjugate roots will yield a quadratic factor that is either strictly positive or strictly negative over the entire real line. Any factors for roots outside the interval from
4224:
3515:
1517:
1409:
17185:
17086:
3941:
15783:
12331:
8783:
17869:
16884:
16695:
1267:
1187:
12020:
4561:
790:
10395:
5980:
6176:
17311:
9715:
2790:
1556:
13139:
11453:
4000:
2858:
2824:
2753:
2291:
16912:
16849:
16063:
13870:
11111:
17575:
17239:
16966:
12643:
11063:
11005:
8723:
2257:
16777:
16723:
16639:
16585:
14915:
14283:
10051:
1298:
1229:
831:
13895:
12595:
18767:
10002:
14952:
7770:
11608:
10888:
10813:
10751:
10531:
9828:
9436:
8984:
7264:
4185:
3242:{\displaystyle {\begin{aligned}11\int _{-1}^{1}{f\left(11x+19\right){dx}}&\approx 11\left\\&=11\left\\&=11\left\\&=11\left\\&=11058.44\end{aligned}}}
12813:
11102:
16099:
15064:
13589:
10489:
10343:
9786:
8441:
7586:
7222:
6567:
6516:
4050:
2051:
12053:
11031:
8319:
2400:
Change the limits so that one can use the weights and abscissae given in Table 1. Also, find the absolute relative true error. The true value is given as 11061.34 m.
13418:{\displaystyle a_{r}:={\frac {(xp_{r},p_{r})}{(p_{r},p_{r})}},\qquad b_{r}:={\frac {(xp_{r},p_{r-1})}{(p_{r-1},p_{r-1})}}={\frac {(p_{r},p_{r})}{(p_{r-1},p_{r-1})}}}
3717:, the problem is the same as that considered above. Other choices lead to other integration rules. Some of these are tabulated below. Equation numbers are given for
17547:
17211:
16938:
16749:
16611:
16018:
14244:
13616:
12677:
12326:
11694:
11662:
11635:
10953:
9750:
7617:
7550:
7523:
7186:
6480:
6453:
3785:
2479:
2440:
10780:
10453:
10424:
10149:
9394:
8931:
741:{\displaystyle \int _{-1}^{1}f(x)\,dx=\int _{-1}^{1}\left(1-x\right)^{\alpha }\left(1+x\right)^{\beta }g(x)\,dx\approx \sum _{i=1}^{n}w_{i}'g\left(x_{i}'\right).}
18480:
1785:
An integral over must be changed into an integral over before applying the
Gaussian quadrature rule. This change of interval can be done in the following way:
229:
17283:
17262:
16989:
16821:
16800:
16662:
16557:
16536:
16111:
13890:
2019:
931:
8324:
4891:
2296:
994:
17424:
17319:
12829:
10295:{\displaystyle f(x)=\prod _{\begin{smallmatrix}1\leq j\leq n\\j\neq i\end{smallmatrix}}{\frac {\left(x-x_{j}\right)^{2}}{\left(x_{i}-x_{j}\right)^{2}}}}
5468:
11462:
4838:{\displaystyle \int _{a}^{b}\omega (x)\,h(x)\,dx=\int _{a}^{b}\omega (x)\,{\big (}\,p_{n}(x)q(x)+r(x)\,{\big )}\,dx=\int _{a}^{b}\omega (x)\,r(x)\,dx.}
1420:
1312:
18426:
Piessens, R. (1971). "Gaussian quadrature formulas for the numerical integration of
Bromwich's integral and the inversion of the laplace transform".
8087:{\displaystyle {\frac {1}{x-x_{i}}}={\frac {1-\left({\frac {x}{x_{i}}}\right)^{k}}{x-x_{i}}}+\left({\frac {x}{x_{i}}}\right)^{k}{\frac {1}{x-x_{i}}}}
14964:
11296:
8795:
12973:
9656:{\displaystyle \int _{a}^{b}\omega (x){\frac {p_{n}(x)}{x-x_{i}}}dx={\frac {a_{n}}{a_{n-1}p_{n-1}(x_{i})}}\int _{a}^{b}\omega (x)p_{n-1}(x)^{2}dx}
8659:{\displaystyle \int _{a}^{b}\omega (x){\frac {p_{n}(x)}{x-x_{i}}}dx={\frac {1}{q(x_{i})}}\int _{a}^{b}\omega (x){\frac {q(x)p_{n}(x)}{x-x_{i}}}dx}
5819:{\displaystyle \int _{a}^{b}\omega (x)\,h(x)\,dx=\int _{a}^{b}\omega (x)\,r(x)\,dx=\sum _{i=1}^{n}w_{i}\,r(x_{i})=\sum _{i=1}^{n}w_{i}\,h(x_{i}).}
18817:
15389:
7152:{\displaystyle w_{i}=\int _{a}^{b}\omega (x)\prod _{\begin{smallmatrix}1\leq j\leq n\\j\neq i\end{smallmatrix}}{\frac {x-x_{j}}{x_{i}-x_{j}}}dx}
9197:{\displaystyle \int _{a}^{b}\omega (x){\frac {p_{n}(x)}{x-x_{i}}}dx={\frac {a_{n}}{p_{n-1}(x_{i})}}\int _{a}^{b}\omega (x)p_{n-1}(x)x^{n-1}dx}
17094:
16995:
6724:{\displaystyle r(x)=\sum _{i=1}^{n}r(x_{i})\prod _{\begin{smallmatrix}1\leq j\leq n\\j\neq i\end{smallmatrix}}{\frac {x-x_{j}}{x_{i}-x_{j}}}}
4998:
859:(the class orthogonal with respect to a weighted inner-product). This is a key observation for computing Gauss quadrature nodes and weights.
1674:
1567:
14250:, which are used as nodes for the Gaussian quadrature can be found by computing the eigenvalues of this matrix. This procedure is known as
13633:
7753:{\displaystyle \prod _{\begin{smallmatrix}1\leq j\leq n\\j\neq i\end{smallmatrix}}\left(x_{i}-x_{j}\right)={\frac {p'_{n}(x_{i})}{a_{n}}}}
18972:
15069:
2214:{\displaystyle \int _{a}^{b}f(x)\,dx\approx {\frac {b-a}{2}}\sum _{i=1}^{n}w_{i}f\left({\frac {b-a}{2}}\xi _{i}+{\frac {a+b}{2}}\right).}
8286:{\displaystyle \int _{a}^{b}\omega (x){\frac {x^{k}p_{n}(x)}{x-x_{i}}}dx=x_{i}^{k}\int _{a}^{b}\omega (x){\frac {p_{n}(x)}{x-x_{i}}}dx}
3619:
9833:
18605:
13430:
5898:
1932:{\displaystyle \int _{a}^{b}f(x)\,dx=\int _{-1}^{1}f\left({\frac {b-a}{2}}\xi +{\frac {a+b}{2}}\right)\,{\frac {dx}{d\xi }}d\xi }
16478:{\displaystyle R_{n}={\frac {-n\left(n-1\right)^{3}2^{2n-1}\left^{4}}{(2n-1)\left^{3}}}f^{(2n-2)}(\xi ),\qquad -1<\xi <1.}
15684:
Stoer and
Bulirsch remark that this error estimate is inconvenient in practice, since it may be difficult to estimate the order
15348:{\displaystyle \int _{a}^{b}\omega (x)\,f(x)\,dx-\sum _{i=1}^{n}w_{i}\,f(x_{i})={\frac {f^{(2n)}(\xi )}{(2n)!}}\,(p_{n},p_{n})}
4484:
18786:
18614:
18586:
18467:
18329:
18257:
18206:
18187:
17899:
17767:
6405:{\displaystyle w_{i}={\frac {a_{n}}{a_{n-1}}}{\frac {\int _{a}^{b}\omega (x)p_{n-1}(x)^{2}dx}{p'_{n}(x_{i})p_{n-1}(x_{i})}}}
18895:
18778:
15701:
11954:{\displaystyle (p_{1},p_{0})=(x-a_{0,0})(p_{0},p_{0})=(xp_{0},p_{0})-a_{0,0}(p_{0},p_{0})=(xp_{0},p_{0})-(xp_{0},p_{0})=0.}
4076:
15675:{\displaystyle {\frac {\left(b-a\right)^{2n+1}\left(n!\right)^{4}}{(2n+1)\left^{3}}}f^{(2n)}(\xi ),\qquad a<\xi <b.}
1944:
18810:
14921:
and therefore have the same eigenvalues (the nodes). The weights can be computed from the corresponding eigenvectors: If
13836:{\displaystyle {\tilde {P}}={\begin{bmatrix}p_{0}(x)&p_{1}(x)&\cdots &p_{n-1}(x)\end{bmatrix}}^{\mathsf {T}}}
10056:
6132:
12294:{\displaystyle (p_{r+1},p_{s})=(xp_{r},p_{s})-a_{r,r}(p_{r},p_{s})-a_{r,r-1}(p_{r-1},p_{s})\cdots -a_{r,0}(p_{0},p_{s})}
17977:
10715:{\displaystyle \int _{a}^{b}\omega (x)f(x)dx=\sum _{j=1}^{n}w_{j}f(x_{j})=\sum _{j=1}^{n}\delta _{ij}w_{j}=w_{i}>0.}
1733:
1626:
17484:
17379:
12597:(which is the case for Gaussian quadrature), the recurrence relation reduces to a three-term recurrence relation: For
18502:
12689:
3592:{\displaystyle \left|\varepsilon _{t}\right|=\left|{\frac {11061.34-11058.44}{11061.34}}\right|\times 100\%=0.0262\%}
75:
53:
94:
The blue curve shows the function whose definite integral on the interval is to be calculated (the integrand). The
46:
18773:
14954:
is a normalized eigenvector (i.e., an eigenvector with euclidean norm equal to one) associated with the eigenvalue
10891:
3485:
1483:
1375:
17151:
17052:
3904:
8728:
1112:
Some low-order quadrature rules are tabulated below (over interval , see the section below for other intervals).
18770:
for abscissas and weights generation for arbitrary weighting functions W(x), integration domains and precisions.
17806:
16857:
16668:
1240:
1160:
18936:
18910:
18803:
18747:
18279:
11966:
4010:
3954:
793:
18669:
753:
18905:
18900:
18679:
4139:
4062:
868:
834:
347:
18852:
18885:
18857:
10348:
4197:
838:
18689:
contains a collection of algorithms for numerical integration (in C# / C++ / Delphi / Visual Basic / etc.)
18544:
Sagar, Robin P. (1991). "A Gaussian quadrature for the calculation of generalized Fermi-Dirac integrals".
3473:{\displaystyle {\begin{aligned}f(25.35085)&=2000\ln \left-9.8(25.35085)\\&=708.4811\end{aligned}}}
3358:{\displaystyle {\begin{aligned}f(12.64915)&=2000\ln \left-9.8(12.64915)\\&=296.8317\end{aligned}}}
126:). To obtain a more accurate result, the interval must be partitioned into many subintervals and then the
18890:
18674:
18366:
Laurie, Dirk P. (1999), "Accurate recovery of recursion coefficients from
Gaussian quadrature formulas",
18069:
3890:
533:
17289:
9361:{\displaystyle x^{n-1}=\left(x^{n-1}-{\frac {p_{n-1}(x)}{a_{n-1}}}\right)+{\frac {p_{n-1}(x)}{a_{n-1}}}}
2759:
1534:
496:{\displaystyle f(x)=\left(1-x\right)^{\alpha }\left(1+x\right)^{\beta }g(x),\quad \alpha ,\beta >-1,}
18632:
Yakimiw, E. (1996). "Accurate computation of weights in classical Gauss–Christoffel quadrature rules".
13099:
11408:
3968:
2830:
2796:
2725:
2262:
2228:
Use the two-point Gauss quadrature rule to approximate the distance in meters covered by a rocket from
16890:
16827:
13960:{\displaystyle \mathbf {e} _{n}={\begin{bmatrix}0&\cdots &0&1\end{bmatrix}}^{\mathsf {T}}}
13846:
4343:{\displaystyle \int _{a}^{b}\omega (x)\,x^{k}p_{n}(x)\,dx=0,\quad {\text{for all }}k=0,1,\ldots ,n-1.}
226:
in 1826. The most common domain of integration for such a rule is taken as , so the rule is stated as
18926:
18834:
18756:, high precision (16 and 256 decimal places) Legendre-Gaussian quadrature weights and abscissas, for
18710:
17553:
17217:
16944:
12600:
11036:
10958:
8671:
2231:
16755:
16701:
16617:
16563:
15972:{\displaystyle \int _{-1}^{1}{f(x)\,dx}={\frac {2}{n(n-1)}}+\sum _{i=2}^{n-1}{w_{i}f(x_{i})}+R_{n}.}
14896:
14264:
12532:{\displaystyle (p_{r+1},p_{s})=(xp_{r},p_{s})-a_{r,s}(p_{s},p_{s})=(xp_{r},p_{s})-(xp_{r},p_{s})=0.}
10007:
7904:{\displaystyle w_{i}={\frac {1}{p'_{n}(x_{i})}}\int _{a}^{b}\omega (x){\frac {p_{n}(x)}{x-x_{i}}}dx}
1276:
1207:
18715:
18696:
18253:
18153:
Gautschi, Walter (1970). "On the construction of
Gaussian quadrature rules from modified moments".
12544:
9674:
799:
40:
9964:
3609:
The integration problem can be expressed in a slightly more general way by introducing a positive
2719:
Next, get the weighting factors and function argument values from Table 1 for the two-point rule,
16023:
14961:, the corresponding weight can be computed from the first component of this eigenvector, namely:
14924:
11571:
10851:
10785:
10727:
10494:
9791:
9399:
8947:
7227:
4153:
18704:
18113:
12779:
11068:
6519:
3616:
into the integrand, and allowing an interval other than . That is, the problem is to calculate
856:
57:
18729:
16068:
15042:
13558:
10458:
10312:
9755:
8410:
7555:
7191:
6536:
6485:
6093:{\displaystyle \int _{a}^{b}p_{n}(x)\,\left(\prod _{i}(x-x_{i})\right)\,\omega (x)\,dx\neq 0,}
4024:
3866:{\displaystyle \left(1-x\right)^{\alpha }\left(1+x\right)^{\beta },\quad \alpha ,\beta >-1}
2024:
18826:
18273:
17757:
15146:
The error of a
Gaussian quadrature rule can be stated as follows. For an integrand which has
12025:
11010:
8298:
4606:
or less (because the degree of the remainder is always less than that of the divisor). Since
3945:
3718:
373:
219:
18951:
18641:
18624:
18553:
18431:
18397:
18170:
18145:
18109:
18060:
17785:
17529:
17193:
16920:
16731:
16593:
15991:
14222:
13594:
12655:
12304:
11667:
11640:
11613:
10931:
9728:
7595:
7528:
7501:
7164:
6458:
6431:
4870:
4132:
4054:
3763:
2445:
2409:
934:
172:
18418:
10756:
10429:
10400:
10125:
350:
rule. The quadrature rule will only be an accurate approximation to the integral above if
8:
18931:
18867:
18753:
18716:
Gaussian Quadrature Rule of Integration – Notes, PPT, Matlab, Mathematica, Maple, Mathcad
11284:{\displaystyle p_{r+1}(x)=(x-a_{r,r})p_{r}(x)-a_{r,r-1}p_{r-1}(x)\cdots -a_{r,0}p_{0}(x)}
9373:
8910:
4189:
18743:
18645:
18557:
18435:
18401:
10837:
of Gaussian quadrature rules. The most popular are the Golub-Welsch algorithm requiring
18946:
18522:
18314:
18241:
17268:
17247:
16974:
16806:
16785:
16647:
16542:
16521:
13972:
13875:
3875:
2004:
852:
377:
223:
160:
18409:
18380:
18162:
18137:
18052:
17887:
17802:
901:
18847:
18726:
18610:
18600:
18582:
18565:
18498:
18463:
18325:
18202:
18183:
17789:
17773:
17763:
17749:
5869:
will not change sign over that interval. Finally, for factors corresponding to roots
20:
17896:
Danloy, Bernard (1973). "Numerical construction of Gaussian quadrature formulas for
5850:
To prove the second part of the claim, consider the factored form of the polynomial
844:
It can be shown (see Press et al., or Stoer and Bulirsch) that the quadrature nodes
19:"Gaussian integration" redirects here. For the integral of a Gaussian function, see
18941:
18842:
18649:
18561:
18532:
18490:
18476:
18455:
18451:
18439:
18405:
18375:
18352:
18296:
18231:
18158:
18133:
18097:
18048:
17882:
17473:{\displaystyle \pm {\sqrt {{\frac {5}{11}}+{\frac {2}{11}}{\sqrt {\frac {5}{3}}}}}}
17368:{\displaystyle \pm {\sqrt {{\frac {5}{11}}-{\frac {2}{11}}{\sqrt {\frac {5}{3}}}}}}
14918:
14258:
11105:
4352:
Note that this will be true for all the orthogonal polynomials above, because each
100:
95:
88:
18783:
18774:
Gaussian Quadrature in Boost.Math, for arbitrary precision and approximation order
18236:
18219:
18128:
Gautschi, Walter (1968). "Construction of Gauss–Christoffel Quadrature Formulas".
11561:{\textstyle a_{r,s}={\frac {\left(xp_{r},p_{s}\right)}{\left(p_{s},p_{s}\right)}}}
18862:
18790:
18620:
18578:
18494:
18486:
18300:
18166:
18141:
18056:
17781:
17759:
Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables
11568:. First of all, the polynomials defined by the recurrence relation starting with
3610:
1469:{\displaystyle \pm {\sqrt {{\frac {3}{7}}+{\frac {2}{7}}{\sqrt {\frac {6}{5}}}}}}
1361:{\displaystyle \pm {\sqrt {{\frac {3}{7}}-{\frac {2}{7}}{\sqrt {\frac {6}{5}}}}}}
176:
10455:, the Gaussian quadrature formula involving the weights and nodes obtained from
10178:
7632:
7282:
7059:
6897:
6637:
18357:
18340:
15743:
4414:
which make the Gaussian quadrature computed integral exact for all polynomials
18536:
18197:
Gil, Amparo; Segura, Javier; Temme, Nico M. (2007), "§5.3: Gauss quadrature",
18101:
18966:
18321:
18258:"Ueber Gauß' neue Methode, die Werthe der Integrale näherungsweise zu finden"
11610:
have leading coefficient one and correct degree. Given the starting point by
376:
quadrature rule is not typically used for integrable functions with endpoint
18653:
18085:
17753:
3482:
Given that the true value is 11061.34 m, the absolute relative true error,
329:{\displaystyle \int _{-1}^{1}f(x)\,dx\approx \sum _{i=1}^{n}w_{i}f(x_{i}),}
16232:{\displaystyle w_{i}={\frac {2}{n(n-1)\left^{2}}},\qquad x_{i}\neq \pm 1.}
15750:
The integration points include the end points of the integration interval.
5977:
because all its roots there are now of even multiplicity. So the integral
18388:
Laurie, Dirk P. (2001). "Computation of Gauss-type quadrature formulas".
18309:
18795:
18180:
A Software Repository for Gaussian Quadratures and Christoffel Functions
17582:
An adaptive variant of this algorithm with 2 interior nodes is found in
8393:{\displaystyle {\frac {1-\left({\frac {x}{x_{i}}}\right)^{k}}{x-x_{i}}}}
4986:{\displaystyle l_{i}(x)=\prod _{j\neq i}{\frac {x-x_{j}}{x_{i}-x_{j}}}.}
516:
is well-approximated by a low-degree polynomial, then alternative nodes
18443:
18245:
18215:
17583:
17140:{\displaystyle \pm {\sqrt {{\frac {1}{3}}+{\frac {2{\sqrt {7}}}{21}}}}}
17041:{\displaystyle \pm {\sqrt {{\frac {1}{3}}-{\frac {2{\sqrt {7}}}{21}}}}}
15746:. It is similar to Gaussian quadrature with the following differences:
5847:
or less, its integral is given exactly by the Gaussian quadrature sum.
873:
180:
18068:
Eaton, John W.; Bateman, David; Hauberg, Søren; Wehbring, Rik (2018).
14257:
For computing the weights and nodes, it is preferable to consider the
12301:
all scalar products vanish except for the first one and the one where
1719:{\displaystyle \pm {\frac {1}{3}}{\sqrt {5+2{\sqrt {\frac {10}{7}}}}}}
1612:{\displaystyle \pm {\frac {1}{3}}{\sqrt {5-2{\sqrt {\frac {10}{7}}}}}}
18734:
18686:
2391:{\displaystyle s=\int _{8}^{30}{\left(2000\ln \left-9.8t\right){dt}}}
1103:{\displaystyle w_{i}={\frac {2}{\left(1-x_{i}^{2}\right)\left^{2}}}.}
532:
will usually give more accurate quadrature rules. These are known as
155:), the 2-point Gaussian quadrature rule even returns an exact result.
130:
trapezoidal rule must be used, which requires many more calculations.
18692:
18603:; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W. (eds.),
17793:
13707:{\displaystyle J{\tilde {P}}=x{\tilde {P}}-p_{n}(x)\mathbf {e} _{n}}
18700:
18527:
11459:
is the maximal degree which can be taken to be infinity, and where
18724:
12961:{\displaystyle p_{r+1}(x)=(x-a_{r,r})p_{r}(x)-a_{r,r-1}p_{r-1}(x)}
12679:
is orthogonal to every polynomial of degree less than or equal to
5590:{\displaystyle h(x_{i})=p_{n}(x_{i})\,q(x_{i})+r(x_{i})=r(x_{i}),}
18341:"Computing Gaussian quadrature rules with high relative accuracy"
13630:
The three-term recurrence relation can be written in matrix form
4613:
is by assumption orthogonal to all monomials of degree less than
18450:
18339:
Laudadio, Teresa; Mastronardi, Nicola; Van Dooren, Paul (2023).
18115:
Methodus nova integralium valores per approximationem inveniendi
17777:
9954:{\displaystyle p_{n+1}(x_{i})=(a)p_{n}(x_{i})+(b)p_{n-1}(x_{i})}
17587:
13548:{\displaystyle (xp_{r},p_{r-1})=(p_{r},xp_{r-1})=(p_{r},p_{r})}
4581:
or less (because the sum of its degree and that of the divisor
92:
Comparison between 2-point Gaussian and trapezoidal quadrature.
18456:"Section 4.6. Gaussian Quadratures and Orthogonal Polynomials"
15722:-point rule in such a way that the resulting rule is of order
15030:{\displaystyle w_{j}=\mu _{0}\left(\phi _{1}^{(j)}\right)^{2}}
18596:
18338:
15711:
are extensions of Gauss quadrature rules generated by adding
11396:{\displaystyle (f(x),g(x))=\int _{a}^{b}\omega (x)f(x)g(x)dx}
9752:
can also be expressed in terms of the orthogonal polynomials
8889:{\displaystyle {\frac {p_{n}(x)}{x-x_{i}}}=a_{n}x^{n-1}+s(x)}
7762:
We can thus write the integral expression for the weights as
2223:
18754:
Tabulated weights and abscissae with Mathematica source code
13087:{\displaystyle p_{r+1}(x)=(x-a_{r})p_{r}(x)-b_{r}p_{r-1}(x)}
5451:(see below for other formulas for the weights). But all the
18067:
17733:
10818:
18289:
Journal of the Serbian Society for Computational Mechanics
15475:{\displaystyle (f,g)=\int _{a}^{b}\omega (x)f(x)g(x)\,dx.}
5969:
This polynomial cannot change sign over the interval from
4359:
is constructed to be orthogonal to the other polynomials
887:
For the simplest integration problem stated above, i.e.,
18118:. Comm. Soc. Sci. Göttingen Math. Vol. 3. S. 29–76.
11108:
orthogonal polynomials) satisfy the recurrence relation
8668:
We can evaluate the integral on the right hand side for
7188:
can be expressed in terms of the orthogonal polynomials
18120:
datiert 1814, auch in Werke, Band 3, 1876, S. 163–196.
17967:{\displaystyle \int _{0}^{1}(-\log x)x^{\alpha }f(x)dx}
5083:{\displaystyle r(x)=\sum _{i=1}^{n}l_{i}(x)\,r(x_{i}).}
18768:
Mathematica source code distributed under the GNU LGPL
18462:(3rd ed.), New York: Cambridge University Press,
18454:; Teukolsky, SA; Vetterling, WT; Flannery, BP (2007),
14613:
13995:
13920:
13742:
11465:
10117:
10059:
802:
756:
18475:
17980:
17902:
17809:
17700:
17556:
17532:
17487:
17427:
17382:
17322:
17292:
17271:
17250:
17220:
17196:
17154:
17097:
17055:
16998:
16977:
16947:
16923:
16893:
16860:
16830:
16809:
16788:
16758:
16734:
16704:
16671:
16650:
16620:
16596:
16566:
16545:
16524:
16247:
16114:
16071:
16026:
15994:
15786:
15501:
15392:
15159:
15072:
15045:
14967:
14927:
14899:
14597:
14291:
14267:
14225:
13981:
13898:
13878:
13849:
13720:
13636:
13597:
13561:
13433:
13147:
13102:
12976:
12832:
12782:
12692:
12658:
12603:
12547:
12334:
12307:
12061:
12028:
11969:
11702:
11670:
11643:
11616:
11574:
11411:
11299:
11114:
11071:
11039:
11013:
10961:
10934:
10854:
10788:
10759:
10730:
10547:
10497:
10461:
10432:
10403:
10351:
10315:
10156:
10128:
10010:
9967:
9836:
9794:
9758:
9731:
9677:
9444:
9402:
9376:
9212:
8992:
8950:
8913:
8798:
8731:
8674:
8460:
8413:
8327:
8301:
8102:
7933:
7773:
7625:
7598:
7558:
7531:
7504:
7275:
7230:
7194:
7167:
7012:
6775:
6575:
6539:
6488:
6461:
6434:
6224:
6135:
5983:
5901:
5607:
5471:
5098:
5001:
4894:
4634:
4487:
4227:
4156:
4122:{\displaystyle x^{\alpha }e^{-x},\quad \alpha >-1}
4079:
4027:
3971:
3907:
3788:
3622:
3523:
3488:
3370:
3255:
2869:
2833:
2799:
2762:
2728:
2489:
2448:
2412:
2299:
2265:
2234:
2059:
2027:
2007:
1947:
1791:
1736:
1677:
1629:
1570:
1537:
1486:
1423:
1378:
1315:
1279:
1243:
1210:
1163:
997:
904:
833:. One may also want to integrate over semi-infinite (
544:
388:
232:
10107:{\textstyle {\frac {1}{b}}p_{n+1}\left(x_{i}\right)}
1992:{\displaystyle {\frac {dx}{d\xi }}={\frac {b-a}{2}}}
16:
Approximation of the definite integral of a function
15125:{\displaystyle \mu _{0}=\int _{a}^{b}\omega (x)dx.}
9370:The term in the brackets is a polynomial of degree
18511:
18460:Numerical Recipes: The Art of Scientific Computing
18313:
18035:
17966:
17863:
17569:
17541:
17515:
17472:
17410:
17367:
17305:
17277:
17256:
17233:
17205:
17179:
17139:
17080:
17040:
16983:
16960:
16932:
16906:
16878:
16843:
16815:
16794:
16771:
16743:
16717:
16689:
16656:
16633:
16605:
16579:
16551:
16530:
16477:
16231:
16093:
16057:
16012:
15971:
15674:
15474:
15347:
15124:
15058:
15029:
14946:
14909:
14878:
14579:
14277:
14238:
14209:
13959:
13884:
13864:
13835:
13706:
13610:
13583:
13547:
13417:
13133:
13086:
12960:
12807:
12768:
12671:
12637:
12589:
12531:
12320:
12293:
12047:
12014:
11953:
11688:
11656:
11629:
11602:
11560:
11447:
11395:
11283:
11096:
11057:
11025:
10999:
10947:
10882:
10823:There are many algorithms for computing the nodes
10807:
10774:
10745:
10714:
10525:
10483:
10447:
10418:
10389:
10337:
10294:
10143:
10106:
10045:
9996:
9953:
9822:
9780:
9744:
9709:
9655:
9430:
9388:
9360:
9196:
8978:
8925:
8888:
8777:
8717:
8658:
8435:
8392:
8313:
8285:
8086:
7903:
7752:
7611:
7580:
7544:
7517:
7488:
7258:
7216:
7180:
7151:
6989:
6723:
6561:
6510:
6474:
6447:
6404:
6208:
6170:
6092:
5959:
5818:
5589:
5403:
5082:
4985:
4837:
4555:
4342:
4179:
4121:
4044:
3994:
3935:
3865:
3672:
3591:
3509:
3472:
3357:
3241:
2852:
2818:
2784:
2747:
2709:
2473:
2434:
2390:
2285:
2251:
2213:
2053:rule then results in the following approximation:
2045:
2013:
1991:
1931:
1764:
1718:
1657:
1611:
1550:
1511:
1468:
1403:
1360:
1292:
1261:
1223:
1181:
1102:
925:
825:
784:
740:
495:
328:
18036:{\displaystyle \int _{0}^{\infty }E_{m}(x)f(x)dx}
12328:meets the same orthogonal polynomial. Therefore,
3673:{\displaystyle \int _{a}^{b}\omega (x)\,f(x)\,dx}
1765:{\displaystyle {\frac {322-13{\sqrt {70}}}{900}}}
1658:{\displaystyle {\frac {322+13{\sqrt {70}}}{900}}}
18964:
18307:
18132:. Vol. 22, no. 102. pp. 251–270.
18047:. Vol. 27, no. 124. pp. 861–869.
17748:
17688:
17676:
17640:
17516:{\displaystyle {\frac {124-7{\sqrt {15}}}{350}}}
17411:{\displaystyle {\frac {124+7{\sqrt {15}}}{350}}}
16105:-th degree and the dash denotes the derivative.
12645:is a polynomial of degree less than or equal to
9671:), the weights are obtained by dividing this by
4474:or less. Divide it by the orthogonal polynomial
103:gives an approximation of the integral equal to
18711:From Lobatto Quadrature to the Euler constant e
18286:
18262:Journal für die Reine und Angewandte Mathematik
12769:{\displaystyle (xp_{r},p_{s})=(p_{r},xp_{s})=0}
6197:distinct roots, all real, in the interval from
5435:, is defined to equal the weighted integral of
2406:First, changing the limits of integration from
361:is well-approximated by a polynomial of degree
18196:
18083:
17756:, eds. (1983) . "Chapter 25.4, Integration".
17711:
15378:is the monic (i.e. the leading coefficient is
15135:
10905:operations, and asymptotic formulas for large
8454:is a polynomial of at most nth degree we have
6746:and is thus fixed by the values it attains at
5960:{\displaystyle p_{n}(x)\,\prod _{i}(x-x_{i}).}
5465:, so the division formula above tells us that
3771:
380:. Instead, if the integrand can be written as
18811:
18572:
17664:
17652:
13625:
12826:. The recurrence relation then simplifies to
4865:or less, we can interpolate it exactly using
4770:
4715:
3510:{\displaystyle \left|\varepsilon _{t}\right|}
1512:{\displaystyle {\frac {18-{\sqrt {30}}}{36}}}
1404:{\displaystyle {\frac {18+{\sqrt {30}}}{36}}}
17180:{\displaystyle {\frac {14-{\sqrt {7}}}{30}}}
17081:{\displaystyle {\frac {14+{\sqrt {7}}}{30}}}
16101:denotes the standard Legendre polynomial of
15753:It is accurate for polynomials up to degree
10782:are non-negative functions, it follows that
10122:Consider the following polynomial of degree
9717:and that yields the expression in equation (
6750:different points. Multiplying both sides by
5895:by one more factor to make a new polynomial
3936:{\displaystyle {\frac {1}{\sqrt {1-x^{2}}}}}
2863:Now we can use the Gauss quadrature formula
933:, the associated orthogonal polynomials are
862:
18214:
8778:{\displaystyle {\frac {p_{n}(x)}{x-x_{i}}}}
6122:is orthogonal to all polynomials of degree
4455:To prove the first part of this claim, let
18818:
18804:
18479:; Sacco, Riccardo; Saleri, Fausto (2000).
17864:{\displaystyle \int _{0}^{1}-\ln(x)f(x)dx}
16879:{\displaystyle \pm {\sqrt {\frac {3}{7}}}}
16690:{\displaystyle \pm {\sqrt {\frac {1}{5}}}}
2224:Example of two-point Gauss quadrature rule
1262:{\displaystyle \pm {\sqrt {\frac {3}{5}}}}
1182:{\displaystyle \pm {\frac {1}{\sqrt {3}}}}
191:or less by a suitable choice of the nodes
18825:
18526:
18430:. Vol. 5, no. 1. pp. 1–9.
18419:"Numerical integration - MATLAB integral"
18416:
18379:
18356:
18235:
18084:Gander, Walter; Gautschi, Walter (2000).
17886:
17722:
15818:
15462:
15315:
15241:
15200:
15187:
14550:
12541:However, if the scalar product satisfies
12015:{\displaystyle p_{0},p_{1},\ldots ,p_{r}}
6074:
6061:
6018:
5921:
5793:
5739:
5698:
5685:
5648:
5635:
5520:
5387:
5367:
5336:
5316:
5288:
5268:
5237:
5217:
5176:
5139:
5126:
5057:
4825:
4812:
4775:
4767:
4720:
4712:
4675:
4662:
4556:{\displaystyle h(x)=p_{n}(x)\,q(x)+r(x).}
4522:
4285:
4255:
4041:
3663:
3650:
2087:
1902:
1819:
991:and the weights are given by the formula
785:{\textstyle {\frac {1}{\sqrt {1-x^{2}}}}}
667:
575:
346:or less. This exact rule is known as the
338:which is exact for polynomials of degree
263:
179:constructed to yield an exact result for
76:Learn how and when to remove this message
18784:Nodes and Weights of Gaussian quadrature
18425:
18177:
18152:
18127:
18070:"Functions of One Variable (GNU Octave)"
17800:
10848:operations, Newton's method for solving
10819:Computation of Gaussian quadrature rules
4617:, it must be orthogonal to the quotient
872:
855:of a polynomial belonging to a class of
87:
39:This article includes a list of general
18631:
18606:NIST Handbook of Mathematical Functions
18573:Stoer, Josef; Bulirsch, Roland (2002),
18220:"Calculation of Gauss Quadrature Rules"
18199:Numerical Methods for Special Functions
15066:is the integral of the weight function
14542:
5884:that are of odd multiplicity, multiply
898:is well-approximated by polynomials on
18965:
18779:Gauss–Kronrod Quadrature in Boost.Math
18387:
18365:
18252:
17895:
17629:
15733:
15695:
13951:
13827:
10923:
9438:. The integral can thus be written as
6129:or less, so the degree of the product
5424:, the weight associated with the node
4588:must equal that of the dividend), and
4433:or less. Furthermore, all these nodes
4205:
877:Graphs of Legendre polynomials (up to
18799:
18725:
18594:
18543:
18108:
17618:
10390:{\displaystyle f(x_{j})=\delta _{ij}}
4217:be a nontrivial polynomial of degree
1780:
18720:Holistic Numerical Methods Institute
15765:is the number of integration points.
9830:. In the 3-term recurrence relation
7764:
6215:
5828:This proves that for any polynomial
25:
17701:Quarteroni, Sacco & Saleri 2000
10118:Proof that the weights are positive
9396:, which is therefore orthogonal to
6171:{\displaystyle \prod _{i}(x-x_{i})}
4004:Chebyshev polynomials (second kind)
13:
18973:Numerical integration (quadrature)
18764:=64, with Mathematica source code.
18575:Introduction to Numerical Analysis
18157:. Vol. 24. pp. 245–260.
17991:
17803:"Gaussian quadrature formulae for
17306:{\displaystyle {\frac {256}{525}}}
15742:, named after Dutch mathematician
15382:) orthogonal polynomial of degree
15141:
14902:
14600:
14427:
14398:
14299:
14270:
11104:and leading coefficient one (i.e.
3586:
3577:
2785:{\displaystyle x_{1}=-0.577350269}
2276:
2245:
1551:{\displaystyle {\frac {128}{225}}}
957:-th polynomial normalized to give
45:it lacks sufficient corresponding
14:
18984:
18746:by Chris Maes and Anton Antonov,
18662:
18163:10.1090/S0025-5718-1970-0285117-6
18138:10.1090/S0025-5718-1968-0228171-0
18086:"Adaptive Quadrature - Revisited"
18053:10.1090/S0025-5718-1973-0331730-X
17888:10.1090/s0025-5718-1965-0178569-1
15484:In the important special case of
13134:{\displaystyle p_{-1}(x)\equiv 0}
11448:{\displaystyle r=0,1,\ldots ,n-1}
10177:
7631:
7281:
7058:
6896:
6636:
6518:. To prove this, note that using
4599:is the remainder, also of degree
3995:{\displaystyle {\sqrt {1-x^{2}}}}
2853:{\displaystyle x_{2}=0.577350269}
2819:{\displaystyle c_{2}=1.000000000}
2748:{\displaystyle c_{1}=1.000000000}
2286:{\displaystyle t=30\mathrm {s} ,}
18896:Gauss–Kronrod quadrature formula
16907:{\displaystyle {\frac {49}{90}}}
16844:{\displaystyle {\frac {32}{45}}}
15702:Gauss–Kronrod quadrature formula
14246:of the polynomials up to degree
13983:
13901:
13892:th standard basis vector, i.e.,
13865:{\displaystyle \mathbf {e} _{n}}
13852:
13694:
10309:are the roots of the polynomial
6213:The weights can be expressed as
3772:§ Gauss–Legendre quadrature
30:
17741:
17727:
17716:
17570:{\displaystyle {\frac {1}{21}}}
17234:{\displaystyle {\frac {1}{15}}}
16961:{\displaystyle {\frac {1}{10}}}
16456:
16209:
15769:Lobatto quadrature of function
15653:
13233:
12638:{\displaystyle s<r-1,xp_{s}}
11664:can be shown by induction. For
11058:{\displaystyle (\cdot ,\cdot )}
11000:{\displaystyle (p_{r},p_{s})=0}
8718:{\displaystyle q(x)=p_{n-1}(x)}
6209:General formula for the weights
4301:
4106:
3844:
2252:{\displaystyle t=8\mathrm {s} }
471:
18748:Wolfram Demonstrations Project
18609:, Cambridge University Press,
18316:Numerical Methods and Software
18268:. S. 301–308und Werke, Band 6.
18024:
18018:
18012:
18006:
17955:
17949:
17933:
17918:
17852:
17846:
17840:
17834:
17705:
17694:
17682:
17677:Kahaner, Moler & Nash 1989
17670:
17658:
17646:
17634:
17623:
17612:
16772:{\displaystyle {\frac {1}{6}}}
16718:{\displaystyle {\frac {5}{6}}}
16634:{\displaystyle {\frac {1}{3}}}
16580:{\displaystyle {\frac {4}{3}}}
16450:
16444:
16439:
16424:
16371:
16356:
16149:
16137:
16088:
16082:
16052:
16046:
16007:
15995:
15949:
15936:
15889:
15886:
15877:
15868:
15862:
15856:
15850:
15838:
15815:
15809:
15647:
15641:
15636:
15627:
15580:
15565:
15459:
15453:
15447:
15441:
15435:
15429:
15405:
15393:
15342:
15316:
15306:
15297:
15292:
15286:
15281:
15272:
15258:
15245:
15197:
15191:
15184:
15178:
15110:
15104:
15012:
15006:
14939:
14933:
14910:{\displaystyle {\mathcal {J}}}
14278:{\displaystyle {\mathcal {J}}}
13814:
13808:
13782:
13776:
13761:
13755:
13727:
13689:
13683:
13664:
13646:
13542:
13516:
13510:
13475:
13469:
13434:
13409:
13371:
13366:
13340:
13328:
13290:
13285:
13250:
13224:
13198:
13193:
13164:
13122:
13116:
13081:
13075:
13043:
13037:
13024:
13005:
12999:
12993:
12955:
12949:
12905:
12899:
12886:
12861:
12855:
12849:
12757:
12728:
12722:
12693:
12584:
12569:
12563:
12548:
12520:
12491:
12485:
12456:
12450:
12424:
12402:
12373:
12367:
12335:
12288:
12262:
12237:
12205:
12177:
12151:
12129:
12100:
12094:
12062:
11942:
11913:
11907:
11878:
11872:
11846:
11824:
11795:
11789:
11763:
11760:
11735:
11729:
11703:
11591:
11585:
11384:
11378:
11372:
11366:
11360:
11354:
11330:
11327:
11321:
11312:
11306:
11300:
11278:
11272:
11237:
11231:
11187:
11181:
11168:
11143:
11137:
11131:
11085:
11072:
11052:
11040:
10988:
10962:
10871:
10865:
10769:
10763:
10740:
10734:
10643:
10630:
10584:
10578:
10572:
10566:
10514:
10501:
10478:
10472:
10413:
10407:
10368:
10355:
10332:
10326:
10166:
10160:
10053:in Eq. (1) can be replaced by
10046:{\displaystyle p_{n-1}(x_{i})}
10040:
10027:
9991:
9978:
9948:
9935:
9916:
9910:
9904:
9891:
9878:
9872:
9866:
9853:
9817:
9811:
9775:
9769:
9704:
9691:
9638:
9631:
9612:
9606:
9582:
9569:
9491:
9485:
9469:
9463:
9425:
9419:
9336:
9330:
9281:
9275:
9169:
9163:
9144:
9138:
9114:
9101:
9039:
9033:
9017:
9011:
8973:
8967:
8883:
8877:
8818:
8812:
8751:
8745:
8712:
8706:
8684:
8678:
8626:
8620:
8607:
8601:
8592:
8586:
8562:
8549:
8507:
8501:
8485:
8479:
8430:
8424:
8253:
8247:
8231:
8225:
8159:
8153:
8127:
8121:
7871:
7865:
7849:
7843:
7819:
7806:
7734:
7721:
7619:yields using L'Hôpital's rule
7575:
7569:
7442:
7436:
7253:
7247:
7211:
7205:
7050:
7044:
6888:
6882:
6861:
6848:
6812:
6806:
6800:
6794:
6628:
6615:
6585:
6579:
6556:
6543:
6505:
6499:
6396:
6383:
6364:
6351:
6321:
6314:
6295:
6289:
6165:
6146:
6071:
6065:
6053:
6034:
6015:
6009:
5951:
5932:
5918:
5912:
5810:
5797:
5756:
5743:
5695:
5689:
5682:
5676:
5645:
5639:
5632:
5626:
5581:
5568:
5559:
5546:
5537:
5524:
5517:
5504:
5488:
5475:
5384:
5371:
5333:
5327:
5313:
5307:
5285:
5272:
5234:
5221:
5214:
5208:
5173:
5167:
5136:
5130:
5123:
5117:
5074:
5061:
5054:
5048:
5011:
5005:
4911:
4905:
4822:
4816:
4809:
4803:
4764:
4758:
4749:
4743:
4737:
4731:
4709:
4703:
4672:
4666:
4659:
4653:
4547:
4541:
4532:
4526:
4519:
4513:
4497:
4491:
4440:will lie in the open interval
4282:
4276:
4252:
4246:
3748:For more information, see ...
3660:
3654:
3647:
3641:
3604:
3450:
3444:
3428:
3422:
3384:
3378:
3335:
3329:
3313:
3307:
3269:
3263:
3214:
3208:
3202:
3196:
3170:
3164:
3155:
3149:
3109:
3103:
3075:
3066:
2519:
2513:
2084:
2078:
2040:
2028:
1816:
1810:
1293:{\displaystyle {\frac {5}{9}}}
1224:{\displaystyle {\frac {8}{9}}}
1080:
1067:
920:
905:
826:{\textstyle {\sqrt {1-x^{2}}}}
664:
658:
572:
566:
465:
459:
398:
392:
320:
307:
260:
254:
1:
18410:10.1016/S0377-0427(00)00506-9
18381:10.1016/S0377-0427(99)00228-9
18237:10.1090/S0025-5718-69-99647-1
17601:
15495:, we have the error estimate
12590:{\displaystyle (xf,g)=(f,xg)}
9710:{\displaystyle p'_{n}(x_{i})}
7161:This integral expression for
5092:Then its integral will equal
218:The modern formulation using
114:, while the correct value is
18566:10.1016/0010-4655(91)90076-W
18495:10.1007/978-3-540-49809-4_10
18301:10.24874/jsscm.2017.11.01.02
17801:Anderson, Donald G. (1965).
17689:Abramowitz & Stegun 1983
17641:Abramowitz & Stegun 1983
17606:
15136:Gil, Segura & Temme 2007
13975:, called the Jacobi matrix:
9997:{\displaystyle p_{n}(x_{i})}
8407:which is then orthogonal to
6111:is always non-negative. But
4466:be any polynomial of degree
4382:is in the span of that set.
7:
18730:"Legendre-Gauss Quadrature"
18675:Encyclopedia of Mathematics
18597:"§3.5(v): Gauss Quadrature"
16058:{\displaystyle P'_{n-1}(x)}
14947:{\displaystyle \phi ^{(j)}}
11293:and scalar product defined
9719:
9667:
7917:
6418:
4574:is the quotient, of degree
10:
18989:
18937:Clenshaw–Curtis quadrature
18911:Chebyshev–Gauss quadrature
18670:"Gauss quadrature formula"
18358:10.1007/s11075-022-01297-9
18278:: CS1 maint: postscript (
18224:Mathematics of Computation
18218:; Welsch, John H. (1969).
17712:Gander & Gautschi 2000
15699:
13626:The Golub-Welsch algorithm
12022:are orthogonal, then also
11603:{\displaystyle p_{0}(x)=1}
10883:{\displaystyle p_{n}(x)=0}
10808:{\displaystyle w_{i}>0}
10746:{\displaystyle \omega (x)}
10526:{\displaystyle f(x_{j})=0}
9823:{\displaystyle p_{n+1}(x)}
9431:{\displaystyle p_{n-1}(x)}
8979:{\displaystyle p_{n-1}(x)}
8907:is a polynomial of degree
8785:is a polynomial of degree
8400:is a polynomial of degree
7927:In the integrand, writing
7259:{\displaystyle p_{n-1}(x)}
6100:since the weight function
4869:interpolation points with
4180:{\displaystyle e^{-x^{2}}}
4011:Chebyshev–Gauss quadrature
3955:Chebyshev–Gauss quadrature
2021:point Gaussian quadrature
866:
837:) and infinite intervals (
18:
18919:
18906:Gauss–Legendre quadrature
18901:Gauss–Laguerre quadrature
18876:
18858:Adaptive Simpson's method
18833:
18537:10.1016/j.jco.2017.10.002
18178:Gautschi, Walter (2020).
18090:BIT Numerical Mathematics
17665:Stoer & Bulirsch 2002
17653:Stoer & Bulirsch 2002
17244:
16971:
16782:
16644:
16518:
16487:Some of the weights are:
12808:{\displaystyle a_{r,s}=0}
11097:{\displaystyle (p_{r})=r}
10894:for evaluation requiring
5876:inside the interval from
4140:Gauss–Laguerre quadrature
4063:Gauss–Laguerre quadrature
1528:
1525:
1306:
1201:
1198:
1193:
1149:
1146:
1132:
1123:
869:Gauss–Legendre quadrature
863:Gauss–Legendre quadrature
835:Gauss–Laguerre quadrature
348:Gauss–Legendre quadrature
18886:Gauss–Hermite quadrature
18312:; Nash, Stephen (1989).
16094:{\displaystyle P_{m}(x)}
15153:continuous derivatives,
15059:{\displaystyle \mu _{0}}
13584:{\displaystyle xp_{r-1}}
10484:{\displaystyle p_{n}(x)}
10338:{\displaystyle p_{n}(x)}
9781:{\displaystyle p_{n}(x)}
8436:{\displaystyle p_{n}(x)}
7581:{\displaystyle p_{n}(x)}
7217:{\displaystyle p_{n}(x)}
6562:{\displaystyle r(x_{i})}
6511:{\displaystyle p_{k}(x)}
4198:Gauss–Hermite quadrature
4045:{\displaystyle e^{-x}\,}
2046:{\displaystyle (\xi ,w)}
839:Gauss–Hermite quadrature
169:Gaussian quadrature rule
18891:Gauss–Jacobi quadrature
18595:Temme, Nico M. (2010),
18102:10.1023/A:1022318402393
15138:) for further details.
12048:{\displaystyle p_{r+1}}
11637:, the orthogonality of
11026:{\displaystyle r\neq s}
10928:Orthogonal polynomials
9665:According to equation (
8314:{\displaystyle k\leq n}
5601:. Thus we finally have
3891:Gauss–Jacobi quadrature
3742:Orthogonal polynomials
750:Common weights include
534:Gauss–Jacobi quadrature
60:more precise citations.
18705:GNU Scientific Library
18693:GNU Scientific Library
18654:10.1006/jcph.1996.0258
18368:J. Comput. Appl. Math.
18037:
17968:
17865:
17571:
17543:
17517:
17474:
17412:
17369:
17307:
17279:
17258:
17235:
17207:
17181:
17141:
17082:
17042:
16985:
16962:
16934:
16908:
16880:
16845:
16817:
16796:
16773:
16745:
16719:
16691:
16658:
16635:
16607:
16581:
16553:
16532:
16479:
16233:
16095:
16059:
16014:
15973:
15921:
15676:
15476:
15349:
15230:
15126:
15060:
15031:
14948:
14911:
14880:
14581:
14279:
14252:Golub–Welsch algorithm
14240:
14211:
13961:
13886:
13866:
13837:
13708:
13618:by a degree less than
13612:
13585:
13549:
13419:
13135:
13088:
12962:
12809:
12770:
12673:
12639:
12591:
12533:
12322:
12295:
12049:
12016:
11955:
11690:
11658:
11631:
11604:
11562:
11449:
11397:
11285:
11098:
11059:
11027:
11001:
10949:
10884:
10809:
10776:
10747:
10716:
10669:
10616:
10527:
10485:
10449:
10420:
10397:. Since the degree of
10391:
10339:
10296:
10145:
10108:
10047:
9998:
9955:
9824:
9782:
9746:
9711:
9657:
9432:
9390:
9362:
9198:
8980:
8927:
8890:
8779:
8719:
8660:
8437:
8394:
8315:
8287:
8088:
7905:
7754:
7613:
7588:. Taking the limit of
7582:
7546:
7525:is the coefficient of
7519:
7490:
7260:
7218:
7182:
7153:
6991:
6844:
6725:
6611:
6563:
6520:Lagrange interpolation
6512:
6476:
6455:is the coefficient of
6449:
6406:
6172:
6094:
5961:
5820:
5782:
5728:
5591:
5405:
5366:
5267:
5197:
5084:
5037:
4987:
4839:
4557:
4344:
4181:
4123:
4046:
3996:
3937:
3867:
3674:
3593:
3511:
3474:
3359:
3243:
2854:
2820:
2786:
2749:
2711:
2475:
2436:
2392:
2287:
2253:
2215:
2135:
2047:
2015:
1993:
1933:
1766:
1720:
1659:
1613:
1552:
1513:
1470:
1405:
1362:
1294:
1263:
1225:
1183:
1104:
927:
884:
857:orthogonal polynomials
827:
786:
742:
697:
497:
330:
293:
220:orthogonal polynomials
156:
18927:Barnes–Hut simulation
18835:Newton–Cotes formulas
18827:Numerical integration
18703:algorithms (see also
18515:Journal of Complexity
18482:Numerical Mathematics
18390:J. Comput. Appl. Math
18110:Gauss, Carl Friedrich
18038:
17969:
17866:
17572:
17544:
17542:{\displaystyle \pm 1}
17518:
17475:
17413:
17370:
17308:
17280:
17259:
17236:
17208:
17206:{\displaystyle \pm 1}
17182:
17142:
17083:
17043:
16986:
16963:
16935:
16933:{\displaystyle \pm 1}
16909:
16881:
16846:
16818:
16797:
16774:
16746:
16744:{\displaystyle \pm 1}
16720:
16692:
16659:
16636:
16608:
16606:{\displaystyle \pm 1}
16582:
16554:
16533:
16480:
16234:
16096:
16060:
16015:
16013:{\displaystyle (i-1)}
15974:
15895:
15677:
15477:
15350:
15210:
15127:
15061:
15032:
14949:
14912:
14881:
14582:
14280:
14241:
14239:{\displaystyle x_{j}}
14212:
13962:
13887:
13867:
13838:
13709:
13613:
13611:{\displaystyle p_{r}}
13586:
13550:
13427:(the last because of
13420:
13136:
13096:(with the convention
13089:
12963:
12810:
12771:
12686:. Therefore, one has
12674:
12672:{\displaystyle p_{r}}
12652:. On the other hand,
12640:
12592:
12534:
12323:
12321:{\displaystyle p_{s}}
12296:
12050:
12017:
11956:
11691:
11689:{\displaystyle r=s=0}
11659:
11657:{\displaystyle p_{r}}
11632:
11630:{\displaystyle p_{0}}
11605:
11563:
11450:
11398:
11286:
11099:
11060:
11033:for a scalar product
11028:
11002:
10950:
10948:{\displaystyle p_{r}}
10892:three-term recurrence
10885:
10810:
10777:
10748:
10717:
10649:
10596:
10528:
10486:
10450:
10421:
10392:
10340:
10302:where, as above, the
10297:
10146:
10109:
10048:
9999:
9956:
9825:
9783:
9747:
9745:{\displaystyle w_{i}}
9712:
9658:
9433:
9391:
9363:
9199:
8981:
8928:
8891:
8780:
8720:
8661:
8438:
8395:
8316:
8288:
8089:
7906:
7755:
7614:
7612:{\displaystyle x_{i}}
7583:
7547:
7545:{\displaystyle x^{n}}
7520:
7518:{\displaystyle a_{n}}
7491:
7261:
7219:
7183:
7181:{\displaystyle w_{i}}
7154:
6992:
6824:
6761:and integrating from
6742:has degree less than
6726:
6591:
6564:
6513:
6477:
6475:{\displaystyle x^{k}}
6450:
6448:{\displaystyle a_{k}}
6407:
6173:
6095:
5962:
5821:
5762:
5708:
5592:
5406:
5346:
5247:
5177:
5085:
5017:
4988:
4840:
4558:
4345:
4182:
4124:
4047:
3997:
3946:Chebyshev polynomials
3938:
3868:
3719:Abramowitz and Stegun
3675:
3594:
3512:
3475:
3360:
3244:
2855:
2821:
2787:
2750:
2712:
2476:
2474:{\displaystyle \left}
2437:
2435:{\displaystyle \left}
2393:
2288:
2254:
2216:
2115:
2048:
2016:
1994:
1934:
1767:
1721:
1660:
1614:
1553:
1514:
1471:
1406:
1363:
1295:
1264:
1226:
1184:
1105:
928:
876:
867:Further information:
828:
787:
743:
677:
498:
331:
273:
91:
18952:Tanh-sinh quadrature
18546:Comput. Phys. Commun
18489:. pp. 425–478.
18345:Numerical Algorithms
17978:
17900:
17807:
17554:
17530:
17485:
17425:
17380:
17320:
17290:
17269:
17248:
17218:
17194:
17152:
17095:
17053:
16996:
16975:
16945:
16921:
16891:
16858:
16828:
16807:
16786:
16756:
16732:
16702:
16669:
16648:
16618:
16594:
16564:
16543:
16522:
16245:
16112:
16069:
16024:
15992:
15784:
15499:
15390:
15157:
15134:See, for instance, (
15070:
15043:
14965:
14925:
14897:
14595:
14551:
14289:
14265:
14223:
13979:
13896:
13876:
13847:
13718:
13634:
13595:
13559:
13431:
13145:
13100:
12974:
12830:
12780:
12690:
12656:
12601:
12545:
12332:
12305:
12059:
12026:
11967:
11700:
11668:
11641:
11614:
11572:
11463:
11409:
11297:
11112:
11069:
11037:
11011:
10959:
10932:
10852:
10786:
10775:{\displaystyle f(x)}
10757:
10728:
10545:
10495:
10459:
10448:{\displaystyle 2n-1}
10430:
10419:{\displaystyle f(x)}
10401:
10349:
10313:
10154:
10144:{\displaystyle 2n-2}
10126:
10057:
10008:
9965:
9834:
9792:
9756:
9729:
9675:
9442:
9400:
9374:
9210:
8990:
8948:
8911:
8796:
8729:
8725:as follows. Because
8672:
8458:
8411:
8325:
8299:
8100:
7931:
7771:
7623:
7596:
7556:
7529:
7502:
7273:
7228:
7192:
7165:
7010:
6773:
6573:
6537:
6486:
6459:
6432:
6222:
6133:
5981:
5899:
5605:
5469:
5096:
4999:
4892:
4871:Lagrange polynomials
4847:Since the remainder
4632:
4485:
4225:
4154:
4133:Laguerre polynomials
4077:
4055:Laguerre polynomials
4025:
3969:
3905:
3786:
3764:Legendre polynomials
3680:for some choices of
3620:
3521:
3486:
3368:
3253:
2867:
2831:
2797:
2760:
2726:
2487:
2446:
2410:
2297:
2263:
2232:
2057:
2025:
2005:
1945:
1789:
1734:
1675:
1627:
1568:
1535:
1484:
1421:
1376:
1313:
1277:
1241:
1208:
1161:
995:
935:Legendre polynomials
902:
800:
754:
542:
386:
230:
173:Carl Friedrich Gauss
18932:Bayesian quadrature
18878:Gaussian quadrature
18744:Gaussian Quadrature
18646:1996JCoPh.129..406Y
18558:1991CoPhC..66..271S
18436:1971JEnMa...5....1P
18402:2001JCoAM.127..201L
18122:English Translation
17995:
17917:
17824:
16045:
15804:
15734:Gauss–Lobatto rules
15709:Gauss–Kronrod rules
15696:Gauss–Kronrod rules
15425:
15174:
15100:
15016:
14552:
14543:
14261:tridiagonal matrix
11350:
10924:Recurrence relation
10562:
9690:
9602:
9459:
9389:{\displaystyle n-2}
9134:
9007:
8926:{\displaystyle n-2}
8582:
8475:
8221:
8206:
8117:
7839:
7805:
7720:
7040:
6878:
6790:
6350:
6285:
5998:
5672:
5622:
5303:
5163:
5113:
4799:
4699:
4649:
4403:, then there exist
4396:to be the zeros of
4242:
4206:Fundamental theorem
4190:Hermite polynomials
3637:
2894:
2667:
2571:
2508:
2320:
2074:
1846:
1806:
1066:
1042:
730:
710:
602:
562:
250:
18947:Lebedev quadrature
18853:Simpson's 3/8 rule
18789:2021-04-14 at the
18727:Weisstein, Eric W.
18601:Olver, Frank W. J.
18444:10.1007/BF01535429
18417:MathWorks (2012).
18033:
17981:
17964:
17903:
17861:
17810:
17750:Abramowitz, Milton
17655:, pp. 172–175
17567:
17539:
17513:
17470:
17408:
17365:
17303:
17275:
17254:
17231:
17203:
17177:
17137:
17078:
17038:
16981:
16958:
16930:
16904:
16876:
16841:
16813:
16792:
16769:
16741:
16715:
16687:
16654:
16631:
16603:
16577:
16549:
16528:
16493:Number of points,
16475:
16229:
16091:
16055:
16027:
16010:
15969:
15787:
15740:Lobatto quadrature
15672:
15472:
15411:
15345:
15160:
15122:
15086:
15056:
15027:
14996:
14944:
14907:
14876:
14867:
14577:
14575:
14275:
14236:
14207:
14198:
13973:tridiagonal matrix
13957:
13943:
13882:
13862:
13833:
13819:
13704:
13608:
13581:
13545:
13415:
13131:
13084:
12958:
12805:
12766:
12669:
12635:
12587:
12529:
12318:
12291:
12045:
12012:
11951:
11686:
11654:
11627:
11600:
11558:
11445:
11393:
11336:
11281:
11094:
11055:
11023:
10997:
10945:
10880:
10805:
10772:
10743:
10712:
10548:
10523:
10481:
10445:
10416:
10387:
10335:
10292:
10214:
10212:
10211:
10141:
10104:
10043:
9994:
9951:
9820:
9778:
9742:
9707:
9678:
9653:
9588:
9445:
9428:
9386:
9358:
9206:We can then write
9194:
9120:
8993:
8976:
8923:
8886:
8775:
8715:
8656:
8568:
8461:
8433:
8390:
8311:
8283:
8207:
8192:
8103:
8084:
7901:
7825:
7793:
7750:
7708:
7668:
7666:
7665:
7609:
7578:
7542:
7515:
7486:
7372:
7318:
7316:
7315:
7256:
7214:
7178:
7149:
7095:
7093:
7092:
7026:
7006:are thus given by
6987:
6933:
6931:
6930:
6864:
6776:
6721:
6673:
6671:
6670:
6559:
6508:
6472:
6445:
6402:
6338:
6271:
6168:
6145:
6090:
6033:
5984:
5957:
5931:
5816:
5658:
5608:
5587:
5401:
5289:
5149:
5099:
5080:
4983:
4932:
4835:
4785:
4685:
4635:
4553:
4340:
4228:
4177:
4119:
4042:
3992:
3933:
3876:Jacobi polynomials
3863:
3670:
3623:
3589:
3507:
3470:
3468:
3355:
3353:
3239:
3237:
2877:
2850:
2816:
2782:
2745:
2707:
2705:
2650:
2554:
2494:
2471:
2432:
2388:
2306:
2283:
2249:
2211:
2060:
2043:
2011:
1989:
1929:
1829:
1792:
1781:Change of interval
1762:
1716:
1655:
1609:
1548:
1509:
1466:
1401:
1358:
1290:
1259:
1221:
1179:
1118:Number of points,
1100:
1054:
1028:
923:
885:
823:
782:
738:
718:
698:
585:
545:
493:
326:
233:
224:Carl Gustav Jacobi
161:numerical analysis
157:
18960:
18959:
18616:978-0-521-19225-5
18588:978-0-387-95452-3
18477:Quarteroni, Alfio
18469:978-0-521-88068-8
18331:978-0-13-627258-8
18208:978-0-89871-634-4
18189:978-1-611976-34-2
17769:978-0-486-61272-0
17754:Stegun, Irene Ann
17734:Eaton et al. 2018
17667:, Thm 3.6.24
17580:
17579:
17565:
17511:
17505:
17468:
17466:
17465:
17454:
17441:
17406:
17400:
17363:
17361:
17360:
17349:
17336:
17301:
17278:{\displaystyle 0}
17257:{\displaystyle 7}
17229:
17175:
17169:
17135:
17133:
17127:
17111:
17076:
17070:
17036:
17034:
17028:
17012:
16984:{\displaystyle 6}
16956:
16902:
16874:
16873:
16839:
16816:{\displaystyle 0}
16795:{\displaystyle 5}
16767:
16713:
16685:
16684:
16657:{\displaystyle 4}
16629:
16575:
16552:{\displaystyle 0}
16531:{\displaystyle 3}
16417:
16204:
15854:
15620:
15313:
14845:
14806:
14734:
14701:
14673:
14640:
14526:
14499:
13971:is the following
13885:{\displaystyle n}
13730:
13667:
13649:
13413:
13332:
13228:
11556:
10290:
10172:
10068:
9586:
9513:
9356:
9301:
9118:
9061:
8944:is orthogonal to
8840:
8773:
8648:
8566:
8529:
8388:
8357:
8275:
8181:
8082:
8047:
8022:
7991:
7957:
7925:
7924:
7893:
7823:
7748:
7626:
7484:
7418:
7351:
7276:
7141:
7053:
6979:
6891:
6719:
6631:
6426:
6425:
6400:
6266:
6178:must be at least
6136:
6024:
5922:
4978:
4917:
4305:
4203:
4202:
3990:
3931:
3930:
3565:
3432:
3317:
2621:
2597:
2552:
2359:
2201:
2170:
2113:
2014:{\displaystyle n}
1987:
1966:
1921:
1895:
1871:
1778:
1777:
1760:
1754:
1714:
1712:
1711:
1689:
1653:
1647:
1607:
1605:
1604:
1582:
1546:
1507:
1501:
1464:
1462:
1461:
1450:
1437:
1399:
1393:
1356:
1354:
1353:
1342:
1329:
1288:
1257:
1256:
1219:
1177:
1176:
1095:
821:
780:
779:
222:was developed by
86:
85:
78:
21:Gaussian integral
18980:
18942:Filon quadrature
18868:Romberg's method
18843:Trapezoidal rule
18820:
18813:
18806:
18797:
18796:
18740:
18739:
18683:
18657:
18627:
18591:
18577:(3rd ed.),
18569:
18552:(2–3): 271–275.
18540:
18530:
18508:
18472:
18447:
18422:
18413:
18396:(1–2): 201–217.
18384:
18383:
18374:(1–2): 165–180,
18362:
18360:
18335:
18319:
18308:Kahaner, David;
18304:
18283:
18277:
18269:
18254:Jacobi, C. G. J.
18249:
18239:
18230:(106): 221–230.
18211:
18193:
18174:
18149:
18119:
18105:
18080:
18078:
18076:
18064:
18042:
18040:
18039:
18034:
18005:
18004:
17994:
17989:
17973:
17971:
17970:
17965:
17945:
17944:
17916:
17911:
17892:
17890:
17870:
17868:
17867:
17862:
17823:
17818:
17797:
17736:
17731:
17725:
17720:
17714:
17709:
17703:
17698:
17692:
17686:
17680:
17674:
17668:
17662:
17656:
17650:
17644:
17638:
17632:
17627:
17621:
17616:
17597:
17593:
17576:
17574:
17573:
17568:
17566:
17558:
17548:
17546:
17545:
17540:
17522:
17520:
17519:
17514:
17512:
17507:
17506:
17501:
17489:
17479:
17477:
17476:
17471:
17469:
17467:
17458:
17457:
17455:
17447:
17442:
17434:
17432:
17417:
17415:
17414:
17409:
17407:
17402:
17401:
17396:
17384:
17374:
17372:
17371:
17366:
17364:
17362:
17353:
17352:
17350:
17342:
17337:
17329:
17327:
17312:
17310:
17309:
17304:
17302:
17294:
17284:
17282:
17281:
17276:
17263:
17261:
17260:
17255:
17240:
17238:
17237:
17232:
17230:
17222:
17212:
17210:
17209:
17204:
17186:
17184:
17183:
17178:
17176:
17171:
17170:
17165:
17156:
17146:
17144:
17143:
17138:
17136:
17134:
17129:
17128:
17123:
17117:
17112:
17104:
17102:
17087:
17085:
17084:
17079:
17077:
17072:
17071:
17066:
17057:
17047:
17045:
17044:
17039:
17037:
17035:
17030:
17029:
17024:
17018:
17013:
17005:
17003:
16990:
16988:
16987:
16982:
16967:
16965:
16964:
16959:
16957:
16949:
16939:
16937:
16936:
16931:
16913:
16911:
16910:
16905:
16903:
16895:
16885:
16883:
16882:
16877:
16875:
16866:
16865:
16850:
16848:
16847:
16842:
16840:
16832:
16822:
16820:
16819:
16814:
16801:
16799:
16798:
16793:
16778:
16776:
16775:
16770:
16768:
16760:
16750:
16748:
16747:
16742:
16724:
16722:
16721:
16716:
16714:
16706:
16696:
16694:
16693:
16688:
16686:
16677:
16676:
16663:
16661:
16660:
16655:
16640:
16638:
16637:
16632:
16630:
16622:
16612:
16610:
16609:
16604:
16586:
16584:
16583:
16578:
16576:
16568:
16558:
16556:
16555:
16550:
16537:
16535:
16534:
16529:
16514:
16505:
16490:
16489:
16484:
16482:
16481:
16476:
16443:
16442:
16418:
16416:
16415:
16414:
16409:
16405:
16401:
16397:
16354:
16353:
16352:
16347:
16343:
16339:
16335:
16314:
16313:
16295:
16294:
16289:
16285:
16262:
16257:
16256:
16238:
16236:
16235:
16230:
16219:
16218:
16205:
16203:
16202:
16201:
16196:
16192:
16191:
16187:
16186:
16173:
16172:
16129:
16124:
16123:
16104:
16100:
16098:
16097:
16092:
16081:
16080:
16064:
16062:
16061:
16056:
16041:
16019:
16017:
16016:
16011:
15987:
15978:
15976:
15975:
15970:
15965:
15964:
15952:
15948:
15947:
15932:
15931:
15920:
15909:
15855:
15853:
15830:
15825:
15803:
15798:
15779:
15764:
15760:
15729:
15721:
15717:
15690:
15681:
15679:
15678:
15673:
15640:
15639:
15621:
15619:
15618:
15617:
15612:
15608:
15604:
15600:
15563:
15562:
15561:
15556:
15552:
15539:
15538:
15524:
15520:
15503:
15494:
15481:
15479:
15478:
15473:
15424:
15419:
15385:
15381:
15377:
15370:
15358:
15354:
15352:
15351:
15346:
15341:
15340:
15328:
15327:
15314:
15312:
15295:
15285:
15284:
15265:
15257:
15256:
15240:
15239:
15229:
15224:
15173:
15168:
15152:
15131:
15129:
15128:
15123:
15099:
15094:
15082:
15081:
15065:
15063:
15062:
15057:
15055:
15054:
15036:
15034:
15033:
15028:
15026:
15025:
15020:
15015:
15004:
14990:
14989:
14977:
14976:
14960:
14953:
14951:
14950:
14945:
14943:
14942:
14919:similar matrices
14916:
14914:
14913:
14908:
14906:
14905:
14892:
14885:
14883:
14882:
14877:
14872:
14871:
14864:
14863:
14846:
14844:
14843:
14828:
14807:
14805:
14804:
14789:
14785:
14784:
14735:
14733:
14732:
14723:
14702:
14700:
14699:
14690:
14686:
14685:
14674:
14672:
14671:
14662:
14641:
14639:
14638:
14629:
14625:
14624:
14604:
14603:
14586:
14584:
14583:
14578:
14576:
14554:
14553:
14527:
14525:
14524:
14509:
14500:
14498:
14497:
14476:
14475:
14454:
14449:
14448:
14431:
14430:
14420:
14419:
14402:
14401:
14357:
14356:
14334:
14333:
14315:
14314:
14303:
14302:
14284:
14282:
14281:
14276:
14274:
14273:
14249:
14245:
14243:
14242:
14237:
14235:
14234:
14216:
14214:
14213:
14208:
14203:
14202:
14195:
14194:
14177:
14176:
14137:
14136:
14087:
14086:
14053:
14052:
14041:
14040:
14007:
14006:
13986:
13970:
13966:
13964:
13963:
13958:
13956:
13955:
13954:
13948:
13947:
13910:
13909:
13904:
13891:
13889:
13888:
13883:
13871:
13869:
13868:
13863:
13861:
13860:
13855:
13842:
13840:
13839:
13834:
13832:
13831:
13830:
13824:
13823:
13807:
13806:
13775:
13774:
13754:
13753:
13732:
13731:
13723:
13713:
13711:
13710:
13705:
13703:
13702:
13697:
13682:
13681:
13669:
13668:
13660:
13651:
13650:
13642:
13621:
13617:
13615:
13614:
13609:
13607:
13606:
13590:
13588:
13587:
13582:
13580:
13579:
13554:
13552:
13551:
13546:
13541:
13540:
13528:
13527:
13509:
13508:
13487:
13486:
13468:
13467:
13449:
13448:
13424:
13422:
13421:
13416:
13414:
13412:
13408:
13407:
13389:
13388:
13369:
13365:
13364:
13352:
13351:
13338:
13333:
13331:
13327:
13326:
13308:
13307:
13288:
13284:
13283:
13265:
13264:
13248:
13243:
13242:
13229:
13227:
13223:
13222:
13210:
13209:
13196:
13192:
13191:
13179:
13178:
13162:
13157:
13156:
13140:
13138:
13137:
13132:
13115:
13114:
13093:
13091:
13090:
13085:
13074:
13073:
13058:
13057:
13036:
13035:
13023:
13022:
12992:
12991:
12967:
12965:
12964:
12959:
12948:
12947:
12932:
12931:
12898:
12897:
12885:
12884:
12848:
12847:
12825:
12814:
12812:
12811:
12806:
12798:
12797:
12775:
12773:
12772:
12767:
12756:
12755:
12740:
12739:
12721:
12720:
12708:
12707:
12685:
12678:
12676:
12675:
12670:
12668:
12667:
12651:
12644:
12642:
12641:
12636:
12634:
12633:
12596:
12594:
12593:
12588:
12538:
12536:
12535:
12530:
12519:
12518:
12506:
12505:
12484:
12483:
12471:
12470:
12449:
12448:
12436:
12435:
12423:
12422:
12401:
12400:
12388:
12387:
12366:
12365:
12353:
12352:
12327:
12325:
12324:
12319:
12317:
12316:
12300:
12298:
12297:
12292:
12287:
12286:
12274:
12273:
12261:
12260:
12236:
12235:
12223:
12222:
12204:
12203:
12176:
12175:
12163:
12162:
12150:
12149:
12128:
12127:
12115:
12114:
12093:
12092:
12080:
12079:
12054:
12052:
12051:
12046:
12044:
12043:
12021:
12019:
12018:
12013:
12011:
12010:
11992:
11991:
11979:
11978:
11960:
11958:
11957:
11952:
11941:
11940:
11928:
11927:
11906:
11905:
11893:
11892:
11871:
11870:
11858:
11857:
11845:
11844:
11823:
11822:
11810:
11809:
11788:
11787:
11775:
11774:
11759:
11758:
11728:
11727:
11715:
11714:
11695:
11693:
11692:
11687:
11663:
11661:
11660:
11655:
11653:
11652:
11636:
11634:
11633:
11628:
11626:
11625:
11609:
11607:
11606:
11601:
11584:
11583:
11567:
11565:
11564:
11559:
11557:
11555:
11551:
11550:
11549:
11537:
11536:
11522:
11518:
11517:
11516:
11504:
11503:
11486:
11481:
11480:
11458:
11454:
11452:
11451:
11446:
11402:
11400:
11399:
11394:
11349:
11344:
11290:
11288:
11287:
11282:
11271:
11270:
11261:
11260:
11230:
11229:
11214:
11213:
11180:
11179:
11167:
11166:
11130:
11129:
11103:
11101:
11100:
11095:
11084:
11083:
11064:
11062:
11061:
11056:
11032:
11030:
11029:
11024:
11006:
11004:
11003:
10998:
10987:
10986:
10974:
10973:
10954:
10952:
10951:
10946:
10944:
10943:
10919:
10904:
10889:
10887:
10886:
10881:
10864:
10863:
10847:
10836:
10829:
10814:
10812:
10811:
10806:
10798:
10797:
10781:
10779:
10778:
10773:
10752:
10750:
10749:
10744:
10721:
10719:
10718:
10713:
10705:
10704:
10692:
10691:
10682:
10681:
10668:
10663:
10642:
10641:
10626:
10625:
10615:
10610:
10561:
10556:
10540:
10536:
10532:
10530:
10529:
10524:
10513:
10512:
10490:
10488:
10487:
10482:
10471:
10470:
10454:
10452:
10451:
10446:
10425:
10423:
10422:
10417:
10396:
10394:
10393:
10388:
10386:
10385:
10367:
10366:
10344:
10342:
10341:
10336:
10325:
10324:
10308:
10301:
10299:
10298:
10293:
10291:
10289:
10288:
10283:
10279:
10278:
10277:
10265:
10264:
10249:
10248:
10243:
10239:
10238:
10237:
10216:
10213:
10150:
10148:
10147:
10142:
10113:
10111:
10110:
10105:
10103:
10099:
10098:
10085:
10084:
10069:
10061:
10052:
10050:
10049:
10044:
10039:
10038:
10026:
10025:
10003:
10001:
10000:
9995:
9990:
9989:
9977:
9976:
9960:
9958:
9957:
9952:
9947:
9946:
9934:
9933:
9903:
9902:
9890:
9889:
9865:
9864:
9852:
9851:
9829:
9827:
9826:
9821:
9810:
9809:
9787:
9785:
9784:
9779:
9768:
9767:
9751:
9749:
9748:
9743:
9741:
9740:
9716:
9714:
9713:
9708:
9703:
9702:
9686:
9662:
9660:
9659:
9654:
9646:
9645:
9630:
9629:
9601:
9596:
9587:
9585:
9581:
9580:
9568:
9567:
9552:
9551:
9535:
9534:
9525:
9514:
9512:
9511:
9510:
9494:
9484:
9483:
9473:
9458:
9453:
9437:
9435:
9434:
9429:
9418:
9417:
9395:
9393:
9392:
9387:
9367:
9365:
9364:
9359:
9357:
9355:
9354:
9339:
9329:
9328:
9312:
9307:
9303:
9302:
9300:
9299:
9284:
9274:
9273:
9257:
9252:
9251:
9228:
9227:
9203:
9201:
9200:
9195:
9187:
9186:
9162:
9161:
9133:
9128:
9119:
9117:
9113:
9112:
9100:
9099:
9083:
9082:
9073:
9062:
9060:
9059:
9058:
9042:
9032:
9031:
9021:
9006:
9001:
8985:
8983:
8982:
8977:
8966:
8965:
8943:
8932:
8930:
8929:
8924:
8906:
8895:
8893:
8892:
8887:
8870:
8869:
8854:
8853:
8841:
8839:
8838:
8837:
8821:
8811:
8810:
8800:
8791:
8784:
8782:
8781:
8776:
8774:
8772:
8771:
8770:
8754:
8744:
8743:
8733:
8724:
8722:
8721:
8716:
8705:
8704:
8665:
8663:
8662:
8657:
8649:
8647:
8646:
8645:
8629:
8619:
8618:
8596:
8581:
8576:
8567:
8565:
8561:
8560:
8541:
8530:
8528:
8527:
8526:
8510:
8500:
8499:
8489:
8474:
8469:
8453:
8442:
8440:
8439:
8434:
8423:
8422:
8406:
8399:
8397:
8396:
8391:
8389:
8387:
8386:
8385:
8369:
8368:
8367:
8362:
8358:
8356:
8355:
8343:
8329:
8320:
8318:
8317:
8312:
8292:
8290:
8289:
8284:
8276:
8274:
8273:
8272:
8256:
8246:
8245:
8235:
8220:
8215:
8205:
8200:
8182:
8180:
8179:
8178:
8162:
8152:
8151:
8142:
8141:
8131:
8116:
8111:
8093:
8091:
8090:
8085:
8083:
8081:
8080:
8079:
8060:
8058:
8057:
8052:
8048:
8046:
8045:
8033:
8023:
8021:
8020:
8019:
8003:
8002:
8001:
7996:
7992:
7990:
7989:
7977:
7963:
7958:
7956:
7955:
7954:
7935:
7919:
7910:
7908:
7907:
7902:
7894:
7892:
7891:
7890:
7874:
7864:
7863:
7853:
7838:
7833:
7824:
7822:
7818:
7817:
7801:
7788:
7783:
7782:
7765:
7759:
7757:
7756:
7751:
7749:
7747:
7746:
7737:
7733:
7732:
7716:
7706:
7701:
7697:
7696:
7695:
7683:
7682:
7667:
7618:
7616:
7615:
7610:
7608:
7607:
7591:
7587:
7585:
7584:
7579:
7568:
7567:
7551:
7549:
7548:
7543:
7541:
7540:
7524:
7522:
7521:
7516:
7514:
7513:
7495:
7493:
7492:
7487:
7485:
7483:
7482:
7478:
7477:
7476:
7456:
7455:
7445:
7435:
7434:
7424:
7419:
7417:
7416:
7415:
7399:
7398:
7394:
7393:
7392:
7371:
7349:
7344:
7340:
7339:
7338:
7317:
7265:
7263:
7262:
7257:
7246:
7245:
7223:
7221:
7220:
7215:
7204:
7203:
7187:
7185:
7184:
7179:
7177:
7176:
7158:
7156:
7155:
7150:
7142:
7140:
7139:
7138:
7126:
7125:
7115:
7114:
7113:
7097:
7094:
7039:
7034:
7022:
7021:
7005:
6996:
6994:
6993:
6988:
6980:
6978:
6977:
6976:
6964:
6963:
6953:
6952:
6951:
6935:
6932:
6877:
6872:
6860:
6859:
6843:
6838:
6789:
6784:
6768:
6764:
6760:
6749:
6745:
6741:
6730:
6728:
6727:
6722:
6720:
6718:
6717:
6716:
6704:
6703:
6693:
6692:
6691:
6675:
6672:
6627:
6626:
6610:
6605:
6568:
6566:
6565:
6560:
6555:
6554:
6532:
6522:one can express
6517:
6515:
6514:
6509:
6498:
6497:
6481:
6479:
6478:
6473:
6471:
6470:
6454:
6452:
6451:
6446:
6444:
6443:
6420:
6411:
6409:
6408:
6403:
6401:
6399:
6395:
6394:
6382:
6381:
6363:
6362:
6346:
6336:
6329:
6328:
6313:
6312:
6284:
6279:
6269:
6267:
6265:
6264:
6249:
6248:
6239:
6234:
6233:
6216:
6204:
6200:
6196:
6192:
6181:
6177:
6175:
6174:
6169:
6164:
6163:
6144:
6128:
6121:
6110:
6099:
6097:
6096:
6091:
6060:
6056:
6052:
6051:
6032:
6008:
6007:
5997:
5992:
5976:
5972:
5966:
5964:
5963:
5958:
5950:
5949:
5930:
5911:
5910:
5894:
5883:
5879:
5875:
5868:
5864:
5860:
5846:
5838:
5825:
5823:
5822:
5817:
5809:
5808:
5792:
5791:
5781:
5776:
5755:
5754:
5738:
5737:
5727:
5722:
5671:
5666:
5621:
5616:
5600:
5596:
5594:
5593:
5588:
5580:
5579:
5558:
5557:
5536:
5535:
5516:
5515:
5503:
5502:
5487:
5486:
5464:
5457:
5450:
5434:
5423:
5410:
5408:
5407:
5402:
5397:
5396:
5383:
5382:
5365:
5360:
5326:
5325:
5302:
5297:
5284:
5283:
5266:
5261:
5233:
5232:
5207:
5206:
5196:
5191:
5162:
5157:
5112:
5107:
5089:
5087:
5086:
5081:
5073:
5072:
5047:
5046:
5036:
5031:
4992:
4990:
4989:
4984:
4979:
4977:
4976:
4975:
4963:
4962:
4952:
4951:
4950:
4934:
4931:
4904:
4903:
4887:
4868:
4864:
4857:
4844:
4842:
4841:
4836:
4798:
4793:
4774:
4773:
4730:
4729:
4719:
4718:
4698:
4693:
4648:
4643:
4627:
4616:
4612:
4605:
4598:
4587:
4580:
4573:
4562:
4560:
4559:
4554:
4512:
4511:
4480:
4473:
4465:
4451:
4439:
4432:
4424:
4413:
4406:
4402:
4395:
4388:
4381:
4375:
4365:
4358:
4349:
4347:
4346:
4341:
4306:
4303:
4275:
4274:
4265:
4264:
4241:
4236:
4220:
4216:
4186:
4184:
4183:
4178:
4176:
4175:
4174:
4173:
4148:
4128:
4126:
4125:
4120:
4102:
4101:
4089:
4088:
4071:
4051:
4049:
4048:
4043:
4040:
4039:
4019:
4001:
3999:
3998:
3993:
3991:
3989:
3988:
3973:
3963:
3942:
3940:
3939:
3934:
3932:
3929:
3928:
3913:
3909:
3899:
3886:
3872:
3870:
3869:
3864:
3840:
3839:
3834:
3830:
3814:
3813:
3808:
3804:
3780:
3760:
3755:
3739:
3724:
3723:
3716:
3705:
3698:
3691:
3687:
3683:
3679:
3677:
3676:
3671:
3636:
3631:
3615:
3598:
3596:
3595:
3590:
3570:
3566:
3561:
3550:
3541:
3537:
3536:
3516:
3514:
3513:
3508:
3506:
3502:
3501:
3479:
3477:
3476:
3471:
3469:
3456:
3437:
3433:
3431:
3408:
3364:
3362:
3361:
3356:
3354:
3341:
3322:
3318:
3316:
3293:
3248:
3246:
3245:
3240:
3238:
3225:
3221:
3217:
3181:
3177:
3173:
3131:
3127:
3123:
3122:
3118:
3088:
3084:
3040:
3036:
3032:
3031:
3027:
3020:
3019:
2999:
2998:
2986:
2982:
2975:
2974:
2954:
2953:
2929:
2928:
2920:
2916:
2893:
2888:
2859:
2857:
2856:
2851:
2843:
2842:
2825:
2823:
2822:
2817:
2809:
2808:
2791:
2789:
2788:
2783:
2772:
2771:
2754:
2752:
2751:
2746:
2738:
2737:
2716:
2714:
2713:
2708:
2706:
2702:
2701:
2693:
2689:
2666:
2661:
2640:
2636:
2635:
2627:
2623:
2622:
2617:
2606:
2598:
2593:
2582:
2570:
2565:
2553:
2548:
2537:
2528:
2507:
2502:
2480:
2478:
2477:
2472:
2470:
2466:
2441:
2439:
2438:
2433:
2431:
2427:
2397:
2395:
2394:
2389:
2387:
2386:
2378:
2374:
2364:
2360:
2358:
2341:
2319:
2314:
2292:
2290:
2289:
2284:
2279:
2258:
2256:
2255:
2250:
2248:
2220:
2218:
2217:
2212:
2207:
2203:
2202:
2197:
2186:
2181:
2180:
2171:
2166:
2155:
2145:
2144:
2134:
2129:
2114:
2109:
2098:
2073:
2068:
2052:
2050:
2049:
2044:
2020:
2018:
2017:
2012:
1998:
1996:
1995:
1990:
1988:
1983:
1972:
1967:
1965:
1957:
1949:
1938:
1936:
1935:
1930:
1922:
1920:
1912:
1904:
1901:
1897:
1896:
1891:
1880:
1872:
1867:
1856:
1845:
1840:
1805:
1800:
1771:
1769:
1768:
1763:
1761:
1756:
1755:
1750:
1738:
1725:
1723:
1722:
1717:
1715:
1713:
1704:
1703:
1692:
1690:
1682:
1664:
1662:
1661:
1656:
1654:
1649:
1648:
1643:
1631:
1618:
1616:
1615:
1610:
1608:
1606:
1597:
1596:
1585:
1583:
1575:
1557:
1555:
1554:
1549:
1547:
1539:
1518:
1516:
1515:
1510:
1508:
1503:
1502:
1497:
1488:
1475:
1473:
1472:
1467:
1465:
1463:
1454:
1453:
1451:
1443:
1438:
1430:
1428:
1410:
1408:
1407:
1402:
1400:
1395:
1394:
1389:
1380:
1367:
1365:
1364:
1359:
1357:
1355:
1346:
1345:
1343:
1335:
1330:
1322:
1320:
1299:
1297:
1296:
1291:
1289:
1281:
1268:
1266:
1265:
1260:
1258:
1249:
1248:
1230:
1228:
1227:
1222:
1220:
1212:
1188:
1186:
1185:
1180:
1178:
1172:
1168:
1139:
1130:
1121:
1115:
1114:
1109:
1107:
1106:
1101:
1096:
1094:
1093:
1092:
1087:
1083:
1079:
1078:
1062:
1047:
1043:
1041:
1036:
1012:
1007:
1006:
990:
983:
979:
973:-th Gauss node,
972:
968:
956:
952:
932:
930:
929:
926:{\displaystyle }
924:
897:
883:
850:
832:
830:
829:
824:
822:
820:
819:
804:
791:
789:
788:
783:
781:
778:
777:
762:
758:
747:
745:
744:
739:
734:
726:
706:
696:
691:
654:
653:
648:
644:
628:
627:
622:
618:
601:
596:
561:
556:
531:
523:
515:
502:
500:
499:
494:
455:
454:
449:
445:
429:
428:
423:
419:
368:
360:
345:
335:
333:
332:
327:
319:
318:
303:
302:
292:
287:
249:
244:
214:
204:
197:
190:
166:
154:
125:
124:
123:
119:
113:
101:trapezoidal rule
96:trapezoidal rule
81:
74:
70:
67:
61:
56:this article by
47:inline citations
34:
33:
26:
18988:
18987:
18983:
18982:
18981:
18979:
18978:
18977:
18963:
18962:
18961:
18956:
18915:
18872:
18829:
18824:
18791:Wayback Machine
18668:
18665:
18660:
18634:J. Comput. Phys
18617:
18589:
18505:
18487:Springer-Verlag
18470:
18332:
18271:
18270:
18209:
18190:
18074:
18072:
18000:
17996:
17990:
17985:
17979:
17976:
17975:
17940:
17936:
17912:
17907:
17901:
17898:
17897:
17881:(91): 477–481.
17819:
17814:
17808:
17805:
17804:
17770:
17744:
17739:
17732:
17728:
17721:
17717:
17710:
17706:
17699:
17695:
17687:
17683:
17675:
17671:
17663:
17659:
17651:
17647:
17639:
17635:
17628:
17624:
17617:
17613:
17609:
17604:
17595:
17591:
17557:
17555:
17552:
17551:
17531:
17528:
17527:
17500:
17490:
17488:
17486:
17483:
17482:
17456:
17446:
17433:
17431:
17426:
17423:
17422:
17395:
17385:
17383:
17381:
17378:
17377:
17351:
17341:
17328:
17326:
17321:
17318:
17317:
17293:
17291:
17288:
17287:
17270:
17267:
17266:
17249:
17246:
17245:
17221:
17219:
17216:
17215:
17195:
17192:
17191:
17164:
17157:
17155:
17153:
17150:
17149:
17122:
17118:
17116:
17103:
17101:
17096:
17093:
17092:
17065:
17058:
17056:
17054:
17051:
17050:
17023:
17019:
17017:
17004:
17002:
16997:
16994:
16993:
16976:
16973:
16972:
16948:
16946:
16943:
16942:
16922:
16919:
16918:
16894:
16892:
16889:
16888:
16864:
16859:
16856:
16855:
16831:
16829:
16826:
16825:
16808:
16805:
16804:
16787:
16784:
16783:
16759:
16757:
16754:
16753:
16733:
16730:
16729:
16705:
16703:
16700:
16699:
16675:
16670:
16667:
16666:
16649:
16646:
16645:
16621:
16619:
16616:
16615:
16595:
16592:
16591:
16567:
16565:
16562:
16561:
16544:
16541:
16540:
16523:
16520:
16519:
16513:
16509:
16504:
16500:
16423:
16419:
16410:
16384:
16380:
16379:
16375:
16374:
16355:
16348:
16325:
16321:
16320:
16316:
16315:
16300:
16296:
16290:
16275:
16271:
16270:
16263:
16261:
16252:
16248:
16246:
16243:
16242:
16214:
16210:
16197:
16182:
16178:
16174:
16162:
16158:
16157:
16153:
16152:
16133:
16128:
16119:
16115:
16113:
16110:
16109:
16102:
16076:
16072:
16070:
16067:
16066:
16031:
16025:
16022:
16021:
15993:
15990:
15989:
15986:
15982:
15960:
15956:
15943:
15939:
15927:
15923:
15922:
15910:
15899:
15834:
15829:
15805:
15799:
15791:
15785:
15782:
15781:
15770:
15762:
15754:
15736:
15723:
15719:
15712:
15704:
15698:
15692:can be useful.
15685:
15626:
15622:
15613:
15593:
15589:
15588:
15584:
15583:
15564:
15557:
15545:
15541:
15540:
15525:
15510:
15506:
15505:
15504:
15502:
15500:
15497:
15496:
15485:
15420:
15415:
15391:
15388:
15387:
15383:
15379:
15376:
15372:
15360:
15356:
15336:
15332:
15323:
15319:
15296:
15271:
15267:
15266:
15264:
15252:
15248:
15235:
15231:
15225:
15214:
15169:
15164:
15158:
15155:
15154:
15147:
15144:
15142:Error estimates
15095:
15090:
15077:
15073:
15071:
15068:
15067:
15050:
15046:
15044:
15041:
15040:
15021:
15005:
15000:
14992:
14991:
14985:
14981:
14972:
14968:
14966:
14963:
14962:
14959:
14955:
14932:
14928:
14926:
14923:
14922:
14901:
14900:
14898:
14895:
14894:
14888:
14866:
14865:
14853:
14849:
14847:
14833:
14829:
14827:
14825:
14820:
14815:
14809:
14808:
14794:
14790:
14788:
14786:
14774:
14770:
14768:
14763:
14758:
14752:
14751:
14746:
14741:
14736:
14728:
14724:
14722:
14720:
14714:
14713:
14708:
14703:
14695:
14691:
14689:
14687:
14681:
14677:
14675:
14667:
14663:
14661:
14658:
14657:
14652:
14647:
14642:
14634:
14630:
14628:
14626:
14620:
14616:
14609:
14608:
14599:
14598:
14596:
14593:
14592:
14574:
14573:
14541:
14540:
14533:
14528:
14514:
14510:
14508:
14501:
14481:
14477:
14459:
14455:
14453:
14432:
14426:
14425:
14424:
14403:
14397:
14396:
14395:
14392:
14391:
14363:
14358:
14346:
14342:
14335:
14323:
14319:
14304:
14298:
14297:
14296:
14292:
14290:
14287:
14286:
14269:
14268:
14266:
14263:
14262:
14247:
14230:
14226:
14224:
14221:
14220:
14197:
14196:
14184:
14180:
14178:
14166:
14162:
14160:
14155:
14150:
14144:
14143:
14138:
14126:
14122:
14120:
14115:
14110:
14104:
14103:
14098:
14093:
14088:
14082:
14078:
14076:
14070:
14069:
14064:
14059:
14054:
14048:
14044:
14042:
14036:
14032:
14029:
14028:
14023:
14018:
14013:
14008:
14002:
13998:
13991:
13990:
13982:
13980:
13977:
13976:
13968:
13950:
13949:
13942:
13941:
13936:
13931:
13926:
13916:
13915:
13914:
13905:
13900:
13899:
13897:
13894:
13893:
13877:
13874:
13873:
13856:
13851:
13850:
13848:
13845:
13844:
13826:
13825:
13818:
13817:
13796:
13792:
13790:
13785:
13770:
13766:
13764:
13749:
13745:
13738:
13737:
13736:
13722:
13721:
13719:
13716:
13715:
13698:
13693:
13692:
13677:
13673:
13659:
13658:
13641:
13640:
13635:
13632:
13631:
13628:
13619:
13602:
13598:
13596:
13593:
13592:
13569:
13565:
13560:
13557:
13556:
13536:
13532:
13523:
13519:
13498:
13494:
13482:
13478:
13457:
13453:
13444:
13440:
13432:
13429:
13428:
13397:
13393:
13378:
13374:
13370:
13360:
13356:
13347:
13343:
13339:
13337:
13316:
13312:
13297:
13293:
13289:
13273:
13269:
13260:
13256:
13249:
13247:
13238:
13234:
13218:
13214:
13205:
13201:
13197:
13187:
13183:
13174:
13170:
13163:
13161:
13152:
13148:
13146:
13143:
13142:
13107:
13103:
13101:
13098:
13097:
13063:
13059:
13053:
13049:
13031:
13027:
13018:
13014:
12981:
12977:
12975:
12972:
12971:
12937:
12933:
12915:
12911:
12893:
12889:
12874:
12870:
12837:
12833:
12831:
12828:
12827:
12816:
12787:
12783:
12781:
12778:
12777:
12751:
12747:
12735:
12731:
12716:
12712:
12703:
12699:
12691:
12688:
12687:
12680:
12663:
12659:
12657:
12654:
12653:
12646:
12629:
12625:
12602:
12599:
12598:
12546:
12543:
12542:
12514:
12510:
12501:
12497:
12479:
12475:
12466:
12462:
12444:
12440:
12431:
12427:
12412:
12408:
12396:
12392:
12383:
12379:
12361:
12357:
12342:
12338:
12333:
12330:
12329:
12312:
12308:
12306:
12303:
12302:
12282:
12278:
12269:
12265:
12250:
12246:
12231:
12227:
12212:
12208:
12187:
12183:
12171:
12167:
12158:
12154:
12139:
12135:
12123:
12119:
12110:
12106:
12088:
12084:
12069:
12065:
12060:
12057:
12056:
12033:
12029:
12027:
12024:
12023:
12006:
12002:
11987:
11983:
11974:
11970:
11968:
11965:
11964:
11936:
11932:
11923:
11919:
11901:
11897:
11888:
11884:
11866:
11862:
11853:
11849:
11834:
11830:
11818:
11814:
11805:
11801:
11783:
11779:
11770:
11766:
11748:
11744:
11723:
11719:
11710:
11706:
11701:
11698:
11697:
11669:
11666:
11665:
11648:
11644:
11642:
11639:
11638:
11621:
11617:
11615:
11612:
11611:
11579:
11575:
11573:
11570:
11569:
11545:
11541:
11532:
11528:
11527:
11523:
11512:
11508:
11499:
11495:
11491:
11487:
11485:
11470:
11466:
11464:
11461:
11460:
11456:
11410:
11407:
11406:
11345:
11340:
11298:
11295:
11294:
11266:
11262:
11250:
11246:
11219:
11215:
11197:
11193:
11175:
11171:
11156:
11152:
11119:
11115:
11113:
11110:
11109:
11079:
11075:
11070:
11067:
11066:
11038:
11035:
11034:
11012:
11009:
11008:
10982:
10978:
10969:
10965:
10960:
10957:
10956:
10939:
10935:
10933:
10930:
10929:
10926:
10910:
10895:
10859:
10855:
10853:
10850:
10849:
10838:
10835:
10831:
10828:
10824:
10821:
10793:
10789:
10787:
10784:
10783:
10758:
10755:
10754:
10729:
10726:
10725:
10700:
10696:
10687:
10683:
10674:
10670:
10664:
10653:
10637:
10633:
10621:
10617:
10611:
10600:
10557:
10552:
10546:
10543:
10542:
10538:
10534:
10508:
10504:
10496:
10493:
10492:
10491:applies. Since
10466:
10462:
10460:
10457:
10456:
10431:
10428:
10427:
10402:
10399:
10398:
10378:
10374:
10362:
10358:
10350:
10347:
10346:
10320:
10316:
10314:
10311:
10310:
10307:
10303:
10284:
10273:
10269:
10260:
10256:
10255:
10251:
10250:
10244:
10233:
10229:
10222:
10218:
10217:
10215:
10210:
10209:
10197:
10196:
10176:
10155:
10152:
10151:
10127:
10124:
10123:
10120:
10094:
10090:
10086:
10074:
10070:
10060:
10058:
10055:
10054:
10034:
10030:
10015:
10011:
10009:
10006:
10005:
9985:
9981:
9972:
9968:
9966:
9963:
9962:
9942:
9938:
9923:
9919:
9898:
9894:
9885:
9881:
9860:
9856:
9841:
9837:
9835:
9832:
9831:
9799:
9795:
9793:
9790:
9789:
9763:
9759:
9757:
9754:
9753:
9736:
9732:
9730:
9727:
9726:
9698:
9694:
9682:
9676:
9673:
9672:
9641:
9637:
9619:
9615:
9597:
9592:
9576:
9572:
9557:
9553:
9541:
9537:
9536:
9530:
9526:
9524:
9506:
9502:
9495:
9479:
9475:
9474:
9472:
9454:
9449:
9443:
9440:
9439:
9407:
9403:
9401:
9398:
9397:
9375:
9372:
9371:
9344:
9340:
9318:
9314:
9313:
9311:
9289:
9285:
9263:
9259:
9258:
9256:
9241:
9237:
9236:
9232:
9217:
9213:
9211:
9208:
9207:
9176:
9172:
9151:
9147:
9129:
9124:
9108:
9104:
9089:
9085:
9084:
9078:
9074:
9072:
9054:
9050:
9043:
9027:
9023:
9022:
9020:
9002:
8997:
8991:
8988:
8987:
8955:
8951:
8949:
8946:
8945:
8934:
8912:
8909:
8908:
8897:
8859:
8855:
8849:
8845:
8833:
8829:
8822:
8806:
8802:
8801:
8799:
8797:
8794:
8793:
8786:
8766:
8762:
8755:
8739:
8735:
8734:
8732:
8730:
8727:
8726:
8694:
8690:
8673:
8670:
8669:
8641:
8637:
8630:
8614:
8610:
8597:
8595:
8577:
8572:
8556:
8552:
8545:
8540:
8522:
8518:
8511:
8495:
8491:
8490:
8488:
8470:
8465:
8459:
8456:
8455:
8444:
8418:
8414:
8412:
8409:
8408:
8401:
8381:
8377:
8370:
8363:
8351:
8347:
8342:
8338:
8337:
8330:
8328:
8326:
8323:
8322:
8300:
8297:
8296:
8268:
8264:
8257:
8241:
8237:
8236:
8234:
8216:
8211:
8201:
8196:
8174:
8170:
8163:
8147:
8143:
8137:
8133:
8132:
8130:
8112:
8107:
8101:
8098:
8097:
8075:
8071:
8064:
8059:
8053:
8041:
8037:
8032:
8028:
8027:
8015:
8011:
8004:
7997:
7985:
7981:
7976:
7972:
7971:
7964:
7962:
7950:
7946:
7939:
7934:
7932:
7929:
7928:
7886:
7882:
7875:
7859:
7855:
7854:
7852:
7834:
7829:
7813:
7809:
7797:
7792:
7787:
7778:
7774:
7772:
7769:
7768:
7742:
7738:
7728:
7724:
7712:
7707:
7705:
7691:
7687:
7678:
7674:
7673:
7669:
7664:
7663:
7651:
7650:
7630:
7624:
7621:
7620:
7603:
7599:
7597:
7594:
7593:
7589:
7563:
7559:
7557:
7554:
7553:
7536:
7532:
7530:
7527:
7526:
7509:
7505:
7503:
7500:
7499:
7472:
7468:
7461:
7457:
7451:
7447:
7446:
7430:
7426:
7425:
7423:
7411:
7407:
7400:
7388:
7384:
7377:
7373:
7355:
7350:
7348:
7334:
7330:
7323:
7319:
7314:
7313:
7301:
7300:
7280:
7274:
7271:
7270:
7235:
7231:
7229:
7226:
7225:
7199:
7195:
7193:
7190:
7189:
7172:
7168:
7166:
7163:
7162:
7134:
7130:
7121:
7117:
7116:
7109:
7105:
7098:
7096:
7091:
7090:
7078:
7077:
7057:
7035:
7030:
7017:
7013:
7011:
7008:
7007:
7004:
7000:
6972:
6968:
6959:
6955:
6954:
6947:
6943:
6936:
6934:
6929:
6928:
6916:
6915:
6895:
6873:
6868:
6855:
6851:
6839:
6828:
6785:
6780:
6774:
6771:
6770:
6766:
6762:
6751:
6747:
6743:
6732:
6712:
6708:
6699:
6695:
6694:
6687:
6683:
6676:
6674:
6669:
6668:
6656:
6655:
6635:
6622:
6618:
6606:
6595:
6574:
6571:
6570:
6550:
6546:
6538:
6535:
6534:
6523:
6493:
6489:
6487:
6484:
6483:
6466:
6462:
6460:
6457:
6456:
6439:
6435:
6433:
6430:
6429:
6390:
6386:
6371:
6367:
6358:
6354:
6342:
6337:
6324:
6320:
6302:
6298:
6280:
6275:
6270:
6268:
6254:
6250:
6244:
6240:
6238:
6229:
6225:
6223:
6220:
6219:
6211:
6202:
6198:
6194:
6191:
6183:
6179:
6159:
6155:
6140:
6134:
6131:
6130:
6123:
6120:
6112:
6101:
6047:
6043:
6028:
6023:
6019:
6003:
5999:
5993:
5988:
5982:
5979:
5978:
5974:
5970:
5945:
5941:
5926:
5906:
5902:
5900:
5897:
5896:
5893:
5885:
5881:
5877:
5874:
5870:
5866:
5862:
5859:
5851:
5840:
5829:
5804:
5800:
5787:
5783:
5777:
5766:
5750:
5746:
5733:
5729:
5723:
5712:
5667:
5662:
5617:
5612:
5606:
5603:
5602:
5598:
5575:
5571:
5553:
5549:
5531:
5527:
5511:
5507:
5498:
5494:
5482:
5478:
5470:
5467:
5466:
5463:
5459:
5456:
5452:
5444:
5436:
5433:
5425:
5422:
5414:
5392:
5388:
5378:
5374:
5361:
5350:
5321:
5317:
5298:
5293:
5279:
5275:
5262:
5251:
5228:
5224:
5202:
5198:
5192:
5181:
5158:
5153:
5108:
5103:
5097:
5094:
5093:
5068:
5064:
5042:
5038:
5032:
5021:
5000:
4997:
4996:
4971:
4967:
4958:
4954:
4953:
4946:
4942:
4935:
4933:
4921:
4899:
4895:
4893:
4890:
4889:
4881:
4873:
4866:
4859:
4848:
4794:
4789:
4769:
4768:
4725:
4721:
4714:
4713:
4694:
4689:
4644:
4639:
4633:
4630:
4629:
4618:
4614:
4611:
4607:
4600:
4589:
4586:
4582:
4575:
4564:
4507:
4503:
4486:
4483:
4482:
4479:
4475:
4467:
4456:
4441:
4438:
4434:
4426:
4415:
4412:
4408:
4404:
4401:
4397:
4394:
4390:
4386:
4385:If we pick the
4377:
4367:
4364:
4360:
4357:
4353:
4302:
4270:
4266:
4260:
4256:
4237:
4232:
4226:
4223:
4222:
4218:
4215:
4211:
4208:
4169:
4165:
4161:
4157:
4155:
4152:
4151:
4146:
4094:
4090:
4084:
4080:
4078:
4075:
4074:
4069:
4032:
4028:
4026:
4023:
4022:
4017:
3984:
3980:
3972:
3970:
3967:
3966:
3962:[−1, 1]
3961:
3924:
3920:
3908:
3906:
3903:
3902:
3897:
3881:
3835:
3820:
3816:
3815:
3809:
3794:
3790:
3789:
3787:
3784:
3783:
3778:
3758:
3754:[−1, 1]
3753:
3730:
3707:
3700:
3693:
3689:
3685:
3681:
3632:
3627:
3621:
3618:
3617:
3613:
3611:weight function
3607:
3601:
3551:
3549:
3545:
3532:
3528:
3524:
3522:
3519:
3518:
3497:
3493:
3489:
3487:
3484:
3483:
3467:
3466:
3454:
3453:
3412:
3407:
3403:
3387:
3371:
3369:
3366:
3365:
3352:
3351:
3339:
3338:
3297:
3292:
3288:
3272:
3256:
3254:
3251:
3250:
3236:
3235:
3223:
3222:
3195:
3191:
3179:
3178:
3145:
3141:
3129:
3128:
3099:
3095:
3062:
3058:
3054:
3050:
3038:
3037:
3015:
3011:
3007:
3003:
2994:
2990:
2970:
2966:
2962:
2958:
2949:
2945:
2944:
2940:
2930:
2921:
2903:
2899:
2895:
2889:
2881:
2870:
2868:
2865:
2864:
2838:
2834:
2832:
2829:
2828:
2804:
2800:
2798:
2795:
2794:
2767:
2763:
2761:
2758:
2757:
2733:
2729:
2727:
2724:
2723:
2704:
2703:
2694:
2676:
2672:
2668:
2662:
2654:
2638:
2637:
2628:
2607:
2605:
2583:
2581:
2580:
2576:
2572:
2566:
2558:
2538:
2536:
2529:
2509:
2503:
2498:
2490:
2488:
2485:
2484:
2453:
2449:
2447:
2444:
2443:
2417:
2413:
2411:
2408:
2407:
2379:
2345:
2340:
2336:
2326:
2322:
2321:
2315:
2310:
2298:
2295:
2294:
2275:
2264:
2261:
2260:
2244:
2233:
2230:
2229:
2226:
2187:
2185:
2176:
2172:
2156:
2154:
2153:
2149:
2140:
2136:
2130:
2119:
2099:
2097:
2069:
2064:
2058:
2055:
2054:
2026:
2023:
2022:
2006:
2003:
2002:
1973:
1971:
1958:
1950:
1948:
1946:
1943:
1942:
1913:
1905:
1903:
1881:
1879:
1857:
1855:
1854:
1850:
1841:
1833:
1801:
1796:
1790:
1787:
1786:
1783:
1749:
1739:
1737:
1735:
1732:
1731:
1702:
1691:
1681:
1676:
1673:
1672:
1642:
1632:
1630:
1628:
1625:
1624:
1595:
1584:
1574:
1569:
1566:
1565:
1538:
1536:
1533:
1532:
1496:
1489:
1487:
1485:
1482:
1481:
1452:
1442:
1429:
1427:
1422:
1419:
1418:
1388:
1381:
1379:
1377:
1374:
1373:
1344:
1334:
1321:
1319:
1314:
1311:
1310:
1280:
1278:
1275:
1274:
1247:
1242:
1239:
1238:
1211:
1209:
1206:
1205:
1167:
1162:
1159:
1158:
1138:
1134:
1129:
1125:
1119:
1088:
1074:
1070:
1058:
1053:
1049:
1048:
1037:
1032:
1021:
1017:
1016:
1011:
1002:
998:
996:
993:
992:
989:
985:
981:
978:
974:
970:
966:
958:
954:
946:
938:
903:
900:
899:
888:
878:
871:
865:
849:
845:
815:
811:
803:
801:
798:
797:
794:Chebyshev–Gauss
773:
769:
757:
755:
752:
751:
722:
714:
702:
692:
681:
649:
634:
630:
629:
623:
608:
604:
603:
597:
589:
557:
549:
543:
540:
539:
529:
525:
521:
517:
506:
450:
435:
431:
430:
424:
409:
405:
404:
387:
384:
383:
362:
351:
339:
314:
310:
298:
294:
288:
277:
245:
237:
231:
228:
227:
206:
203:
199:
196:
192:
184:
177:quadrature rule
164:
133:
131:
121:
117:
116:
115:
104:
93:
82:
71:
65:
62:
52:Please help to
51:
35:
31:
24:
17:
12:
11:
5:
18986:
18976:
18975:
18958:
18957:
18955:
18954:
18949:
18944:
18939:
18934:
18929:
18923:
18921:
18917:
18916:
18914:
18913:
18908:
18903:
18898:
18893:
18888:
18882:
18880:
18874:
18873:
18871:
18870:
18865:
18860:
18855:
18850:
18848:Simpson's rule
18845:
18839:
18837:
18831:
18830:
18823:
18822:
18815:
18808:
18800:
18794:
18793:
18781:
18776:
18771:
18765:
18751:
18741:
18722:
18713:
18708:
18690:
18684:
18664:
18663:External links
18661:
18659:
18658:
18640:(2): 406–430.
18629:
18615:
18592:
18587:
18570:
18541:
18509:
18503:
18473:
18468:
18448:
18423:
18414:
18385:
18363:
18336:
18330:
18305:
18284:
18250:
18216:Golub, Gene H.
18212:
18207:
18194:
18188:
18175:
18150:
18125:
18124:by Wikisource.
18106:
18081:
18065:
18032:
18029:
18026:
18023:
18020:
18017:
18014:
18011:
18008:
18003:
17999:
17993:
17988:
17984:
17963:
17960:
17957:
17954:
17951:
17948:
17943:
17939:
17935:
17932:
17929:
17926:
17923:
17920:
17915:
17910:
17906:
17893:
17860:
17857:
17854:
17851:
17848:
17845:
17842:
17839:
17836:
17833:
17830:
17827:
17822:
17817:
17813:
17798:
17768:
17745:
17743:
17740:
17738:
17737:
17726:
17723:MathWorks 2012
17715:
17704:
17693:
17681:
17669:
17657:
17645:
17633:
17622:
17610:
17608:
17605:
17603:
17600:
17578:
17577:
17564:
17561:
17549:
17538:
17535:
17524:
17523:
17510:
17504:
17499:
17496:
17493:
17480:
17464:
17461:
17453:
17450:
17445:
17440:
17437:
17430:
17419:
17418:
17405:
17399:
17394:
17391:
17388:
17375:
17359:
17356:
17348:
17345:
17340:
17335:
17332:
17325:
17314:
17313:
17300:
17297:
17285:
17274:
17264:
17253:
17242:
17241:
17228:
17225:
17213:
17202:
17199:
17188:
17187:
17174:
17168:
17163:
17160:
17147:
17132:
17126:
17121:
17115:
17110:
17107:
17100:
17089:
17088:
17075:
17069:
17064:
17061:
17048:
17033:
17027:
17022:
17016:
17011:
17008:
17001:
16991:
16980:
16969:
16968:
16955:
16952:
16940:
16929:
16926:
16915:
16914:
16901:
16898:
16886:
16872:
16869:
16863:
16852:
16851:
16838:
16835:
16823:
16812:
16802:
16791:
16780:
16779:
16766:
16763:
16751:
16740:
16737:
16726:
16725:
16712:
16709:
16697:
16683:
16680:
16674:
16664:
16653:
16642:
16641:
16628:
16625:
16613:
16602:
16599:
16588:
16587:
16574:
16571:
16559:
16548:
16538:
16527:
16516:
16515:
16511:
16506:
16502:
16497:
16474:
16471:
16468:
16465:
16462:
16459:
16455:
16452:
16449:
16446:
16441:
16438:
16435:
16432:
16429:
16426:
16422:
16413:
16408:
16404:
16400:
16396:
16393:
16390:
16387:
16383:
16378:
16373:
16370:
16367:
16364:
16361:
16358:
16351:
16346:
16342:
16338:
16334:
16331:
16328:
16324:
16319:
16312:
16309:
16306:
16303:
16299:
16293:
16288:
16284:
16281:
16278:
16274:
16269:
16266:
16260:
16255:
16251:
16228:
16225:
16222:
16217:
16213:
16208:
16200:
16195:
16190:
16185:
16181:
16177:
16171:
16168:
16165:
16161:
16156:
16151:
16148:
16145:
16142:
16139:
16136:
16132:
16127:
16122:
16118:
16090:
16087:
16084:
16079:
16075:
16054:
16051:
16048:
16044:
16040:
16037:
16034:
16030:
16009:
16006:
16003:
16000:
15997:
15984:
15968:
15963:
15959:
15955:
15951:
15946:
15942:
15938:
15935:
15930:
15926:
15919:
15916:
15913:
15908:
15905:
15902:
15898:
15894:
15891:
15888:
15885:
15882:
15879:
15876:
15873:
15870:
15867:
15864:
15861:
15858:
15852:
15849:
15846:
15843:
15840:
15837:
15833:
15828:
15824:
15821:
15817:
15814:
15811:
15808:
15802:
15797:
15794:
15790:
15780:on interval :
15767:
15766:
15751:
15744:Rehuel Lobatto
15738:Also known as
15735:
15732:
15700:Main article:
15697:
15694:
15671:
15668:
15665:
15662:
15659:
15656:
15652:
15649:
15646:
15643:
15638:
15635:
15632:
15629:
15625:
15616:
15611:
15607:
15603:
15599:
15596:
15592:
15587:
15582:
15579:
15576:
15573:
15570:
15567:
15560:
15555:
15551:
15548:
15544:
15537:
15534:
15531:
15528:
15523:
15519:
15516:
15513:
15509:
15471:
15468:
15465:
15461:
15458:
15455:
15452:
15449:
15446:
15443:
15440:
15437:
15434:
15431:
15428:
15423:
15418:
15414:
15410:
15407:
15404:
15401:
15398:
15395:
15374:
15344:
15339:
15335:
15331:
15326:
15322:
15318:
15311:
15308:
15305:
15302:
15299:
15294:
15291:
15288:
15283:
15280:
15277:
15274:
15270:
15263:
15260:
15255:
15251:
15247:
15244:
15238:
15234:
15228:
15223:
15220:
15217:
15213:
15209:
15206:
15203:
15199:
15196:
15193:
15190:
15186:
15183:
15180:
15177:
15172:
15167:
15163:
15143:
15140:
15121:
15118:
15115:
15112:
15109:
15106:
15103:
15098:
15093:
15089:
15085:
15080:
15076:
15053:
15049:
15024:
15019:
15014:
15011:
15008:
15003:
14999:
14995:
14988:
14984:
14980:
14975:
14971:
14957:
14941:
14938:
14935:
14931:
14904:
14875:
14870:
14862:
14859:
14856:
14852:
14848:
14842:
14839:
14836:
14832:
14826:
14824:
14821:
14819:
14816:
14814:
14811:
14810:
14803:
14800:
14797:
14793:
14787:
14783:
14780:
14777:
14773:
14769:
14767:
14764:
14762:
14759:
14757:
14754:
14753:
14750:
14747:
14745:
14742:
14740:
14737:
14731:
14727:
14721:
14719:
14716:
14715:
14712:
14709:
14707:
14704:
14698:
14694:
14688:
14684:
14680:
14676:
14670:
14666:
14660:
14659:
14656:
14653:
14651:
14648:
14646:
14643:
14637:
14633:
14627:
14623:
14619:
14615:
14614:
14612:
14607:
14602:
14572:
14569:
14566:
14563:
14560:
14557:
14549:
14546:
14539:
14536:
14534:
14532:
14529:
14523:
14520:
14517:
14513:
14507:
14504:
14502:
14496:
14493:
14490:
14487:
14484:
14480:
14474:
14471:
14468:
14465:
14462:
14458:
14452:
14447:
14444:
14441:
14438:
14435:
14429:
14423:
14418:
14415:
14412:
14409:
14406:
14400:
14394:
14393:
14390:
14387:
14384:
14381:
14378:
14375:
14372:
14369:
14366:
14364:
14362:
14359:
14355:
14352:
14349:
14345:
14341:
14338:
14336:
14332:
14329:
14326:
14322:
14318:
14313:
14310:
14307:
14301:
14295:
14294:
14285:with elements
14272:
14233:
14229:
14206:
14201:
14193:
14190:
14187:
14183:
14179:
14175:
14172:
14169:
14165:
14161:
14159:
14156:
14154:
14151:
14149:
14146:
14145:
14142:
14139:
14135:
14132:
14129:
14125:
14121:
14119:
14116:
14114:
14111:
14109:
14106:
14105:
14102:
14099:
14097:
14094:
14092:
14089:
14085:
14081:
14077:
14075:
14072:
14071:
14068:
14065:
14063:
14060:
14058:
14055:
14051:
14047:
14043:
14039:
14035:
14031:
14030:
14027:
14024:
14022:
14019:
14017:
14014:
14012:
14009:
14005:
14001:
13997:
13996:
13994:
13989:
13985:
13953:
13946:
13940:
13937:
13935:
13932:
13930:
13927:
13925:
13922:
13921:
13919:
13913:
13908:
13903:
13881:
13859:
13854:
13829:
13822:
13816:
13813:
13810:
13805:
13802:
13799:
13795:
13791:
13789:
13786:
13784:
13781:
13778:
13773:
13769:
13765:
13763:
13760:
13757:
13752:
13748:
13744:
13743:
13741:
13735:
13729:
13726:
13701:
13696:
13691:
13688:
13685:
13680:
13676:
13672:
13666:
13663:
13657:
13654:
13648:
13645:
13639:
13627:
13624:
13605:
13601:
13578:
13575:
13572:
13568:
13564:
13544:
13539:
13535:
13531:
13526:
13522:
13518:
13515:
13512:
13507:
13504:
13501:
13497:
13493:
13490:
13485:
13481:
13477:
13474:
13471:
13466:
13463:
13460:
13456:
13452:
13447:
13443:
13439:
13436:
13411:
13406:
13403:
13400:
13396:
13392:
13387:
13384:
13381:
13377:
13373:
13368:
13363:
13359:
13355:
13350:
13346:
13342:
13336:
13330:
13325:
13322:
13319:
13315:
13311:
13306:
13303:
13300:
13296:
13292:
13287:
13282:
13279:
13276:
13272:
13268:
13263:
13259:
13255:
13252:
13246:
13241:
13237:
13232:
13226:
13221:
13217:
13213:
13208:
13204:
13200:
13195:
13190:
13186:
13182:
13177:
13173:
13169:
13166:
13160:
13155:
13151:
13130:
13127:
13124:
13121:
13118:
13113:
13110:
13106:
13083:
13080:
13077:
13072:
13069:
13066:
13062:
13056:
13052:
13048:
13045:
13042:
13039:
13034:
13030:
13026:
13021:
13017:
13013:
13010:
13007:
13004:
13001:
12998:
12995:
12990:
12987:
12984:
12980:
12957:
12954:
12951:
12946:
12943:
12940:
12936:
12930:
12927:
12924:
12921:
12918:
12914:
12910:
12907:
12904:
12901:
12896:
12892:
12888:
12883:
12880:
12877:
12873:
12869:
12866:
12863:
12860:
12857:
12854:
12851:
12846:
12843:
12840:
12836:
12804:
12801:
12796:
12793:
12790:
12786:
12765:
12762:
12759:
12754:
12750:
12746:
12743:
12738:
12734:
12730:
12727:
12724:
12719:
12715:
12711:
12706:
12702:
12698:
12695:
12666:
12662:
12632:
12628:
12624:
12621:
12618:
12615:
12612:
12609:
12606:
12586:
12583:
12580:
12577:
12574:
12571:
12568:
12565:
12562:
12559:
12556:
12553:
12550:
12528:
12525:
12522:
12517:
12513:
12509:
12504:
12500:
12496:
12493:
12490:
12487:
12482:
12478:
12474:
12469:
12465:
12461:
12458:
12455:
12452:
12447:
12443:
12439:
12434:
12430:
12426:
12421:
12418:
12415:
12411:
12407:
12404:
12399:
12395:
12391:
12386:
12382:
12378:
12375:
12372:
12369:
12364:
12360:
12356:
12351:
12348:
12345:
12341:
12337:
12315:
12311:
12290:
12285:
12281:
12277:
12272:
12268:
12264:
12259:
12256:
12253:
12249:
12245:
12242:
12239:
12234:
12230:
12226:
12221:
12218:
12215:
12211:
12207:
12202:
12199:
12196:
12193:
12190:
12186:
12182:
12179:
12174:
12170:
12166:
12161:
12157:
12153:
12148:
12145:
12142:
12138:
12134:
12131:
12126:
12122:
12118:
12113:
12109:
12105:
12102:
12099:
12096:
12091:
12087:
12083:
12078:
12075:
12072:
12068:
12064:
12042:
12039:
12036:
12032:
12009:
12005:
12001:
11998:
11995:
11990:
11986:
11982:
11977:
11973:
11950:
11947:
11944:
11939:
11935:
11931:
11926:
11922:
11918:
11915:
11912:
11909:
11904:
11900:
11896:
11891:
11887:
11883:
11880:
11877:
11874:
11869:
11865:
11861:
11856:
11852:
11848:
11843:
11840:
11837:
11833:
11829:
11826:
11821:
11817:
11813:
11808:
11804:
11800:
11797:
11794:
11791:
11786:
11782:
11778:
11773:
11769:
11765:
11762:
11757:
11754:
11751:
11747:
11743:
11740:
11737:
11734:
11731:
11726:
11722:
11718:
11713:
11709:
11705:
11685:
11682:
11679:
11676:
11673:
11651:
11647:
11624:
11620:
11599:
11596:
11593:
11590:
11587:
11582:
11578:
11554:
11548:
11544:
11540:
11535:
11531:
11526:
11521:
11515:
11511:
11507:
11502:
11498:
11494:
11490:
11484:
11479:
11476:
11473:
11469:
11444:
11441:
11438:
11435:
11432:
11429:
11426:
11423:
11420:
11417:
11414:
11392:
11389:
11386:
11383:
11380:
11377:
11374:
11371:
11368:
11365:
11362:
11359:
11356:
11353:
11348:
11343:
11339:
11335:
11332:
11329:
11326:
11323:
11320:
11317:
11314:
11311:
11308:
11305:
11302:
11280:
11277:
11274:
11269:
11265:
11259:
11256:
11253:
11249:
11245:
11242:
11239:
11236:
11233:
11228:
11225:
11222:
11218:
11212:
11209:
11206:
11203:
11200:
11196:
11192:
11189:
11186:
11183:
11178:
11174:
11170:
11165:
11162:
11159:
11155:
11151:
11148:
11145:
11142:
11139:
11136:
11133:
11128:
11125:
11122:
11118:
11093:
11090:
11087:
11082:
11078:
11074:
11054:
11051:
11048:
11045:
11042:
11022:
11019:
11016:
10996:
10993:
10990:
10985:
10981:
10977:
10972:
10968:
10964:
10942:
10938:
10925:
10922:
10879:
10876:
10873:
10870:
10867:
10862:
10858:
10833:
10826:
10820:
10817:
10804:
10801:
10796:
10792:
10771:
10768:
10765:
10762:
10742:
10739:
10736:
10733:
10711:
10708:
10703:
10699:
10695:
10690:
10686:
10680:
10677:
10673:
10667:
10662:
10659:
10656:
10652:
10648:
10645:
10640:
10636:
10632:
10629:
10624:
10620:
10614:
10609:
10606:
10603:
10599:
10595:
10592:
10589:
10586:
10583:
10580:
10577:
10574:
10571:
10568:
10565:
10560:
10555:
10551:
10522:
10519:
10516:
10511:
10507:
10503:
10500:
10480:
10477:
10474:
10469:
10465:
10444:
10441:
10438:
10435:
10415:
10412:
10409:
10406:
10384:
10381:
10377:
10373:
10370:
10365:
10361:
10357:
10354:
10334:
10331:
10328:
10323:
10319:
10305:
10287:
10282:
10276:
10272:
10268:
10263:
10259:
10254:
10247:
10242:
10236:
10232:
10228:
10225:
10221:
10208:
10205:
10202:
10199:
10198:
10195:
10192:
10189:
10186:
10183:
10180:
10179:
10175:
10171:
10168:
10165:
10162:
10159:
10140:
10137:
10134:
10131:
10119:
10116:
10102:
10097:
10093:
10089:
10083:
10080:
10077:
10073:
10067:
10064:
10042:
10037:
10033:
10029:
10024:
10021:
10018:
10014:
9993:
9988:
9984:
9980:
9975:
9971:
9961:the term with
9950:
9945:
9941:
9937:
9932:
9929:
9926:
9922:
9918:
9915:
9912:
9909:
9906:
9901:
9897:
9893:
9888:
9884:
9880:
9877:
9874:
9871:
9868:
9863:
9859:
9855:
9850:
9847:
9844:
9840:
9819:
9816:
9813:
9808:
9805:
9802:
9798:
9777:
9774:
9771:
9766:
9762:
9739:
9735:
9706:
9701:
9697:
9693:
9689:
9685:
9681:
9652:
9649:
9644:
9640:
9636:
9633:
9628:
9625:
9622:
9618:
9614:
9611:
9608:
9605:
9600:
9595:
9591:
9584:
9579:
9575:
9571:
9566:
9563:
9560:
9556:
9550:
9547:
9544:
9540:
9533:
9529:
9523:
9520:
9517:
9509:
9505:
9501:
9498:
9493:
9490:
9487:
9482:
9478:
9471:
9468:
9465:
9462:
9457:
9452:
9448:
9427:
9424:
9421:
9416:
9413:
9410:
9406:
9385:
9382:
9379:
9353:
9350:
9347:
9343:
9338:
9335:
9332:
9327:
9324:
9321:
9317:
9310:
9306:
9298:
9295:
9292:
9288:
9283:
9280:
9277:
9272:
9269:
9266:
9262:
9255:
9250:
9247:
9244:
9240:
9235:
9231:
9226:
9223:
9220:
9216:
9193:
9190:
9185:
9182:
9179:
9175:
9171:
9168:
9165:
9160:
9157:
9154:
9150:
9146:
9143:
9140:
9137:
9132:
9127:
9123:
9116:
9111:
9107:
9103:
9098:
9095:
9092:
9088:
9081:
9077:
9071:
9068:
9065:
9057:
9053:
9049:
9046:
9041:
9038:
9035:
9030:
9026:
9019:
9016:
9013:
9010:
9005:
9000:
8996:
8975:
8972:
8969:
8964:
8961:
8958:
8954:
8922:
8919:
8916:
8885:
8882:
8879:
8876:
8873:
8868:
8865:
8862:
8858:
8852:
8848:
8844:
8836:
8832:
8828:
8825:
8820:
8817:
8814:
8809:
8805:
8769:
8765:
8761:
8758:
8753:
8750:
8747:
8742:
8738:
8714:
8711:
8708:
8703:
8700:
8697:
8693:
8689:
8686:
8683:
8680:
8677:
8655:
8652:
8644:
8640:
8636:
8633:
8628:
8625:
8622:
8617:
8613:
8609:
8606:
8603:
8600:
8594:
8591:
8588:
8585:
8580:
8575:
8571:
8564:
8559:
8555:
8551:
8548:
8544:
8539:
8536:
8533:
8525:
8521:
8517:
8514:
8509:
8506:
8503:
8498:
8494:
8487:
8484:
8481:
8478:
8473:
8468:
8464:
8432:
8429:
8426:
8421:
8417:
8384:
8380:
8376:
8373:
8366:
8361:
8354:
8350:
8346:
8341:
8336:
8333:
8310:
8307:
8304:
8282:
8279:
8271:
8267:
8263:
8260:
8255:
8252:
8249:
8244:
8240:
8233:
8230:
8227:
8224:
8219:
8214:
8210:
8204:
8199:
8195:
8191:
8188:
8185:
8177:
8173:
8169:
8166:
8161:
8158:
8155:
8150:
8146:
8140:
8136:
8129:
8126:
8123:
8120:
8115:
8110:
8106:
8078:
8074:
8070:
8067:
8063:
8056:
8051:
8044:
8040:
8036:
8031:
8026:
8018:
8014:
8010:
8007:
8000:
7995:
7988:
7984:
7980:
7975:
7970:
7967:
7961:
7953:
7949:
7945:
7942:
7938:
7923:
7922:
7913:
7911:
7900:
7897:
7889:
7885:
7881:
7878:
7873:
7870:
7867:
7862:
7858:
7851:
7848:
7845:
7842:
7837:
7832:
7828:
7821:
7816:
7812:
7808:
7804:
7800:
7796:
7791:
7786:
7781:
7777:
7745:
7741:
7736:
7731:
7727:
7723:
7719:
7715:
7711:
7704:
7700:
7694:
7690:
7686:
7681:
7677:
7672:
7662:
7659:
7656:
7653:
7652:
7649:
7646:
7643:
7640:
7637:
7634:
7633:
7629:
7606:
7602:
7577:
7574:
7571:
7566:
7562:
7539:
7535:
7512:
7508:
7481:
7475:
7471:
7467:
7464:
7460:
7454:
7450:
7444:
7441:
7438:
7433:
7429:
7422:
7414:
7410:
7406:
7403:
7397:
7391:
7387:
7383:
7380:
7376:
7370:
7367:
7364:
7361:
7358:
7354:
7347:
7343:
7337:
7333:
7329:
7326:
7322:
7312:
7309:
7306:
7303:
7302:
7299:
7296:
7293:
7290:
7287:
7284:
7283:
7279:
7255:
7252:
7249:
7244:
7241:
7238:
7234:
7213:
7210:
7207:
7202:
7198:
7175:
7171:
7148:
7145:
7137:
7133:
7129:
7124:
7120:
7112:
7108:
7104:
7101:
7089:
7086:
7083:
7080:
7079:
7076:
7073:
7070:
7067:
7064:
7061:
7060:
7056:
7052:
7049:
7046:
7043:
7038:
7033:
7029:
7025:
7020:
7016:
7002:
6986:
6983:
6975:
6971:
6967:
6962:
6958:
6950:
6946:
6942:
6939:
6927:
6924:
6921:
6918:
6917:
6914:
6911:
6908:
6905:
6902:
6899:
6898:
6894:
6890:
6887:
6884:
6881:
6876:
6871:
6867:
6863:
6858:
6854:
6850:
6847:
6842:
6837:
6834:
6831:
6827:
6823:
6820:
6817:
6814:
6811:
6808:
6805:
6802:
6799:
6796:
6793:
6788:
6783:
6779:
6715:
6711:
6707:
6702:
6698:
6690:
6686:
6682:
6679:
6667:
6664:
6661:
6658:
6657:
6654:
6651:
6648:
6645:
6642:
6639:
6638:
6634:
6630:
6625:
6621:
6617:
6614:
6609:
6604:
6601:
6598:
6594:
6590:
6587:
6584:
6581:
6578:
6558:
6553:
6549:
6545:
6542:
6507:
6504:
6501:
6496:
6492:
6469:
6465:
6442:
6438:
6424:
6423:
6414:
6412:
6398:
6393:
6389:
6385:
6380:
6377:
6374:
6370:
6366:
6361:
6357:
6353:
6349:
6345:
6341:
6335:
6332:
6327:
6323:
6319:
6316:
6311:
6308:
6305:
6301:
6297:
6294:
6291:
6288:
6283:
6278:
6274:
6263:
6260:
6257:
6253:
6247:
6243:
6237:
6232:
6228:
6210:
6207:
6187:
6167:
6162:
6158:
6154:
6151:
6148:
6143:
6139:
6116:
6089:
6086:
6083:
6080:
6077:
6073:
6070:
6067:
6064:
6059:
6055:
6050:
6046:
6042:
6039:
6036:
6031:
6027:
6022:
6017:
6014:
6011:
6006:
6002:
5996:
5991:
5987:
5956:
5953:
5948:
5944:
5940:
5937:
5934:
5929:
5925:
5920:
5917:
5914:
5909:
5905:
5889:
5872:
5855:
5815:
5812:
5807:
5803:
5799:
5796:
5790:
5786:
5780:
5775:
5772:
5769:
5765:
5761:
5758:
5753:
5749:
5745:
5742:
5736:
5732:
5726:
5721:
5718:
5715:
5711:
5707:
5704:
5701:
5697:
5694:
5691:
5688:
5684:
5681:
5678:
5675:
5670:
5665:
5661:
5657:
5654:
5651:
5647:
5644:
5641:
5638:
5634:
5631:
5628:
5625:
5620:
5615:
5611:
5586:
5583:
5578:
5574:
5570:
5567:
5564:
5561:
5556:
5552:
5548:
5545:
5542:
5539:
5534:
5530:
5526:
5523:
5519:
5514:
5510:
5506:
5501:
5497:
5493:
5490:
5485:
5481:
5477:
5474:
5461:
5454:
5440:
5429:
5418:
5400:
5395:
5391:
5386:
5381:
5377:
5373:
5370:
5364:
5359:
5356:
5353:
5349:
5345:
5342:
5339:
5335:
5332:
5329:
5324:
5320:
5315:
5312:
5309:
5306:
5301:
5296:
5292:
5287:
5282:
5278:
5274:
5271:
5265:
5260:
5257:
5254:
5250:
5246:
5243:
5240:
5236:
5231:
5227:
5223:
5220:
5216:
5213:
5210:
5205:
5201:
5195:
5190:
5187:
5184:
5180:
5175:
5172:
5169:
5166:
5161:
5156:
5152:
5148:
5145:
5142:
5138:
5135:
5132:
5129:
5125:
5122:
5119:
5116:
5111:
5106:
5102:
5079:
5076:
5071:
5067:
5063:
5060:
5056:
5053:
5050:
5045:
5041:
5035:
5030:
5027:
5024:
5020:
5016:
5013:
5010:
5007:
5004:
4982:
4974:
4970:
4966:
4961:
4957:
4949:
4945:
4941:
4938:
4930:
4927:
4924:
4920:
4916:
4913:
4910:
4907:
4902:
4898:
4877:
4834:
4831:
4828:
4824:
4821:
4818:
4815:
4811:
4808:
4805:
4802:
4797:
4792:
4788:
4784:
4781:
4778:
4772:
4766:
4763:
4760:
4757:
4754:
4751:
4748:
4745:
4742:
4739:
4736:
4733:
4728:
4724:
4717:
4711:
4708:
4705:
4702:
4697:
4692:
4688:
4684:
4681:
4678:
4674:
4671:
4668:
4665:
4661:
4658:
4655:
4652:
4647:
4642:
4638:
4609:
4584:
4552:
4549:
4546:
4543:
4540:
4537:
4534:
4531:
4528:
4525:
4521:
4518:
4515:
4510:
4506:
4502:
4499:
4496:
4493:
4490:
4477:
4436:
4410:
4399:
4392:
4362:
4355:
4339:
4336:
4333:
4330:
4327:
4324:
4321:
4318:
4315:
4312:
4309:
4300:
4297:
4294:
4291:
4288:
4284:
4281:
4278:
4273:
4269:
4263:
4259:
4254:
4251:
4248:
4245:
4240:
4235:
4231:
4213:
4207:
4204:
4201:
4200:
4195:
4192:
4187:
4172:
4168:
4164:
4160:
4149:
4143:
4142:
4137:
4135:
4129:
4118:
4115:
4112:
4109:
4105:
4100:
4097:
4093:
4087:
4083:
4072:
4066:
4065:
4060:
4057:
4052:
4038:
4035:
4031:
4020:
4014:
4013:
4008:
4005:
4002:
3987:
3983:
3979:
3976:
3964:
3958:
3957:
3952:
3949:
3943:
3927:
3923:
3919:
3916:
3912:
3900:
3894:
3893:
3888:
3878:
3873:
3862:
3859:
3856:
3853:
3850:
3847:
3843:
3838:
3833:
3829:
3826:
3823:
3819:
3812:
3807:
3803:
3800:
3797:
3793:
3781:
3775:
3774:
3769:
3766:
3761:
3756:
3750:
3749:
3746:
3743:
3740:
3728:
3669:
3666:
3662:
3659:
3656:
3653:
3649:
3646:
3643:
3640:
3635:
3630:
3626:
3606:
3603:
3588:
3585:
3582:
3579:
3576:
3573:
3569:
3564:
3560:
3557:
3554:
3548:
3544:
3540:
3535:
3531:
3527:
3505:
3500:
3496:
3492:
3465:
3462:
3459:
3457:
3455:
3452:
3449:
3446:
3443:
3440:
3436:
3430:
3427:
3424:
3421:
3418:
3415:
3411:
3406:
3402:
3399:
3396:
3393:
3390:
3388:
3386:
3383:
3380:
3377:
3374:
3373:
3350:
3347:
3344:
3342:
3340:
3337:
3334:
3331:
3328:
3325:
3321:
3315:
3312:
3309:
3306:
3303:
3300:
3296:
3291:
3287:
3284:
3281:
3278:
3275:
3273:
3271:
3268:
3265:
3262:
3259:
3258:
3234:
3231:
3228:
3226:
3224:
3220:
3216:
3213:
3210:
3207:
3204:
3201:
3198:
3194:
3190:
3187:
3184:
3182:
3180:
3176:
3172:
3169:
3166:
3163:
3160:
3157:
3154:
3151:
3148:
3144:
3140:
3137:
3134:
3132:
3130:
3126:
3121:
3117:
3114:
3111:
3108:
3105:
3102:
3098:
3094:
3091:
3087:
3083:
3080:
3077:
3074:
3071:
3068:
3065:
3061:
3057:
3053:
3049:
3046:
3043:
3041:
3039:
3035:
3030:
3026:
3023:
3018:
3014:
3010:
3006:
3002:
2997:
2993:
2989:
2985:
2981:
2978:
2973:
2969:
2965:
2961:
2957:
2952:
2948:
2943:
2939:
2936:
2933:
2931:
2927:
2924:
2919:
2915:
2912:
2909:
2906:
2902:
2898:
2892:
2887:
2884:
2880:
2876:
2873:
2872:
2861:
2860:
2849:
2846:
2841:
2837:
2826:
2815:
2812:
2807:
2803:
2792:
2781:
2778:
2775:
2770:
2766:
2755:
2744:
2741:
2736:
2732:
2700:
2697:
2692:
2688:
2685:
2682:
2679:
2675:
2671:
2665:
2660:
2657:
2653:
2649:
2646:
2643:
2641:
2639:
2634:
2631:
2626:
2620:
2616:
2613:
2610:
2604:
2601:
2596:
2592:
2589:
2586:
2579:
2575:
2569:
2564:
2561:
2557:
2551:
2547:
2544:
2541:
2535:
2532:
2530:
2527:
2524:
2521:
2518:
2515:
2512:
2506:
2501:
2497:
2493:
2492:
2469:
2465:
2462:
2459:
2456:
2452:
2430:
2426:
2423:
2420:
2416:
2385:
2382:
2377:
2373:
2370:
2367:
2363:
2357:
2354:
2351:
2348:
2344:
2339:
2335:
2332:
2329:
2325:
2318:
2313:
2309:
2305:
2302:
2282:
2278:
2274:
2271:
2268:
2247:
2243:
2240:
2237:
2225:
2222:
2210:
2206:
2200:
2196:
2193:
2190:
2184:
2179:
2175:
2169:
2165:
2162:
2159:
2152:
2148:
2143:
2139:
2133:
2128:
2125:
2122:
2118:
2112:
2108:
2105:
2102:
2096:
2093:
2090:
2086:
2083:
2080:
2077:
2072:
2067:
2063:
2042:
2039:
2036:
2033:
2030:
2010:
1986:
1982:
1979:
1976:
1970:
1964:
1961:
1956:
1953:
1928:
1925:
1919:
1916:
1911:
1908:
1900:
1894:
1890:
1887:
1884:
1878:
1875:
1870:
1866:
1863:
1860:
1853:
1849:
1844:
1839:
1836:
1832:
1828:
1825:
1822:
1818:
1815:
1812:
1809:
1804:
1799:
1795:
1782:
1779:
1776:
1775:
1772:
1759:
1753:
1748:
1745:
1742:
1729:
1726:
1710:
1707:
1701:
1698:
1695:
1688:
1685:
1680:
1669:
1668:
1665:
1652:
1646:
1641:
1638:
1635:
1622:
1619:
1603:
1600:
1594:
1591:
1588:
1581:
1578:
1573:
1562:
1561:
1558:
1545:
1542:
1530:
1527:
1523:
1522:
1519:
1506:
1500:
1495:
1492:
1479:
1476:
1460:
1457:
1449:
1446:
1441:
1436:
1433:
1426:
1415:
1414:
1411:
1398:
1392:
1387:
1384:
1371:
1368:
1352:
1349:
1341:
1338:
1333:
1328:
1325:
1318:
1308:
1304:
1303:
1300:
1287:
1284:
1272:
1269:
1255:
1252:
1246:
1235:
1234:
1231:
1218:
1215:
1203:
1200:
1196:
1195:
1192:
1189:
1175:
1171:
1166:
1156:
1152:
1151:
1148:
1145:
1141:
1140:
1136:
1131:
1127:
1122:
1099:
1091:
1086:
1082:
1077:
1073:
1069:
1065:
1061:
1057:
1052:
1046:
1040:
1035:
1031:
1027:
1024:
1020:
1015:
1010:
1005:
1001:
987:
976:
962:
942:
922:
919:
916:
913:
910:
907:
864:
861:
847:
818:
814:
810:
807:
776:
772:
768:
765:
761:
737:
733:
729:
725:
721:
717:
713:
709:
705:
701:
695:
690:
687:
684:
680:
676:
673:
670:
666:
663:
660:
657:
652:
647:
643:
640:
637:
633:
626:
621:
617:
614:
611:
607:
600:
595:
592:
588:
584:
581:
578:
574:
571:
568:
565:
560:
555:
552:
548:
527:
519:
492:
489:
486:
483:
480:
477:
474:
470:
467:
464:
461:
458:
453:
448:
444:
441:
438:
434:
427:
422:
418:
415:
412:
408:
403:
400:
397:
394:
391:
325:
322:
317:
313:
309:
306:
301:
297:
291:
286:
283:
280:
276:
272:
269:
266:
262:
259:
256:
253:
248:
243:
240:
236:
201:
194:
171:, named after
84:
83:
66:September 2018
38:
36:
29:
15:
9:
6:
4:
3:
2:
18985:
18974:
18971:
18970:
18968:
18953:
18950:
18948:
18945:
18943:
18940:
18938:
18935:
18933:
18930:
18928:
18925:
18924:
18922:
18918:
18912:
18909:
18907:
18904:
18902:
18899:
18897:
18894:
18892:
18889:
18887:
18884:
18883:
18881:
18879:
18875:
18869:
18866:
18864:
18861:
18859:
18856:
18854:
18851:
18849:
18846:
18844:
18841:
18840:
18838:
18836:
18832:
18828:
18821:
18816:
18814:
18809:
18807:
18802:
18801:
18798:
18792:
18788:
18785:
18782:
18780:
18777:
18775:
18772:
18769:
18766:
18763:
18759:
18755:
18752:
18749:
18745:
18742:
18737:
18736:
18731:
18728:
18723:
18721:
18717:
18714:
18712:
18709:
18706:
18702:
18698:
18694:
18691:
18688:
18685:
18681:
18677:
18676:
18671:
18667:
18666:
18655:
18651:
18647:
18643:
18639:
18635:
18630:
18626:
18622:
18618:
18612:
18608:
18607:
18602:
18598:
18593:
18590:
18584:
18580:
18576:
18571:
18567:
18563:
18559:
18555:
18551:
18547:
18542:
18538:
18534:
18529:
18524:
18520:
18516:
18510:
18506:
18504:0-387-98959-5
18500:
18496:
18492:
18488:
18484:
18483:
18478:
18474:
18471:
18465:
18461:
18457:
18453:
18449:
18445:
18441:
18437:
18433:
18429:
18424:
18420:
18415:
18411:
18407:
18403:
18399:
18395:
18391:
18386:
18382:
18377:
18373:
18369:
18364:
18359:
18354:
18350:
18346:
18342:
18337:
18333:
18327:
18323:
18322:Prentice-Hall
18318:
18317:
18311:
18306:
18302:
18298:
18294:
18290:
18285:
18281:
18275:
18267:
18263:
18259:
18255:
18251:
18247:
18243:
18238:
18233:
18229:
18225:
18221:
18217:
18213:
18210:
18204:
18200:
18195:
18191:
18185:
18181:
18176:
18172:
18168:
18164:
18160:
18156:
18151:
18147:
18143:
18139:
18135:
18131:
18126:
18123:
18117:
18116:
18111:
18107:
18103:
18099:
18096:(1): 84–101.
18095:
18091:
18087:
18082:
18071:
18066:
18062:
18058:
18054:
18050:
18046:
18030:
18027:
18021:
18015:
18009:
18001:
17997:
17986:
17982:
17961:
17958:
17952:
17946:
17941:
17937:
17930:
17927:
17924:
17921:
17913:
17908:
17904:
17894:
17889:
17884:
17880:
17876:
17872:
17858:
17855:
17849:
17843:
17837:
17831:
17828:
17825:
17820:
17815:
17811:
17799:
17795:
17791:
17787:
17783:
17779:
17775:
17771:
17765:
17761:
17760:
17755:
17751:
17747:
17746:
17735:
17730:
17724:
17719:
17713:
17708:
17702:
17697:
17691:, p. 888
17690:
17685:
17678:
17673:
17666:
17661:
17654:
17649:
17643:, p. 887
17642:
17637:
17631:
17626:
17620:
17615:
17611:
17599:
17589:
17585:
17562:
17559:
17550:
17536:
17533:
17526:
17525:
17508:
17502:
17497:
17494:
17491:
17481:
17462:
17459:
17451:
17448:
17443:
17438:
17435:
17428:
17421:
17420:
17403:
17397:
17392:
17389:
17386:
17376:
17357:
17354:
17346:
17343:
17338:
17333:
17330:
17323:
17316:
17315:
17298:
17295:
17286:
17272:
17265:
17251:
17243:
17226:
17223:
17214:
17200:
17197:
17190:
17189:
17172:
17166:
17161:
17158:
17148:
17130:
17124:
17119:
17113:
17108:
17105:
17098:
17091:
17090:
17073:
17067:
17062:
17059:
17049:
17031:
17025:
17020:
17014:
17009:
17006:
16999:
16992:
16978:
16970:
16953:
16950:
16941:
16927:
16924:
16917:
16916:
16899:
16896:
16887:
16870:
16867:
16861:
16854:
16853:
16836:
16833:
16824:
16810:
16803:
16789:
16781:
16764:
16761:
16752:
16738:
16735:
16728:
16727:
16710:
16707:
16698:
16681:
16678:
16672:
16665:
16651:
16643:
16626:
16623:
16614:
16600:
16597:
16590:
16589:
16572:
16569:
16560:
16546:
16539:
16525:
16517:
16507:
16498:
16496:
16492:
16491:
16488:
16485:
16472:
16469:
16466:
16463:
16460:
16457:
16453:
16447:
16436:
16433:
16430:
16427:
16420:
16411:
16406:
16402:
16398:
16394:
16391:
16388:
16385:
16381:
16376:
16368:
16365:
16362:
16359:
16349:
16344:
16340:
16336:
16332:
16329:
16326:
16322:
16317:
16310:
16307:
16304:
16301:
16297:
16291:
16286:
16282:
16279:
16276:
16272:
16267:
16264:
16258:
16253:
16249:
16239:
16226:
16223:
16220:
16215:
16211:
16206:
16198:
16193:
16188:
16183:
16179:
16175:
16169:
16166:
16163:
16159:
16154:
16146:
16143:
16140:
16134:
16130:
16125:
16120:
16116:
16106:
16085:
16077:
16073:
16049:
16042:
16038:
16035:
16032:
16028:
16004:
16001:
15998:
15979:
15966:
15961:
15957:
15953:
15944:
15940:
15933:
15928:
15924:
15917:
15914:
15911:
15906:
15903:
15900:
15896:
15892:
15883:
15880:
15874:
15871:
15865:
15859:
15847:
15844:
15841:
15835:
15831:
15826:
15822:
15819:
15812:
15806:
15800:
15795:
15792:
15788:
15777:
15773:
15758:
15752:
15749:
15748:
15747:
15745:
15741:
15731:
15727:
15718:points to an
15715:
15710:
15703:
15693:
15689:
15682:
15669:
15666:
15663:
15660:
15657:
15654:
15650:
15644:
15633:
15630:
15623:
15614:
15609:
15605:
15601:
15597:
15594:
15590:
15585:
15577:
15574:
15571:
15568:
15558:
15553:
15549:
15546:
15542:
15535:
15532:
15529:
15526:
15521:
15517:
15514:
15511:
15507:
15492:
15488:
15482:
15469:
15466:
15463:
15456:
15450:
15444:
15438:
15432:
15426:
15421:
15416:
15412:
15408:
15402:
15399:
15396:
15368:
15364:
15337:
15333:
15329:
15324:
15320:
15309:
15303:
15300:
15289:
15278:
15275:
15268:
15261:
15253:
15249:
15242:
15236:
15232:
15226:
15221:
15218:
15215:
15211:
15207:
15204:
15201:
15194:
15188:
15181:
15175:
15170:
15165:
15161:
15151:
15139:
15137:
15132:
15119:
15116:
15113:
15107:
15101:
15096:
15091:
15087:
15083:
15078:
15074:
15051:
15047:
15037:
15022:
15017:
15009:
15001:
14997:
14993:
14986:
14982:
14978:
14973:
14969:
14936:
14929:
14920:
14891:
14886:
14873:
14868:
14860:
14857:
14854:
14850:
14840:
14837:
14834:
14830:
14822:
14817:
14812:
14801:
14798:
14795:
14791:
14781:
14778:
14775:
14771:
14765:
14760:
14755:
14748:
14743:
14738:
14729:
14725:
14717:
14710:
14705:
14696:
14692:
14682:
14678:
14668:
14664:
14654:
14649:
14644:
14635:
14631:
14621:
14617:
14610:
14605:
14590:
14587:
14570:
14567:
14564:
14561:
14558:
14555:
14547:
14544:
14537:
14535:
14530:
14521:
14518:
14515:
14511:
14505:
14503:
14494:
14491:
14488:
14485:
14482:
14478:
14472:
14469:
14466:
14463:
14460:
14456:
14450:
14445:
14442:
14439:
14436:
14433:
14421:
14416:
14413:
14410:
14407:
14404:
14388:
14385:
14382:
14379:
14376:
14373:
14370:
14367:
14365:
14360:
14353:
14350:
14347:
14343:
14339:
14337:
14330:
14327:
14324:
14320:
14316:
14311:
14308:
14305:
14260:
14255:
14253:
14231:
14227:
14217:
14204:
14199:
14191:
14188:
14185:
14181:
14173:
14170:
14167:
14163:
14157:
14152:
14147:
14140:
14133:
14130:
14127:
14123:
14117:
14112:
14107:
14100:
14095:
14090:
14083:
14079:
14073:
14066:
14061:
14056:
14049:
14045:
14037:
14033:
14025:
14020:
14015:
14010:
14003:
13999:
13992:
13987:
13974:
13944:
13938:
13933:
13928:
13923:
13917:
13911:
13906:
13879:
13857:
13820:
13811:
13803:
13800:
13797:
13793:
13787:
13779:
13771:
13767:
13758:
13750:
13746:
13739:
13733:
13724:
13699:
13686:
13678:
13674:
13670:
13661:
13655:
13652:
13643:
13637:
13623:
13603:
13599:
13591:differs from
13576:
13573:
13570:
13566:
13562:
13537:
13533:
13529:
13524:
13520:
13513:
13505:
13502:
13499:
13495:
13491:
13488:
13483:
13479:
13472:
13464:
13461:
13458:
13454:
13450:
13445:
13441:
13437:
13425:
13404:
13401:
13398:
13394:
13390:
13385:
13382:
13379:
13375:
13361:
13357:
13353:
13348:
13344:
13334:
13323:
13320:
13317:
13313:
13309:
13304:
13301:
13298:
13294:
13280:
13277:
13274:
13270:
13266:
13261:
13257:
13253:
13244:
13239:
13235:
13230:
13219:
13215:
13211:
13206:
13202:
13188:
13184:
13180:
13175:
13171:
13167:
13158:
13153:
13149:
13128:
13125:
13119:
13111:
13108:
13104:
13094:
13078:
13070:
13067:
13064:
13060:
13054:
13050:
13046:
13040:
13032:
13028:
13019:
13015:
13011:
13008:
13002:
12996:
12988:
12985:
12982:
12978:
12968:
12952:
12944:
12941:
12938:
12934:
12928:
12925:
12922:
12919:
12916:
12912:
12908:
12902:
12894:
12890:
12881:
12878:
12875:
12871:
12867:
12864:
12858:
12852:
12844:
12841:
12838:
12834:
12823:
12819:
12802:
12799:
12794:
12791:
12788:
12784:
12763:
12760:
12752:
12748:
12744:
12741:
12736:
12732:
12725:
12717:
12713:
12709:
12704:
12700:
12696:
12683:
12664:
12660:
12649:
12630:
12626:
12622:
12619:
12616:
12613:
12610:
12607:
12604:
12581:
12578:
12575:
12572:
12566:
12560:
12557:
12554:
12551:
12539:
12526:
12523:
12515:
12511:
12507:
12502:
12498:
12494:
12488:
12480:
12476:
12472:
12467:
12463:
12459:
12453:
12445:
12441:
12437:
12432:
12428:
12419:
12416:
12413:
12409:
12405:
12397:
12393:
12389:
12384:
12380:
12376:
12370:
12362:
12358:
12354:
12349:
12346:
12343:
12339:
12313:
12309:
12283:
12279:
12275:
12270:
12266:
12257:
12254:
12251:
12247:
12243:
12240:
12232:
12228:
12224:
12219:
12216:
12213:
12209:
12200:
12197:
12194:
12191:
12188:
12184:
12180:
12172:
12168:
12164:
12159:
12155:
12146:
12143:
12140:
12136:
12132:
12124:
12120:
12116:
12111:
12107:
12103:
12097:
12089:
12085:
12081:
12076:
12073:
12070:
12066:
12055:, because in
12040:
12037:
12034:
12030:
12007:
12003:
11999:
11996:
11993:
11988:
11984:
11980:
11975:
11971:
11961:
11948:
11945:
11937:
11933:
11929:
11924:
11920:
11916:
11910:
11902:
11898:
11894:
11889:
11885:
11881:
11875:
11867:
11863:
11859:
11854:
11850:
11841:
11838:
11835:
11831:
11827:
11819:
11815:
11811:
11806:
11802:
11798:
11792:
11784:
11780:
11776:
11771:
11767:
11755:
11752:
11749:
11745:
11741:
11738:
11732:
11724:
11720:
11716:
11711:
11707:
11683:
11680:
11677:
11674:
11671:
11649:
11645:
11622:
11618:
11597:
11594:
11588:
11580:
11576:
11552:
11546:
11542:
11538:
11533:
11529:
11524:
11519:
11513:
11509:
11505:
11500:
11496:
11492:
11488:
11482:
11477:
11474:
11471:
11467:
11442:
11439:
11436:
11433:
11430:
11427:
11424:
11421:
11418:
11415:
11412:
11403:
11390:
11387:
11381:
11375:
11369:
11363:
11357:
11351:
11346:
11341:
11337:
11333:
11324:
11318:
11315:
11309:
11303:
11291:
11275:
11267:
11263:
11257:
11254:
11251:
11247:
11243:
11240:
11234:
11226:
11223:
11220:
11216:
11210:
11207:
11204:
11201:
11198:
11194:
11190:
11184:
11176:
11172:
11163:
11160:
11157:
11153:
11149:
11146:
11140:
11134:
11126:
11123:
11120:
11116:
11107:
11091:
11088:
11080:
11076:
11049:
11046:
11043:
11020:
11017:
11014:
10994:
10991:
10983:
10979:
10975:
10970:
10966:
10940:
10936:
10921:
10917:
10913:
10908:
10902:
10898:
10893:
10877:
10874:
10868:
10860:
10856:
10845:
10841:
10816:
10802:
10799:
10794:
10790:
10766:
10760:
10737:
10731:
10722:
10709:
10706:
10701:
10697:
10693:
10688:
10684:
10678:
10675:
10671:
10665:
10660:
10657:
10654:
10650:
10646:
10638:
10634:
10627:
10622:
10618:
10612:
10607:
10604:
10601:
10597:
10593:
10590:
10587:
10581:
10575:
10569:
10563:
10558:
10553:
10549:
10537:not equal to
10520:
10517:
10509:
10505:
10498:
10475:
10467:
10463:
10442:
10439:
10436:
10433:
10426:is less than
10410:
10404:
10382:
10379:
10375:
10371:
10363:
10359:
10352:
10329:
10321:
10317:
10285:
10280:
10274:
10270:
10266:
10261:
10257:
10252:
10245:
10240:
10234:
10230:
10226:
10223:
10219:
10206:
10203:
10200:
10193:
10190:
10187:
10184:
10181:
10173:
10169:
10163:
10157:
10138:
10135:
10132:
10129:
10115:
10100:
10095:
10091:
10087:
10081:
10078:
10075:
10071:
10065:
10062:
10035:
10031:
10022:
10019:
10016:
10012:
10004:vanishes, so
9986:
9982:
9973:
9969:
9943:
9939:
9930:
9927:
9924:
9920:
9913:
9907:
9899:
9895:
9886:
9882:
9875:
9869:
9861:
9857:
9848:
9845:
9842:
9838:
9814:
9806:
9803:
9800:
9796:
9772:
9764:
9760:
9737:
9733:
9724:
9722:
9721:
9699:
9695:
9687:
9683:
9679:
9670:
9669:
9663:
9650:
9647:
9642:
9634:
9626:
9623:
9620:
9616:
9609:
9603:
9598:
9593:
9589:
9577:
9573:
9564:
9561:
9558:
9554:
9548:
9545:
9542:
9538:
9531:
9527:
9521:
9518:
9515:
9507:
9503:
9499:
9496:
9488:
9480:
9476:
9466:
9460:
9455:
9450:
9446:
9422:
9414:
9411:
9408:
9404:
9383:
9380:
9377:
9368:
9351:
9348:
9345:
9341:
9333:
9325:
9322:
9319:
9315:
9308:
9304:
9296:
9293:
9290:
9286:
9278:
9270:
9267:
9264:
9260:
9253:
9248:
9245:
9242:
9238:
9233:
9229:
9224:
9221:
9218:
9214:
9204:
9191:
9188:
9183:
9180:
9177:
9173:
9166:
9158:
9155:
9152:
9148:
9141:
9135:
9130:
9125:
9121:
9109:
9105:
9096:
9093:
9090:
9086:
9079:
9075:
9069:
9066:
9063:
9055:
9051:
9047:
9044:
9036:
9028:
9024:
9014:
9008:
9003:
8998:
8994:
8970:
8962:
8959:
8956:
8952:
8941:
8937:
8920:
8917:
8914:
8904:
8900:
8880:
8874:
8871:
8866:
8863:
8860:
8856:
8850:
8846:
8842:
8834:
8830:
8826:
8823:
8815:
8807:
8803:
8789:
8767:
8763:
8759:
8756:
8748:
8740:
8736:
8709:
8701:
8698:
8695:
8691:
8687:
8681:
8675:
8666:
8653:
8650:
8642:
8638:
8634:
8631:
8623:
8615:
8611:
8604:
8598:
8589:
8583:
8578:
8573:
8569:
8557:
8553:
8546:
8542:
8537:
8534:
8531:
8523:
8519:
8515:
8512:
8504:
8496:
8492:
8482:
8476:
8471:
8466:
8462:
8451:
8447:
8427:
8419:
8415:
8404:
8382:
8378:
8374:
8371:
8364:
8359:
8352:
8348:
8344:
8339:
8334:
8331:
8308:
8305:
8302:
8293:
8280:
8277:
8269:
8265:
8261:
8258:
8250:
8242:
8238:
8228:
8222:
8217:
8212:
8208:
8202:
8197:
8193:
8189:
8186:
8183:
8175:
8171:
8167:
8164:
8156:
8148:
8144:
8138:
8134:
8124:
8118:
8113:
8108:
8104:
8094:
8076:
8072:
8068:
8065:
8061:
8054:
8049:
8042:
8038:
8034:
8029:
8024:
8016:
8012:
8008:
8005:
7998:
7993:
7986:
7982:
7978:
7973:
7968:
7965:
7959:
7951:
7947:
7943:
7940:
7936:
7921:
7914:
7912:
7898:
7895:
7887:
7883:
7879:
7876:
7868:
7860:
7856:
7846:
7840:
7835:
7830:
7826:
7814:
7810:
7802:
7798:
7794:
7789:
7784:
7779:
7775:
7767:
7766:
7763:
7760:
7743:
7739:
7729:
7725:
7717:
7713:
7709:
7702:
7698:
7692:
7688:
7684:
7679:
7675:
7670:
7660:
7657:
7654:
7647:
7644:
7641:
7638:
7635:
7627:
7604:
7600:
7572:
7564:
7560:
7537:
7533:
7510:
7506:
7496:
7479:
7473:
7469:
7465:
7462:
7458:
7452:
7448:
7439:
7431:
7427:
7420:
7412:
7408:
7404:
7401:
7395:
7389:
7385:
7381:
7378:
7374:
7368:
7365:
7362:
7359:
7356:
7352:
7345:
7341:
7335:
7331:
7327:
7324:
7320:
7310:
7307:
7304:
7297:
7294:
7291:
7288:
7285:
7277:
7269:We can write
7267:
7250:
7242:
7239:
7236:
7232:
7208:
7200:
7196:
7173:
7169:
7159:
7146:
7143:
7135:
7131:
7127:
7122:
7118:
7110:
7106:
7102:
7099:
7087:
7084:
7081:
7074:
7071:
7068:
7065:
7062:
7054:
7047:
7041:
7036:
7031:
7027:
7023:
7018:
7014:
6997:
6984:
6981:
6973:
6969:
6965:
6960:
6956:
6948:
6944:
6940:
6937:
6925:
6922:
6919:
6912:
6909:
6906:
6903:
6900:
6892:
6885:
6879:
6874:
6869:
6865:
6856:
6852:
6845:
6840:
6835:
6832:
6829:
6825:
6821:
6818:
6815:
6809:
6803:
6797:
6791:
6786:
6781:
6777:
6758:
6754:
6739:
6735:
6713:
6709:
6705:
6700:
6696:
6688:
6684:
6680:
6677:
6665:
6662:
6659:
6652:
6649:
6646:
6643:
6640:
6632:
6623:
6619:
6612:
6607:
6602:
6599:
6596:
6592:
6588:
6582:
6576:
6551:
6547:
6540:
6530:
6526:
6521:
6502:
6494:
6490:
6467:
6463:
6440:
6436:
6422:
6415:
6413:
6391:
6387:
6378:
6375:
6372:
6368:
6359:
6355:
6347:
6343:
6339:
6333:
6330:
6325:
6317:
6309:
6306:
6303:
6299:
6292:
6286:
6281:
6276:
6272:
6261:
6258:
6255:
6251:
6245:
6241:
6235:
6230:
6226:
6218:
6217:
6214:
6206:
6190:
6186:
6160:
6156:
6152:
6149:
6141:
6137:
6126:
6119:
6115:
6108:
6104:
6087:
6084:
6081:
6078:
6075:
6068:
6062:
6057:
6048:
6044:
6040:
6037:
6029:
6025:
6020:
6012:
6004:
6000:
5994:
5989:
5985:
5967:
5954:
5946:
5942:
5938:
5935:
5927:
5923:
5915:
5907:
5903:
5892:
5888:
5858:
5854:
5848:
5844:
5836:
5832:
5826:
5813:
5805:
5801:
5794:
5788:
5784:
5778:
5773:
5770:
5767:
5763:
5759:
5751:
5747:
5740:
5734:
5730:
5724:
5719:
5716:
5713:
5709:
5705:
5702:
5699:
5692:
5686:
5679:
5673:
5668:
5663:
5659:
5655:
5652:
5649:
5642:
5636:
5629:
5623:
5618:
5613:
5609:
5584:
5576:
5572:
5565:
5562:
5554:
5550:
5543:
5540:
5532:
5528:
5521:
5512:
5508:
5499:
5495:
5491:
5483:
5479:
5472:
5458:are roots of
5448:
5443:
5439:
5432:
5428:
5421:
5417:
5411:
5398:
5393:
5389:
5379:
5375:
5368:
5362:
5357:
5354:
5351:
5347:
5343:
5340:
5337:
5330:
5322:
5318:
5310:
5304:
5299:
5294:
5290:
5280:
5276:
5269:
5263:
5258:
5255:
5252:
5248:
5244:
5241:
5238:
5229:
5225:
5218:
5211:
5203:
5199:
5193:
5188:
5185:
5182:
5178:
5170:
5164:
5159:
5154:
5150:
5146:
5143:
5140:
5133:
5127:
5120:
5114:
5109:
5104:
5100:
5090:
5077:
5069:
5065:
5058:
5051:
5043:
5039:
5033:
5028:
5025:
5022:
5018:
5014:
5008:
5002:
4993:
4980:
4972:
4968:
4964:
4959:
4955:
4947:
4943:
4939:
4936:
4928:
4925:
4922:
4918:
4914:
4908:
4900:
4896:
4885:
4880:
4876:
4872:
4862:
4858:is of degree
4855:
4851:
4845:
4832:
4829:
4826:
4819:
4813:
4806:
4800:
4795:
4790:
4786:
4782:
4779:
4776:
4761:
4755:
4752:
4746:
4740:
4734:
4726:
4722:
4706:
4700:
4695:
4690:
4686:
4682:
4679:
4676:
4669:
4663:
4656:
4650:
4645:
4640:
4636:
4625:
4621:
4603:
4596:
4592:
4578:
4571:
4567:
4550:
4544:
4538:
4535:
4529:
4523:
4516:
4508:
4504:
4500:
4494:
4488:
4471:
4463:
4459:
4453:
4449:
4445:
4430:
4422:
4418:
4383:
4380:
4374:
4370:
4350:
4337:
4334:
4331:
4328:
4325:
4322:
4319:
4316:
4313:
4310:
4307:
4304:for all
4298:
4295:
4292:
4289:
4286:
4279:
4271:
4267:
4261:
4257:
4249:
4243:
4238:
4233:
4229:
4199:
4196:
4193:
4191:
4188:
4170:
4166:
4162:
4158:
4150:
4145:
4144:
4141:
4138:
4136:
4134:
4130:
4116:
4113:
4110:
4107:
4103:
4098:
4095:
4091:
4085:
4081:
4073:
4068:
4067:
4064:
4061:
4058:
4056:
4053:
4036:
4033:
4029:
4021:
4016:
4015:
4012:
4009:
4006:
4003:
3985:
3981:
3977:
3974:
3965:
3960:
3959:
3956:
3953:
3950:
3947:
3944:
3925:
3921:
3917:
3914:
3910:
3901:
3896:
3895:
3892:
3889:
3884:
3879:
3877:
3874:
3860:
3857:
3854:
3851:
3848:
3845:
3841:
3836:
3831:
3827:
3824:
3821:
3817:
3810:
3805:
3801:
3798:
3795:
3791:
3782:
3777:
3776:
3773:
3770:
3767:
3765:
3762:
3757:
3752:
3751:
3747:
3744:
3741:
3737:
3733:
3729:
3726:
3725:
3722:
3721:(A & S).
3720:
3714:
3710:
3703:
3696:
3667:
3664:
3657:
3651:
3644:
3638:
3633:
3628:
3624:
3612:
3602:
3599:
3583:
3580:
3574:
3571:
3567:
3562:
3558:
3555:
3552:
3546:
3542:
3538:
3533:
3529:
3525:
3503:
3498:
3494:
3490:
3480:
3463:
3460:
3458:
3447:
3441:
3438:
3434:
3425:
3419:
3416:
3413:
3409:
3404:
3400:
3397:
3394:
3391:
3389:
3381:
3375:
3348:
3345:
3343:
3332:
3326:
3323:
3319:
3310:
3304:
3301:
3298:
3294:
3289:
3285:
3282:
3279:
3276:
3274:
3266:
3260:
3232:
3229:
3227:
3218:
3211:
3205:
3199:
3192:
3188:
3185:
3183:
3174:
3167:
3161:
3158:
3152:
3146:
3142:
3138:
3135:
3133:
3124:
3119:
3115:
3112:
3106:
3100:
3096:
3092:
3089:
3085:
3081:
3078:
3072:
3069:
3063:
3059:
3055:
3051:
3047:
3044:
3042:
3033:
3028:
3024:
3021:
3016:
3012:
3008:
3004:
3000:
2995:
2991:
2987:
2983:
2979:
2976:
2971:
2967:
2963:
2959:
2955:
2950:
2946:
2941:
2937:
2934:
2932:
2925:
2922:
2917:
2913:
2910:
2907:
2904:
2900:
2896:
2890:
2885:
2882:
2878:
2874:
2847:
2844:
2839:
2835:
2827:
2813:
2810:
2805:
2801:
2793:
2779:
2776:
2773:
2768:
2764:
2756:
2742:
2739:
2734:
2730:
2722:
2721:
2720:
2717:
2698:
2695:
2690:
2686:
2683:
2680:
2677:
2673:
2669:
2663:
2658:
2655:
2651:
2647:
2644:
2642:
2632:
2629:
2624:
2618:
2614:
2611:
2608:
2602:
2599:
2594:
2590:
2587:
2584:
2577:
2573:
2567:
2562:
2559:
2555:
2549:
2545:
2542:
2539:
2533:
2531:
2525:
2522:
2516:
2510:
2504:
2499:
2495:
2482:
2467:
2463:
2460:
2457:
2454:
2450:
2428:
2424:
2421:
2418:
2414:
2404:
2401:
2398:
2383:
2380:
2375:
2371:
2368:
2365:
2361:
2355:
2352:
2349:
2346:
2342:
2337:
2333:
2330:
2327:
2323:
2316:
2311:
2307:
2303:
2300:
2280:
2272:
2269:
2266:
2241:
2238:
2235:
2221:
2208:
2204:
2198:
2194:
2191:
2188:
2182:
2177:
2173:
2167:
2163:
2160:
2157:
2150:
2146:
2141:
2137:
2131:
2126:
2123:
2120:
2116:
2110:
2106:
2103:
2100:
2094:
2091:
2088:
2081:
2075:
2070:
2065:
2061:
2037:
2034:
2031:
2008:
2001:Applying the
1999:
1984:
1980:
1977:
1974:
1968:
1962:
1959:
1954:
1951:
1939:
1926:
1923:
1917:
1914:
1909:
1906:
1898:
1892:
1888:
1885:
1882:
1876:
1873:
1868:
1864:
1861:
1858:
1851:
1847:
1842:
1837:
1834:
1830:
1826:
1823:
1820:
1813:
1807:
1802:
1797:
1793:
1773:
1757:
1751:
1746:
1743:
1740:
1730:
1727:
1708:
1705:
1699:
1696:
1693:
1686:
1683:
1678:
1671:
1670:
1666:
1650:
1644:
1639:
1636:
1633:
1623:
1621:±0.538469...
1620:
1601:
1598:
1592:
1589:
1586:
1579:
1576:
1571:
1564:
1563:
1559:
1543:
1540:
1531:
1524:
1520:
1504:
1498:
1493:
1490:
1480:
1478:±0.861136...
1477:
1458:
1455:
1447:
1444:
1439:
1434:
1431:
1424:
1417:
1416:
1412:
1396:
1390:
1385:
1382:
1372:
1370:±0.339981...
1369:
1350:
1347:
1339:
1336:
1331:
1326:
1323:
1316:
1309:
1305:
1301:
1285:
1282:
1273:
1271:±0.774597...
1270:
1253:
1250:
1244:
1237:
1236:
1232:
1216:
1213:
1204:
1197:
1190:
1173:
1169:
1164:
1157:
1154:
1153:
1143:
1142:
1117:
1116:
1113:
1110:
1097:
1089:
1084:
1075:
1071:
1063:
1059:
1055:
1050:
1044:
1038:
1033:
1029:
1025:
1022:
1018:
1013:
1008:
1003:
999:
965:
961:
950:
945:
941:
937:, denoted by
936:
917:
914:
911:
908:
895:
891:
881:
875:
870:
860:
858:
854:
842:
840:
836:
816:
812:
808:
805:
795:
774:
770:
766:
763:
759:
748:
735:
731:
727:
723:
719:
715:
711:
707:
703:
699:
693:
688:
685:
682:
678:
674:
671:
668:
661:
655:
650:
645:
641:
638:
635:
631:
624:
619:
615:
612:
609:
605:
598:
593:
590:
586:
582:
579:
576:
569:
563:
558:
553:
550:
546:
537:
536:rules, i.e.,
535:
513:
509:
503:
490:
487:
484:
481:
478:
475:
472:
468:
462:
456:
451:
446:
442:
439:
436:
432:
425:
420:
416:
413:
410:
406:
401:
395:
389:
381:
379:
378:singularities
375:
370:
369:or less on .
366:
358:
354:
349:
343:
336:
323:
315:
311:
304:
299:
295:
289:
284:
281:
278:
274:
270:
267:
264:
257:
251:
246:
241:
238:
234:
225:
221:
216:
213:
209:
188:
182:
178:
174:
170:
162:
152:
148:
144:
140:
136:
129:
111:
107:
102:
97:
90:
80:
77:
69:
59:
55:
49:
48:
42:
37:
28:
27:
22:
18877:
18863:Boole's rule
18761:
18757:
18733:
18719:
18673:
18637:
18633:
18604:
18574:
18549:
18545:
18518:
18514:
18485:. New York:
18481:
18459:
18428:J. Eng. Math
18427:
18393:
18389:
18371:
18367:
18348:
18344:
18315:
18310:Moler, Cleve
18295:(1): 11–19.
18292:
18288:
18274:cite journal
18265:
18261:
18227:
18223:
18198:
18179:
18154:
18129:
18114:
18093:
18089:
18075:28 September
18073:. Retrieved
18044:
17878:
17874:
17758:
17742:Bibliography
17729:
17718:
17707:
17696:
17684:
17672:
17660:
17648:
17636:
17625:
17614:
17581:
16494:
16486:
16240:
16107:
15980:
15775:
15771:
15768:
15756:
15739:
15737:
15725:
15713:
15708:
15705:
15687:
15683:
15490:
15486:
15483:
15366:
15362:
15149:
15145:
15133:
15038:
14889:
14887:
14591:
14588:
14256:
14251:
14218:
13629:
13426:
13095:
12969:
12821:
12817:
12681:
12647:
12540:
11962:
11404:
11292:
10927:
10920:operations.
10915:
10911:
10906:
10900:
10896:
10843:
10839:
10830:and weights
10822:
10723:
10121:
9725:
9718:
9666:
9664:
9369:
9205:
8939:
8935:
8902:
8898:
8787:
8667:
8449:
8445:
8402:
8294:
8095:
7926:
7915:
7761:
7497:
7268:
7266:as follows.
7160:
6999:The weights
6998:
6756:
6752:
6737:
6733:
6533:in terms of
6528:
6524:
6427:
6416:
6212:
6188:
6184:
6182:. Therefore
6124:
6117:
6113:
6106:
6102:
5968:
5890:
5886:
5856:
5852:
5849:
5842:
5834:
5830:
5827:
5446:
5441:
5437:
5430:
5426:
5419:
5415:
5412:
5091:
4994:
4883:
4878:
4874:
4860:
4853:
4849:
4846:
4628:. Therefore
4623:
4619:
4601:
4594:
4590:
4576:
4569:
4565:
4469:
4461:
4457:
4454:
4447:
4443:
4428:
4420:
4416:
4384:
4378:
4372:
4368:
4351:
4209:
4131:Generalized
3948:(first kind)
3882:
3735:
3731:
3712:
3708:
3701:
3694:
3608:
3600:
3481:
2862:
2718:
2483:
2405:
2402:
2399:
2293:as given by
2227:
2000:
1940:
1784:
1774:0.236927...
1728:±0.90618...
1667:0.478629...
1560:0.568889...
1521:0.347855...
1413:0.652145...
1302:0.555556...
1233:0.888889...
1191:±0.57735...
1111:
984:-th root of
963:
959:
948:
943:
939:
893:
889:
886:
879:
843:
749:
538:
524:and weights
511:
507:
504:
382:
371:
364:
356:
352:
341:
337:
217:
211:
207:
198:and weights
186:
168:
158:
150:
146:
142:
138:
134:
127:
109:
105:
72:
63:
44:
18760:=2 through
18699:version of
18695:— includes
18351:: 767–793.
17630:Jacobi 1826
16241:Remainder:
16020:st zero of
15981:Abscissas:
10724:Since both
10345:. Clearly
3605:Other forms
2848:0.577350269
2814:1.000000000
2780:0.577350269
2743:1.000000000
953:. With the
181:polynomials
58:introducing
18528:1607.08404
18155:Math. Comp
18130:Math. Comp
18045:Math. Comp
17875:Math. Comp
17619:Gauss 1815
17602:References
17584:GNU Octave
15386:and where
14219:The zeros
10909:requiring
10890:using the
10541:, we have
8792:, we have
8321:, because
5839:of degree
4425:of degree
4221:such that
4070:[0, ∞)
4018:[0, ∞)
3745:A & S
372:The Gauss–
210:= 1, ...,
183:of degree
41:references
18735:MathWorld
18680:EMS Press
18521:: 22–54.
18452:Press, WH
17992:∞
17983:∫
17942:α
17928:
17922:−
17905:∫
17832:
17826:−
17812:∫
17607:Citations
17596:integrate
17534:±
17495:−
17429:±
17339:−
17324:±
17198:±
17162:−
17099:±
17015:−
17000:±
16925:±
16862:±
16736:±
16673:±
16598:±
16508:Weights,
16467:ξ
16458:−
16448:ξ
16434:−
16392:−
16366:−
16330:−
16308:−
16280:−
16265:−
16224:±
16221:≠
16167:−
16144:−
16108:Weights:
16036:−
16002:−
15915:−
15897:∑
15881:−
15845:−
15793:−
15789:∫
15661:ξ
15645:ξ
15515:−
15427:ω
15413:∫
15355:for some
15290:ξ
15212:∑
15208:−
15176:ω
15162:∫
15102:ω
15088:∫
15075:μ
15048:μ
14998:ϕ
14983:μ
14930:ϕ
14858:−
14838:−
14818:⋯
14799:−
14779:−
14766:⋱
14761:⋱
14756:⋮
14744:⋱
14739:⋱
14711:⋮
14706:⋱
14650:⋯
14589:That is,
14562:…
14519:−
14486:−
14470:−
14443:−
14408:−
14383:…
14351:−
14259:symmetric
14189:−
14171:−
14153:⋯
14131:−
14118:⋱
14113:⋱
14108:⋮
14096:⋱
14091:⋱
14067:⋮
14062:⋱
14021:⋯
13929:⋯
13801:−
13788:⋯
13728:~
13671:−
13665:~
13647:~
13574:−
13503:−
13462:−
13402:−
13383:−
13321:−
13302:−
13278:−
13126:≡
13109:−
13068:−
13047:−
13012:−
12942:−
12926:−
12909:−
12868:−
12614:−
12489:−
12406:−
12244:−
12241:⋯
12217:−
12198:−
12181:−
12133:−
11997:…
11911:−
11828:−
11742:−
11440:−
11431:…
11352:ω
11338:∫
11244:−
11241:⋯
11224:−
11208:−
11191:−
11150:−
11065:, degree
11050:⋅
11044:⋅
11018:≠
10732:ω
10672:δ
10651:∑
10598:∑
10564:ω
10550:∫
10440:−
10376:δ
10267:−
10227:−
10204:≠
10191:≤
10185:≤
10174:∏
10136:−
10020:−
9928:−
9624:−
9604:ω
9590:∫
9562:−
9546:−
9500:−
9461:ω
9447:∫
9412:−
9381:−
9349:−
9323:−
9294:−
9268:−
9254:−
9246:−
9222:−
9181:−
9156:−
9136:ω
9122:∫
9094:−
9048:−
9009:ω
8995:∫
8960:−
8918:−
8864:−
8827:−
8760:−
8699:−
8635:−
8584:ω
8570:∫
8516:−
8477:ω
8463:∫
8443:. So, if
8375:−
8335:−
8306:≤
8295:provided
8262:−
8223:ω
8209:∫
8168:−
8119:ω
8105:∫
8069:−
8009:−
7969:−
7944:−
7880:−
7841:ω
7827:∫
7685:−
7658:≠
7645:≤
7639:≤
7628:∏
7466:−
7405:−
7382:−
7366:≤
7360:≤
7353:∏
7328:−
7308:≠
7295:≤
7289:≤
7278:∏
7240:−
7128:−
7103:−
7085:≠
7072:≤
7066:≤
7055:∏
7042:ω
7028:∫
6966:−
6941:−
6923:≠
6910:≤
6904:≤
6893:∏
6880:ω
6866:∫
6826:∑
6792:ω
6778:∫
6706:−
6681:−
6663:≠
6650:≤
6644:≤
6633:∏
6593:∑
6376:−
6307:−
6287:ω
6273:∫
6259:−
6153:−
6138:∏
6082:≠
6063:ω
6041:−
6026:∏
5986:∫
5939:−
5924:∏
5764:∑
5710:∑
5674:ω
5660:∫
5624:ω
5610:∫
5348:∑
5305:ω
5291:∫
5249:∑
5179:∑
5165:ω
5151:∫
5115:ω
5101:∫
5019:∑
4965:−
4940:−
4926:≠
4919:∏
4801:ω
4787:∫
4701:ω
4687:∫
4651:ω
4637:∫
4335:−
4326:…
4244:ω
4230:∫
4163:−
4114:−
4108:α
4096:−
4086:α
4034:−
3978:−
3918:−
3880:25.4.33 (
3858:−
3852:β
3846:α
3837:β
3811:α
3799:−
3727:Interval
3639:ω
3625:∫
3587:%
3578:%
3572:×
3556:−
3530:ε
3495:ε
3439:−
3417:−
3401:
3324:−
3302:−
3286:
3107:0.5773503
3073:0.5773503
3070:−
2935:≈
2883:−
2879:∫
2777:−
2656:−
2652:∫
2588:−
2560:−
2556:∫
2543:−
2496:∫
2455:−
2403:Solution
2366:−
2350:−
2334:
2308:∫
2174:ξ
2161:−
2117:∑
2104:−
2095:≈
2062:∫
2032:ξ
1978:−
1963:ξ
1927:ξ
1918:ξ
1874:ξ
1862:−
1835:−
1831:∫
1794:∫
1744:−
1679:±
1590:−
1572:±
1494:−
1425:±
1332:−
1317:±
1245:±
1165:±
1133:Weights,
1026:−
980:, is the
909:−
809:−
767:−
679:∑
675:≈
651:β
625:α
613:−
591:−
587:∫
551:−
547:∫
485:−
479:β
473:α
452:β
426:α
414:−
275:∑
271:≈
239:−
235:∫
128:composite
112:(1) = –10
18967:Category
18787:Archived
18701:QUADPACK
18579:Springer
18256:(1826).
18201:, SIAM,
18182:. SIAM.
18112:(1815).
17794:65-12253
17778:64-60036
16499:Points,
16043:′
15761:, where
15371:, where
13555:, since
13141:) where
11696:one has
9788:and now
9688:′
8986:we have
8933:. Since
7803:′
7718:′
6731:because
6348:′
5597:for all
4995:We have
4888:, where
4407:weights
3563:11061.34
3559:11058.44
3553:11061.34
3464:708.4811
3448:25.35085
3426:25.35085
3382:25.35085
3349:296.8317
3333:12.64915
3311:12.64915
3267:12.64915
3233:11058.44
3212:708.4811
3200:296.8317
3168:25.35085
3153:12.64915
1124:Points,
1064:′
851:are the
728:′
708:′
374:Legendre
18682:, 2001
18642:Bibcode
18625:2723248
18554:Bibcode
18432:Bibcode
18398:Bibcode
18246:2004418
18171:0285177
18146:0228171
18061:0331730
17786:0167642
16065:, here
15988:is the
13872:is the
11963:Now if
8096:yields
6769:yields
4481:to get
4194:25.4.46
4147:(−∞, ∞)
4059:25.4.45
4007:25.4.40
3951:25.4.38
3898:(−1, 1)
3779:(−1, 1)
3768:25.4.29
967:(1) = 1
175:, is a
167:-point
120:⁄
108:(–1) +
54:improve
18687:ALGLIB
18623:
18613:
18585:
18501:
18466:
18328:
18244:
18205:
18186:
18169:
18144:
18059:
17792:
17784:
17776:
17766:
17679:, §5.2
17588:MATLAB
15039:where
13967:, and
13714:where
11455:where
8896:where
7498:where
6428:where
5413:where
4563:where
4389:nodes
4376:, and
3706:, and
3692:. For
3688:, and
3584:0.0262
3414:140000
3410:140000
3299:140000
3295:140000
3249:since
2481:gives
2347:140000
2343:140000
969:, the
796:) and
505:where
43:, but
18920:Other
18599:, in
18523:arXiv
18242:JSTOR
17592:quadl
15493:) = 1
12820:<
11106:monic
10955:with
3715:) = 1
1941:with
853:roots
163:, an
141:) = 7
18611:ISBN
18583:ISBN
18499:ISBN
18464:ISBN
18326:ISBN
18280:link
18203:ISBN
18184:ISBN
18077:2018
17974:and
17790:LCCN
17774:LCCN
17764:ISBN
17594:and
17586:and
16470:<
16464:<
15664:<
15658:<
14917:are
14893:and
12815:for
12776:and
12608:<
11405:for
11007:for
10800:>
10753:and
10707:>
10533:for
7224:and
6193:has
4371:<
4366:for
4210:Let
4111:>
3855:>
3697:= −1
3420:2100
3395:2000
3305:2100
3280:2000
2353:2100
2328:2000
882:= 5)
482:>
205:for
18718:at
18650:doi
18638:129
18562:doi
18533:doi
18491:doi
18440:doi
18406:doi
18394:127
18376:doi
18372:112
18353:doi
18297:doi
18232:doi
18159:doi
18134:doi
18098:doi
18049:doi
18043:".
17925:log
17883:doi
17590:as
17509:350
17492:124
17404:350
17387:124
17299:525
17296:256
15759:– 3
15728:+ 1
15716:+ 1
15359:in
13622:).
12970:or
12824:− 1
12684:− 1
12650:− 1
9723:).
8790:− 1
8405:− 1
7592:to
7552:in
6765:to
6569:as
6482:in
6201:to
5973:to
5880:to
5865:to
5845:− 1
4863:− 1
4604:− 1
4579:− 1
4472:− 1
4431:− 1
3885:= 0
3704:= 1
3575:100
3517:is
3442:9.8
3327:9.8
2442:to
2369:9.8
2259:to
1758:900
1741:322
1651:900
1634:322
1544:225
1541:128
841:).
367:− 1
344:− 1
215:.
189:− 1
159:In
153:+ 3
149:– 3
145:– 8
18969::
18732:.
18678:,
18672:,
18648:.
18636:.
18621:MR
18619:,
18581:,
18560:.
18550:66
18548:.
18531:.
18519:45
18517:.
18497:.
18458:,
18438:.
18404:.
18392:.
18370:,
18349:92
18347:.
18343:.
18324:.
18320:.
18293:11
18291:.
18276:}}
18272:{{
18264:.
18260:.
18240:.
18228:23
18226:.
18222:.
18167:MR
18165:.
18142:MR
18140:.
18094:40
18092:.
18088:.
18057:MR
18055:.
17879:19
17877:.
17873:.
17829:ln
17788:.
17782:MR
17780:.
17772:.
17752:;
17598:.
17563:21
17503:15
17452:11
17439:11
17398:15
17347:11
17334:11
17227:15
17173:30
17159:14
17131:21
17074:30
17060:14
17032:21
16954:10
16900:90
16897:49
16837:45
16834:32
16473:1.
16227:1.
15365:,
14254:.
13843:,
13245::=
13159::=
12527:0.
11949:0.
10815:.
10710:0.
10114:.
6205:.
6127:-1
4452:.
4446:,
4338:1.
3699:,
3684:,
3398:ln
3283:ln
3189:11
3139:11
3116:19
3101:11
3082:19
3064:11
3048:11
3025:19
3009:11
2980:19
2964:11
2938:11
2914:19
2905:11
2875:11
2687:19
2678:11
2648:11
2609:30
2585:30
2540:30
2505:30
2425:30
2331:ln
2317:30
2273:30
1752:70
1747:13
1706:10
1645:70
1640:13
1599:10
1529:0
1526:5
1505:36
1499:30
1491:18
1397:36
1391:30
1383:18
1307:4
1202:0
1199:3
1194:1
1155:2
1150:2
1147:0
1144:1
18819:e
18812:t
18805:v
18762:n
18758:n
18750:.
18738:.
18707:)
18697:C
18656:.
18652::
18644::
18628:.
18568:.
18564::
18556::
18539:.
18535::
18525::
18507:.
18493::
18446:.
18442::
18434::
18421:.
18412:.
18408::
18400::
18378::
18361:.
18355::
18334:.
18303:.
18299::
18282:)
18266:1
18248:.
18234::
18192:.
18173:.
18161::
18148:.
18136::
18104:.
18100::
18079:.
18063:.
18051::
18031:x
18028:d
18025:)
18022:x
18019:(
18016:f
18013:)
18010:x
18007:(
18002:m
17998:E
17987:0
17962:x
17959:d
17956:)
17953:x
17950:(
17947:f
17938:x
17934:)
17931:x
17919:(
17914:1
17909:0
17891:.
17885::
17871:"
17859:x
17856:d
17853:)
17850:x
17847:(
17844:f
17841:)
17838:x
17835:(
17821:1
17816:0
17796:.
17560:1
17537:1
17498:7
17463:3
17460:5
17449:2
17444:+
17436:5
17393:7
17390:+
17358:3
17355:5
17344:2
17331:5
17273:0
17252:7
17224:1
17201:1
17167:7
17125:7
17120:2
17114:+
17109:3
17106:1
17068:7
17063:+
17026:7
17021:2
17010:3
17007:1
16979:6
16951:1
16928:1
16871:7
16868:3
16811:0
16790:5
16765:6
16762:1
16739:1
16711:6
16708:5
16682:5
16679:1
16652:4
16627:3
16624:1
16601:1
16573:3
16570:4
16547:0
16526:3
16512:i
16510:w
16503:i
16501:x
16495:n
16461:1
16454:,
16451:)
16445:(
16440:)
16437:2
16431:n
16428:2
16425:(
16421:f
16412:3
16407:]
16403:!
16399:)
16395:2
16389:n
16386:2
16382:(
16377:[
16372:)
16369:1
16363:n
16360:2
16357:(
16350:4
16345:]
16341:!
16337:)
16333:2
16327:n
16323:(
16318:[
16311:1
16305:n
16302:2
16298:2
16292:3
16287:)
16283:1
16277:n
16273:(
16268:n
16259:=
16254:n
16250:R
16216:i
16212:x
16207:,
16199:2
16194:]
16189:)
16184:i
16180:x
16176:(
16170:1
16164:n
16160:P
16155:[
16150:)
16147:1
16141:n
16138:(
16135:n
16131:2
16126:=
16121:i
16117:w
16103:m
16089:)
16086:x
16083:(
16078:m
16074:P
16053:)
16050:x
16047:(
16039:1
16033:n
16029:P
16008:)
16005:1
15999:i
15996:(
15985:i
15983:x
15967:.
15962:n
15958:R
15954:+
15950:)
15945:i
15941:x
15937:(
15934:f
15929:i
15925:w
15918:1
15912:n
15907:2
15904:=
15901:i
15893:+
15890:]
15887:)
15884:1
15878:(
15875:f
15872:+
15869:)
15866:1
15863:(
15860:f
15857:[
15851:)
15848:1
15842:n
15839:(
15836:n
15832:2
15827:=
15823:x
15820:d
15816:)
15813:x
15810:(
15807:f
15801:1
15796:1
15778:)
15776:x
15774:(
15772:f
15763:n
15757:n
15755:2
15726:n
15724:2
15720:n
15714:n
15688:n
15686:2
15670:.
15667:b
15655:a
15651:,
15648:)
15642:(
15637:)
15634:n
15631:2
15628:(
15624:f
15615:3
15610:]
15606:!
15602:)
15598:n
15595:2
15591:(
15586:[
15581:)
15578:1
15575:+
15572:n
15569:2
15566:(
15559:4
15554:)
15550:!
15547:n
15543:(
15536:1
15533:+
15530:n
15527:2
15522:)
15518:a
15512:b
15508:(
15491:x
15489:(
15487:ω
15470:.
15467:x
15464:d
15460:)
15457:x
15454:(
15451:g
15448:)
15445:x
15442:(
15439:f
15436:)
15433:x
15430:(
15422:b
15417:a
15409:=
15406:)
15403:g
15400:,
15397:f
15394:(
15384:n
15380:1
15375:n
15373:p
15369:)
15367:b
15363:a
15361:(
15357:ξ
15343:)
15338:n
15334:p
15330:,
15325:n
15321:p
15317:(
15310:!
15307:)
15304:n
15301:2
15298:(
15293:)
15287:(
15282:)
15279:n
15276:2
15273:(
15269:f
15262:=
15259:)
15254:i
15250:x
15246:(
15243:f
15237:i
15233:w
15227:n
15222:1
15219:=
15216:i
15205:x
15202:d
15198:)
15195:x
15192:(
15189:f
15185:)
15182:x
15179:(
15171:b
15166:a
15150:n
15148:2
15120:.
15117:x
15114:d
15111:)
15108:x
15105:(
15097:b
15092:a
15084:=
15079:0
15052:0
15023:2
15018:)
15013:)
15010:j
15007:(
15002:1
14994:(
14987:0
14979:=
14974:j
14970:w
14958:j
14956:x
14940:)
14937:j
14934:(
14903:J
14890:J
14874:.
14869:]
14861:1
14855:n
14851:a
14841:1
14835:n
14831:b
14823:0
14813:0
14802:1
14796:n
14792:b
14782:2
14776:n
14772:a
14749:0
14730:2
14726:b
14718:0
14697:2
14693:b
14683:1
14679:a
14669:1
14665:b
14655:0
14645:0
14636:1
14632:b
14622:0
14618:a
14611:[
14606:=
14601:J
14571:.
14568:n
14565:,
14559:,
14556:2
14548:,
14545:1
14538:=
14531:k
14522:1
14516:k
14512:b
14506:=
14495:k
14492:,
14489:1
14483:k
14479:J
14473:1
14467:k
14464:,
14461:k
14457:J
14451:=
14446:1
14440:k
14437:,
14434:k
14428:J
14422:=
14417:i
14414:,
14411:1
14405:k
14399:J
14389:n
14386:,
14380:,
14377:2
14374:,
14371:1
14368:=
14361:k
14354:1
14348:k
14344:a
14340:=
14331:i
14328:,
14325:k
14321:J
14317:=
14312:i
14309:,
14306:k
14300:J
14271:J
14248:n
14232:j
14228:x
14205:.
14200:]
14192:1
14186:n
14182:a
14174:1
14168:n
14164:b
14158:0
14148:0
14141:1
14134:2
14128:n
14124:a
14101:0
14084:2
14080:b
14074:0
14057:1
14050:1
14046:a
14038:1
14034:b
14026:0
14016:0
14011:1
14004:0
14000:a
13993:[
13988:=
13984:J
13969:J
13952:T
13945:]
13939:1
13934:0
13924:0
13918:[
13912:=
13907:n
13902:e
13880:n
13858:n
13853:e
13828:T
13821:]
13815:)
13812:x
13809:(
13804:1
13798:n
13794:p
13783:)
13780:x
13777:(
13772:1
13768:p
13762:)
13759:x
13756:(
13751:0
13747:p
13740:[
13734:=
13725:P
13700:n
13695:e
13690:)
13687:x
13684:(
13679:n
13675:p
13662:P
13656:x
13653:=
13644:P
13638:J
13620:r
13604:r
13600:p
13577:1
13571:r
13567:p
13563:x
13543:)
13538:r
13534:p
13530:,
13525:r
13521:p
13517:(
13514:=
13511:)
13506:1
13500:r
13496:p
13492:x
13489:,
13484:r
13480:p
13476:(
13473:=
13470:)
13465:1
13459:r
13455:p
13451:,
13446:r
13442:p
13438:x
13435:(
13410:)
13405:1
13399:r
13395:p
13391:,
13386:1
13380:r
13376:p
13372:(
13367:)
13362:r
13358:p
13354:,
13349:r
13345:p
13341:(
13335:=
13329:)
13324:1
13318:r
13314:p
13310:,
13305:1
13299:r
13295:p
13291:(
13286:)
13281:1
13275:r
13271:p
13267:,
13262:r
13258:p
13254:x
13251:(
13240:r
13236:b
13231:,
13225:)
13220:r
13216:p
13212:,
13207:r
13203:p
13199:(
13194:)
13189:r
13185:p
13181:,
13176:r
13172:p
13168:x
13165:(
13154:r
13150:a
13129:0
13123:)
13120:x
13117:(
13112:1
13105:p
13082:)
13079:x
13076:(
13071:1
13065:r
13061:p
13055:r
13051:b
13044:)
13041:x
13038:(
13033:r
13029:p
13025:)
13020:r
13016:a
13009:x
13006:(
13003:=
13000:)
12997:x
12994:(
12989:1
12986:+
12983:r
12979:p
12956:)
12953:x
12950:(
12945:1
12939:r
12935:p
12929:1
12923:r
12920:,
12917:r
12913:a
12906:)
12903:x
12900:(
12895:r
12891:p
12887:)
12882:r
12879:,
12876:r
12872:a
12865:x
12862:(
12859:=
12856:)
12853:x
12850:(
12845:1
12842:+
12839:r
12835:p
12822:r
12818:s
12803:0
12800:=
12795:s
12792:,
12789:r
12785:a
12764:0
12761:=
12758:)
12753:s
12749:p
12745:x
12742:,
12737:r
12733:p
12729:(
12726:=
12723:)
12718:s
12714:p
12710:,
12705:r
12701:p
12697:x
12694:(
12682:r
12665:r
12661:p
12648:r
12631:s
12627:p
12623:x
12620:,
12617:1
12611:r
12605:s
12585:)
12582:g
12579:x
12576:,
12573:f
12570:(
12567:=
12564:)
12561:g
12558:,
12555:f
12552:x
12549:(
12524:=
12521:)
12516:s
12512:p
12508:,
12503:r
12499:p
12495:x
12492:(
12486:)
12481:s
12477:p
12473:,
12468:r
12464:p
12460:x
12457:(
12454:=
12451:)
12446:s
12442:p
12438:,
12433:s
12429:p
12425:(
12420:s
12417:,
12414:r
12410:a
12403:)
12398:s
12394:p
12390:,
12385:r
12381:p
12377:x
12374:(
12371:=
12368:)
12363:s
12359:p
12355:,
12350:1
12347:+
12344:r
12340:p
12336:(
12314:s
12310:p
12289:)
12284:s
12280:p
12276:,
12271:0
12267:p
12263:(
12258:0
12255:,
12252:r
12248:a
12238:)
12233:s
12229:p
12225:,
12220:1
12214:r
12210:p
12206:(
12201:1
12195:r
12192:,
12189:r
12185:a
12178:)
12173:s
12169:p
12165:,
12160:r
12156:p
12152:(
12147:r
12144:,
12141:r
12137:a
12130:)
12125:s
12121:p
12117:,
12112:r
12108:p
12104:x
12101:(
12098:=
12095:)
12090:s
12086:p
12082:,
12077:1
12074:+
12071:r
12067:p
12063:(
12041:1
12038:+
12035:r
12031:p
12008:r
12004:p
12000:,
11994:,
11989:1
11985:p
11981:,
11976:0
11972:p
11946:=
11943:)
11938:0
11934:p
11930:,
11925:0
11921:p
11917:x
11914:(
11908:)
11903:0
11899:p
11895:,
11890:0
11886:p
11882:x
11879:(
11876:=
11873:)
11868:0
11864:p
11860:,
11855:0
11851:p
11847:(
11842:0
11839:,
11836:0
11832:a
11825:)
11820:0
11816:p
11812:,
11807:0
11803:p
11799:x
11796:(
11793:=
11790:)
11785:0
11781:p
11777:,
11772:0
11768:p
11764:(
11761:)
11756:0
11753:,
11750:0
11746:a
11739:x
11736:(
11733:=
11730:)
11725:0
11721:p
11717:,
11712:1
11708:p
11704:(
11684:0
11681:=
11678:s
11675:=
11672:r
11650:r
11646:p
11623:0
11619:p
11598:1
11595:=
11592:)
11589:x
11586:(
11581:0
11577:p
11553:)
11547:s
11543:p
11539:,
11534:s
11530:p
11525:(
11520:)
11514:s
11510:p
11506:,
11501:r
11497:p
11493:x
11489:(
11483:=
11478:s
11475:,
11472:r
11468:a
11457:n
11443:1
11437:n
11434:,
11428:,
11425:1
11422:,
11419:0
11416:=
11413:r
11391:x
11388:d
11385:)
11382:x
11379:(
11376:g
11373:)
11370:x
11367:(
11364:f
11361:)
11358:x
11355:(
11347:b
11342:a
11334:=
11331:)
11328:)
11325:x
11322:(
11319:g
11316:,
11313:)
11310:x
11307:(
11304:f
11301:(
11279:)
11276:x
11273:(
11268:0
11264:p
11258:0
11255:,
11252:r
11248:a
11238:)
11235:x
11232:(
11227:1
11221:r
11217:p
11211:1
11205:r
11202:,
11199:r
11195:a
11188:)
11185:x
11182:(
11177:r
11173:p
11169:)
11164:r
11161:,
11158:r
11154:a
11147:x
11144:(
11141:=
11138:)
11135:x
11132:(
11127:1
11124:+
11121:r
11117:p
11092:r
11089:=
11086:)
11081:r
11077:p
11073:(
11053:)
11047:,
11041:(
11021:s
11015:r
10995:0
10992:=
10989:)
10984:s
10980:p
10976:,
10971:r
10967:p
10963:(
10941:r
10937:p
10918:)
10916:n
10914:(
10912:O
10907:n
10903:)
10901:n
10899:(
10897:O
10878:0
10875:=
10872:)
10869:x
10866:(
10861:n
10857:p
10846:)
10844:n
10842:(
10840:O
10834:i
10832:w
10827:i
10825:x
10803:0
10795:i
10791:w
10770:)
10767:x
10764:(
10761:f
10741:)
10738:x
10735:(
10702:i
10698:w
10694:=
10689:j
10685:w
10679:j
10676:i
10666:n
10661:1
10658:=
10655:j
10647:=
10644:)
10639:j
10635:x
10631:(
10628:f
10623:j
10619:w
10613:n
10608:1
10605:=
10602:j
10594:=
10591:x
10588:d
10585:)
10582:x
10579:(
10576:f
10573:)
10570:x
10567:(
10559:b
10554:a
10539:i
10535:j
10521:0
10518:=
10515:)
10510:j
10506:x
10502:(
10499:f
10479:)
10476:x
10473:(
10468:n
10464:p
10443:1
10437:n
10434:2
10414:)
10411:x
10408:(
10405:f
10383:j
10380:i
10372:=
10369:)
10364:j
10360:x
10356:(
10353:f
10333:)
10330:x
10327:(
10322:n
10318:p
10306:j
10304:x
10286:2
10281:)
10275:j
10271:x
10262:i
10258:x
10253:(
10246:2
10241:)
10235:j
10231:x
10224:x
10220:(
10207:i
10201:j
10194:n
10188:j
10182:1
10170:=
10167:)
10164:x
10161:(
10158:f
10139:2
10133:n
10130:2
10101:)
10096:i
10092:x
10088:(
10082:1
10079:+
10076:n
10072:p
10066:b
10063:1
10041:)
10036:i
10032:x
10028:(
10023:1
10017:n
10013:p
9992:)
9987:i
9983:x
9979:(
9974:n
9970:p
9949:)
9944:i
9940:x
9936:(
9931:1
9925:n
9921:p
9917:)
9914:b
9911:(
9908:+
9905:)
9900:i
9896:x
9892:(
9887:n
9883:p
9879:)
9876:a
9873:(
9870:=
9867:)
9862:i
9858:x
9854:(
9849:1
9846:+
9843:n
9839:p
9818:)
9815:x
9812:(
9807:1
9804:+
9801:n
9797:p
9776:)
9773:x
9770:(
9765:n
9761:p
9738:i
9734:w
9720:1
9705:)
9700:i
9696:x
9692:(
9684:n
9680:p
9668:2
9651:x
9648:d
9643:2
9639:)
9635:x
9632:(
9627:1
9621:n
9617:p
9613:)
9610:x
9607:(
9599:b
9594:a
9583:)
9578:i
9574:x
9570:(
9565:1
9559:n
9555:p
9549:1
9543:n
9539:a
9532:n
9528:a
9522:=
9519:x
9516:d
9508:i
9504:x
9497:x
9492:)
9489:x
9486:(
9481:n
9477:p
9470:)
9467:x
9464:(
9456:b
9451:a
9426:)
9423:x
9420:(
9415:1
9409:n
9405:p
9384:2
9378:n
9352:1
9346:n
9342:a
9337:)
9334:x
9331:(
9326:1
9320:n
9316:p
9309:+
9305:)
9297:1
9291:n
9287:a
9282:)
9279:x
9276:(
9271:1
9265:n
9261:p
9249:1
9243:n
9239:x
9234:(
9230:=
9225:1
9219:n
9215:x
9192:x
9189:d
9184:1
9178:n
9174:x
9170:)
9167:x
9164:(
9159:1
9153:n
9149:p
9145:)
9142:x
9139:(
9131:b
9126:a
9115:)
9110:i
9106:x
9102:(
9097:1
9091:n
9087:p
9080:n
9076:a
9070:=
9067:x
9064:d
9056:i
9052:x
9045:x
9040:)
9037:x
9034:(
9029:n
9025:p
9018:)
9015:x
9012:(
9004:b
8999:a
8974:)
8971:x
8968:(
8963:1
8957:n
8953:p
8942:)
8940:x
8938:(
8936:s
8921:2
8915:n
8905:)
8903:x
8901:(
8899:s
8884:)
8881:x
8878:(
8875:s
8872:+
8867:1
8861:n
8857:x
8851:n
8847:a
8843:=
8835:i
8831:x
8824:x
8819:)
8816:x
8813:(
8808:n
8804:p
8788:n
8768:i
8764:x
8757:x
8752:)
8749:x
8746:(
8741:n
8737:p
8713:)
8710:x
8707:(
8702:1
8696:n
8692:p
8688:=
8685:)
8682:x
8679:(
8676:q
8654:x
8651:d
8643:i
8639:x
8632:x
8627:)
8624:x
8621:(
8616:n
8612:p
8608:)
8605:x
8602:(
8599:q
8593:)
8590:x
8587:(
8579:b
8574:a
8563:)
8558:i
8554:x
8550:(
8547:q
8543:1
8538:=
8535:x
8532:d
8524:i
8520:x
8513:x
8508:)
8505:x
8502:(
8497:n
8493:p
8486:)
8483:x
8480:(
8472:b
8467:a
8452:)
8450:x
8448:(
8446:q
8431:)
8428:x
8425:(
8420:n
8416:p
8403:k
8383:i
8379:x
8372:x
8365:k
8360:)
8353:i
8349:x
8345:x
8340:(
8332:1
8309:n
8303:k
8281:x
8278:d
8270:i
8266:x
8259:x
8254:)
8251:x
8248:(
8243:n
8239:p
8232:)
8229:x
8226:(
8218:b
8213:a
8203:k
8198:i
8194:x
8190:=
8187:x
8184:d
8176:i
8172:x
8165:x
8160:)
8157:x
8154:(
8149:n
8145:p
8139:k
8135:x
8128:)
8125:x
8122:(
8114:b
8109:a
8077:i
8073:x
8066:x
8062:1
8055:k
8050:)
8043:i
8039:x
8035:x
8030:(
8025:+
8017:i
8013:x
8006:x
7999:k
7994:)
7987:i
7983:x
7979:x
7974:(
7966:1
7960:=
7952:i
7948:x
7941:x
7937:1
7920:)
7918:2
7916:(
7899:x
7896:d
7888:i
7884:x
7877:x
7872:)
7869:x
7866:(
7861:n
7857:p
7850:)
7847:x
7844:(
7836:b
7831:a
7820:)
7815:i
7811:x
7807:(
7799:n
7795:p
7790:1
7785:=
7780:i
7776:w
7744:n
7740:a
7735:)
7730:i
7726:x
7722:(
7714:n
7710:p
7703:=
7699:)
7693:j
7689:x
7680:i
7676:x
7671:(
7661:i
7655:j
7648:n
7642:j
7636:1
7605:i
7601:x
7590:x
7576:)
7573:x
7570:(
7565:n
7561:p
7538:n
7534:x
7511:n
7507:a
7480:)
7474:i
7470:x
7463:x
7459:(
7453:n
7449:a
7443:)
7440:x
7437:(
7432:n
7428:p
7421:=
7413:i
7409:x
7402:x
7396:)
7390:j
7386:x
7379:x
7375:(
7369:n
7363:j
7357:1
7346:=
7342:)
7336:j
7332:x
7325:x
7321:(
7311:i
7305:j
7298:n
7292:j
7286:1
7254:)
7251:x
7248:(
7243:1
7237:n
7233:p
7212:)
7209:x
7206:(
7201:n
7197:p
7174:i
7170:w
7147:x
7144:d
7136:j
7132:x
7123:i
7119:x
7111:j
7107:x
7100:x
7088:i
7082:j
7075:n
7069:j
7063:1
7051:)
7048:x
7045:(
7037:b
7032:a
7024:=
7019:i
7015:w
7003:i
7001:w
6985:x
6982:d
6974:j
6970:x
6961:i
6957:x
6949:j
6945:x
6938:x
6926:i
6920:j
6913:n
6907:j
6901:1
6889:)
6886:x
6883:(
6875:b
6870:a
6862:)
6857:i
6853:x
6849:(
6846:r
6841:n
6836:1
6833:=
6830:i
6822:=
6819:x
6816:d
6813:)
6810:x
6807:(
6804:r
6801:)
6798:x
6795:(
6787:b
6782:a
6767:b
6763:a
6759:)
6757:x
6755:(
6753:ω
6748:n
6744:n
6740:)
6738:x
6736:(
6734:r
6714:j
6710:x
6701:i
6697:x
6689:j
6685:x
6678:x
6666:i
6660:j
6653:n
6647:j
6641:1
6629:)
6624:i
6620:x
6616:(
6613:r
6608:n
6603:1
6600:=
6597:i
6589:=
6586:)
6583:x
6580:(
6577:r
6557:)
6552:i
6548:x
6544:(
6541:r
6531:)
6529:x
6527:(
6525:r
6506:)
6503:x
6500:(
6495:k
6491:p
6468:k
6464:x
6441:k
6437:a
6421:)
6419:1
6417:(
6397:)
6392:i
6388:x
6384:(
6379:1
6373:n
6369:p
6365:)
6360:i
6356:x
6352:(
6344:n
6340:p
6334:x
6331:d
6326:2
6322:)
6318:x
6315:(
6310:1
6304:n
6300:p
6296:)
6293:x
6290:(
6282:b
6277:a
6262:1
6256:n
6252:a
6246:n
6242:a
6236:=
6231:i
6227:w
6203:b
6199:a
6195:n
6189:n
6185:p
6180:n
6166:)
6161:i
6157:x
6150:x
6147:(
6142:i
6125:n
6118:n
6114:p
6109:)
6107:x
6105:(
6103:ω
6088:,
6085:0
6079:x
6076:d
6072:)
6069:x
6066:(
6058:)
6054:)
6049:i
6045:x
6038:x
6035:(
6030:i
6021:(
6016:)
6013:x
6010:(
6005:n
6001:p
5995:b
5990:a
5975:b
5971:a
5955:.
5952:)
5947:i
5943:x
5936:x
5933:(
5928:i
5919:)
5916:x
5913:(
5908:n
5904:p
5891:n
5887:p
5882:b
5878:a
5873:i
5871:x
5867:b
5863:a
5857:n
5853:p
5843:n
5841:2
5837:)
5835:x
5833:(
5831:h
5814:.
5811:)
5806:i
5802:x
5798:(
5795:h
5789:i
5785:w
5779:n
5774:1
5771:=
5768:i
5760:=
5757:)
5752:i
5748:x
5744:(
5741:r
5735:i
5731:w
5725:n
5720:1
5717:=
5714:i
5706:=
5703:x
5700:d
5696:)
5693:x
5690:(
5687:r
5683:)
5680:x
5677:(
5669:b
5664:a
5656:=
5653:x
5650:d
5646:)
5643:x
5640:(
5637:h
5633:)
5630:x
5627:(
5619:b
5614:a
5599:i
5585:,
5582:)
5577:i
5573:x
5569:(
5566:r
5563:=
5560:)
5555:i
5551:x
5547:(
5544:r
5541:+
5538:)
5533:i
5529:x
5525:(
5522:q
5518:)
5513:i
5509:x
5505:(
5500:n
5496:p
5492:=
5489:)
5484:i
5480:x
5476:(
5473:h
5462:n
5460:p
5455:i
5453:x
5449:)
5447:x
5445:(
5442:i
5438:l
5431:i
5427:x
5420:i
5416:w
5399:,
5394:i
5390:w
5385:)
5380:i
5376:x
5372:(
5369:r
5363:n
5358:1
5355:=
5352:i
5344:=
5341:x
5338:d
5334:)
5331:x
5328:(
5323:i
5319:l
5314:)
5311:x
5308:(
5300:b
5295:a
5286:)
5281:i
5277:x
5273:(
5270:r
5264:n
5259:1
5256:=
5253:i
5245:=
5242:x
5239:d
5235:)
5230:i
5226:x
5222:(
5219:r
5215:)
5212:x
5209:(
5204:i
5200:l
5194:n
5189:1
5186:=
5183:i
5174:)
5171:x
5168:(
5160:b
5155:a
5147:=
5144:x
5141:d
5137:)
5134:x
5131:(
5128:r
5124:)
5121:x
5118:(
5110:b
5105:a
5078:.
5075:)
5070:i
5066:x
5062:(
5059:r
5055:)
5052:x
5049:(
5044:i
5040:l
5034:n
5029:1
5026:=
5023:i
5015:=
5012:)
5009:x
5006:(
5003:r
4981:.
4973:j
4969:x
4960:i
4956:x
4948:j
4944:x
4937:x
4929:i
4923:j
4915:=
4912:)
4909:x
4906:(
4901:i
4897:l
4886:)
4884:x
4882:(
4879:i
4875:l
4867:n
4861:n
4856:)
4854:x
4852:(
4850:r
4833:.
4830:x
4827:d
4823:)
4820:x
4817:(
4814:r
4810:)
4807:x
4804:(
4796:b
4791:a
4783:=
4780:x
4777:d
4771:)
4765:)
4762:x
4759:(
4756:r
4753:+
4750:)
4747:x
4744:(
4741:q
4738:)
4735:x
4732:(
4727:n
4723:p
4716:(
4710:)
4707:x
4704:(
4696:b
4691:a
4683:=
4680:x
4677:d
4673:)
4670:x
4667:(
4664:h
4660:)
4657:x
4654:(
4646:b
4641:a
4626:)
4624:x
4622:(
4620:q
4615:n
4610:n
4608:p
4602:n
4597:)
4595:x
4593:(
4591:r
4585:n
4583:p
4577:n
4572:)
4570:x
4568:(
4566:q
4551:.
4548:)
4545:x
4542:(
4539:r
4536:+
4533:)
4530:x
4527:(
4524:q
4520:)
4517:x
4514:(
4509:n
4505:p
4501:=
4498:)
4495:x
4492:(
4489:h
4478:n
4476:p
4470:n
4468:2
4464:)
4462:x
4460:(
4458:h
4450:)
4448:b
4444:a
4442:(
4437:i
4435:x
4429:n
4427:2
4423:)
4421:x
4419:(
4417:h
4411:i
4409:w
4405:n
4400:n
4398:p
4393:i
4391:x
4387:n
4379:x
4373:n
4369:j
4363:j
4361:p
4356:n
4354:p
4332:n
4329:,
4323:,
4320:1
4317:,
4314:0
4311:=
4308:k
4299:,
4296:0
4293:=
4290:x
4287:d
4283:)
4280:x
4277:(
4272:n
4268:p
4262:k
4258:x
4253:)
4250:x
4247:(
4239:b
4234:a
4219:n
4214:n
4212:p
4171:2
4167:x
4159:e
4117:1
4104:,
4099:x
4092:e
4082:x
4037:x
4030:e
3986:2
3982:x
3975:1
3926:2
3922:x
3915:1
3911:1
3887:)
3883:β
3861:1
3849:,
3842:,
3832:)
3828:x
3825:+
3822:1
3818:(
3806:)
3802:x
3796:1
3792:(
3759:1
3738:)
3736:x
3734:(
3732:ω
3713:x
3711:(
3709:ω
3702:b
3695:a
3690:ω
3686:b
3682:a
3668:x
3665:d
3661:)
3658:x
3655:(
3652:f
3648:)
3645:x
3642:(
3634:b
3629:a
3614:ω
3581:=
3568:|
3547:|
3543:=
3539:|
3534:t
3526:|
3504:|
3499:t
3491:|
3461:=
3451:)
3445:(
3435:]
3429:)
3423:(
3405:[
3392:=
3385:)
3379:(
3376:f
3346:=
3336:)
3330:(
3320:]
3314:)
3308:(
3290:[
3277:=
3270:)
3264:(
3261:f
3230:=
3219:]
3215:)
3209:(
3206:+
3203:)
3197:(
3193:[
3186:=
3175:]
3171:)
3165:(
3162:f
3159:+
3156:)
3150:(
3147:f
3143:[
3136:=
3125:]
3120:)
3113:+
3110:)
3104:(
3097:(
3093:f
3090:+
3086:)
3079:+
3076:)
3067:(
3060:(
3056:f
3052:[
3045:=
3034:]
3029:)
3022:+
3017:2
3013:x
3005:(
3001:f
2996:2
2992:c
2988:+
2984:)
2977:+
2972:1
2968:x
2960:(
2956:f
2951:1
2947:c
2942:[
2926:x
2923:d
2918:)
2911:+
2908:x
2901:(
2897:f
2891:1
2886:1
2845:=
2840:2
2836:x
2811:=
2806:2
2802:c
2774:=
2769:1
2765:x
2740:=
2735:1
2731:c
2699:x
2696:d
2691:)
2684:+
2681:x
2674:(
2670:f
2664:1
2659:1
2645:=
2633:x
2630:d
2625:)
2619:2
2615:8
2612:+
2603:+
2600:x
2595:2
2591:8
2578:(
2574:f
2568:1
2563:1
2550:2
2546:8
2534:=
2526:t
2523:d
2520:)
2517:t
2514:(
2511:f
2500:8
2468:]
2464:1
2461:,
2458:1
2451:[
2429:]
2422:,
2419:8
2415:[
2384:t
2381:d
2376:)
2372:t
2362:]
2356:t
2338:[
2324:(
2312:8
2304:=
2301:s
2281:,
2277:s
2270:=
2267:t
2246:s
2242:8
2239:=
2236:t
2209:.
2205:)
2199:2
2195:b
2192:+
2189:a
2183:+
2178:i
2168:2
2164:a
2158:b
2151:(
2147:f
2142:i
2138:w
2132:n
2127:1
2124:=
2121:i
2111:2
2107:a
2101:b
2092:x
2089:d
2085:)
2082:x
2079:(
2076:f
2071:b
2066:a
2041:)
2038:w
2035:,
2029:(
2009:n
1985:2
1981:a
1975:b
1969:=
1960:d
1955:x
1952:d
1924:d
1915:d
1910:x
1907:d
1899:)
1893:2
1889:b
1886:+
1883:a
1877:+
1869:2
1865:a
1859:b
1852:(
1848:f
1843:1
1838:1
1827:=
1824:x
1821:d
1817:)
1814:x
1811:(
1808:f
1803:b
1798:a
1709:7
1700:2
1697:+
1694:5
1687:3
1684:1
1637:+
1602:7
1593:2
1587:5
1580:3
1577:1
1459:5
1456:6
1448:7
1445:2
1440:+
1435:7
1432:3
1386:+
1351:5
1348:6
1340:7
1337:2
1327:7
1324:3
1286:9
1283:5
1254:5
1251:3
1217:9
1214:8
1174:3
1170:1
1137:i
1135:w
1128:i
1126:x
1120:n
1098:.
1090:2
1085:]
1081:)
1076:i
1072:x
1068:(
1060:n
1056:P
1051:[
1045:)
1039:2
1034:i
1030:x
1023:1
1019:(
1014:2
1009:=
1004:i
1000:w
988:n
986:P
982:i
977:i
975:x
971:i
964:n
960:P
955:n
951:)
949:x
947:(
944:n
940:P
921:]
918:1
915:,
912:1
906:[
896:)
894:x
892:(
890:f
880:n
848:i
846:x
817:2
813:x
806:1
792:(
775:2
771:x
764:1
760:1
736:.
732:)
724:i
720:x
716:(
712:g
704:i
700:w
694:n
689:1
686:=
683:i
672:x
669:d
665:)
662:x
659:(
656:g
646:)
642:x
639:+
636:1
632:(
620:)
616:x
610:1
606:(
599:1
594:1
583:=
580:x
577:d
573:)
570:x
567:(
564:f
559:1
554:1
530:'
528:i
526:w
522:'
520:i
518:x
514:)
512:x
510:(
508:g
491:,
488:1
476:,
469:,
466:)
463:x
460:(
457:g
447:)
443:x
440:+
437:1
433:(
421:)
417:x
411:1
407:(
402:=
399:)
396:x
393:(
390:f
365:n
363:2
359:)
357:x
355:(
353:f
342:n
340:2
324:,
321:)
316:i
312:x
308:(
305:f
300:i
296:w
290:n
285:1
282:=
279:i
268:x
265:d
261:)
258:x
255:(
252:f
247:1
242:1
212:n
208:i
202:i
200:w
195:i
193:x
187:n
185:2
165:n
151:x
147:x
143:x
139:x
137:(
135:y
122:3
118:2
110:y
106:y
79:)
73:(
68:)
64:(
50:.
23:.
Text is available under the Creative Commons Attribution-ShareAlike License. Additional terms may apply.