dim = len(dist)
//sample from the inumpyut dist
samples = dist.sample(sample, **kws)
qoi_dists = []
for i in range(0, len(poly)):
//sample the polynomial solution
if dim == 1:
dataset = poly[i](samples)
else:
dataset = poly[i](*samples)
lo = dataset.min()
After Change
dim = len(dist)
//sample from the input dist
samples = numpy.atleast_2d(dist.sample(sample, **kws))
qoi_dist = chaospy.GaussianKDE(poly(*samples))
return qoi_dist