def arc(r, center, R):
            return center+np.sqrt(R**2-r**2)