I’m new to using Julia and I’m attempting to simulate a reaction network while introducing a disturbance that changes the concentration of one species over time. I’ve learned that I can use a callback to alter the value of a specific species at a given time to simulate the disturbance. Although I’ve read the documentation, all the descriptions and examples are for cases using differential equations, not reaction networks. Here’s the code I’ve used to create the reaction network and necessary parameters: