functional residual capacity , frc , is the lung volume at the end of a relaxed expiration to atmospheric pressure . the frc is affected by lung diseases and thus is a central measure for diagnosis , treatment and monitoring in anesthesia and intensive care . the calculation of frc from measurement of multiple breath nitrogen washout ( mbnw ) when abruptly changing the fraction of inspired oxygen ( fio2 ) was introduced by darling . the method entailed the measurement of the mixed expiratory fraction of nitrogen before and after a step increase in oxygen and can be explained with the following simple example . the inspired fraction of oxygen , fio2 , is 0.21 and the balance is n2 , i.e. fin2 = 0.79 . assume dry alveolar fractions of o2 = 0.16 , fao2 , of n2 = 0.79 , fan2 , and of co2 = 0.05 , faco2 . if the equilibrium between fin2 and fan2 is disrupted by having the person inhale at fio2 of 1 while exhaling into a bag , eventually all alveolar nitrogen will be transferred to the bag . frc equals vn2 divided by the change in fraction of nitrogen that was in the lung : 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { \rm{frc } } = { { { \rm{v}}{{\rm{n}}_{\rm{2}}}{\rm{bag } } } \over { \delta { { \rm{f}}_{\rm{a}}}{{\rm{n}}_{\rm{2}}}}},\delta { { \rm{f}}_{\rm{a}}}{{\rm{n}}_{\rm{2 } } } = 0.79$$ 
 \end{document}. vn2bag can be found by multiplying the bag nitrogen fraction with the total volume accumulated in the bag . with the advent of the nitrogen analyzer , the mass spectrometer   and the use of high precision spirometers the washout could be followed on a breath - by - breath basis . this refinement eliminated the need for the bag and replaced the numerator in equation ( 1 ) with the expired volumes of nitrogen summed over the number of breaths . the modification is heavily reliant on the exact synchronization of analysis of fn2 and expired volume for the calculation of breath - by - breath volume of nitrogen . the example above utilized fio2 to step fan2 from 0.79 to 0 . instead , any step change in fio2 could be used to create a corresponding step change in nitrogen . established the ratio of functional residual air , r , divided by the sum of functional residual air and effective tidal volume , t , as  dilution rate , r , describing the nitrogen clearance curve : 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ 
 r = \frac{r}{{r   +   t\prime } } . this factor multiplied by initial alveolar fraction of nitrogen and raised to the power of breath number allowed for the calculation of fraction of nitrogen in the expired volume of air . fowler   expanded the concept of the alveolar dilution factor to encompass three lung phases and introduced the term  pulmonary nitrogen clearance delay , adding another characterization of the nitrogen clearance curve . the list has since been extended by investigators in pulmonary physiology , see [ 8 , 9 ] . recently , a modification of the mbnw method has been introduced , implicitly utilizing a first order difference equation for the calculation of frc based on mass conservation and a change in fio2 of 0.1 . in the present study , the difference equation is solved for the step response using basic digital signal processing ( dsp ) techniques   where the step is interpreted as the nitrogen fraction change initiating the measurement . the resulting step response equation by comparing calculated clearance data to experimental and clinical washout data . finally , we discuss technical and clinical limitations in addition to characteristics of the measurement . the first order difference equation based on mass conservation of nitrogen over one breath after a step change in inhaled n2 fraction is described in equation ( 3 ) . \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{v}}_{{\rm{ti}}}}{\rm { \times } } { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{_{{\rm{low}}}}{\rm { + frc } } \times { { \rm{f}}_{\rm{a}}}{{\rm{n}}_{\rm{2}}}_{_{n{\rm { - 1}}}}{\rm { = } } { { \rm{v}}_{{\rm{te}}}}{\rm { \times } } { { \rm{f}}_{{\rm{\bar e}}}}{{\rm{n}}_{\rm{2}}}_{_n}{\rm { + frc \times } } { { \rm{f}}_{\rm{a}}}{{\rm{n}}_{\rm{2}}}_{_n}$$ 
 \end{document } the step change in nitrogen occurs at breath index n = 0 . the step is from \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{_{{\rm{high}}}}$$ 
 \end{document } to \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{_{{\rm{low}}}}$$ 
 \end{document } which will lead to n2 washing out . \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{{\rm{\bar e}}}}{{\rm{n}}_{\rm{2}}}$$ 
 \end{document } is an average breath fraction of expired gas which includes gas from anatomical ( series ) dead space . \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{\rm{a}}}{{\rm{n}}_{\rm{2}}}_{_n}$$ 
 \end{document } depends on the fraction of the immediately preceding breath \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{\rm{a}}}{{\rm{n}}_{\rm{2}}}_{_{n - 1}}$$ 
 \end{document}. vti and vte designate inspiratory and expiratory tidal volumes . in all equations a refinement is to measure the lung volume to the exclusion of tidal anatomical dead space , vdtanat . thus , the vte and vti of equation ( 3 ) must be replaced by the alveolar tidal volumes , vate and vati . also , using endtidal co2 fraction , fetco2 , as approximating alveolar fraction , faco2 , the following equations can be used to calculate the alveolar tidal volumes : 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { \rm{v}}{{\rm{a}}_{{\rm{te}}}}{\rm { = } } { { { { \rm{v}}_{\rm{t}}}{\rm{c}}{{\rm{o}}_{\rm{2 } } } } \over { { { \rm{f}}_{{\rm{et}}}}{\rm{c}}{{\rm{o}}_{\rm{2}}}}}$$ 
 \end{document } , where vtco2 is the tidal excretion of co2 . \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { \rm{v}}{{\rm{a}}_{{\rm{ti}}}}{\rm { = v}}{{\rm{a}}_{{\rm{te}}}}{\rm { + } } \left ( { { { { { \rm{v}}_{\rm{t}}}{\rm{c}}{{\rm{o}}_{\rm{2 } } } } \over { { \rm{rq}}}}{\rm { - } } { { \rm{v}}_{\rm{t}}}{\rm{c}}{{\rm{o}}_{\rm{2 } } } } \right)$$ 
 \end{document } rq is the respiratory quotient defined as \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { \rm{rq = } } { { { { \rm{v}}_{\rm{t}}}{\rm{c}}{{\rm{o}}_{\rm{2 } } } } \over { { { \rm{v}}_{\rm{t}}}{{\rm{o}}_{\rm{2}}}}}$$ 
 \end{document}. the value of rq depends on the composition of substrate being metabolized . for example , rq = 1 when carbohydrates are being metabolized and rq = 0.7 when fat is being metabolized . the range of rq is typically from 0.7 to 1.0 with reservation for gluconeo- and lipogenesis . a fixed value of rq of 0.8 is used in equation ( 5 ) . for the purpose of accumulating the in- and expired volumes of nitrogen , the fractions of nitrogen are calculated as residue . \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}{\rm { = 1 - } } { { \rm{f}}_{\rm{i}}}{{\rm{o}}_{\rm{2}}}$$ 
 \end{document}\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{\rm{2}}}{\rm { = 1 - } } { { \rm{f}}_{{\rm{et}}}}{\rm{c}}{{\rm{o}}_{\rm{2}}}{\rm { - } } { { \rm{f}}_{{\rm{et}}}}{{\rm{o}}_{\rm{2}}}$$ 
 \end{document } substituting for alveolar volumes and endtidal fractions , equation ( 3 ) becomes : 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { \rm{v}}{{\rm{a}}_{{\rm{ti}}}}{\rm { \times } } { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{_{{\rm{low}}}}{\rm { + frc \times } } { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{\rm{2}}}_{_{n - 1 } } = { \rm{v}}{{\rm{a}}_{{\rm{te}}}}{\rm { \times } } { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{\rm{2}}}_n + { \rm{frc \times } } { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{\rm{2}}}_{_n}$$ 
 \end{document } this is a first order difference equation . the fin2 and fetn2 are based on samples from the continuous fco2 and fo2 signals according to equations 6 and 7 , see figure 1 . vati and vate are calculated according to equations ( 4 ) and ( 5 ) . 1sampling at two frequencies , the waveform sampling rate and two samples per breath sampling rate . lower frame : the continuous sampling of fo2 ( left y - axis ) and fco2 ( right y - axis ) during a n2 washout procedure . from these tracings ( enlarged in upper frame ) two samples of fo2 and fco2 are sampled once during inspiration and as endtidal values for entering into equations 6 and 7 . sampling at two frequencies , the waveform sampling rate and two samples per breath sampling rate . lower frame : the continuous sampling of fo2 ( left y - axis ) and fco2 ( right y - axis ) during a n2 washout procedure . from these tracings ( enlarged in upper frame ) two samples of fo2 and fco2 are sampled once during inspiration and as endtidal values for entering into equations 6 and 7 . in the presence of alveolar ( parallel ) dead space , fetco2 < faco2 and , as a result , the equations include alveolar dead space . first , the frc is found by summing both sides of equation 8 over the total number of breaths , 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { \rm{frc } } = { { \sum ( { \rm{v}}{{\rm{a}}_{{\rm{te}}}}{\rm { \times } } { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{\rm{2}}}{\rm { - v}}{{\rm{a}}_{{\rm{ti}}}}{\rm { \times } } { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{_{{\rm{low } } } } ) } \over { { { \rm{f}}_{{\rm{etinit}}}}{{\rm{n}}_{\rm{2}}}{\rm { - } } { { \rm{f}}_{{\rm{etlatest}}}}{{\rm{n}}_{\rm{2 } } } } } = { { \sum ( { \rm{v}}{{\rm{a}}_{\rm{e}}}{{\rm{n}}_{\rm{2}}}{\rm { - v}}{{\rm{a}}_{\rm{i}}}{{\rm{n}}_{\rm{2 } } } ) } \over { { { \rm{f}}_{{\rm{etinit}}}}{{\rm{n}}_{\rm{2}}}{\rm { - } } { { \rm{f}}_{{\rm{etlatest}}}}{{\rm{n}}_{\rm{2 } } } } } = { { \sum { \rm{\delta v}}{{\rm{a}}_{\rm{n}}}_{_{\rm{2 } } } } \over { { \rm{\delta } } { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{\rm{2}}}}}$$ 
 \end{document}. fetinitn2 denotes the initial condition at time n equals 1 which is the time sample before the step ( \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{{\rm{etinit}}}}{{\rm{n}}_{\rm{2}}}{\rm { = } } { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{\rm{2}}}_{_{{\rm { - 1}}}}$$ 
 \end{document } ) . fetlatestn2 represents \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{\rm{2}}}_{_n}$$ 
 \end{document } when n equals the most recent breath . in general , vate , vati and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{_{{\rm{low}}}}$$ 
 \end{document } can vary with n and may be replaced by \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { \rm{v}}{{\rm{a}}_{{\rm{te}}}}_n$$ 
 \end{document } , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } \begin{document } 
  $ $ { \rm{v}}{{\rm{a}}_{{\rm{ti}}}}_{_n}$$ 
 \end{document } and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{_n}$$ 
 \end{document } respectively . second , equation 8 is solved for the step response using the frc calculated by equation ( 9 ) 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ \eqalign { & { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{{2_n } } } = { { { \rm{v}}{{\rm{a}}_{{\rm{ti}}}}{\rm { \times } } { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{_{{\rm{low } } } } } \over { { \rm{v}}{{\rm{a}}_{{\rm{te}}}}{\rm { + frc } } } } \times \left ( { { { 1 - { \lambda^{n + 1 } } } \over { 1 - \lambda } } } \right ) + { \lambda^{n + 1 } } \times { { \rm{f}}_{{\rm{etinit}}}}{{\rm{n}}_2 } \cr & \lambda = { { { \rm{frc } } } \over { { \rm{frc + v}}{{\rm{a}}_{{\rm{te } } } } } } \cr } $ $ \end{document } in this solution , the values , vati , vate and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{_{{\rm{low}}}}$$ 
 \end{document } are assumed to remain constant during the washout . for a full derivation , see appendix a. this is an equation that gives the fetn2 for a particular breath without iteration and can be used to determine the applicability of a first order model . some implications of equation 10 are immediately apparent : 
  is related to the time constant of the n2 clearance.as the vate increases , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{\rm{2}}}_{_n}$$ 
 \end{document } decays faster consistent with a faster lung washout at larger tidal breaths.if the frc becomes large compared to the vate , then  approaches 1 and there is very slow decay consistent with a small tidal breath requiring more breaths to change the concentration of a large lung volume ( frc).as n   , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_2}_n \to { { { \rm{v}}{{\rm{a}}_{{\rm{ti } } } } \cdot { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{{\rm{low } } } } \over { { \rm{v}}{{\rm{a}}_{{\rm{te}}}}}}$$ 
 \end{document}. the \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{\rm{2}}}_{_n}$$ 
 \end{document } will not reach \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{_{{\rm{low}}}}$$ 
 \end{document } when rq  1 since vati  vate . one would expect a good match between the measured fetn2 sequence compared to the calculated values if the system being measured is indeed first order . is related to the time constant of the n2 clearance . as the vate increases , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{\rm{2}}}_{_n}$$ 
 \end{document } decays faster consistent with a faster lung washout at larger tidal breaths . if the frc becomes large compared to the vate , then  approaches 1 and there is very slow decay consistent with a small tidal breath requiring more breaths to change the concentration of a large lung volume ( frc ) . as n   , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_2}_n \to { { { \rm{v}}{{\rm{a}}_{{\rm{ti } } } } \cdot { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{{\rm{low } } } } \over { { \rm{v}}{{\rm{a}}_{{\rm{te}}}}}}$$ 
 \end{document}. the \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } $ $ { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{\rm{2}}}_{_n}$$ 
 \end{document } will not reach \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{_{{\rm{low}}}}$$ 
 \end{document } when rq  1 since vati  vate . one would expect a good match between the measured fetn2 sequence compared to the calculated values if the system being measured is indeed first order . the modified mbnw has been extensively validated in a metabolic lung model . in this paper , the prediction of fetn2 by the frc step response equation was validated by comparing calculated to measured values of fetn2 from washout procedures in one experimental and two clinical settings . the clinical measurements were made as part of two other studies for which the ethical committee of the medical faculty at gteborg university approved the studies and informed consent was obtained from the patients , next of kin or volunteers , respectively . thus , they represent examples and are not the result of a selection powered to demonstrate differences between the normal and disease categories . model . an experimental setup using an o2-consuming / co2-producing metabolic lung model previously described [ 10 , 13 ] . the model was ventilated in volume controlled ventilation with a ge engstrm carestation ventilator ( madison , wi , usa ) . settings comprised two frcs ( 1.8 , 2.9 l ) , two tidal volumes ( 500 , 750 ml ) and two steps of fio2 ( 0.1 , 0.3 ) at \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { \rm{\dot vc}}{{\rm{o}}_2}$$ 
 \end{document } 200 ml / min and rq 1 . frc measurements were computed with the commercial software algorithm based on equation 9 detailed above . patient recordings were performed either in the operating theater using ge adu/5 ventilator or in the intensive care unit ( icu ) using a servo 900c or servo 300 ventilator ( siemens elema , solna , sweden ) . subjects with assumed normal lungs were measured peroperatively and patients with diseased lungs in the intensive care unit , where positive pressure ventilation was instituted due to respiratory failure . the eigenvalues ,  , were calculated for one patient with acute respiratory distress syndrome ( ards ) , one with heart failure ( hf ) , two patients with chronic obstructive pulmonary disease ( copd ) , and one obese patient with normal lungs ( perioperative registration ) . mono- vs biexponential fits were examined in one of the icu patients ( copd ) and in the bariatric patient . recordings in healthy volunteers ( two female and one male ) spontaneously breathing while attached to a high flow cpap system via a snorkeling mouthpiece and equipped with a nose clip in supine , sitting and standing position . in patients and volunteers , respiratory monitoring and gas sampling were performed with datex - ohmeda spirometry module in an as/3 intensive care monitor . data were collected using datex - ohmeda s/5 collect version 4 software ( datex - ohmeda , helsinki , finland ) and analyzed using lung monitor ( plug - in to datex - ohmeda s/5 collect ) or proprietary software programmed in test- point ( measurement computing corporation , norton , ma , usa ) . agreement between calculated and measured values of fetn2 was assessed by bland & altman analysis   in terms of bias and limits of agreement.1 an error percentage , calculated according to   below 30 was accepted as sign of good agreement . calculated and measured values were analyzed by means of regression and r. patient and volunteer characteristics and results are reported in table 1 . table 1patients and volunteers entered in evaluation of step response derivationpatients i d , sex and agedescriptionpeep ( cm h20)vat ( ml)eigenvaluesbl f102 ; fn2frc ( ml)#biasupper & lower limit of agreementslope;constant030930 , m , 76icu , sepsis , copd53750.790.4148070.822.13 & 3.761.15;6.49100.840.12090150.880.22346200.913045031110 , m , 65icu , sepsis , copd122350.890.4185060.212.32 & 2.71.13;5.35180.900.122000.20.3031120 , m , 48icu , hf102160.850.4122860.042.1 & 2.21.09;4.46150.860.113360.20.3040216 , f , 42icu , ards , sle102800.720.560080080.040.7 & 0.630.93;2.7150.770.11040040415 , f , 35laparoscopic , bariatric surgery04350.640.479860.081.42 & 1.260.93;3.34100.670.5950200.740.11200all patients350.2971.989 & 2.581.07;3.05volunteers i d , sex and ageco , f , 40sitting , standing0variable0.830.2250060.011.12 & 1.20.99;1.19100.880.20.32900coh , m , 29sitting0variable0.740.2140041.467.6 & 4.70.98;0.4550.830.22400100.850.452668ab , m , 38supine0variable0.770.20.3175060.351.7 & 2.41.03;1.1850.790.22184100.770.352125all volunteers160.163.5 & 3.21.02;0.97age , sex , location of measurement , peep and fn2 changes , baseline f1o2 and frc . patients and volunteers entered in evaluation of step response derivation age , sex , location of measurement , peep and fn2 changes , baseline f1o2 and frc . the first order difference equation based on mass conservation of nitrogen over one breath after a step change in inhaled n2 fraction is described in equation ( 3 ) . \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{v}}_{{\rm{ti}}}}{\rm { \times } } { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{_{{\rm{low}}}}{\rm { + frc } } \times { { \rm{f}}_{\rm{a}}}{{\rm{n}}_{\rm{2}}}_{_{n{\rm { - 1}}}}{\rm { = } } { { \rm{v}}_{{\rm{te}}}}{\rm { \times } } { { \rm{f}}_{{\rm{\bar e}}}}{{\rm{n}}_{\rm{2}}}_{_n}{\rm { + frc \times } } { { \rm{f}}_{\rm{a}}}{{\rm{n}}_{\rm{2}}}_{_n}$$ 
 \end{document } the step change in nitrogen occurs at breath index n = 0 . the step is from \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{_{{\rm{high}}}}$$ 
 \end{document } to \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{_{{\rm{low}}}}$$ 
 \end{document } which will lead to n2 washing out . \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{{\rm{\bar e}}}}{{\rm{n}}_{\rm{2}}}$$ 
 \end{document } is an average breath fraction of expired gas which includes gas from anatomical ( series ) dead space . \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{\rm{a}}}{{\rm{n}}_{\rm{2}}}_{_n}$$ 
 \end{document } depends on the fraction of the immediately preceding breath \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{\rm{a}}}{{\rm{n}}_{\rm{2}}}_{_{n - 1}}$$ 
 \end{document}. vti and vte designate inspiratory and expiratory tidal volumes . in all equations a refinement is to measure the lung volume to the exclusion of tidal anatomical dead space , vdtanat . thus , the vte and vti of equation ( 3 ) must be replaced by the alveolar tidal volumes , vate and vati . also , using endtidal co2 fraction , fetco2 , as approximating alveolar fraction , faco2 , the following equations can be used to calculate the alveolar tidal volumes : 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { \rm{v}}{{\rm{a}}_{{\rm{te}}}}{\rm { = } } { { { { \rm{v}}_{\rm{t}}}{\rm{c}}{{\rm{o}}_{\rm{2 } } } } \over { { { \rm{f}}_{{\rm{et}}}}{\rm{c}}{{\rm{o}}_{\rm{2}}}}}$$ 
 \end{document } , where vtco2 is the tidal excretion of co2 . \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { \rm{v}}{{\rm{a}}_{{\rm{ti}}}}{\rm { = v}}{{\rm{a}}_{{\rm{te}}}}{\rm { + } } \left ( { { { { { \rm{v}}_{\rm{t}}}{\rm{c}}{{\rm{o}}_{\rm{2 } } } } \over { { \rm{rq}}}}{\rm { - } } { { \rm{v}}_{\rm{t}}}{\rm{c}}{{\rm{o}}_{\rm{2 } } } } \right)$$ 
 \end{document } rq is the respiratory quotient defined as \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { \rm{rq = } } { { { { \rm{v}}_{\rm{t}}}{\rm{c}}{{\rm{o}}_{\rm{2 } } } } \over { { { \rm{v}}_{\rm{t}}}{{\rm{o}}_{\rm{2}}}}}$$ 
 \end{document}. the value of rq depends on the composition of substrate being metabolized . for example , rq = 1 when carbohydrates are being metabolized and rq = 0.7 when fat is being metabolized . the range of rq is typically from 0.7 to 1.0 with reservation for gluconeo- and lipogenesis . a fixed value of rq of 0.8 is used in equation ( 5 ) . for the purpose of accumulating the in- and expired volumes of nitrogen , the fractions of nitrogen are calculated as residue . \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}{\rm { = 1 - } } { { \rm{f}}_{\rm{i}}}{{\rm{o}}_{\rm{2}}}$$ 
 \end{document}\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{\rm{2}}}{\rm { = 1 - } } { { \rm{f}}_{{\rm{et}}}}{\rm{c}}{{\rm{o}}_{\rm{2}}}{\rm { - } } { { \rm{f}}_{{\rm{et}}}}{{\rm{o}}_{\rm{2}}}$$ 
 \end{document } substituting for alveolar volumes and endtidal fractions , equation ( 3 ) becomes : 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { \rm{v}}{{\rm{a}}_{{\rm{ti}}}}{\rm { \times } } { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{_{{\rm{low}}}}{\rm { + frc \times } } { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{\rm{2}}}_{_{n - 1 } } = { \rm{v}}{{\rm{a}}_{{\rm{te}}}}{\rm { \times } } { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{\rm{2}}}_n + { \rm{frc \times } } { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{\rm{2}}}_{_n}$$ 
 \end{document } this is a first order difference equation . the fin2 and fetn2 are based on samples from the continuous fco2 and fo2 signals according to equations 6 and 7 , see figure 1 . vati and vate are calculated according to equations ( 4 ) and ( 5 ) . 1sampling at two frequencies , the waveform sampling rate and two samples per breath sampling rate . lower frame : the continuous sampling of fo2 ( left y - axis ) and fco2 ( right y - axis ) during a n2 washout procedure . from these tracings ( enlarged in upper frame ) two samples of fo2 and fco2 are sampled once during inspiration and as endtidal values for entering into equations 6 and 7 . sampling at two frequencies , the waveform sampling rate and two samples per breath sampling rate . lower frame : the continuous sampling of fo2 ( left y - axis ) and fco2 ( right y - axis ) during a n2 washout procedure . from these tracings ( enlarged in upper frame ) two samples of fo2 and fco2 are sampled once during inspiration and as endtidal values for entering into equations 6 and 7 . in the presence of alveolar ( parallel ) dead space , fetco2 < faco2 and , as a result , the equations include alveolar dead space . first , the frc is found by summing both sides of equation 8 over the total number of breaths , 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { \rm{frc } } = { { \sum ( { \rm{v}}{{\rm{a}}_{{\rm{te}}}}{\rm { \times } } { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{\rm{2}}}{\rm { - v}}{{\rm{a}}_{{\rm{ti}}}}{\rm { \times } } { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{_{{\rm{low } } } } ) } \over { { { \rm{f}}_{{\rm{etinit}}}}{{\rm{n}}_{\rm{2}}}{\rm { - } } { { \rm{f}}_{{\rm{etlatest}}}}{{\rm{n}}_{\rm{2 } } } } } = { { \sum ( { \rm{v}}{{\rm{a}}_{\rm{e}}}{{\rm{n}}_{\rm{2}}}{\rm { - v}}{{\rm{a}}_{\rm{i}}}{{\rm{n}}_{\rm{2 } } } ) } \over { { { \rm{f}}_{{\rm{etinit}}}}{{\rm{n}}_{\rm{2}}}{\rm { - } } { { \rm{f}}_{{\rm{etlatest}}}}{{\rm{n}}_{\rm{2 } } } } } = { { \sum { \rm{\delta v}}{{\rm{a}}_{\rm{n}}}_{_{\rm{2 } } } } \over { { \rm{\delta } } { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{\rm{2}}}}}$$ 
 \end{document}. fetinitn2 denotes the initial condition at time n equals 1 which is the time sample before the step ( \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{{\rm{etinit}}}}{{\rm{n}}_{\rm{2}}}{\rm { = } } { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{\rm{2}}}_{_{{\rm { - 1}}}}$$ 
 \end{document } ) . fetlatestn2 represents \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{\rm{2}}}_{_n}$$ 
 \end{document } when n equals the most recent breath . in general , vate , vati and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{_{{\rm{low}}}}$$ 
 \end{document } can vary with n and may be replaced by \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { \rm{v}}{{\rm{a}}_{{\rm{te}}}}_n$$ 
 \end{document } , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } \begin{document } 
  $ $ { \rm{v}}{{\rm{a}}_{{\rm{ti}}}}_{_n}$$ 
 \end{document } and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{_n}$$ 
 \end{document } respectively . second , equation 8 is solved for the step response using the frc calculated by equation ( 9 ) 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ \eqalign { & { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{{2_n } } } = { { { \rm{v}}{{\rm{a}}_{{\rm{ti}}}}{\rm { \times } } { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{_{{\rm{low } } } } } \over { { \rm{v}}{{\rm{a}}_{{\rm{te}}}}{\rm { + frc } } } } \times \left ( { { { 1 - { \lambda^{n + 1 } } } \over { 1 - \lambda } } } \right ) + { \lambda^{n + 1 } } \times { { \rm{f}}_{{\rm{etinit}}}}{{\rm{n}}_2 } \cr & \lambda = { { { \rm{frc } } } \over { { \rm{frc + v}}{{\rm{a}}_{{\rm{te } } } } } } \cr } $ $ \end{document } in this solution , the values , vati , vate and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{_{{\rm{low}}}}$$ 
 \end{document } are assumed to remain constant during the washout . for a full derivation , see appendix a. this is an equation that gives the fetn2 for a particular breath without iteration and can be used to determine the applicability of a first order model . some implications of equation 10 are immediately apparent : 
  is related to the time constant of the n2 clearance.as the vate increases , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{\rm{2}}}_{_n}$$ 
 \end{document } decays faster consistent with a faster lung washout at larger tidal breaths.if the frc becomes large compared to the vate , then  approaches 1 and there is very slow decay consistent with a small tidal breath requiring more breaths to change the concentration of a large lung volume ( frc).as n   , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_2}_n \to { { { \rm{v}}{{\rm{a}}_{{\rm{ti } } } } \cdot { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{{\rm{low } } } } \over { { \rm{v}}{{\rm{a}}_{{\rm{te}}}}}}$$ 
 \end{document}. the \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{\rm{2}}}_{_n}$$ 
 \end{document } will not reach \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{_{{\rm{low}}}}$$ 
 \end{document } when rq  1 since vati  vate . one would expect a good match between the measured fetn2 sequence compared to the calculated values if the system being measured is indeed first order . is related to the time constant of the n2 clearance . as the vate increases , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{\rm{2}}}_{_n}$$ 
 \end{document } decays faster consistent with a faster lung washout at larger tidal breaths . if the frc becomes large compared to the vate , then  approaches 1 and there is very slow decay consistent with a small tidal breath requiring more breaths to change the concentration of a large lung volume ( frc ) . as n   , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_2}_n \to { { { \rm{v}}{{\rm{a}}_{{\rm{ti } } } } \cdot { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{{\rm{low } } } } \over { { \rm{v}}{{\rm{a}}_{{\rm{te}}}}}}$$ 
 \end{document}. the \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } $ $ { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_{\rm{2}}}_{_n}$$ 
 \end{document } will not reach \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{_{{\rm{low}}}}$$ 
 \end{document } when rq  1 since vati  vate . one would expect a good match between the measured fetn2 sequence compared to the calculated values if the system being measured is indeed first order . the modified mbnw has been extensively validated in a metabolic lung model . in this paper , the prediction of fetn2 by the frc step response equation was validated by comparing calculated to measured values of fetn2 from washout procedures in one experimental and two clinical settings . the clinical measurements were made as part of two other studies for which the ethical committee of the medical faculty at gteborg university approved the studies and informed consent was obtained from the patients , next of kin or volunteers , respectively . thus , they represent examples and are not the result of a selection powered to demonstrate differences between the normal and disease categories . model . an experimental setup using an o2-consuming / co2-producing metabolic lung model previously described [ 10 , 13 ] . the model was ventilated in volume controlled ventilation with a ge engstrm carestation ventilator ( madison , wi , usa ) . settings comprised two frcs ( 1.8 , 2.9 l ) , two tidal volumes ( 500 , 750 ml ) and two steps of fio2 ( 0.1 , 0.3 ) at \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { \rm{\dot vc}}{{\rm{o}}_2}$$ 
 \end{document } 200 ml / min and rq 1 . frc measurements were computed with the commercial software algorithm based on equation 9 detailed above . patient recordings were performed either in the operating theater using ge adu/5 ventilator or in the intensive care unit ( icu ) using a servo 900c or servo 300 ventilator ( siemens elema , solna , sweden ) . subjects with assumed normal lungs were measured peroperatively and patients with diseased lungs in the intensive care unit , where positive pressure ventilation was instituted due to respiratory failure . the eigenvalues ,  , were calculated for one patient with acute respiratory distress syndrome ( ards ) , one with heart failure ( hf ) , two patients with chronic obstructive pulmonary disease ( copd ) , and one obese patient with normal lungs ( perioperative registration ) . mono- vs biexponential fits were examined in one of the icu patients ( copd ) and in the bariatric patient . recordings in healthy volunteers ( two female and one male ) spontaneously breathing while attached to a high flow cpap system via a snorkeling mouthpiece and equipped with a nose clip in supine , sitting and standing position . in patients and volunteers , respiratory monitoring and gas sampling were performed with datex - ohmeda spirometry module in an as/3 intensive care monitor . data were collected using datex - ohmeda s/5 collect version 4 software ( datex - ohmeda , helsinki , finland ) and analyzed using lung monitor ( plug - in to datex - ohmeda s/5 collect ) or proprietary software programmed in test- point ( measurement computing corporation , norton , ma , usa ) . agreement between calculated and measured values of fetn2 was assessed by bland & altman analysis   in terms of bias and limits of agreement.1 an error percentage , calculated according to   below 30 was accepted as sign of good agreement . calculated and measured values were analyzed by means of regression and r. patient and volunteer characteristics and results are reported in table 1 . table 1patients and volunteers entered in evaluation of step response derivationpatients i d , sex and agedescriptionpeep ( cm h20)vat ( ml)eigenvaluesbl f102 ; fn2frc ( ml)#biasupper & lower limit of agreementslope;constant030930 , m , 76icu , sepsis , copd53750.790.4148070.822.13 & 3.761.15;6.49100.840.12090150.880.22346200.913045031110 , m , 65icu , sepsis , copd122350.890.4185060.212.32 & 2.71.13;5.35180.900.122000.20.3031120 , m , 48icu , hf102160.850.4122860.042.1 & 2.21.09;4.46150.860.113360.20.3040216 , f , 42icu , ards , sle102800.720.560080080.040.7 & 0.630.93;2.7150.770.11040040415 , f , 35laparoscopic , bariatric surgery04350.640.479860.081.42 & 1.260.93;3.34100.670.5950200.740.11200all patients350.2971.989 & 2.581.07;3.05volunteers i d , sex and ageco , f , 40sitting , standing0variable0.830.2250060.011.12 & 1.20.99;1.19100.880.20.32900coh , m , 29sitting0variable0.740.2140041.467.6 & 4.70.98;0.4550.830.22400100.850.452668ab , m , 38supine0variable0.770.20.3175060.351.7 & 2.41.03;1.1850.790.22184100.770.352125all volunteers160.163.5 & 3.21.02;0.97age , sex , location of measurement , peep and fn2 changes , baseline f1o2 and frc . results of bland & altman analysis of agreement and regression analysis . patients and volunteers entered in evaluation of step response derivation age , sex , location of measurement , peep and fn2 changes , baseline f1o2 and frc . model . there was an excellent correspondence between 178 pairs of measured and calculated values of fetn2 in six washout procedures using the step response equation , see figure 2 . bland & altman analysis showed bias of 0 , upper and lower limits of agreement of 0.49 and 0.5 and error 1.3% , see figure 3 . the regression equation showed a slope of 1.0 , intercept 0.12 and r 0.98 with a p < 0.001 . fn2 0.3 , frc 1800 ml , vat 270 ml , eigenvalue 0.86 , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { \rm{\dot vc}}{{\rm{o}}_2}$$ 
 \end{document } 187 ml / min . 3bland & altman plot of agreement between pairs of measured and calculated fetn2 according to step response derivation in metabolic lung model . bias 0 , upper and lower limits of agreement of 0.25 and 0.5 . calculated and measured fetn2 during n2 washout in metabolic lung model . fn2 0.3 , frc 1800 ml , vat 270 ml , eigenvalue 0.86 , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { \rm{\dot vc}}{{\rm{o}}_2}$$ 
 \end{document } 187 ml / min . the regression equation has a slope of 1.02 and a r of 0.999 . bland & altman plot of agreement between pairs of measured and calculated fetn2 according to step response derivation in metabolic lung model . bias 0 , upper and lower limits of agreement of 0.25 and 0.5 . patients . peep levels varied from 0 to 20 cm h2o , frc ranged from 600 to 3000 ml , fetn2 from 0.1 to 0.3 and baseline fio2 from 0.2 to 0.5 . there was good agreement between 1222 pairs of measured and the calculated values of fetn2 . bland & altman analysis showed bias of 0.3 , upper and lower limits of agreement of 2 and 2.6 , error 4.9% , see figure 4 . the regression equation showed a slope of 1.07 , intercept 3.0 , r 0.98 with a p < 0.001 . at comparable peep levels  in ards ( ards ) patient was 0.734 , in patient with heart failure ( hf ) 0.846 , in copd patients 0.842 and 0.889 ( copd ) and in perioperative normal adult lung ( bariatric patient ) 0.654 ( periop ) . 4bland & altman plot of agreement between pairs of measured and calculated fetn2 according to step response derivation in 35 n2 washout procedures in five patients . the regression equation showed a slope of 1.07 , intercept 3 and r 0.98 with a p < 0.001.fig . 5 as a function of peep in five patients . bland & altman plot of agreement between pairs of measured and calculated fetn2 according to step response derivation in 35 n2 washout procedures in five patients . the regression equation showed a slope of 1.07 , intercept 3 and r 0.98 with a p < 0.001 . as a function of peep in five patients . in the analysis of mono- vs. bicompartmental analysis we found a difference between mono- and biexponential fit to measured fetn2 in the icu ( copd ) patient whereas there was no notable difference in the bariatric patient . the number of patients and measurements preclude consideration statistical significance . volunteers . peep levels varied from 0 , 5 to 10 cm h2o , frc ranged from 1400 to 2900 ml , fetn2 0.2 to 0.45 and baseline fio2 from 0.21 to 0.3 . there was good agreement between 524 pairs of measured and the calculated values of fetn2 using the step response equation . bland & altman analysis showed bias of 0.16 , upper and lower limits of agreement at 3.2 and 3.5 and an error 6.4% . the regression equation showed a slope of 0.97 , intercept 1.9 and r 0.98 with a p < 0.001 ( figure 6 ) . 6bland & altman plot of agreement between pairs of measured and calculated fetn2 according to step response derivation in 16 n2 washout procedures in three healthy volunteers . bland & altman analysis showed bias of 0.16 , lower and upper limits of agreement at 3.2 and 3.5 . the regression equation showed a slope of 0.97 , intercept at 1.9 and r 0.98 with a p < 0.001 . bland & altman plot of agreement between pairs of measured and calculated fetn2 according to step response derivation in 16 n2 washout procedures in three healthy volunteers . bland & altman analysis showed bias of 0.16 , lower and upper limits of agreement at 3.2 and 3.5 . the regression equation showed a slope of 0.97 , intercept at 1.9 and r 0.98 with a p < 0.001 . we take as a starting point the mathematical derivation of a modified version of the multiple breath nitrogen washout method for calculation of frc . using the result of the frc calculation , a step response equation was derived based on techniques from dsp and used to calculate the course of fetn2 during washout . these calculated values were then compared to measured data from an experimental setup , patients and volunteers . we found excellent alignment between measured and calculated nitrogen clearance using a single compartment , first order metabolic lung model . the first order model also fitted measurements in patients with various respiratory conditions and in healthy , spontaneously breathing volunteers . the eigenvalue of the step response was shown to be dependent on peep level and to be different among healthy and diseased states of the lung . the sample size , however , does not allow for any statistically significant inferences , but the results may indicate a diagnostic and prognostic method to be verified in further studies . technical and clinical characteristics and limitations of the modified mbnw method and first order step response equation are addressed in the following . the present study is  to our knowledge  the first attempt at deriving the equivalent of the alveolar dilution ratio first introduced by darling   by solving the underlying first order difference equation based on mass conservation of alveolar contents of n2 . the eigenvalue ,  , is equal to darling s r in equation ( 2 ) . this ratio was revived by wagner   as a term in linear programming of ventilation / volume ratios and by rossing . wagner and rossing published their research before oppenheim and schaffer published their work   on dsp , which forms the basis of this study . using techniques made popular by the field of dsp , enables one to frame the problem in the form of a difference equation and solve for the step response . dsp has entered the mainstream subsequent to many of the early works on mbnw . reformulating the subject in a dsp framework the present study details nitrogen clearance as a first order process characterized by an eigenvalue ,  and compares actual values of fetn2 to values predicted by the step response equation . the degree of agreement confirms that the nitrogen washout is dominated by a first order process and the tools of dsp are applicable in this setting . clinical impact and uses for . this paper provides an alternative derivation for  in the more general framework of dsp . results in the present study showed higher  in healthy adult lungs and in patients with copd and lower values in patients with ards . in addition ,  is reminiscent of the strain relationship ( vt / frc ) which according to the  baby lung2 the bariatric patient ( see figure 5 ) the strain is perilously close to this limit . it may not be clinically justifiable to measure  by setting peep at zero ( zeep ) , but as the value is linearly related to peep by virtue of the lung compliance , a series measurement at peep levels above zero can be used to extrapolate to a value of  corresponding to zeep . the  or the strain relationship could be considered a discriminator for classifying lung states but may not be sufficient in isolation . since  has a dependency on the ventilator settings through frc ( and as a corollary peep ) and vate , calculation of  at several tidal volumes would likely be a better discriminate . further , these ratios in conjunction with other features such as hemodynamics would be an excellent start for a feature vector that could be used for pattern classification . the pattern classifier would be used to distinguish between different lung states and act as a guide for optimal ventilation therapy . the washout procedure posits that the o2-consumption and co2-excretion , the end - expiratory lung volume as well as hemodynamics are stable at baseline before washout and that no interventions are undertaken during the washout period which will alter this , e.g. changes in inotropic support with inherent thermogenic effects or changes in ventilatory settings . in our formulation , rq is fixed at 0.8 , although in reality rq can be between 0.7 and 1.0 depending on the metabolic state of the patient . in this context it can be shown that for a fio2 of 0.5 , a fetn2 of 0.1 , a rq of 0.1 and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { \rm{\dot vc}}{{\rm{o}}_2}$$ 
 \end{document } 200 ml / min the error in frc amounts to approximately 180 ml bias . equality in equation 3 assumes that n2 is neither taken up nor excreted from blood and tissues . estimates of tissue nitrogen uptake and excretion do , however , exist [ 18 , 19 ] . the modified mbnw utilizes a fn2 of 0.1 and accordingly the volume excreted will be diminished to approximately 10 ml . however , repeating the fn2 with an identical change in opposite direction and averaging the results for frc compensates for this error . the fin2 and fetn2 are calculated as residue after subtraction of appropriate fractions of o2 and co2 and entered into the calculation of vn2 by multiplication with va . this assumes that the dead spaces of o2 , co2 and n2 are identical within the clinical range of expired volume , which is supported by the findings in [ 20 , 21 ] . furthermore , feto2 and fetco2 measurements are used to approximate fao2 and faco2 although there is evidence to the contrary in the case of fetco2 [ 22 , 23 ] . typically , the fetco2 will be lower than faco2 and feto2 the modified mbnw method incorporates breath - by - breath measurement and calculation of vati , and vate allowing for variation of tidal volumes during the washout . in the step response derivation the discrepancy becomes obvious when examining frc measurements performed during spontaneous ventilation where breath - to - breath variation of inspiratory and expiratory tidal volumes in quiet breathing typically amounted to 50100 ml during the washout procedure . the potential of using  from the step response solution of the difference equation to classify lung pathologies applies more to mechanical ventilation where the volumes are constant . washout was clearly not monoexponential in the case of emphysematous ( copd ) lung with compartments of differing mechanical properties . a better curve fit was accomplished by a biexponential equation , which in the framework of dsp translates into a second order difference equation . the modified mbnw method based on a first order difference equation , however , will still be able to calculate the frc given the measurement period is extended to allow for slow compartments to empty in the case of copd lungs . the first order difference equation solution to a step input assumes that washout is monoexponential and deviations from this will cause increased lack of agreement between measured and calculated values . this was demonstrated in the fitting of washout data to a mono- vs. a biexponential equation in the case of the copd / icu and the perioperative bariatric patient . the present study is  to our knowledge  the first attempt at deriving the equivalent of the alveolar dilution ratio first introduced by darling   by solving the underlying first order difference equation based on mass conservation of alveolar contents of n2 . the eigenvalue ,  , is equal to darling s r in equation ( 2 ) . this ratio was revived by wagner   as a term in linear programming of ventilation / volume ratios and by rossing . wagner and rossing published their research before oppenheim and schaffer published their work   on dsp , which forms the basis of this study . using techniques made popular by the field of dsp , enables one to frame the problem in the form of a difference equation and solve for the step response . dsp has entered the mainstream subsequent to many of the early works on mbnw . reformulating the subject in a dsp framework should allow algorithms from that field to improve or extend mbnw algorithms . the present study details nitrogen clearance as a first order process characterized by an eigenvalue ,  and compares actual values of fetn2 to values predicted by the step response equation . the degree of agreement confirms that the nitrogen washout is dominated by a first order process and the tools of dsp are applicable in this setting . clinical impact and uses for . the term  = frc/(frc+vate ) this paper provides an alternative derivation for  in the more general framework of dsp . results in the present study showed higher  in healthy adult lungs and in patients with copd and lower values in patients with ards . in addition ,  is reminiscent of the strain relationship ( vt / frc ) which according to the  baby lung2 the bariatric patient ( see figure 5 ) the strain is perilously close to this limit . it may not be clinically justifiable to measure  by setting peep at zero ( zeep ) , but as the value is linearly related to peep by virtue of the lung compliance , a series measurement at peep levels above zero can be used to extrapolate to a value of  corresponding to zeep . the  or the strain relationship could be considered a discriminator for classifying lung states but may not be sufficient in isolation . since  has a dependency on the ventilator settings through frc ( and as a corollary peep ) and vate , calculation of  at several tidal volumes would likely be a better discriminate . further , these ratios in conjunction with other features such as hemodynamics would be an excellent start for a feature vector that could be used for pattern classification . the pattern classifier would be used to distinguish between different lung states and act as a guide for optimal ventilation therapy . the washout procedure posits that the o2-consumption and co2-excretion , the end - expiratory lung volume as well as hemodynamics are stable at baseline before washout and that no interventions are undertaken during the washout period which will alter this , e.g. changes in inotropic support with inherent thermogenic effects or changes in ventilatory settings . in our formulation , rq is fixed at 0.8 , although in reality rq can be between 0.7 and 1.0 depending on the metabolic state of the patient . in this context it can be shown that for a fio2 of 0.5 , a fetn2 of 0.1 , a rq of 0.1 and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { \rm{\dot vc}}{{\rm{o}}_2}$$ 
 \end{document } 200 ml / min the error in frc amounts to approximately 180 ml bias . equality in equation 3 assumes that n2 is neither taken up nor excreted from blood and tissues . estimates of tissue nitrogen uptake and excretion do , however , exist [ 18 , 19 ] . the modified mbnw utilizes a fn2 of 0.1 and accordingly the volume excreted will be diminished to approximately 10 ml . however , repeating the fn2 with an identical change in opposite direction and averaging the results for frc compensates for this error . the fin2 and fetn2 are calculated as residue after subtraction of appropriate fractions of o2 and co2 and entered into the calculation of vn2 by multiplication with va . this assumes that the dead spaces of o2 , co2 and n2 are identical within the clinical range of expired volume , which is supported by the findings in [ 20 , 21 ] . furthermore , feto2 and fetco2 measurements are used to approximate fao2 and faco2 although there is evidence to the contrary in the case of fetco2 [ 22 , 23 ] . typically , the fetco2 will be lower than faco2 and feto2 will be higher than fao2 due to imperfect mixing and alveolar dead space . the modified mbnw method incorporates breath - by - breath measurement and calculation of vati , and vate allowing for variation of tidal volumes during the washout . in the step response derivation the discrepancy becomes obvious when examining frc measurements performed during spontaneous ventilation where breath - to - breath variation of inspiratory and expiratory tidal volumes in quiet breathing typically amounted to 50100 ml during the washout procedure . the potential of using  from the step response solution of the difference equation to classify lung pathologies applies more to mechanical ventilation where the volumes are constant . washout was clearly not monoexponential in the case of emphysematous ( copd ) lung with compartments of differing mechanical properties . a better curve fit was accomplished by a biexponential equation , which in the framework of dsp translates into a second order difference equation . the modified mbnw method based on a first order difference equation , however , will still be able to calculate the frc given the measurement period is extended to allow for slow compartments to empty in the case of copd lungs . the first order difference equation solution to a step input assumes that washout is monoexponential and deviations from this will cause increased lack of agreement between measured and calculated values . this was demonstrated in the fitting of washout data to a mono- vs. a biexponential equation in the case of the copd / icu and the perioperative bariatric patient . we have presented a modified mbnw method in terms of a first order difference equation which was solved for the step response of fetn2 to a step change in fio2 . a central concept of the step response equation , the eigenvalue , showed different values in healthy and diseased lungs . the mbnw difference equation has some interesting implications for the roles of serial and parallel dead spaces . the step response equation opens up research avenues within digital signal processing and the characterization of ventilatory states . the derivation of the step response solution is as follows . in equation ( 8) , define the step response as \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ a = { \rm{v}}{{\rm{a}}_{{\rm{ti}}}}{\rm { \times f}}{f_{\rm{i}}}{{\rm{n}}_2}_{_{{\rm{low}}}}$$ 
 \end{document } for n  0 , e vate + frc , b = frc , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_n } = { { \rm{f}}_{{\rm{et}}}}{{\rm{n}}_2}_{_n}$$ 
 \end{document } and the initial condition as fetinitn2 : utilizing a substitution technique : 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ n = 0:{{\rm{f}}_0 } = \left ( { a + b \times { { \rm{f}}_{{\rm{etinit}}}}{\rm{n}}2 } \right)/e$$ 
 \end{document}\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ n = 1:{{\rm{f}}_{\rm{1 } } } = a / e + b \times \left ( { \left ( { a + b \times { { \rm{f}}_{{\rm{etinit}}}}{{\rm{n}}_2 } } \right)/{e^2 } } \right)$$ 
 \end{document}\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ n = 2:{{\rm{f}}_2 } = a / e + b \times a/{e^2 } + a \times { b^2}/{e^3 } \times + { { \rm{f}}_{{\rm{etinit}}}}{{\rm{n}}_2 } \times { b^3}/{e^3}$$ 
 \end{document } in general , 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_n } = { a \over e}\mathop \sum \limits_{i = 0}^n { \left ( { { b \over e } } \right)^i } + { \left ( { { b \over e } } \right)^{n + 1 } } \times { { \rm{f}}_{{\rm{etinit}}}}{{\rm{n}}_2}$$ 
 \end{document } let  = b / e and equation 13 simplifies into : 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_n } = { a \over e}\left ( { { { 1 - { \lambda^{n + 1 } } } \over { 1 - \lambda } } } \right ) + { \lambda^{n + 1 } } \times { { \rm{f}}_{{\rm{etinit}}}}{{\rm{n}}_2}$$ 
 \end{document } substitute back to the original variables : 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document } 
 $ $ { { \rm{f}}_n } = { { { \rm{v}}{{\rm{a}}_{{\rm{ti}}}}{\rm { \times } } { { \rm{f}}_{\rm{i}}}{{\rm{n}}_{\rm{2}}}_{_{{\rm{low } } } } } \over { { \rm{v}}{{\rm{a}}_{{\rm{te } } } } + { \rm{frc } } } } \times \left ( { { { 1 - { \lambda^{n + 1 } } } \over { 1 - \lambda } } } \right ) + { \lambda^{n + 1 } } \times { { \rm{f}}_{{\rm{etinit}}}}{{\rm{n}}_2};\lambda = { { { \rm{frc } } } \over { { \rm{frc + v}}{{\rm{a}}_{{\rm{te}}}}}}$$ 
 \end{document }
objectiveto derive a difference equation based on mass conservation and on alveolar tidal volumes for the calculation of functional residual capacity . 
 derive an equation for the frc from the difference equation . 
 furthermore , to derive and validate a step response equation as a solution of the difference equation within the framework of digital signal processing where the frc is known a priori.methodsa difference equation for the calculation of functional residual capacity is derived and solved as step response of a first order system . 
 the step response equation calculates endtidal fractions of nitrogen during multiple breath nitrogen clearance . 
 the step response equation contains the eigenvalue defined as the ratio of frc to the sum of frc and alveolar tidal ventilation . 
 agreement of calculated nitrogen fractions with measured fractions is demonstrated with data from a metabolic lung model , measurements from patients in positive pressure ventilation and volunteers breathing spontaneously . 
 examples of eigenvalue are given and compared between diseased and healthy lungs and between ventilatory settings.resultscomparison of calculated and measured fractions of endtidal nitrogen demonstrates a high degree of agreement in terms of regression and bias and limits of agreement ( precision ) in bland & altman analysis . 
 examples illustrate the use of the eigenvalue as a possible discriminator between disease states.conclusionthe first order step response equation reliably calculates endtidal fractions of nitrogen during washout based on a functional residual capacity . 
 the eigenvalue may be a clinically valuable index alone or in conjunction with other indices in the analysis of respiratory states and may aid in the setting of the ventilator .