Shooting methods for a PT-symmetric periodic eigenvalue problem

We present a rigorous analysis of the performance of some one-step discretization schemes for a class of PT-symmetric singular boundary eigenvalue problem which encompasses a number of different problems whose investigation has been inspired by the 2003 article of Benilov et al. (J Fluid Mech 497:201-224, 2003). These discretization schemes are analyzed as initial value problems rather than as discrete boundary problems, since this is the setting which ties in most naturally with the formulation of the problem which one is forced to adopt due to the presence of an interior singularity. We also devise and analyze a variable step scheme for dealing with the singular points. Numerical results show better agreement between our results and those obtained from small-epsilon asymptotics than has been shown in results presented hitherto.