We combined the Gay-Berne potential (with parameters optimized for typical liquid crystals) and the technique for simulation of mixtures. Using the Monte Carlo method, we studied spatial and orientational correlations for equimolar concentration. Preliminary two dimensional (2D) molecular dynamics results on time correlation functions are also presented.