Update new_phylo.R

added error check
correctly can specify number of root nodes