i hate bugs

I just spent about an hour scrutinizing a few hundred lines of code before determining that

bottom = (2 * pi * Sigmas_new(j)^2)^2;

should have been

bottom = (2 * pi * Sigmas_new(j)^2)^(dim/2);

At least it works now. (I was writing an expectation maximization algorithm).