Abstract: Some important physical effects appear due to short range order in nematic liquid crystals. On the basis of the cell theory, new numerical method is presented, which are used to calculate these effects.We do not use ann functional expansions but we obtain accurate orientational distribution functions by solving the equation of equilibrium state, using an iterative method. Then various physical quantities can be calculated. The accurate numerical results of order parameter at phase transition and entropy change etc. are given.

Key words: nematic liquid crystals, short range order, cell theory, orientational distribution function