Abstract: Constraint-based metabolic flux analysis describes the space of viable flux configurations for a metabolic network as a high-dimensional polytope defined by the linear constraints that enforce the balancing of production and consumption fluxes for each chemical species in the system. Here we compute the distribution of viable fluxes with a method that scales linearly with system size.