c4-6 (779483), страница 2
Текст из файла (страница 2)
[1]Dahlquist, G., and Bjorck, A. 1974, Numerical Methods (Englewood Cliffs, NJ: Prentice-Hall),§7.7, p. 318.Johnson, L.W., and Riess, R.D. 1982, Numerical Analysis, 2nd ed. (Reading, MA: AddisonWesley), §6.2.5, p. 307.Abramowitz, M., and Stegun, I.A. 1964, Handbook of Mathematical Functions, Applied Mathematics Series, Volume 55 (Washington: National Bureau of Standards; reprinted 1968 byDover Publications, New York), equations 25.4.58ff.Sample page from NUMERICAL RECIPES IN C: THE ART OF SCIENTIFIC COMPUTING (ISBN 0-521-43108-5)Copyright (C) 1988-1992 by Cambridge University Press.Programs Copyright (C) 1988-1992 by Numerical Recipes Software.Permission is granted for internet users to make one paper copy for their own personal use.
Further reproduction, or any copying of machinereadable files (including this one) to any servercomputer, is strictly prohibited. To order Numerical Recipes books,diskettes, or CDROMsvisit website http://www.nr.com or call 1-800-872-7423 (North America only),or send email to trade@cup.cam.ac.uk (outside North America).float quad3d(float (*func)(float, float, float), float x1, float x2)Returns the integral of a user-supplied function func over a three-dimensional region specifiedby the limits x1, x2, and by the user-supplied functions yy1, yy2, z1, and z2, as defined in(4.6.2).
(The functions y1 and y2 are here called yy1 and yy2 to avoid conflict with the namesof Bessel functions in some C libraries). Integration is performed by calling qgaus recursively.{float qgaus(float (*func)(float), float a, float b);float f1(float x);.















