julia> ex = x^2 + a * y^2
x^2 + a*(y^2)
julia> C = forward_backward_contractor(ex, vars, [a])
#96 (generic function with 1 method)
julia> C(IntervalBox(-10..10, -10..10), 0..1, 17)
(([-1, 1], [-0.242536, 0.242536]), [0, 1800])
The third argument is the value of the parameter.