first_step function

Estimate the first-order inhomogeneity