Using scipy you can create samples of random variable from multivariate normal distribution. See here.
If you already have a sample that was generated using multivariate normal distribution, but you don't know its parameters, you can use maximum likelihood estimator to estimate them (see this example).
If you have a sample and you want to test if it was generated using multivariate normal distribution, you can use goodness of fit test. See this discussion for relevant scipy methods.
