Apply multivariate normal log-likelihood find the posterior distribution

I can use the chaospy to samples 10000 output values. but can not use chaospy that take 20*3 input.