parameter beta = 1/5 parameter gamma = 1/9 parameter population = 1000 compartment S = population-1 compartment I = 1 compartment R = 0 S+I->I+I @beta I->R @gamma