analyse_n.m

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

% [ K, L, NR, null_l, pinv_NR, independent_metabolites, N1] = analyse_N(N)
%
% analyse stoichiometric matrix 
% 
% K                        right kernel                        N * K    = 0
% L                        link matrix                           L * NR = N
% NR                       reduced stoichiometric matrix (independent metabolites)
% null_l                   left kernel                         null_l*N = 0
% pinv_NR                  pseudoinverse of NR
% independent_metabolites  indices of independent metabolites
% N1                       matrix for dependent metabolites

function [ K, L, NR, null_l, pinv_NR, independent_metabolites, N1] = analyse_N(N)

 [n_metab,n_react] = size(N);

 K=null(full(N),'r');

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

if nargout>3,
 null_l =  null(full(N'),'r')';

pinv_NR = pinv(NR);
N1                    = N(dependent_metabolites,:);
end