In this paper, the mechanical properties of erythrocytes were studied numerically based upon the mechanical model originally developed by Pamplona and Calladine (ASME J. Biomech. Eng., 115, p. 149, 1993) for liposomes. The case under study is the erythrocyte stretched by a pair of laser beams in opposite directions within buffer solutions. The study aims to elucidate the effect of radiation pressure from the optical laser because up to now little is known about its influence on the cell deformation. Following an earlier study by Guck (Phys. Rev. Lett., 84, p. 5451, 2000; Biophys. J., 81, p. 767, 2001), the empirical results of the radiation pressure were introduced and imposed on the cell surface to simulate the real experimental situation. In addition, an algorithm is specially designed to implement the simulation. For better understanding of the radiation pressure on the cell deformation, a large number of simulations were conducted for different properties of cell membrane. Results are first discussed parametrically and then evaluated by comparing with the experimental data reported by Guck An optimization approach through minimizing the errors between experimental and numerical data is used to determine the optimal values of membrane properties. The results showed that an average shear stiffness around , when the nondimensional ratio of shear modulus to bending modulus ranges from 10 to 300. These values are in a good agreement with those reported in literature.