function entropy( d1, d2 ) {
	var half;
	var hd1;
	var hd2;
	var out;

	if (
		isnan( d1 ) ||
		isnan( d2 ) ||
		d1 <= 0.0 ||
		d2 <= 0.0
	) {
		return NaN;
	}
	half = ( d1 + d2 ) / 2.0;
	hd1 = d1 / 2.0;
	hd2 = d2 / 2.0;
	out = ln( d2 / d1 ) + gammaln( hd1 ) + gammaln( hd2 ) - gammaln( half );
	out += ( 1.0-hd1 ) * digamma( hd1 );
	out += ( -1.0-hd2 ) * digamma( hd2 );
	out += half * digamma( half );
	return out;
}
Returns the differential entropy of an F distribution.

@param {PositiveNumber} d1 - numerator degrees of freedom
@param {PositiveNumber} d2 - denominator degrees of freedom
@returns {number} entropy

@example
var v = entropy( 3.0, 7.0 );
// returns ~1.298

@example
var v = entropy( 4.0, 12.0 );
// returns ~1.12

@example
var v = entropy( 8.0, 7.0 );
// returns ~1.193

@example
var v = entropy( 1.0, -0.1 );
// returns NaN

@example
var v = entropy( -0.1, 1.0 );
// returns NaN

@example
var v = entropy( 2.0, NaN );
// returns NaN

@example
var v = entropy( NaN, 2.0 );
// returns NaN