Accuracy studies of the numerical method of characteristics for axisymmetric, steady supersonic flows