fit_env function

Maximum likelihood fit of the environmental birth-death model