o Simulation of OU process with rTraitCont() did not work correctly.
This now uses formula from Gillespie (1996) reduced to a BM
- process when alpha = 0 avoiding division by zero.
+ process when alpha = 0 to avoid division by zero. The option
+ 'linear' has been removed.