docs: def points(self, size=1.0, highlight=None, colorlist=None, opacity=1.0):
          """Display the system as points.
  
          :param float size: the size of the points.
  
  
          """
          if colorlist is None:
              colorlist = [get_atom_color(t) for t in self.topology['atom_types']]
          if highlight is not None:
              if isinstance(highlight, int):
                  colorlist[highlight] = 0xff0000
              if isinstance(highlight, (list, np.ndarray)):
                  for i in highlight:
                      colorlist[i] = 0xff0000
  
          sizes = [size] * len(self.topology['atom_types'])
  
          points = self.add_representation('points', {'coordinates': self.coordinates.astype('float32'),
                                                      'colors': colorlist,
                                                      'sizes': sizes,
                                                      'opacity': opacity})
          # Update closure
          def update(self=self, points=points):
              self.update_representation(points, {'coordinates': self.coordinates.astype('float32')})
  
          self.update_callbacks.append(update)
          self.autozoom(self.coordinates)