updated main_longitudinal
correcting longitudinal tracking equations of motion (almost done): note
that you have to pass now phi_s = pi instead of 0 in the main file to
have no acceleration