def factorial(n):
"""Computes the factorial of n"""
if n == 0:
return 1
result = 1
for i in inclusive_range(2, n):
result = result * i
return result
def declination_measurement_points(bandwidth):
"""This function returns a list of declination angles in
radians where measurements should be made if the
transform has a certain bandwidth"""
thetas=[pi*(2*j+1)/(4*bandwidth) for j in inclusive_range(0,2*bandwidth-1)]
return thetas
def azimuthal_measurement_points(bandwidth):
"""This function returns a list of azimuthal angles in
radians where measurements should be made if the
transform has a certain bandwidth"""
phis=[2*pi*k/(2*bandwidth) for k in inclusive_range(0,2*bandwidth-1)]
return phis
