Integral representation of the probability Distribution
Generating a PDF based on the radius of the unit circle
create a linspace from -1 to 1based on the no.of Points needed (1-r).^p
If its a uniform pdf then lets say its 1 every where. then sum(pdf(:)) gives row*col.
Now we want to generate a pdf such that sum(pdf(:)) = row * col * sampling rate;
use the distance from the centre to create the pdf.
The code is genpdf1