Monte Carlo simulations of mono– and polydisperse two–dimensional crystals are reported. The particles in the studied system, interacting through hard–core repulsive Yukawa potential, form a solid phase of hexagonal lattice. The elastic properties of crystalline Yukawa systems are determined in the NpT ensemble with variable shape of the periodic box. Effects of the Debye screening length (κ), ...