C语言实现大地坐标系到高斯坐标系的转换

WGS84坐标系的经纬度和各向速度转换为高斯坐标系的三个位置和三个速度，有没有人知道公式，越简明越好，最好能提供C语言版本。

有现成公式可以套用，很方便的，需要时与我联系。