func GetSunEclipticLatitudeForEarth(jd float64) float64 {
t := calendarutil.GetJulianThousandYears(jd)
B0 := GetEarthB0(t)
B1 := GetEarthB1(t)
B2 := GetEarthB2(t)
B3 := GetEarthB3(t)
B4 := GetEarthB4(t)
B := ((((B4*t)+B3)*t+B2)*t+B1)*t + B0
return B
}
/**
* 按照儒略日计算地球和太阳的距离
*
* @param jd
* 儒略日
* @return 地球和太阳的距离，单位是天文单位(au)
*/
func GetSunRadiusForEarth(jd float64) float64 {
t := calendarutil.GetJulianThousandYears(jd)
R0 := GetEarthR0(t)
R1 := GetEarthR1(t)
R2 := GetEarthR2(t)
R3 := GetEarthR3(t)
R4 := GetEarthR4(t)
R5 := GetEarthR5(t)
R := ((((R5*t+R4)*t+R3)*t+R2)*t+R1)*t + R0
return R
}
func GetSunEclipticLongitudeForEarth(jd float64) float64 {
t := calendarutil.GetJulianThousandYears(jd)
L0 := GetEarthL0(t)
L1 := GetEarthL1(t)
L2 := GetEarthL2(t)
L3 := GetEarthL3(t)
L4 := GetEarthL4(t)
L5 := GetEarthL5(t)
L := ((((L5*t+L4)*t+L3)*t+L2)*t+L1)*t + L0
return mathutil.Mod2Pi(L)
}