How to take advantage of vectorization when computing the pdf for a multivariate gaussian?

Question

I\'ve been spending a few hours googling about this problem and it seems I can\'t find any information.

I tried coding a multivariate gaussian pdf as:

def