Readme file ------------------------------- 2/14/12
Numerical Procedures from "Inverse Method for Estimating Respiration
Rates from Decay Time Series"
contact: forneyd@alum.mit.edu

ILaplaceNN.m
[rho0,lk,gp,A] = iLaplaceNN(t,g,kmin,kmax,n)
Calculates non-negatively constrained Inverse Laplace transform of
decay data time series t,g

MultipoolAnalysis.m
[a,k,c,r,type,constpool,initdpool,npoolactive,npooltot]=MultipoolAnalysis(lk,rho0,t1)
Estimates the multi-pool model associated with the discrete
distribution rho0(lk).  rho0, lk are outputs from ILaplaceNN.m

ReguILT.m
[rho_c,lk,A,L,lambda_c,err_c,rgh_c,lambda,err,rgh,flag]=ReguILT(t,g,kmin,kmax,n)
Calculates the regularized non-negative inverse Laplace transform of
the decay data prescribed in t,g.   

FitLognormalToDecay.m
[mu,sigma]=FitLognormalToDecay(t,g)
Estimate the two parameters, mu and sigma, of the lognormal distribution
that best fit the decay data.

LognormalDecay.m
g=LognormalDecay(pars,t)
This function calculates the decay g, at the times t, assuming a
lognormal distribution of rates with parameters [mu,sigma]=pars. 


See individual .m and .readme files for further usage.
