Let's say I generated samples from a distribution $f$ whose functional form I do not know. I want to approximate this distirbution with a multi-variate normal.
We know that if samples are generated using a MVN, then the MLE for mean snd varaince are simply the sample mean and sample co-variance.
My question is that, in my case where I want to approximate with a MVN, would I just use the sample mean and sample co-variance as paramters for the MVN ?