  We demonstrate that a first order isotropic-to-nematic phase transition in
liquid crystals can be succesfully modeled within the generalized Landau-de
Gennes theory by selecting an appropriate combination of elastic constants. The
numerical simulations of the model established in this paper qualitatively
reproduce the experimentally observed configurations that include interfaces
and topological defects in the nematic phase.
