Il faut utiliser la méthode BLS décrite précédemment. Soit vous la codez vous-même, soit vous utilisez une fonction pré-existente qui intègre déjà cet algorithme BLS. Il y a de telles fonctions en python ou fortran par exemple.