An efficient numerical method for the optical fiber with an arbitrary index distribution is given in this paper. The index profile of fiber is simulated by a multi-layer structure, and the vectorial wave equation is replaced with the matrices. The propagation constant is calculated by using iterative multiplication of four by four materices. The further reduced compunction of two by two matrices may be obtained for the fiber with small index difference.