disk01_positive_rule, an Octave code which demonstrates how to compute a quadrature rule of a particular precision to estimate integrals over the interior of the unit positive disk in 2D.
The unit positive disk in 2D is the set of points (X,Y) such that 0 <= X, 0 <= Y, and X^2+Y^2 <= 1.
The program sets up the nonlinear equations that characterize the points and weights of the rule, and then calls fsolve() to solve the nonlinear system.
The computer code and data files described and made available on this web page are distributed under the MIT license
disk01_positive_rule is available in a MATLAB version and an Octave version.
alpert_rule, an Octave code which sets up an Alpert quadrature rule for functions which are regular, log(x) singular, or 1/sqrt(x) singular.
annulus_rule, an Octave code which computes a quadrature rule for estimating integrals of a function over the interior of a circular annulus in 2D.
disk01_positive_monte_carlo, an Octave code which uses the Monte Carlo method to estimate integrals over the unit positive disk in 2D.