Function distanceQuantilesEdsd( plxMas, plxErrorMas, lPc, qpoints, ... )

Description:
Calculates arbitrary quantiles for the distance estimate using the Exponentially Decreasing Space Density prior.

Note this function has to numerically integrate the PDF to determine quantile values, so it is relatively slow.

Parameters:
plxMas (floating point)
parallax in mas
plxErrorMas (floating point)
parallax error in mas
lPc (floating point)
length scale in parsec
qpoints (floating point, one or more)
one or more required quantile cut points, each in the range 0..1
Return Value (array of floating point):
array with one element for each of the supplied qpoints giving the corresponding distance in parsec
Examples:
Signature:
double[] distanceQuantilesEdsd(double, double, double, double...)