network_print_formulae.m

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

% network_print_formulae(network,actions,metabolites)
%
% print the chemical formulae
%
% actions, metabolites (optional): names to be used in the formulae
% (column lists of strings)

function network_print_formulae(network,actions,metabolites)

     if ~exist('actions','var'), actions=[]; end
     if isempty(actions), actions=network.actions; end
     if ~exist('metabolites','var'), metabolites=network.metabolites; end

     for i=1:length(actions)
     fprintf('%i:\t%s:\t',i,actions{i})
     dum=find(network.N(:,i)<0);
     if length(dum), fprintf('%i %s',-network.N(dum(1),i), metabolites{dum(1)}); end
     for j=2:length(dum);  fprintf(' + %i %s',-network.N(dum(j),i), metabolites{dum(j)}); end

 if network.reversible(i),  fprintf(' <=> '); else fprintf(' -> '); end
     dum=find(network.N(:,i)>0);
     if length(dum), fprintf('%i %s',network.N(dum(1),i), metabolites{dum(1)}); end
     for j=2:length(dum);  fprintf(' + %i %s',network.N(dum(j),i), metabolites{dum(j)}); end
     fprintf('\n')
end