[numeric] The first neq rows in A and
b are treated as linear equalities.
nleq
[numeric] The nleq rows after neq are treated as
inequations of the form a.x<=b. All remaining rows are treated as strict inequations
of the form a.x<b.
eps
[numeric] Absolute values < eps are treated as zero.
method
[character] At the moment, only the 'elimination' method is implemented.
Examples
# An infeasible system:
# x + y == 0
# x > 0
# y > 0
A <- matrix(c(1,1,1,0,0,1),byrow=TRUE,nrow=3)
b <- rep(0,3)
is_feasible(A=A,b=b,neq=1,nleq=0)
# A feasible system:
# x + y == 0
# x >= 0
# y >= 0
A <- matrix(c(1,1,1,0,0,1),byrow=TRUE,nrow=3)
b <- rep(0,3)
is_feasible(A=A,b=b,neq=1,nleq=2)