Mu and lambda for inverse Gaussian distributed values.

procedure InverseGaussianFit(const X: TVec; out mu: double; out lambda: double; var muConfInt: TTwoElmReal; var lambdaConfInt: TTwoElmReal; Alpha: double = 0.05); overload;