logo资料库

数学函数手册 Handbook of Mathematical Functions.pdf

第1页 / 共1058页
第2页 / 共1058页
第3页 / 共1058页
第4页 / 共1058页
第5页 / 共1058页
第6页 / 共1058页
第7页 / 共1058页
第8页 / 共1058页
资料共1058页,剩余部分请下载后查看
Preface’ Science Foundation The present volume is an outgrowth held at Cambridge, Mass., on September National purpose of the meeting was to evaluate of the availability opinion tables would continue to exist. that of a Conference on Mathematical 15-16, 1954, under Tables the auspices of the The in the light of Technology. tables and the Massachusetts Institute the need for mathematical of large scale computing machines. in spite of the increasing use of the new machines It was the consensus of for the basic need Numerical tables of mathematical functions are in continual demand by scien- A greater variety of functions and higher accuracy of tabula- in- tables serve tists and engineers. tion are now required as a result of scientific advances and, especially, creasing use of automatic mainly forpreliminarysurveys For those without easy access to machines, for machine operation. such tables are, of course, indispensable. of problems before programming connection, computers. of the latter the the In Consequently, the Conference version of the classical recognized tables of functions of Jahnke-Emde. that there was a pressing need for a the project, the National Science Foundation to prepare such a volume and established modernized ment of Standards mittee, with Professor Philip M. Morse of the Massachusetts as chairman, to advise ~course of its preparation. of A. Erdelyi, M. C. Gray, N. Metropolis, and J. W. Tukey. Todd, C. B. Tompkins, the staff of the National to the Chairman, In addition requested the National To imple- Bureau an Ad Hoc Advisory Com- of Technology Institute Bureau of Standards the Committee J. B. Rosser, H. C. Thacher, during the consisted Jr., John The primary aim has been to include a maximum of useful information within of a moderately in all fields. large volume, with particular attention An attempt has been made to cover the entire to supplement To carry out the goal set forth by tbe Ad Hoc Committee, the in computation work, as well as by providing tables by including the mathematical to the needs of field of special it has been that numerical methods properties the limits scientists functions. necessary are important which demonstrate the use and extension of the tables. the direction of the late Milton Abramowitz, The Handbook was prepared under Irene A. Stegun. Its success has depended greatly upon Their efforts together with appreciated. The particular the cooperation contributions are acknowledged at appropriate places in the text. Science Foundation for the preparation of the material the cooperation of of the Ad HOC of these and The sponsor- is are greatly and many mathematicians. Committee other ship of gratefully the National recognized. that individuals is hoped It in many cases acquaint will Washington, D.C. this volume will not only meet the needs of all table users but its users with new functions. ALLEN V. ASTIN, L?imctor.
to the Ninth Printing Preface The enthusiastic reception accorded the “Handbook of Mathematical Functions” is little short of unprecedented in the long history of mathe- matical tables that began when John Napier published his tables of loga- rithms the first copy came from the Assistant Secretary of Com- merce for Science and Technology, presented the 100,OOOth copy of the Handbook to the President. Today, total distribution the 150,000 mark at a scarcely diminished rate. in 1614. Only four and one-half years after the press in 1964, Myron Tribus, then Science Advisor to Lee A. DuBridge, is approaching The success of the Handbook has not ended our interest in the subject. On the contrary, we continue our close watch over the growing and chang- ing world of computation and to discuss with outside experts and among ourselves the various proposals for possible extension or supplementation of the formulas, methods and tables that make up the Handbook. In keeping with previous policy, a number of errors discovered since the last printing have been corrected. Aside from this, the mathematical tables and accompanying text are unaltered. However, some noteworthy changes have been made in Chapter 2: Physical Constants and Conversion Factors, pp. 6-8. The table on page 7 has been revised to give the values of physical constants obtained in a recent reevaluation; and pages 6 and 8 have been modified to reflect changes in definition and nomenclature of physical units and in the values adopted for the acceleration due to gravity in the revised Potsdam system. The record of continuing acceptance of the Handbook, the praise that has come from all quarters, and the fact that it is one of the most-quoted scientific publications in recent years are evidence that the hope expressed by Dr. Astin in his Preface is being amply fulfilled. LEWIS M. BRANSCOMB, Director National Bureau of Standards November 1970
Foreword This volume is the result of the cooperative The National Bureau tables and has had under of Standards consideration, of a compendium like the present one. During the NBS Applied Mathematics Division long been effort of many persons and a number out the on Tables, on May 15, 19.52, Dr. Abramo- turning IO years, has for at a Conference least for such an undertaking, support. but preliminary plans technical advice and financial Division since 1943 it has published of the National Research Council has also had an the quarterly editorial journal, supervision “Mathe- being (MTAC),, of organizations. mathematical production called by witz indicated of t,hat Division mentioned the need for The Mathematics interest in tables; active matical Tables and Aids exercised by a Committee Subsequent to for “that consensus With Institute using on several in tbe published Report tables of various and engineers National Science Foundation table production. Massachusetts needs scientists reached forth for example, task of table making also agreed Computer, tables for The Report and that participants, J. H. Curtiss, help implement that with interpolation suggested that the NSF contribute the following R. W. Hamming, these and other the advent to Computation” of the Division. the NBS Conference was drawn on Tables to the desirability in 1952 its support of Technology a z-day Conference on September kinds. Twenty-eight persons tables as well as table producers. cpnclusions recomlmendations, of of high-speed the Conference. cornputting and but definitely did not remove the need of financing the attention of activity the in on Tables was called at the 15-16, 1954, the attended, to discuss representing This conference which were set There was general agreement, equipment changed the tables”. It was for the Occasional for of Tables and a set of formulas and functions useful the production to the occasional computer”. of such a Handbook need is for a Handbook “an outstanding tables of usually encountered techniques the NBS undertake and other financial assistance. The Conference Committee: P. M. Morse (Chairman), D. H. Lehmer, C. B. Tompkins, recommendations. elected, from M. Abramowitz, J. W. Tukey, its to undertook Division independent The Bureau of Standards to produce National Science Foundation made funds available. to the Mathematics of the Bureau, which vide the NSF with Committee Tables of the Mathematics some changes of membership, The present and that judgments as Division that Conferences sometimes recommendations the Committee reconstituted is evidence became volume their was on grants ffor the work, on Revision of the National Research Council. the Committee which is signing can sometimes get acted on. the recommended To provide carried out the work, tables and the technical guidance and to pro- the Conference of Mathematical This, after this Foreword. reach conclusions V ,/”
VI FOREWORD Active work was started at the Bureau for the various chapters, and of Dr. Abramowitz. the general direction required in 1956. The overall plan, the enthusiasm Since his untimely of Irene A. Stegun. the selection task the effort has the of have had many discussions about Though many details have had t’o be argued out as they the same as were death, The workers at of the volume have remained the Committee to begin the and the members style and layout. the basic specifications of authors were contributions continued under Bureau content, came up, outlined by of the task carried out by the staff of the NBS Computing the Massachusetts Institute The Committee wishes here to register in planning, quality expert collaborators ciation of the willingness with which the plans. We hope this resulting memorial We regret he did not live the vision and industry to to see its publication. its commendation of Technology Conference of 1954. of the magnitude and Section and their these Tables, and its appre- into volume will be judged by its users to be a worthy its various suggestions were incorporated and editing collecting of its chief architect, Milton Abramowitz. P. M. MORSE, Chairman. A. ERD~LYI M. C. GRAY N. C. METROPOLIS J. B. ROSSER H. C. THACHER. JOHN TODD ‘C. B. TOMPKINS J. W. TUKEY. Jr.
Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables Edited by Milton Abramowitz and Irene A. Stegun 1. Introduction The present Handbook scientific functions problems. the volume ext,ends by E. Jahnke to workers during provide hensive and self-contained matical neering Funct.ions invaluable editions’ present by giving more numerical of mathematical functions. also been increased. The classification the chapters of that of An A. Fletcher, In general, graphs, for automatic principal mathematical lated the chapters or computers, polynomial functions, investigators that arise The well-known with summary in physical has been designed to a compre- of the mathe- and engi- Tables of and F. Emde has been in its many The of these authors accurate larger collections tabulated fields half-century. and more the of in these past the work tables, and by giving extensive properties The number of functions covered has of functions this Handbook of Mathematical in Index J. C. P. Miller, and organization is similar Tables to by and L. Rosenhead. contain numerical rational approximations tables, properties and statements of the tabu- those of computa- the of particularly to illustrate importance. tional are given the computation also range. outside their each chapter there books and papers properties matical found. Also more sive tioned tables Council (formerly to Computation). above, is important lists of listed to be found quarterly Mathematical Many of function At numerical the use of the examples tables and lie values which text in giving proofs of the mathe- in the chapter may be the the is a short bibliography in which stated in the bibliographies the end of are numerical tables. tables are given in and current in information the National Mathematics the Comprehen- Index men- on new Research of Computation Tables and Other Aids notations commonly Higher The ma.thematical are those particularly adopted Transcendental used in this Hand- in standard book Func- texts, 1-3, by A. ErdBlyi, W. Magnus, tions, Volumes F. Oberhettinger and F. G. Tricomi (McGraw- Hill, 1953-55). Some alternative notations have also been listed. The introduction of new symbols has been kept to a minimum, and an effort has been made to avoid the use of conflicting notation. 2. Accuracy of the Tables The number of significant figures given in each table has depended to some extent on the number available in existing tabulations. There has been no attempt the throughout to make it uniform Handbook, which would have been a costly and laborious undertaking. In most tables at least five significant figures have been provided, and the tabular’ intervals have generally been chosen to ensure that linear interpolation will yield. four- or five-figure accuracy, which suffices in most physical applications. Users requiring higher recent, in 1960 by McGraw-Hill, was 1 The most 2 The second edition, with L. J. Comrie added as co-author, was published and Scientific Com- the sixth, with F. Loesch added as cc-author, U.S.A., and Teubner, Germany. published in two volumes puting Service Ltd., Great Britain. in 1962 by Addison-Wesley, U.S.A., precision in their by use of higher-order described below. interpolates may obtain them interpolation procedures, In certain tables many-figured function values are given at irregular in the argument. An example is provided by Table 9.4. The pur- pose of these tables is to furnish “key values” for the checking of programs for automatic computers; no question of interpolation arises. intervals The maximum end-figure error, or “tolerance” is 6/& of 1 unit in the tables in this Handbook everywhere func- tions, and 1 unit in the case of the higher functions except in a few cases where it has been permitted to rise to 2 units. in the case of the elementary IX /-
. X INTRODUCTION 3. Auxiliary Functions and Arguments One of the objects of this Handbook is to pro- tables or computing methods which enable functions over tabulated the ranges of real values of their parameters. to achieve this object, functions the original and auxiliary arguments An example will make functions to remove at frequent use has the their to co e with t fi e pro- vide the user to evaluate complete In order been made of auxiliary infinite singularities, infinite cedure clear. part of ranges. The exponential integral of positive argument is given by recludes direct unctions Ei(x)-In P however, in function; this inter- x are well- region. the latter higher function to nine slightly The tabulated For +/(x1-~0> and jP instance, we have the required interpolate. In the present jo=.89717 4302 ji=.89823 7113 p=.527 The most convenient way to evaluate on a desk calculating machine in turn on the keyboard, plications check is reading unity. We obtain the formula is.to set o and ji and carry out t d e multi- a partial dial by then provided and p cumulatively; the multiplier by l-p j.6z,E.‘;9;72;&39717 4302)+.527(.89823 7113) it Since is known of 3 X 10 -6 in the linear result this answer .89773. The maximum possible error is composed of the error committed there is a possible error formula, we round off this in that to
INTRODUCTION XI is, that by the 3 X lo-‘, roundingJ last and so certainly formula. this example, is the 5-point one, given by .4403X 10m5, plus cannot exceed .8X lo-‘. the (2) Lagrange’s formula relevant f=A-,(p)f_z+A-,(p)f-1+Ao(p>fo+A,(p)fi In the range p=O(.Ol)l. for p=.52, Tables of the coefficients An(p) are given 25 for formula in each evaluation we accumulate multiplier now have register since the following their sum subtable. .53 and .54 in We evaluate +A&)fa in chapter the Again, in the is unity. We turn. the An(p) x 7.952 7.953 7.954 m=&(x) 9757 .89772 .89774 0379 .89775 0999 10622 10620 -2 in the third and fourth columns are first and second differences of the values of (see below) ; the smallness of the second The numbers the xezEl(x) difference provides a check on the three tions. linear is now obtained required : interpolation interpola- value The by fn=.3(.89772 9757)+.7(.89774 0379) = 239773 7192. In cases where the correct order of the Lagrange polynomial interpolations polynomials check on their adequacy. is not known, one of the prelimina may have to be performed wit T of two or more different orders as a (3) Aitken’s method of iterative for carrying The scheme tion. in the present example linear out is as follows: interpola- this process .; 1 2 3 4 5 & 7.9 8.1 7.8 8.2 7.7 : Yn=ze”G@) 7113 89823 89717 4302 7888 89927 8737 : 89608 7306 . 90029 . 89497 9666 Yo. I 89773 :89774 44034 48264 2 90220 4 98773 2 35221 Yo. 1, (I Yo, 1.2. I Yo.1.a.s.n .89773 71499 2394 1216 2706 . 89773 71938 ii 89773 71930 30 X,-X .0473 0527 . 1473 1527 . 2473 2527 -. -. -. Here yo,n=- 1 Yo x.--20 Yn 20-x x,-x Yo.1 Yo. 1 Yo.1 ,n=- G--z1 l/O.” . . ., m--l.m.n-- 1. x,-x x,-x 1 ~n-%n . l/0.1. Yo.1. . . . ., n-1.98 -, m-1.n x,-x x,-x 1 and x~--5 are used as on a the cross-product If the quantities when Z.-X forming their accumulation multipliers desk machine, in the multiplier at that stage. An extra decimal place intermediate carried interpolates of rounding guard against accumulation is the divisor register (~~-2) the in -(x,-x) to be used is usually to safe- errors. in which to some extent, but the tabular rate of convergence and at accumulation this example, with of rounding the the given argument, tabular to the nearest of the remaining values are used to achieve the the same errors, tabular then argu- The order is immaterial maximum time minimize we begin, as in argument take ments, and so on. of The number nearest values required emerges naturally tabular achieve a given precision the course of the iterations. example six values were used, even though known in advance extra row confirms check. a valuable to in in the present it was that five would suffice. The the convergence and provides Thus (4) Difference difference notation formulas. We use the central (chapter 25), S2fl safz wa Here Sf1l2=f1-f0, 8f3/a=fz-f1, . . . ,, a2/1=sf3ia-afiia=fa-2fi+fo ~af3~~=~aja-~aj~=fa-3j2+3fi-k 8'fa=~aj~fsla-6~3~2=f4-~f~+~ja-4f~+fo and so on. In table is as follows, the present example difference written function, high differences provides a check on the function values the relevant part of the the differences being lace of the The sma ness of the in units of the as is customary. last decimal B xe=El(x) 7:9 8.0 .89717 . 89823 4302 7113 SY 2754 2036 -2 -2 S4f -34 -39 Applying, formula for example, Everett’s interpolation j~=(l-P)fo+E2(P)~*jo+E4(P)~4jo+ . +Pfl+F2(P)~afl+F4(P)~4fl+ . . . . * * the numerical and tion from Table 25.1, we find takin toe flf cients Es(p), E4 !l t at values of the interpola- and F,(p) ), F,(p) ,,/
INTRODUCTION XII 10Qf.6,= .473(89717 + .527(89823 7193. = 89773 4302) + 7113) .061196(2 + .063439(2 2754) - 2036) .012(34) - .012(39) We may formula shows is approximately notice that in passing that in a linear Everett’s interpolate the error mPwfo+ F2(P)wl= m(P) + ~2(P)lk?f0+wJ Since the maximum range O
分享到:
收藏