Function ped_repulsion 

pub fn ped_repulsion(p: &Pedestrian, q: &Pedestrian, params: &Params) -> Vec2

Pairwise repulsion A * exp((r_ij - d_ij) / B) * e_ij along the line connecting q to p, where r_ij = p.radius + q.radius.