def unit_vectors(draw, min_angle=0, max_angle=2 * np.pi): phi = draw(st.floats(min_angle, max_angle, False, False)) return np.array((np.cos(phi), np.sin(phi)), dtype=np.float64)