Abstract
We propose a lattice model of elongated molecules interacting via an angular part of the Berne-Pechukas potential. Monte Carlo simulations of this model are performed near the nematic to isotropic phase transition. We found that this transition is a stronger first order in the case of more elongated molecules. The comparison with the other simulations and the experimental data is performed.