Fixed sign error in initial velocity

vy = -2.*kappa*omega_r*ampl*cos(phi) - shear*x0

should be

vy = -2.*kappa*omega_r*ampl*cos(phi) + shear*x0
New input parameters
Added Python script that visualizes particle orbit

Agreement with exact trajectory is perfect