added some commented out ICs for testing dynamics at high concentrations
and %%comments about the results
Part 3/4 is done. Added code to find the lysogeny/lysis transition point
TODO: Tested in Octave, not in Matlab.