**discrete_pdf_sample_2d**,
a MATLAB code which
constructs a Probability Density Function (PDF)
from a table of sample data, and uses that PDF to create new samples.

The program presented here is hard-wired to handle a specific problem. However, the ideas used in the program are easily extended to other regions and other dimensions.

For the problem given here, we assume we have sample values of a function
**F(X,Y)** for each subregion of a region. These values might actually
represent population counts, a density, the integral of some function over
the subregion, or simply an abstract function. We implicitly assumed that
all the values are positive.

The particular region studied here is the unit square, which has been broken down into a 20x20 array of equal subsquares.

If we normalize by the sum of the data values, the result is a PDF associated with each subregion. By assigning an arbitrary order to the subregions, we can add the PDF values up to the given subregion to get a CDF (cumulative density function) for that subregion. Now given an arbitrary random value U, we can locate the subregion whose CDF value just exceeds U. Choosing a random point within this subregion gives us the sample point. If we choose many such sample points, the statistics for this sample will tend to the discrete PDF that we defined from the data we were given.

The computer code and data files described and made available on this web page are distributed under the MIT license

**discrete_pdf_sample_2d** is available in
a C version and
a C++ version and
a FORTRAN90 version and
a MATLAB version.

fem1d_sample, a MATLAB code which samples a scalar or vector finite element function of one variable, defined by FEM files, returning interpolated values at the sample points.

fem2d_sample, a MATLAB code which evaluates a finite element function defined on an order 3 or order 6 triangulation.

prob, a MATLAB code which evaluates and inverts a number of probabilistic distributions.

random_data, a MATLAB code which generates sample points for various probability distributions, spatial dimensions, and geometries;

walker_sample, a MATLAB code which efficiently samples a discrete probability vector using Walker sampling.

- discrete_cdf_to_xy.m, finds XY points corresponding to discrete CDF values.
- r8mat_write.m, writes an R8MAT file.
- set_discrete_cdf.m, sets a CDF from a discrete PDF.