Number of contacts (a) and net economic benefit (b) as functions of p . Black lines indicate the result with control, and light grey indicates the result with no control. Each point is the mean of 20 simulation runs with error bars showing ±1 standard deviation. The parameter values are: q = 0.5, , r
= 2, α = 0.25, and c = 0.05.