The iteration formulas for numerical solution of ray equation have been derived according to transformed ray equation. Two typical numerical examples shoaled that the method is more general and requires much less calculation time to reach a desired accuracy.