How to fix co variance matrix going singular in Gaussian mixture model implementation ??.

I am implementing GMM in matlab without using any machine learning library. I am able to initialize the parameters ,perform expectation and maximization for 1 iteration,but when I put expectation and maximization in loop for several iterations till convergence,the co variance matrices go singular. How can I fix this issue.