Math 如何计算任意给定坐标下的月亮相角？

这几天我遇到了一个问题，我试图精确计算月球在任何给定坐标（lon，lat）下的相位角。如果我知道我当前的坐标，月亮和太阳的地心赤经、赤纬、高度和天顶坐标

鉴于所有这些数据，我如何实现这一点？我想要的结果是以度为单位的，其中0度表示新月直接指向我，180度表示新月指向另一个方向，以此类推

---

我试过一个公式，在这个公式中，我测量了太阳天顶和月亮天顶的某个坐标之间的夹角，在纸上这是可行的，结果是不错的，但没有达到我想要的精度。任何帮助都将不胜感激。

这不是一个真正的软件问题-也许可以在这里发布您的问题？谢谢。我不知道这个问题的存在。我投票结束这个问题，因为，正如已经建议和同意的，这不是一个真正的编程问题。