disk_rule, a C code which computes quadrature rules for approximating integrals over the interior of the general disk in 2D.
The user specifies values NT and NR, where NT is the number of equally spaced angles, and NR controls the number of radial points, the center of the disk (XC,YC), and the radius of the disk RC. The program returns vectors W(NR*NT), X(NR*NT) and Y(NR*NT), which define the rule Q(f).
To use a rule that is equally powerful in R and T, typically, set NT = 2 * NR.
Given NT, NR, and the quadrature vectors W, X, Y, the integral I(f) is estimated by Q(f) as follows:
s = w' * f(x,y);
area = pi * rc ^ 2;
q = area * s;
The information on this web page is distributed under the MIT license.
disk_rule is available in a C version and a C++ version and a Fortran90 version and a MATLAB version and an Octave version and a Python version.
c_rule, a C code which computes a quadrature rule which estimates the integral of a function f(x), which might be defined over a one dimensional region (a line) or more complex shapes such as a circle, a triangle, a quadrilateral, a polygon, or a higher dimensional region, and which might include an associated weight function w(x).
disk01_integrals, a C code which returns the exact value of the integral of any monomial over the interior of the unit disk in 2D.
disk01_monte_carlo, a C code which applies a Monte Carlo method to estimate the integral of a function over the interior of the unit disk in 2D;