From: Analysis of Metabolic Subnetworks by Flux Cone Projection
Input:
   • J (the set of ProCEMs)
Output:
   • E (the set of EFPs)
Initialization:
E := ∅;
Main procedure: