def power_law_integrate(p_min,p_max,slope):
	"""Integrates a power law

	Args:
		p_min (float): The lower bound of the power law
		p_max (float): The upper bound of the power law
		slope (float): The slope of the power law

	Returns:
		(float): The integral of the power law x^slope from p_min
		to p_max
	"""
	upper_bound = 1/(slope+1)*p_max**(slope+1)
	lower_bound = 1/(slope+1)*p_min**(slope+1)
	return upper_bound-lower_bound