# sqrt integral not working

I can't seem to integrate sqrt in SageMath -

Jupyter cells in website CoCalc

```
In:
x1 = var('x1')
g1 = function('g11')
s1 = sqrt(g1(x1)); s1
g0i = integral(s1,x1); g0i
Out:
+Infinity
```

But sqrt works fine otherwise -

```
In:
x1 = var('x1')
g1 = function('g11')
s1 = sqrt(g1(x1)); s1
Out:
sqrt(g11(x1))
```

And sqrt works in other expressions, e.g. inverse -

```
In:
x1 = var('x1')
g1 = function('g11')
s1 = sqrt(g1(x1)); s1
g0i = integral(s1/g1(x1),x1); g0i
Out:
integrate(1/sqrt(g11(x1)), x1) ---- note the processing of
sqrt for the output
```

But when the final expression depends on the first power of sqrt -

```
In:
x1 = var('x1')
g1 = function('g11')
s1 = sqrt(g1(x1)); s1
g0i = integral(g1(x1)/s1,x1); g0i
Out:
+Infinity
```

I'm interested in processing in symbolic math, not in specific numbers. I need to use this for sagemanifolds, and just need for the differentiation that occurs, ultimately, to return the original functions - sqrt and g1.

integration is trying several engines, one of them is wrong:

Thanks for reporting this bug. This is now the issue #36222.