function C=stats_cov(stats)
	m1=mstep(stats);
	C=m1.cov;
end

function M=mstep(S)
	M.mu=S.amu/S.ag;
	M.cov=msym(S.aco/S.ag - M.mu*M.mu'); 
end

function A=msym(X), A=(X+X')/2; end
