To extend the work of Gott, Park, and Lee (1989), statistical properties of gravitational lensing in a wide variety of cosmological models involving non-zero cosmological constant is investigated, using the redshifts of both lens and source and observed angular separation of images for gravitational lens systems. We assume singular isothermal sphere as lensing galaxy in homogenous and isotropic Friedmann­Lemaitre-Robertson- Walker universe, Schechter luminosity function, standard angular diameter distance formula and other galaxy parameters used in Fukugita and Turner (1991). To find the most adequate flat cosmological model and put a limit on the value of dimensionless cosmological constant $lambda_0$, the mean value of the angular separation of images, probability distribution of angular separation and cumulative probability are calculated for given source and lens redshifts and compared with the observed values through several statistical methods. When there is no angular selection effect, models with highest value of $lambda_0$ is preferred generally. When the angular selection effects are considered, the preferred model depends on the shape of the selection functions and statistical methods; yet, models with large $lambda_0$ are preferred in general. However, the present data can not rule out any of the flat universe models with enough confidence. This approach can potentially select out best model. But at the moment, we need more data.