'Compute unit specific parameters in
Laird, Lange, Stram (see help(Unit)).
Displays (3.2)-(3.5).'

def fit(self, a, D, sigma):

    self._compute_S(D, sigma)
    self._compute_W()
    self._compute_r(a)
    self._compute_b(D)
