fastBM(tree,a=0,mu=0,sig2=1,bounds=c(-Inf,Inf),internal=FALSE,nsim=1)"phylo" format.mu!=0.nsim=1) or matrix containing the tip states for the n species in the tree, and (optionally) the ancestral states for internal nodes.branching.diffusion