The parameters for individual transitions were off. They are now fixed.
Now have nonexponential distributions in the main two codes.
Fixed bugs. Coalesced S was done wrong so two S would go to the same I. Rider was using s=0, i=1 except where i=2. Ooops.