## define a model
T = [-0.5 0.5;
    1 -1]
C = [1; -3]

r₁(x) = - (x.>1) + (x.<=1)
r₂(x) = x.*0
R₁(x) = (x.>1).*(2 .- x) + (x.<=1).*x

r = (
    r = function (x)
        [r₁(x) r₂(x)]
    end,
    R = function (x)
        [R₁(x) r₂(x)]
    end
)

bounds = [0 12; -Inf Inf]
model = SFFM.Model( T, C, r, Bounds = bounds)
