lmult -> mult and data y has n rows, time points
data y2 has n rows, time points
increase verbosity requirement for negative states printout
exported sir function