Application of effective index method to higher order modes of photonic crystal fibers