Integral

The INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL) was a space telescope for observing gamma rays of energies up to 8 MeV. It was launched by the European Space Agency (ESA) into Earth orbit in 2002, and is designed to provide imaging and spectroscopy of cosmic sources. In the MeV energy range, it is the most sensitive gamma ray observatory in space. It is sensitive to higher energy photons than X-ray instruments such as NuSTAR, the Neil Gehrels Swift Observatory, XMM-Newton, and lower than other gamma-ray instruments such Fermi and HESS. Photons in INTEGRAL's energy range are emitted by relativistic particles in violent sources, radioactivity from unstable isotopes produced during nucleosynthesis, X-ray binaries, and astronomical transients of all types, including gamma-ray bursts. The spacecraft's instruments have very wide fields of view, which is particularly useful for detecting gamma-ray emission from transient sources as they can continuously monitor large parts of the sky. INTEGRAL was an ESA mission with additional contributions from European member states including Italy, France, Germany, and Spain. Cooperation partners are the Russian Space Agency with IKI (military CP Command Punkt KW) and NASA.…

of the region bounded by its graph and the horizontal axis; in the above graph as an example, the integral of f(x) between a and b is the yellow (−) area subtracted from the blue (+) area|alt=Definite integral example In mathematics, an integral is the continuous analog of a sum, and is used to calculate areas, volumes, and their generalizations. The process of computing an integral, called integration, is one of the two fundamental operations of calculus, along with differentiation. Integration was initially used to solve problems in mathematics and physics, such as finding the area under a curve, or determining displacement from velocity. Usage of integration expanded to a wide variety of scientific fields thereafter. A definite integral computes the signed area of the region in the plane that is bounded by the graph of a given function between two points in the real line. Conventionally, areas above the horizontal axis of the plane are positive while areas below are negative. Integrals also refer to the concept of an antiderivative, a function whose derivative is the given function; in this case, they are also called indefinite integrals. The fundamental theorem of calculus relates definite integration to differentiation and provides a method to compute the definite integral of a function when its antiderivative is known; differentiation and integration are inverse operations. Although methods of calculating areas and volumes dated from ancient Greek mathematics, the principles of integration were formulated independently by Isaac Newton and Gottfried Wilhelm Leibniz in the late 17th century, who thought of the area under a curve as an infinite sum of rectangles of infinitesimal width. Bernhard Riemann later gave a rigorous definition of integrals, which is based on a limiting procedure that approximates the area of a curvilinear region by breaking the region into infinitesimally thin vertical slabs. In the early 20th century, Henri Lebesgue generalized Riemann's formulation by introducing what is now referred to as the Lebesgue integral; it is more general than Riemann's in the sense that a wider class of functions are Lebesgue-integrable. Integrals may be generalized depending on the type of the function as well as the domain over which the integration is performed. For example, a line integral is defined for functions of two or more variables, and the interval of integration is replaced by a curve connecting two points in space. In a surface integral, the curve is replaced by a piece of a surface in three-dimensional space. The next significant advances in integral calculus did not begin to appear until the 17th century. At this time, the work of Cavalieri with his method of indivisibles, and work by Fermat, began to lay the foundations of modern calculus, with Cavalieri computing the integrals of up to degree in Cavalieri's quadrature formula. The case n = −1 required the invention of a function, the hyperbolic logarithm, achieved by quadrature of the hyperbola in 1647. Further steps were made in the early 17th century by Barrow and Torricelli, who provided the first hints of a connection between integration and differentiation. Barrow provided the first proof of the fundamental theorem of calculus. Wallis generalized Cavalieri's method, computing integrals of to a general power, including negative powers and fractional powers. \int_{a}^{b} f(x) \,dx. The integral symbol represents integration. The symbol , called the differential of the variable , indicates that the variable of integration is . The function is called the integrand, the points and are called the limits (or bounds) of integration, and the integral is said to be over the interval , called the interval of integration. A function is said to be if its integral is well-defined and finite. (The type of integral being considered determines the type of integrability this means, thus Riemann integrable means that the upper and lower Riemann sums converge over the domain, Lebesgue integrable means that the Lebesgue integral exists and is finite, and so on.) If limits are specified, the integral is called a definite integral. When the limits are omitted, as in \int f(x) \,dx, the integral is called an indefinite integral, which represents a class of functions (the antiderivative) whose derivative is the integrand. The fundamental theorem of calculus relates the evaluation of definite integrals to indefinite integrals. There are several extensions of the notation for integrals to encompass integration on unbounded domains and/or in multiple dimensions (see later sections of this article). In advanced settings, it is not uncommon to leave out when only the simple Riemann integral is being used, or the exact type of integral is immaterial. For instance, one might write \int_a^b (c_1f+c_2g) = c_1\int_a^b f + c_2\int_a^b g to express the linearity of the integral, a property shared by the Riemann integral and all generalizations thereof. ==Interpretations== Integrals appear in many practical situations. For instance, from the length, width and depth of a swimming pool which is rectangular with a flat bottom, one can determine the volume of water it can contain, the area of its surface, and the length of its edge. But if it is oval with a rounded bottom, integrals are required to find exact and rigorous values for these quantities. In each case, one may divide the sought quantity into infinitely many infinitesimal pieces, then sum the pieces to achieve an accurate approximation. As another example, to find the area of the region bounded by the graph of the function \sqrt{x} between and , one can divide the interval into five pieces ( ), then construct rectangles using the right end height of each piece (thus ) and sum their areas to get the approximation \textstyle \sqrt{\frac{1}{5}}\left(\frac{1}{5}-0\right)+\sqrt{\frac{2}{5}}\left(\frac{2}{5}-\frac{1}{5}\right)+\cdots+\sqrt{\frac{5}{5}}\left(\frac{5}{5}-\frac{4}{5}\right)\approx 0.7497, which is larger than the exact value. Alternatively, when replacing these subintervals by ones with the left end height of each piece, the approximation one gets is too low: with twelve such subintervals the approximated area is only 0.6203. However, when the number of pieces increases to infinity, it will reach a limit which is the exact value of the area sought (in this case, ). One writes \int_{0}^{1} \sqrt{x} \,dx = \frac{2}{3}, which means is the result of a weighted sum of function values, , multiplied by the infinitesimal step widths, denoted by , on the interval . ==Formal definitions== There are many ways of formally defining an integral, not all of which are equivalent. The differences exist mostly to deal with differing special cases which may not be integrable under other definitions, but are also occasionally for pedagogical reasons. The most commonly used definitions are Riemann integrals and Lebesgue integrals. a = x_0 \le t_1 \le x_1 \le t_2 \le x_2 \le \cdots \le x_{n-1} \le t_n \le x_n = b . \,\! This partitions the interval into sub-intervals indexed by , each of which is "tagged" with a specific point . A Riemann sum of a function with respect to such a tagged partition is defined as \sum_{i=1}^n f(t_i) \, \Delta_i ; thus each term of the sum is the area of a rectangle with height equal to the function value at the chosen point of the given sub-interval, and width the same as the width of sub-interval, . The mesh of such a tagged partition is the width of the largest sub-interval formed by the partition, . The Riemann integral of a function over the interval is equal to if: For all \varepsilon > 0 there exists \delta > 0 such that, for any tagged partition [a, b] with mesh less than \delta , \left| S - \sum_{i=1}^n f(t_i) \, \Delta_i \right| When the chosen tags are the maximum (respectively, minimum) value of the function in each interval, the Riemann sum becomes an upper (respectively, lower) Darboux sum, suggesting the close connection between the Riemann integral and the Darboux integral. \int f = \int_0^\infty f^*(t)\,dt where the integral on the right is an ordinary improper Riemann integral ( is a strictly decreasing positive function, and therefore has a well-defined improper Riemann integral). For a suitable class of functions (the measurable functions) this defines the Lebesgue integral. A general measurable function is Lebesgue-integrable if the sum of the absolute values of the areas of the regions between the graph of and the -axis is finite: \int_E |f|\,d\mu In that case, the integral is, as in the Riemannian case, the difference between the area above the -axis and the area below the -axis: \int_E f \,d\mu = \int_E f^+ \,d\mu - \int_E f^- \,d\mu where \begin{alignat}{3} & f^+(x) &&{}={} \max \{f(x),0\} &&{}={} \begin{cases} f(x), & \text{if } f(x) > 0, \\[2pt] 0, & \text{otherwise,} \end{cases}\\[4pt] & f^-(x) &&{}={} \max \{-f(x),0\} &&{}={} \begin{cases} -f(x), & \text{if } f(x) ==Properties== ===Linearity=== The collection of Riemann-integrable functions on a closed interval forms a vector space under the operations of pointwise addition and multiplication by a scalar, and the operation of integration f \mapsto \int_a^b f(x) \; dx is a linear functional on this vector space. Thus, the collection of integrable functions is closed under taking linear combinations, and the integral of a linear combination is the linear combination of the integrals: \int_a^b (\alpha f + \beta g)(x) \, dx = \alpha \int_a^b f(x) \,dx + \beta \int_a^b g(x) \, dx. \, Similarly, the set of real-valued Lebesgue-integrable functions on a given measure space with measure is closed under taking linear combinations and hence form a vector space, and the Lebesgue integral f\mapsto \int_E f \, d\mu is a linear functional on this vector space, so that: \int_E (\alpha f + \beta g) \, d\mu = \alpha \int_E f \, d\mu + \beta \int_E g \, d\mu. More generally, consider the vector space of all measurable functions on a measure space , taking values in a locally compact complete topological vector space over a locally compact topological field . Then one may define an abstract integration map assigning to each function an element of or the symbol , f\mapsto\int_E f \,d\mu, \, that is compatible with linear combinations. In this situation, the linearity holds for the subspace of functions whose integral is an element of (i.e. "finite"). The most important special cases arise when is , , or a finite extension of the field of p-adic numbers, and is a finite-dimensional vector space over , and when and is a complex Hilbert space. Linearity, together with some natural continuity properties and normalization for a certain class of "simple" functions, may be used to give an alternative definition of the integral. This is the approach of Daniell for the case of real-valued functions on a set , generalized by Nicolas Bourbaki to functions with values in a locally compact topological vector space. See for an axiomatic characterization of the integral. Upper and lower bounds. An integrable function on , is necessarily bounded on that interval. Thus there are real numbers and so that for all in . Since the lower and upper sums of over are therefore bounded by, respectively, and , it follows that m(b - a) \leq \int_a^b f(x) \, dx \leq M(b - a). Inequalities between functions. If for each in then each of the upper and lower sums of is bounded above by the upper and lower sums, respectively, of . Thus \int_a^b f(x) \, dx \leq \int_a^b g(x) \, dx. This is a generalization of the above inequalities, as is the integral of the constant function with value over . In addition, if the inequality between functions is strict, then the inequality between integrals is also strict. That is, if for each in , then \int_a^b f(x) \, dx Subintervals. If is a subinterval of and is non-negative for all , then \int_c^d f(x) \, dx \leq \int_a^b f(x) \, dx. Products and absolute values of functions. If and are two functions, then we may consider their pointwise products and powers, and absolute values: (fg)(x)= f(x) g(x), \; f^2 (x) = (f(x))^2, \; |f| (x) = |f(x)|. If is Riemann-integrable on then the same is true for , and \left| \int_a^b f(x) \, dx \right| \leq \int_a^b | f(x) | \, dx. Moreover, if and are both Riemann-integrable then is also Riemann-integrable, and \left( \int_a^b (fg)(x) \, dx \right)^2 \leq \left( \int_a^b f(x)^2 \, dx \right) \left( \int_a^b g(x)^2 \, dx \right). This inequality, known as the Cauchy–Schwarz inequality, plays a prominent role in Hilbert space theory, where the left hand side is interpreted as the inner product of two square-integrable functions and on the interval . Hölder's inequality. Suppose that and are two real numbers, with , and and are two Riemann-integrable functions. Then the functions and are also integrable and the following Hölder's inequality holds: \left|\int f(x)g(x)\,dx\right| \leq \left(\int \left|f(x)\right|^p\,dx \right)^{1/p} \left(\int\left|g(x)\right|^q\,dx\right)^{1/q}. For , Hölder's inequality becomes the Cauchy–Schwarz inequality. Minkowski inequality. Suppose that is a real number and and are Riemann-integrable functions. Then and are also Riemann-integrable and the following Minkowski inequality holds: \left(\int \left|f(x)+g(x)\right|^p\,dx \right)^{1/p} \leq \left(\int \left|f(x)\right|^p\,dx \right)^{1/p} + \left(\int \left|g(x)\right|^p\,dx \right)^{1/p}. An analogue of this inequality for Lebesgue integral is used in construction of L p spaces. \int_a^b f(x) \, dx over an interval is defined if . This means that the upper and lower sums of the function are evaluated on a partition whose values are increasing. Geometrically, this signifies that integration takes place "left to right", evaluating within intervals where an interval with a higher index lies to the right of one with a lower index. The values and , the end-points of the interval, are called the limits of integration of . Integrals can also be defined if : \int_a^b f(x) \, dx = - \int_b^a f(x) \, dx. With , this implies: \int_a^a f(x) \, dx = 0. The first convention is necessary in consideration of taking integrals over subintervals of ; the second says that an integral taken over a degenerate interval, or a point, should be zero. One reason for the first convention is that the integrability of on an interval implies that is integrable on any subinterval , but in particular integrals have the property that if is any element of , then: \int_a^b f(x) \, dx = \int_a^c f(x) \, dx + \int_c^b f(x) \, dx. With the first convention, the resulting relation \begin{align} \int_a^c f(x) \, dx &{}= \int_a^b f(x) \, dx - \int_c^b f(x) \, dx \\ &{} = \int_a^b f(x) \, dx + \int_b^c f(x) \, dx \end{align} is then well-defined for any cyclic permutation of , , and . F(x) = \int_a^x f(t)\, dt. Then, is continuous on , differentiable on the open interval , and F'(x) = f(x) for all in . f(x) = F'(x). If is integrable on then \int_a^b f(x)\,dx = F(b) - F(a). This allows straightforward integration of functions that can be recognized as derivatives of other functions. For example, \int_0^\pi \sin(x) \,dx = -\cos(x) \big|^{x = \pi}_{x = 0} = -\cos(\pi) - \big(-\cos(0)\big) = 2. \int_{0}^{\infty} \frac{dx}{(x+1)\sqrt{x}} = \pi has unbounded intervals for both domain and range. A "proper" Riemann integral assumes the integrand is defined and finite on a closed and bounded interval, bracketed by the limits of integration. An improper integral occurs when one or more of these conditions is not satisfied. In some cases such integrals may be defined by considering the limit of a sequence of proper Riemann integrals on progressively larger intervals. If the interval is unbounded, for instance at its upper end, then the improper integral is the limit as that endpoint goes to infinity: \int_a^\infty f(x)\,dx = \lim_{b \to \infty} \int_a^b f(x)\,dx. If the integrand is only defined or finite on a half-open interval, for instance , then again a limit may provide a finite result: \int_a^b f(x)\,dx = \lim_{\varepsilon \to 0} \int_{a+\epsilon}^{b} f(x)\,dx. That is, the improper integral is the limit of proper integrals as one endpoint of the interval of integration approaches either a specified real number, or , or . In more complicated cases, limits are required at both endpoints, or at interior points. Just as the definite integral of a positive function of one variable represents the area of the region between the graph of the function and the x-axis, the double integral of a positive function of two variables represents the volume of the region between the surface defined by the function and the plane that contains its domain. For example, a function in two dimensions depends on two real variables, x and y, and the integral of a function f over the rectangle R given as the Cartesian product of two intervals R=[a,b]\times [c,d] can be written \int_R f(x,y)\,dA where the differential indicates that integration is taken with respect to area. This double integral can be defined using Riemann sums, and represents the (signed) volume under the graph of over the domain R. Under suitable conditions (e.g., if f is continuous), Fubini's theorem states that this integral can be expressed as an equivalent iterated integral \int_a^b\left[\int_c^d f(x,y)\,dy\right]\,dx. This reduces the problem of computing a double integral to computing one-dimensional integrals. Because of this, another notation for the integral over R uses a double integral sign: \iint_R f(x,y) \, dA. Integration over more general domains is possible. The integral of a function f, with respect to volume, over an n-dimensional region D of \mathbb{R}^n is denoted by symbols such as: \int_D f(\mathbf x) d^n\mathbf x \ = \int_D f\,dV. W=\mathbf F\cdot\mathbf s. For an object moving along a path in a vector field such as an electric field or gravitational field, the total work done by the field on the object is obtained by summing up the differential work done in moving from to . This gives the line integral W=\int_C \mathbf F\cdot d\mathbf s. A surface integral generalizes double integrals to integration over a surface (which may be a curved set in space); it can be thought of as the double integral analog of the line integral. The function to be integrated may be a scalar field or a vector field. The value of the surface integral is the sum of the field at all points on the surface. This can be achieved by splitting the surface into surface elements, which provide the partitioning for Riemann sums. For an example of applications of surface integrals, consider a vector field on a surface ; that is, for each point in , is a vector. Imagine that a fluid flows through , such that determines the velocity of the fluid at . The flux is defined as the quantity of fluid flowing through in unit amount of time. To find the flux, one need to take the dot product of with the unit surface normal to at each point, which will give a scalar field, which is integrated over the surface: \int_S {\mathbf v}\cdot \,d{\mathbf S}. The fluid flux in this example may be from a physical fluid such as water or air, or from electrical or magnetic flux. Thus surface integrals have applications in physics, particularly with the classical theory of electromagnetism. \int_\gamma f(z)\,dz. This is known as a contour integral. E(x,y,z)\,dx + F(x,y,z)\,dy + G(x,y,z)\, dz where E, F, G are functions in three dimensions. A differential one-form can be integrated over an oriented path, and the resulting integral is just another way of writing a line integral. Here the basic differentials dx, dy, dz measure infinitesimal oriented lengths parallel to the three coordinate axes. A differential two-form is a sum of the form G(x,y,z) \, dx\wedge dy + E(x,y,z) \, dy\wedge dz + F(x,y,z) \, dz\wedge dx. Here the basic two-forms dx\wedge dy, dz\wedge dx, dy\wedge dz measure oriented areas parallel to the coordinate two-planes. The symbol \wedge denotes the wedge product, which is similar to the cross product in the sense that the wedge product of two forms representing oriented lengths represents an oriented area. A two-form can be integrated over an oriented surface, and the resulting integral is equivalent to the surface integral giving the flux of E\mathbf i+F\mathbf j+G\mathbf k . Unlike the cross product, and the three-dimensional vector calculus, the wedge product and the calculus of differential forms makes sense in arbitrary dimension and on more general manifolds (curves, surfaces, and their higher-dimensional analogs). The exterior derivative plays the role of the gradient and curl of vector calculus, and Stokes' theorem simultaneously generalizes the three theorems of vector calculus: the divergence theorem, Green's theorem, and the Kelvin-Stokes theorem. Integrals can be used for computing the area of a two-dimensional region that has a curved boundary, as well as computing the volume of a three-dimensional object that has a curved boundary. The area of a two-dimensional region can be calculated using the aforementioned definite integral. The volume of a three-dimensional object such as a disc or washer can be computed by disc integration using the equation for the volume of a cylinder, \pi r^2 h , where r is the radius. In the case of a simple disc created by rotating a curve about the -axis, the radius is given by , and its height is the differential . Using an integral with bounds and , the volume of the disc is equal to: \pi \int_a^b f^2 (x) \, dx. Integrals are also used in physics, in areas like kinematics to find quantities like displacement, time, and velocity. For example, in rectilinear motion, the displacement of an object over the time interval [a,b] is given by x(b)-x(a) = \int_a^b v(t) \,dt, where v(t) is the velocity expressed as a function of time. The work done by a force F(x) (given as a function of position) from an initial position A to a final position B is: W_{A\rightarrow B} = \int_A^B F(x)\,dx. Integrals are also used in thermodynamics, where thermodynamic integration is used to calculate the difference in free energy between two given states. ==Computation== ===Analytical=== The most basic technique for computing definite integrals of one real variable is based on the fundamental theorem of calculus. Let be the function of to be integrated over a given interval . Then, find an antiderivative of ; that is, a function such that on the interval. Provided the integrand and integral have no singularities on the path of integration, by the fundamental theorem of calculus, \int_a^b f(x)\,dx=F(b)-F(a). Sometimes it is necessary to use one of the many techniques that have been developed to evaluate integrals. Most of these techniques rewrite one integral as a different one which is hopefully more tractable. Techniques include integration by substitution, integration by parts, integration by trigonometric substitution, and integration by partial fractions. Alternative methods exist to compute more complex integrals. Many nonelementary integrals can be expanded in a Taylor series and integrated term by term. Occasionally, the resulting infinite series can be summed analytically. The method of convolution using Meijer G-functions can also be used, assuming that the integrand can be written as a product of Meijer G-functions. There are also many less common ways of calculating definite integrals; for instance, Parseval's identity can be used to transform an integral over a rectangular region into an infinite sum. Occasionally, an integral can be evaluated by a trick; for an example of this, see Gaussian integral. Computations of volumes of solids of revolution can usually be done with disk integration or shell integration. Specific results which have been worked out by various techniques are collected in the list of integrals. ===Symbolic=== Many problems in mathematics, physics, and engineering involve integration where an explicit formula for the integral is desired. Extensive tables of integrals have been compiled and published over the years for this purpose. With the spread of computers, many professionals, educators, and students have turned to computer algebra systems that are specifically designed to perform difficult or tedious tasks, including integration. Symbolic integration has been one of the motivations for the development of the first such systems, like Macsyma and Maple. A major mathematical difficulty in symbolic integration is that in many cases, a relatively simple function does not have integrals that can be expressed in closed form involving only elementary functions, include rational and exponential functions, logarithm, trigonometric functions and inverse trigonometric functions, and the operations of multiplication and composition. The Risch algorithm provides a general criterion to determine whether the antiderivative of an elementary function is elementary and to compute the integral if is elementary. However, functions with closed expressions of antiderivatives are the exception, and consequently, computerized algebra systems have no hope of being able to find an antiderivative for a randomly constructed elementary function. On the positive side, if the 'building blocks' for antiderivatives are fixed in advance, it may still be possible to decide whether the antiderivative of a given function can be expressed using these blocks and operations of multiplication and composition and to find the symbolic answer whenever it exists. The Risch algorithm, implemented in Mathematica, Maple and other computer algebra systems, does just that for functions and antiderivatives built from rational functions, radicals, logarithm, and exponential functions. Some special integrands occur often enough to warrant special study. In particular, it may be useful to have, in the set of antiderivatives, the special functions (like the Legendre functions, the hypergeometric function, the gamma function, the incomplete gamma function and so on). Extending Risch's algorithm to include such functions is possible but challenging and has been an active research subject. More recently a new approach has emerged, using D-finite functions, which are the solutions of linear differential equations with polynomial coefficients. Most of the elementary and special functions are D-finite, and the integral of a D-finite function is also a D-finite function. This provides an algorithm to express the antiderivative of a D-finite function as the solution of a differential equation. This theory also allows one to compute the definite integral of a D-function as the sum of a series given by the first coefficients and provides an algorithm to compute any coefficient. Rule-based integration systems facilitate integration. Rubi, a computer algebra system rule-based integrator, pattern matches an extensive system of symbolic integration rules to integrate a wide variety of integrands. This system uses over 6600 integration rules to compute integrals. The method of brackets is a generalization of Ramanujan's master theorem that can be applied to a wide range of univariate and multivariate integrals. A set of rules are applied to the coefficients and exponential terms of the integrand's power series expansion to determine the integral. The method is closely related to the Mellin transform. ===Numerical=== Definite integrals may be approximated using several methods of numerical integration. The rectangle method relies on dividing the region under the function into a series of rectangles corresponding to function values and multiplies by the step width to find the sum. A better approach, the trapezoidal rule, replaces the rectangles used in a Riemann sum with trapezoids. The trapezoidal rule weights the first and last values by one half, then multiplies by the step width to obtain a better approximation. The idea behind the trapezoidal rule, that more accurate approximations to the function yield better approximations to the integral, can be carried further: Simpson's rule approximates the integrand by a piecewise quadratic function. Riemann sums, the trapezoidal rule, and Simpson's rule are examples of a family of quadrature rules called the Newton–Cotes formulas. The degree Newton–Cotes quadrature rule approximates the polynomial on each subinterval by a degree polynomial. This polynomial is chosen to interpolate the values of the function on the interval. Higher degree Newton–Cotes approximations can be more accurate, but they require more function evaluations, and they can suffer from numerical inaccuracy due to Runge's phenomenon. One solution to this problem is Clenshaw–Curtis quadrature, in which the integrand is approximated by expanding it in terms of Chebyshev polynomials. Romberg's method halves the step widths incrementally, giving trapezoid approximations denoted by , , and so on, where is half of . For each new step size, only half the new function values need to be computed; the others carry over from the previous size. It then interpolate a polynomial through the approximations, and extrapolate to . Gaussian quadrature evaluates the function at the roots of a set of orthogonal polynomials. An -point Gaussian method is exact for polynomials of degree up to . The computation of higher-dimensional integrals (for example, volume calculations) makes important use of such alternatives as Monte Carlo integration. ===Mechanical=== The area of an arbitrary two-dimensional shape can be determined using a measuring instrument called planimeter. The volume of irregular objects can be measured with precision by the fluid displaced as the object is submerged. ===Geometrical=== Area can sometimes be found via geometrical compass-and-straightedge constructions of an equivalent square. ==External links== Online Integral Calculator, Wolfram Alpha. ===Online books=== Keisler, H. Jerome, Elementary Calculus: An Approach Using Infinitesimals, University of Wisconsin Stroyan, K. D., A Brief Introduction to Infinitesimal Calculus, University of Iowa Mauch, Sean, Sean's Applied Math Book, CIT, an online textbook that includes a complete introduction to calculus Crowell, Benjamin, Calculus, Fullerton College, an online textbook Garrett, Paul, Notes on First-Year Calculus Hussain, Faraz, Understanding Calculus, an online textbook Johnson, William Woolsey (1909) Elementary Treatise on Integral Calculus, link from HathiTrust. Kowalk, W. P., Integration Theory , University of Oldenburg. A new concept to an old problem. Online textbook Sloughter, Dan, Difference Equations to Differential Equations , an introduction to calculus Numerical Methods of Integration at Holistic Numerical Methods Institute P. S. Wang, Evaluation of Definite Integrals by Symbolic Manipulation (1972) — a cookbook of definite integral techniques

Text derived from https://en.wikipedia.org/w/index.php?oldid=1375552562 under CC BY-SA 4.0; modifications have been made.