“…A number of methods, such as the perturbation method [7,11], the power series expansion method [8], the modified WKB method [9], the direct integration iteration method [10], the matrix method [12] and other kinds of numerical methods [3,[13][14] have been proposed to compute the normalised cutoff frequencies of modes of arbitrarily graded index fibers, each method claiming to be simpler than the other. However, all these methods are time consuming although these are valid for any kind of profile.…”