There you go, being all confrontational and aggressive again without looking into where and why we differ and seeing if in fact our results are in agreement and just expressed in different ways.It does not help that it has taken over 80 posts to discover what you think the results should be.
OK, straight up, I am not familiar with the formalism you are using, as I do not come across it very often. If you care to state the exact titles and ISBN's of the Gron and Rindler books (or the Amazon links), I will order them and at least then we will be using the same references and notation and perhaps we might understand each other better.
Now for your derivation.
The strong field aproximation metric is:
ds^2= \left(e^{\frac{2\Phi}{c^2}}\right)c^2dt^2 - \left(e^{\frac{2\Phi}{c^2 }}\right)^{-1}dr^2 - r^2(d \theta^2 - \sin^2 \theta d \phi^2)
and the Regular Schwarzschild metric is:
ds^2= \left(1-\frac{2GM}{rc^2}}\right)c^2dt^2 - \left(1-\frac{2GM}{rc^2 }}\right)^{-1}dr^2 - r^2(d \theta^2 - \sin^2 \theta d \phi^2)
From the above two forms of the metric it is easy to make the identity:
e^{2\Phi/c^2}=\left(1-\frac{2GM}{rc^2}\right)
Solve the above for \phi:
\phi = \frac{1}{2}c^2 \log \left(1-\frac{2GM}{rc^2}\right)
Now you state:
\vec{F}=-grad\Phi
(as does Wikipedia here
http://en.wikipedia.org/wiki/Force#Potential_energy)
so:
\vec{F}= -grad\Phi = -grad\left(\frac{1}{2}c^2 \log \left(1-\frac{2GM}{rc^2}\right)\right)
\vec{F} = \frac{GM}{r^2}\left(1-\frac{2GM}{rc^2}\right)^{-1}