function bisect( p, mu, s ) {
	var a;
	var b;
	var c;
	var m;
	var n;

	n = 1;
	a = mu - s;
	b = mu + s;
	while ( n < MAX_ITERATIONS ) {
		m = ( a + b ) / 2.0;
		if ( b - a < TOLERANCE ) {
			return m;
		}
		c = cosineCDF( m, mu, s );
		if ( p > c ) {
			a = m;
		} else {
			b = m;
		}
		n += 1;
	}
	return m;
}
Bisection method to find quantile as there is no closed-form expression for the inverse of the CDF.

@private
@param {Probability} p - input value
@param {number} mu - location parameter
@param {NonNegativeNumber} s - scale parameter
@returns {number} evaluated quantile function