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
	xs = [(dx*y + xmin) for y in ns]
	return dx*sum(funct(x)*w for x,w in izip(xs,ws))

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