code: def _hammer_function(self, x,x0,r_x0, s_x0):
          '''
          Creates the function to define the exclusion zones
          '''
          return norm.logcdf((np.sqrt((np.square(np.atleast_2d(x)[:,None,:]-np.atleast_2d(x0)[None,:,:])).sum(-1))- r_x0)/s_x0f _hammer_function(self, x,x0,r_x0, s_x0):
          '''
          Creates the function to define the exclusion zones
          '''
          return norm.logcdf((np.sqrt((np.square(np.atleast_2d(x)[:,None,:]-np.atleast_2d(x0)[None,:,:])).sum(-1))- r_x0)/s_x0)
docs: Creates the function to define the exclusion zones