Ajout d'options, algorithme stable pour dark soliton th(x/sqrt(2))