ES_SAMPLE_MODEL - Sample all model parameters for a network result = es_sample_model(N, W, ind_ext, es_constraints, es_options) This script reconstructs one model instance, including the sampled elasticities. For samplig multiple model instances, see 'es_sample_multiple'. Inputs (with nm: # metabolites; nr: # reactions) N - Stoichiometric matrix (nm x nr) W - Allosteric regulation matrix (nr x nm) ind_ext - indices of external metabolites For the inputs 'es_constraints' and 'es_options', see 'es_default_options' Outputs All output data are stored in a data structure 'result' (see 'es_sample_elasticities') es_options and es_constraints are returned with the (possible) changes made. Sampling consists of two stages: 1. Steady state (calling 'es_sample_steady_state') 2. Elasticities and MCA (calling 'es_sample_elasticities')