x = MultivariateGaussian(np.array([1,2]), np.eye(2))
x.sample(3) # take 3 samples of x
print (x.entropy) # evaluates to the differential entropy
x.pdf(np.random.randn(2)) # evaluate the pdf at this input

x = MultivariateGaussian(np.array([1,2]), np.eye(2))
A = np.random.randn(10,2)
b = np.random.randn(10)
y = A@x+b # gives another MultivariateGaussian with the correct mean and covariance
y = 2*x # broadcasts correctly for multiplication
y = x+3 # and addition
y = x-2 # and subtraction
y = x/2 # and division
y = x[0] # and indexing
y = x.concat(x) # and concatenation (assuming independence)

x = MultivariateGaussian(np.array([1,2]), np.eye(2))
y = x.condition([0],np.array([1]))