clpBNR
?- {X==Y+1}, interval_degree(X,N). N = 1, X::real(-1.0Inf, 1.0Inf), Y::real(-1.0Inf, 1.0Inf). ?- interval_degree(42,N). N = 0.