def METHOD_NAME(self, y, t, gamma, v):
    """
    Calculates the model RHS.
    """
    dS = -gamma * y[0] * y[1]
    dI = gamma * y[0] * y[1] - v * y[1]
    dR = v * y[1]
    return np.array([dS, dI, dR])