In spite of the fact that many simulation programs have already been commercialized, this article provides a reliable methodology for the simulation of the director field of liquid crystal. The present method is easier to be applied than ref. [1,2], especially in various asymmetric pretilt orientations and twists angles. Through this technique, ITRI/OES is developing a series of spreadsheets. Present work is simple and easy to be extended to dynamics and two/three dimensional simulations.The giving grid test also shows that for small discrete numbers the errors will become serious while position is near the rubbing surface. Some of the published simulation works do not mention about this.

No abstract