Suppose A be a point on a geodesically complete Riemannian manifold (M,g),and define a distance function on M, r(x)=d(x,A).One can see intuitively that the gradient vector of the distance function at x is just the targent vector at x of the geodesic connecting x and A, with length 1.
My question is where can I find a formal proof ?

