def unit_sigma_points(dim, degree): # 1D sigma-points (x) and weights (w) x, w = hermegauss(degree) # nD sigma-points by cartesian product return cartesian([x] * dim).T # column/sigma-point