Hi there,

I have noticed the following problem:

```
sage: f = function('f')(x)
sage: var('h')
sage: integrate(exp(h)*exp(x)*f(x),x)
integrate(e^(h + x)*f(x), x)
```

The workaround seems to be using the `sympy`

backend for symbolic integration

```
sage: integrate(exp(h)*exp(x)*f(x),x,algorithm='sympy')
e^h*integrate(e^x*f(x), x)
```

which always seems to be a good idea as I learned from @Emmanuel Charpentier over here.

Now I would like to force the use of `algorith='sympy'`

for simplifying these `integrate(...)`

expressions globally.
Unfortunately, the `simplify()`

command does not allow to set this option.

```
sage: integrate(exp(h)*exp(x)*f(x),x)
integrate(e^(h + x)*f(x), x)
sage: _.simplify()
integrate(e^(h + x)*f(x), x)
```

**TL;DR** How can I force sage to pull out these type of exponential constants from the integral with the `simplify()`

command?