simulation of a probability on Matlab

Hi everyone,

I have this density of proability

p(x1,x2,x3,x4)=log2(1+x1/(1+(x2+x3+x4)))*x1*x2*x3*x4/sigma*exp(-(x1^2+x2^2+x3^2+x4^2)/sigma)

and I want to simulate the multiple integral of p(x1,x2,x3,x4)dx1 dx2 dx3 dx4

Is it possible on Matlab, if yes how can I do it?

Thank you so much for your help