def pink_power(f, alpha=1, scale=1):
	""" power spectral density for pink noise 
	inputs:
		f => frequency
		alpha => exponent characteristic to pink noise
		scale => scaling factor of proportionality
	returns: power spectral density"""
	return scale / f**alpha