sqrt integral not working

edit

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.