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}.