code: def rvs(self, size=1, **kwargs):
          
          size = int(size)
          dtype = [(p, float) for p in self.params]
          arr = numpy.zeros(size, dtype=dtype)
          remaining = size
          keepidx = 0
          while remaining:
              draws = self._draw(size=remaining, **kwargs)
              mask = self._constraints(draws)
              addpts = mask.sum()
              arr[keepidx:keepidx+addpts] = draws[mask]
              keepidx += addpts
              remaining = size - keepidx
          return arr
docs: Returns random values for all of the parameters.