add the mathematica code to calculate three-body bound state of excitons