code: def _mom(self, kloc, cache, **kwargs):
          
          if evaluation.get_dependencies(*list(self.inverse_map)):
              raise StochasticallyDependentError(
                  "Joint distribution with dependencies not supported.")
          output = 1.
          for dist in evaluation.sorted_dependencies(self):
              if dist not in self.inverse_map:
                  continue
              idx = self.inverse_map[dist]
              kloc_ = kloc[idx].reshape(1)
              output *= evaluation.evaluate_moment(dist, kloc_, cache=cache)
          return output
docs: Example:
              >>> dist = chaospy.J(chaospy.Uniform(), chaospy.Normal())
              >>> print(numpy.around(dist.mom([[0, 0, 1], [0, 1, 1]]), 4))
              [1. 0. 0.]
              >>> d0 = chaospy.Uniform()
              >>> dist = chaospy.J(d0, d0+chaospy.Uniform())
              >>> print(numpy.around(dist.mom([1, 1]), 4))
              0.5833