reduce_n.m

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

% [L, NR,independent_metabolites, N1, dependent_metabolites ] = reduce_N(N)
%
% compute reduced stoichiometric matrix NR for independent metabolites
% and link martix L fulfilling L*NR = N

function [L, NR, independent_metabolites, N1,dependent_metabolites] = reduce_N(N)

[n_metab,n_react] = size(N);

[echelon,independent_metabolites] = rref(N');
NR                                = N(independent_metabolites,:);
L                                 = echelon(1:length(independent_metabolites),:)';

dependent_metabolites = setdiff(1:size(N,1),independent_metabolites);
N1                    = N(dependent_metabolites,:);

L=sparse(L);