 Sigma0(:, :, i) = randpd(d) + 2*eye(d);  
end
mix0 = normalize(rand(1, nstates)); 
%% random data
localEv = randn(nobs, d); 
%% fit using mixModelFit
fprintf('Mix Gauss\n');
initParams.mu = mu0;
initParams.Sigma = Sigma0;
initParams.mixWeight = mix0; 
%mixGauss = mixModelFit(localEv, nstates, 'gauss', 'initParams', initParams, 'verbose', true, 'mixPrior', 'none');
mixGauss = mixGaussFit(localEv, nstates, 'initParams', initParams, 'verbose', true, 'mixPrior', 'none');
%% create the initial dgm
G           =