I have made some corrections to acceleration formula

http://www.part.lt/img/562af2c949145...8bca6b6221.png

and have made numerical integration.

Now it looks it fits very well with experimental data

( about 0.103 arc seconds per revolution ).

But I still have no idea how to solve it analytically.

Maybe some advice? Thank you.