submodel_formula_strings.m

％The Metabolic Networks Toolbox contains functions to create, ％modify, display, and simulate bioche

function sub_kinetics = submodel_formula_strings(kinetics,n_met,met_indices,act_indices)

new_met_indices = zeros(n_met,1);
new_met_indices(met_indices) = 1:length(met_indices);

  keep_parameters = [];

for it = 1:length(act_indices),
  index = act_indices(it);
  keep_parameters = [ keep_parameters;  kinetics.reactions{index}.parameters{:}.index];
end

kinetics.reactions = kinetics.reactions(act_indices);

kinetics.parameters = kinetics.parameters(keep_parameters);
kinetics.parameter_values = kinetics.parameter_values(keep_parameters);

sub_kinetics = kinetics;