put constraints on mu in mixture of normals; works seems good for K = 2
1. (important) bug fix, for targetbimix; ans1,2 instead of p1, 2
2. add tolerance instead of .ne. 0
3. expand a,b to -300, 300
1. (important) bug fix, for targetbimix; ans1,2 instead of p1, 2
2. add tolerance instead of .ne. 0
3. expand a,b to -300, 300
1. (important) introduce DP mixture
2. normalize prior
3. introduce constant in DPM
4. allow shared mu, sigma, omega
5. introduce z
6. (important) mu constraints
7. new updating MCMC
8.  new priors
9. latent variable; alpha, theta
10. add z2omega function
(important) bugfix, correct mu constraints
Revert "1. (important) bug fix, for targetbimix; ans1,2 instead of p1, 2"

This reverts commit 36687ec577bf50c9776cb7b4c8d49dcd42aaba1d.
Revert "Revert "1. (important) bug fix, for targetbimix; ans1,2 instead of p1, 2""

This reverts commit 13329b0720762f6a448e473b823306c4e9fab167.
Revert "(important) bugfix, correct mu constraints"

This reverts commit a54aca85785a4353468d507174fd85ed3d4324bd.

and actually the mu constraints in MLE method is right
Merge branch 'mlemuconstraint' into mixture-whole
Merge branch 'bayesmuconstraint' into mixture-whole
change back to 100 instead of 300 for root finding
1. tol change to 0.01 instead of 0.0001, and make mu[K] as 0 when omega[K] below tol
2. add alphasave and thetasave
3. make plot function default without DP parameters