Computer code developed in R to generate the simulation results for the manuscript: 'Spatial heterogeneity generates contrasting dynamical regimes in a tritrophic Bt and non-Bt model'

 Computer code developed in R to investigate how spatial heterogeneity between Bt and Non-Bt patches influences the dynamics and stability of a tritrophic system. We developed a two-patch ordinary differential equation (ODE) model incorporating logistic growth, interspecific competition, Holling type II predation, Bt-induced mortality, and caterpillar-mediated competitive release of aphids.