annulus_rule, a Fortran90 code which computes a quadrature rule over the interior of the annulus in 2D, with center (XC,YC), inner radius R1 and outer radius R2.
The user specifies values NT and NR, where NT is the number of equally spaced angles, and NR controls the number of radial points.
To use a rule that is equally powerful in R and T, typically, set NT = 4 * NR.
The information on this web page is distributed under the MIT license.
annulus_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.
f90_rule, a Fortran90 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).