This function is useful in creating dummy points for quadrature
schemes (see quadscheme). Given the point pattern specified by the vectors
x and y, the function
creates a pattern of nrad * nper * length(x) dummy points.
Around each data point (x[i],y[i]) there are
nrad * nper dummy points, lying on nrad radii
emanating from (x[i],y[i]), with nper dummy points
equally spaced along each radius.
The (equal) spacing of dummy points along each radius is
controlled by the factor fctr.
The distance from a data point to the furthest of its associated
dummy points is fctr * M
where M is the mean nearest neighbour distance
for the data points.
If the point pattern specified by x and y contains
only one point, the nearest neighbour distance is infinite,
so the value Mdefault will be used in place of M.