Tries to test for flux consistent reactions in one solve
[fluxConsistent, sol]=checkFluxConsistency (model, epsilon)
model – cobra model structure
epsilon – flux threshold
fluxConsistent – empty
sol – result of solveCobraLP function