getAreaOverlappingD​isks

The code returns overlapping area of a 2D bumpy disks pair

Created by Philip Wang
% 000 PW Created on Feb. 23, 2020 at Yale University

This program solves the overlapping area of two bumpy disks.
See the following for bumpy disks:
Papanikolaou, S., O'Hern, C. S., & Shattuck, M. D. (2013).
Isostaticity at frictional jamming. Physical review letters, 110(19), 198002.
The journal paper can also be obtained from:
https://jamming.research.yale.edu/files/papers/bump.pdf

The program will be used as part of the method in a paper to be submitted in 2020 on frictional study of jammed particls.

input: position of the particles and the bumps
return: Alist = [label_p1, label_p2, A_analytical, A_mc];