Abscissas and Weight Coefficients for Lobatto Quadrature H. H. Michels Mathematics of Computation, Vol. 17, No. 83. (Jul., 1963), pp. 237-244. Stable URL: http://links.jstor.org/sici?sici=0025-5718%28196307%2917%3a83%3c237%3aaawcfl%3e2.0.co%3b2-7 Mathematics of Computation is currently published by American Mathematical Society. Your use of the JSTOR archive indicates your acceptance of JSTOR's Terms and Conditions of Use, available at http://www.jstor.org/about/terms.html. JSTOR's Terms and Conditions of Use provides, in part, that unless you have obtained prior permission, you may not download an entire issue of a journal or multiple copies of articles, and you may use content in the JSTOR archive only for your personal, non-commercial use. Please contact the publisher regarding any further use of this work. Publisher contact information may be obtained at http://www.jstor.org/journals/ams.html. Each copy of any part of a JSTOR transmission must contain the same copyright notice that appears on the screen or printed page of such transmission. The JSTOR Archive is a trusted digital repository providing for long-term preservation and access to leading academic journals and scholarly literature from around the world. The Archive is supported by libraries, scholarly societies, publishers, and foundations. It is an initiative of JSTOR, a not-for-profit organization with a mission to help the scholarly community take advantage of advances in technology. For more information regarding JSTOR, please contact support@jstor.org. http://www.jstor.org Fri Jan 18 07:34:19 2008
Abscissas and Weight Coefficients for Lobatto Quadrature By H. H. Michels 1. Introduction. Recently, the numerical evaluation of certain collision illtegrals was studied using several different mechanical quadrature formulas, including Gaussian quadrature of high order [I, 21 and various n'ewton-cotes formulas. It was found that high accuracy could not easily be obtained, owing to the particular behavior of the integralid at the end points of integration, aiid it seemed likely that a "closed" type Gaussian formula of high order might be more eficient for this particular application. The existence of Gaussian-type quadrature formulas with one or more prescribed abscissas has beeii investigated by Lobatto [3] and Radau [I]. For the case where both ends of the integration interval are preassigned (Lobatto quadrature), the free abscissas aiid the corresponding weight coefficients have been evaluated by Radau [j]up to order 11. Nore recently abscissas and ~veights for Lobatto quadrature have beeii reported by Rabinowitz [GI for selected odd order up to 65. In some cases, however, an even-order quadrature formula may be desired aiid the results for such formulas of high order are reported in this communication. 2. Method of Computation. We are concerned with the Lobatto quadrature formulas of order n normalized by a change of variables to the interval (- 1, I) Formula (I) is exact for all polynoi~lials f(x) of degree S 2n - 3, whereas Gaussian quadrature rules are exact for degree 62n - 1. However, if the function f(x) is zero at both ends of the integration interval, only i. - 2 ordinates are involved in the calculation and a higher effective degree of precision is obtainable thail if an open Gaussian type formula is used. The free abscissas xi (12 = 2, 3,..., n - 1) are the zeros of the first derivative of the Legendre polynomial of order n - I The corresponding le eight coefficients Hi, can be found froin the expression TI-hereP?I-I(~b) is the normalized Legendre polynoinial of order n - 1. The ~veights corresponding to the fixed abscissas at n: = il are found to be Received Septe~~lber 18, 1962 237
X first approximation to the zeros of the derivative of the Legendre polyiloinial P,'(x) can be obtained in several ways. It can be sho~~n [7]that Since the zeros of P,,'(z) are the same as those of the associated Legendre polynomial P,' (.xi through the relation (6) p,,'(:c) = (2-I)~'~P,~ (e) equation (5) can be used to relate the zeros of P,;(X) to the successive zeros of the Bessel function Jl(z). A better approxinlation to the zeros of P,'(.E) can be obtained by making use of the inequalities derived by Szego [8]for the zeros of the generalized Jacobi polynoillial ~,'",~'(x). An examination of the upper and lower bounds of the zeros of Pi,-,(x) showed that two or three decimal places could be established using the relatioil where jl,a are the successire zeros of the Bessel functioil Jl(r).These initial approxiillations to the roots were iinproved using Xen-ton-Raphson iteration The Legendre polynomials and their derivatives were coillputed using the recursion formulas The eight coefficients were coinputed directly using equations (3) and (4) 3. Results. Abscissas and weights for Lobatto quadrature are presented in Table I for order n = 3(l)16, 24, 32, 40, 48, 64, 80, 96. All computatioiis were perforilled on an IBAI 7090 digital computer using extended precision routines. The tolerance for iteration 011 the roots mas set at 1 X 10-'! Several hand calculations of the roots and weight coefficients were performed. Coinplete agreement
TABLEI-Continucil - - -.-- - - -- -- -- - ----- - - - - -- - Abscissas Weights - -. - -- - -- - - - - -- 0.846:347.5646.7187231687 0.1:3498192G(i 8960834912 0. fi8(i1884g90 817.i742607 0.183611iSA52 03.5.7009201 0.48290',)8210 91:33620175 0.2207(i77935 6611008(i09 0.2492809301 0(i2399!)2e77 0.2410137903 000763504K 0.0000000000 0000000000 0.2.51930849:3 334467:3604 71 = 14 1.0000000000 0000000000 0.0109890109 8901098901 0.93903.504v52 672609013.5 0.0068372844 9768128463 0. 8678O3Oq538 503472qJ100 0.1165860558 9871163154 0.7288Ci85990 91:32(jl-1-059 0. 1600218s517 6295214241 0. 3*506394029 28G47Oti332 0.1948261493 7341611864 0.3427240133 4271284$504 0.2191262.530 0977075487 0.1163318688 8370:386766 0.23161279 t % 684570,7889 12 = 1,; 1.0000000000 0000000000 0.0093238093 2380932381 0. 96*724592(55 038:3837280 0.0380298930 2800121910 0.885082014~297ti2988~3 0.1016600703 2371806700 0.7(i;i;il!)6899 5181520070 0.14OS5116998 0242810946 0.60R2.532031 6984<571I 12 0.1727896472 336009490t5 0,42O(i:380347 1:367218092 0.19698723t59 6461333609 0.215:35339<553 637912:382;3 0.2119735859 26820920133 0.0000~00000 0000000000 0.2170481 ISv{ 48815649<51 // = 16 1.0000000000 0000000000 0.008333:3:33:3 3:333.3333:3.3 0. 969.5680462 702179:3285 0.0308303~510 0391990.540 0.8992005.3:10 9347209299 0.089393(iS7.3 2593080099 0.7920082918 6181.506393 0.1212338821 32311098:2.5 0. 0.52:3887028 8249;308047 0.lt54O269808 0716128081 0 4860594218 8713761178 0.17749191;33 9170412530 0 2998304689 0076~~20810 0.19369002:18 2320338432 0 101326270.5 2194914784 0.2019.583081 7822987119 iz = 21 1.0000000000 0000000000 0.0036231884 037971014,5 0. 986730t5cj3c50,5160883t53 0 0222388.534 6471120899 0.9357482209 2988655503 0.0396316813 3346780947 0. 9O77Ovj6751 13.706t52200 0.0363098187 2464619902 0.8434640701 3487204062 0.0719818620 3329398222 0 7641704824 2049330779 0.0863690209 6792906822 0.67124010,72 6412869984 0.0992148276 8108358741 0. e566:3;313579 7929031219 0.1102900868 9296860tll 0.4.51:3163732 14::22(?182.j 0 1193971937 0249131900 0 3282476133 7q5,51091203 0.126:3736120 280208001:3 0 199:3212e533 908312663721 0 1~10949418 7K60:1942:3.5 0 0668579927 :37228v57811 0.1334768 3-38 66986377(i0 71 = 3" 1.0000000000 0000000000 0.0020161290 32258064.?2 0 9926089339 7276155937 0.012398106e5 0137384379 0.97t52946904 8270932806 0. 02219938?28 892919646? 0.94828 C8L38-2 1723237808 0.0317731354 1091546578 0 9118499390 637t3190407 0.0410342013 8606272:3:3:3 0. 8663,724760 1267e5e51 983 0.049885271:3 362X20701 0.8122447317 77442:344e71 0.0<?82404972 -i80s5e5869t5.? 0.7.500644930 3667479772 0 0660168772 571.543439.3 0.6801297,>36 1.5.5,5081394 0.0731371396 0267903261 0.6040:323871 4812112614 0. 079q5305?~569210623229 0.3216:322628 81.56329061 0.08<51334979 496682:30.52 0.4S40477172 018469:39GO 0 0808903729,i73578:2:307
Abscissas Weights
Abscissas Weights
Xbscie
244 H. H. MICHELS with the machine results was found in all cases to 21 decimals. In addition, the following relation was used as a check on the accuracy of the results Equation (12) was satisfied to within 2 units in the 21st decimal place for all cases reported here. In the table only the positive abscissas are reported since all abscissas and weights are symmetric with z-a = -z,, and H-I, = HA. 4. Acknowledgements. This work was sponsored in part by the Air Force Office of Scientific Research of the Office of Aerosp&ce Research, and in part by the United Aircraft Corporation Research Laboratories. The author wishes to thank Miss Anne Putnam for her assistance with the calculations. Research Laboratories United Aircraft Corporation East Hartford 8, Connecticut 1. P. DAVIS & P. RABINOWITZ, "Abscissas and Weights for Gaussian Quadratures of High Order," htbsj. of Research, v. 56, 1956, p. 35. 2. P. DAVIS & P. RABINOWITZ, "Additional Abscissas and Weights for Gaussian Quadrature of High Order: values for n = 64, 80, and 96," NBS J. of Research, v. 60, 1958, p. 613. 3. R. LOBATTO, Lessen over de Integraal-Rekening, The Hague, 1852, p. 207. 4. M. R. RADAU, "Etude sur les formules d'ap roximation qui servant 3, calculer la valeur numcrique d'une integrale dhfinie." J. de Math., 6),v. 6, 1880, p. %3. 5. Ibid, p. 307. 6. P. RABINOWITZ, "Abscissas and Weights for Lobatto Quadrature of High Order," Math. Comp., v. 14, 1960, p. 47. 7. E. T. WHITAKER & G. N. WATSOS,Modern Analysis, 4th Edition, Cambridge University Press, 1927. 8. H. SZEG~, Orthogonal Polynomials, American Math. Soc. Publ., v. 23, American Math. Soc., N. Y., 1939.
http://www.jstor.org LINKED CITATIONS - Page 1 of 1 - You have printed the following article: Abscissas and Weight Coefficients for Lobatto Quadrature H. H. Michels Mathematics of Computation, Vol. 17, No. 83. (Jul., 1963), pp. 237-244. Stable URL: http://links.jstor.org/sici?sici=0025-5718%28196307%2917%3a83%3c237%3aaawcfl%3e2.0.co%3b2-7 This article references the following linked citations. If you are trying to access articles from an off-campus location, you may be required to first logon via your library web site to access JSTOR. Please visit your library's website or contact a librarian to learn about options for remote access to JSTOR. [Bibliography] 6 Abscissas and Weights for Lobatto Quadrature of High Order Philip Rabinowitz Mathematics of Computation, Vol. 14, No. 69. (Jan., 1960), pp. 47-52. Stable URL: http://links.jstor.org/sici?sici=0025-5718%28196001%2914%3a69%3c47%3aaawflq%3e2.0.co%3b2-z NOTE: The reference numbering from the original has been maintained in this citation list.