Fine tuned Erwig's et al. code as captured from their paper.