docs: def randmatrix(m, n, random_seed=None):
      """Creates an m x n matrix of random values drawn using
      the Xavier Glorot method."""
      val = np.sqrt(6.0 / (m + n))
      np.random.seed(random_seed)
      return np.random.uniform(-val, val, size=(m, n))