code: def add(self, num_particles, D):
          
          self._plist += self._generate(num_particles, D, box=self.box,
                                        rs=self.rs)
docs: Add particles with diffusion coefficient `D` at random positions.