I need code of matlab or python, thanks

[EM] Write a 2-dimensional RNG for a Gaussian mixture model (GMM) pdf with 2 sub- populations. Use any function/sub-routine available in your language of choice. - Implement the expectation maximization (EM) algorithm for estimating the pdf parameters of 2-D GMMs from samples (Refer to the Noisy Clustering Paper linked on blackboard for the relevant update equations). - Compare the quality and speed your GMM-EM estimation on 300 samples of GMM distributions featuring each of the following: a) spherical covariance matrices, b) ellipsoidal covariance matrices, and c) poorly-separated subpopulations. Apply your GMM-EM algorithm to fit the "old faithful" data set to a GMM pdf with two components [EM] Write a 2-dimensional RNG for a Gaussian mixture model (GMM) pdf with 2 sub- populations. Use any function/sub-routine available in your language of choice. - Implement the expectation maximization (EM) algorithm for estimating the pdf parameters of 2-D GMMs from samples (Refer to the Noisy Clustering Paper linked on blackboard for the relevant update equations). - Compare the quality and speed your GMM-EM estimation on 300 samples of GMM distributions featuring each of the following: a) spherical covariance matrices, b) ellipsoidal covariance matrices, and c) poorly-separated subpopulations. Apply your GMM-EM algorithm to fit the "old faithful" data set to a GMM pdf with two components