Here is a Mathematica notebook. I did some corrections and got some equations. But I think there is still something off. There are two problems here. In the last set of equations for "sfecomp", \phi should be function of "r" (\phi[r]). So, I think, I have to define it separately somehow. Another question is These equation shows the covariant derivatives od \phi. They should be evaluated completely and we should get ordinary differential equation. Can anyone help please?
Thank You.