The accuracy of mode field analysis directly affects the performance of devices based on liquid crystal optical waveguides (LCOWs). The following procedures were used to accurately solve the mode field distribution in LCOW, with a dielectric tensor being graded along one dimension. The anisotropic magnetic field coupling equations of liquid crystal are discretized by the finite-difference method, under the case that the gradient characteristic of liquid crystal field-induced reorientation has been considered. Then, the eigenvalue equation that can accurately solve the modes in the LCOW is obtained by improving the traditional difference operator matrix. Finally, a numerical example is given to compare the solution results with the coupling relationship between the field components in the eigenmode, which is to verify the effectiveness of the proposed mode solver.