We have discussed muon production and their role in high energy cosmic showers, simulated them and obtained their distribution on ground. Also Fitting of ellipse on this distributions was made, we calculated main parameters of the ellipse, also calculated standard deviations for each of the main parameters and Covariance matrix of geometric parameters and we can see that for 0 degrees major and minor axes are very close to each other meaning that it can quite easily be approximated to circle where for high zenith angles we have distinct difference between them meaning ellipse like shape. We can also see high fluctuations for high zenith angles because standard deviations for different parameters can increase more than four times for 70 Degree zenith angle in reference to 0 Degree