samer@44: function C=stats_cov(stats)
samer@44: 	m1=mstep(stats);
samer@44: 	C=m1.cov;
samer@44: end
samer@44: 
samer@44: function M=mstep(S)
samer@44: 	M.mu=S.amu/S.ag;
samer@44: 	M.cov=msym(S.aco/S.ag - M.mu*M.mu'); 
samer@44: end
samer@44: 
samer@44: function A=msym(X), A=(X+X')/2; end