function qnew=q_after_abcd(q_start, distance) % calculates q parameter after propagating through free space qnew=q_old+distance; end