non phylo model implementation
figure p330 partiel