Implements the Modified Omori Law (λ(t)=K(t+c)^-p) to model the decay of aftershock activity following the 2016 Central Italy main shock. Calculates key seismological parameters (p, c, K) using non-linear least squares (scipy.optimize.curve_fit) and provides an R² assessment (R²=0.9531) of the model fit. The analysis focuses on M ≥ 2.0 events.

Python code used to carry out the statistical analysis of magnitude correlations published in https://npg.copernicus.org/articles/27/1/2020/npg-27-1-2020.html for the self-similar aftershock rates (SSAR) model presented in https://journals.aps.org/pre/abstract/10.1103/PhysRevE.94.022314, the catalog ("TimeMagLT.txt") used in the analysis, along …