Flux balance analysis