Glik[[2,All,All]] /. t->Pi/2 1 1 Out[19]= {{0, -, 0, 0}, {-, 0, 0, 0}, {0, 0, 0, 0}, {0, 0, 0, 0}} r r Thus in the geodesic equation for theta U={r'[s],theta'[s],p'[s],ct'[s]} Simplify[Evaluate[Sum[Glik[[2,i,j]] U[[i]] U[[j]],{i,4},{j,4} ] /. t->Pi/2] ] 2 r'[s] theta'[s] Out[23]= ----------------- r So if t=Pi/2 and theta'[0]=0, then theta''[s]=0, so theta remains Pi/2