def gaussLog(xmin, xmax, funct, ng=10):
	"""
	Gauss quadature with Log singularity at x=xmin:
	xmin, xmax: boundaries of integration domain
	funct: integrand function
	ng: Gauss integration order
	"""
	ng = max((min((ng, _ngmax)), _ngmin))
	ns = _nodesLog[ng-1]
	ws = _weightsLog[ng-1];
	dx = xmax - xmin
	x = map(lambda y: (dx*y + xmin), ns)
	return dx*sum(mult(funct(x), ws))

	Gauss quadature with Log singularity at x=xmin:
	xmin, xmax: boundaries of integration domain
	funct: integrand function
	ng: Gauss integration order
	
recipes/Python/52292_Gaussiquadrature_or_without_Log/recipe-52292.py