Implementation of Wilkinson formulas
