I have drawn a random bivariate normal distribution.

matrix m1 = (1,0)
matrix c1 = (41, -22 \ -22, 6)
drawnorm x y, n(10000) means(m1) cov(c1)

How can I estimate the parameters of bivariate normal distribution using GMM?

I need to store these parameters and work on the analysis.

Any help would be much appreciated!!!
Thanks & Regards,