how to solve dAdt $=$ Kon $*$ I $-$ Koff $*$ A$,$ dIdt $=$ Koff $*$ A $-$ Kon $*$ I and dXdt $=$ psi $*$ s$*$ A $-$ degradation$\_$rates $*$ X by Gillespie algorithm in python, where s and degradation rates are $20000*1$ array following the log$-$normal distribution $($ I would like to run the loop for $20000$ times$)$