simulate_bm_model function

Simulate a Brownian motion model for multivariate trait co-evolution.