sim_env_bd function

Simulate birth-death tree dependent on an environmental curve