Gielis公式 python代码