The system was initially solved using the parameters in Table 1. The estimates for \(\epsilon\)\(\delta\),  \(\mu\), and \(\mu^*\) were taken from Battista et al. (2017) , while estimates for \(\alpha\)\(\beta_1\)\(\beta_2\) , and \(\zeta\) were initialized with reasonable guesses since no directly applicable data could be found. Finally,  \(\gamma\) was taken from online resources \cite{cause}.