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, the s and degradation rates is a $20000*1$ array with random value