The method was first described by Cunningham & Jones in 1993 for the analysis of dynamic PET data obtained in the brain. It assumes that the tissue impulse response function (IRF) can be described as a combination of many exponentials. Since A tissue TAC can be expressed as a convolution of measured arterial input function with IRF, C(t) can be expressed as:
where,  is a convolution operator, C(t) is the bone tissue activity concentration of tracer (in units: MBq/ml) over a period of time t, C(t) is the plasma concentration of tracer (in units: MBq/ml) over a period of time t, IRF(t) is equal to the sum of exponentials, β values are fixed between 0.0001 sec and 0.1 sec in intervals of 0.0001, n is the number of α components that resulted from the analysis and β, β,..., β corresponds to the respective α, α,..., α components from the resulted spectrum. The values of α are then estimated from the analysis by fitting multi-exponential to the IRF. The intercept of the linear fit to the slow component of this exponential curve is considered the plasma clearance (K) to the bone mineral.