sigma_mle function

Computes the ith iteration for computing the squared imprecision estimates.