long double d;
if(UT>=UTo) d=UT-UTo;
else d=UTo-UT;
if(d>=0.1) {
UTo=UT;
UT=UT_rise(UTo,GHA(UTo,G_sun(t_century(days(year,month,date),UTo)),ecliptic_longitude(L_sun(t_century(days(year,month,date),UTo)),G_sun(t_century(days(year,month,date),UTo)))),glong,e(h,glat,sun_deviation(earth_tilt(t_century(days(year,month,date),UTo)),ecliptic_longitude(L_sun(t_century(days(year,month,date),UTo)),G_sun(t_century(days(year,month,date),UTo))))));
result_rise(UT,UTo,glong,glat,year,month,date);
}
return UT; }
//判断并返回结果(日出)
long double result_set(long double UT, long double UTo, long double glong, long double glat, int year, int month, int date){
long double d;
if(UT>=UTo) d=UT-UTo;
else d=UTo-UT;
if(d>=0.1){
UTo=UT;
UT=UT_set(UTo,GHA(UTo,G_