Cleaned-up a bit the code. Calculation of the best-fitting correction to the bias can now be implemented.