JuMP, one constraint with two inequalities in the expression, is it possible?

Specific Domains Optimization (Mathematical)
2

Happy new year and welcome,
I don’t know if JuMP supports such constraints, but you could split them up into two sets of constraints, e.g.,

@constraint(model, con4[i in Vtotal, j in Vcustomer; i != j], x[i,j] * Menergy[i,j] - (1-x[i,j])*Q
<= y[i] - y[j] )

and

@constraint(model, con5[i in Vtotal, j in Vcustomer; i != j], y[i] - y[j] <= x[i,j] * Menergy[i,j] + (1-x[i,j])*Q)

(Note that I did not check if the above syntax is correct)