Sample code for the two-compartment model simulation in our study from Compartmental structure in the secondary lymphoid tissue can slow down in vivo HIV-1 evolution in the presence of strong CTL responses

This is a sample code for the Gillespie type simulation of the virus dynamics in the two-compartment model. For details and instructions of how to compile and run it, see description within the code (commented lines).