funman update

[mwdchang]

Summary

Notes from dan bryce:

mono-decrease

# Susceptible is decreasing,  dS/dt <= 0, S_{t+1} - S_t <= 0
c1 = LinearConstraint(name="c1", additive_bounds=Interval(ub=0, closed_upper_bounds=True), variables=["Susceptible"], derivative=True)

mono-increase

# R is increasing, dR/dt >= 0,   R_{t+1} - R_t >= 0
c2 = LinearConstraint(name="c2", additive_bounds=Interval(lb=0), variables=["Recovered"], derivative=True)