integrate((tau+t0)^(-1/2)/sqrt(t-tau),tau=0..t)