The system of equations is as follows: ub'[x] == I m/x ub[x] + I (1 + (2 xi)/x) vb[x]; vb'[x] == -I m/x vb[x] + I (1 + (2 xi)/x) ub[x]; with initial conditions ub[xin]==1 and vb[xin]==1. An example set of parameters is xin=600xi where xi=3000 and m=1, and the ODES are to be solved from xin to xout=0.1. In order to get the correct physical solution xin has to be far greater than 2xi at the initial moment. Attempting to solve the equations numerically yields a solution with a very high frequenc...
