Experiments to study high-energy cosmic rays (CRs) employ Monte Carlo codes for extensive air shower (EAS) simulations to figure out the properties of CRs. COSMOS and CORSIKA among EAS simulation codes are currently being used to analyze the data of the Telescope Array experiment. We have generated a library of about 10,000 simulated EASs with the primary energy ranging from $10^{18.5}eV$ to $10^{20}eV$ and the zenith angle of primary particles ranging from 0 to 45 degree for proton and iron primaries. We have compared the results predicted by CORSIKA and COSMOS under the same condition. In this talk, we show the differences in the energy spectra at the ground, the longitudinal shower profile as a function of atmospheric depth, the Calorimetric energy, and the Xmax distribution. We also discuss the lateral distribution function obtained from GEANT4 simulations which is being used to measure the detector response.