integrate sqrt(c^2+k^2*sin^2(c*tau))/v d(tau)