1568-1644 Solar eclipses, make 1617(8) and 1620(12) Lunar eclipses the same as Chinese records , C++ (gcc)

1568-1644 Solar eclipses, make 1617(8) and 1620(12) Lunar eclipses the...

Language: Layout:

// C++ Solar eclipse program with some Lunar eclipse(1568-1644)=>adjusted MinSun(printed) with Zhuanzhong_C = 27.554368
// Author: Mr Lam/ Yee Din(林貽典) email : yeedlam@yahoo.com
// Date : Nov. 4, 2014. Copyright reserved, please contact the author if you want to use it.
 
#include <iostream>
using namespace std;
#include <cmath>
 
const int start_year = 1568; // enter start_year
const int end_year = 1644; // enter end_year
 
const double Tropical_year = 365.2425;
const int Base_year = 1281;
//===================================================================
// 4 Paramters for adjustments
const double Qiying_P = 55.06, Runying_P = 20.2050;
const double Jiaoying_P = 26.0388-0.10243;   
const double Zhuanying_P = 13.0205-0.785192;
//-------------------------------
// Constants for Moon to determine solar eclipse purpose
const double Shuoshi_C = 29.530593, Zhuanzhong_C = 27.554368;      // 27.5546;
const double Jiao_eclipse_C = 2.318369, Jiaozhong_C = 27.212224;
const double half_Zhuanzhong = Zhuanzhong_C/2;
//-------------------------------
// Constant for Moon positions : Lead or Lag
const double MoonRef_C = 1.976225;                                 //  1.975993;
const double MoonMeanSpeed = 13.36878311;                          // 13.36875;
const double MinSun_drift=0.023452;        //  MinSun=MinSun_drift*(calculated_year-Base_year);
//  RefSun_ = half_year-LeapDay;
//  RefSun=RefSun_- MinSun;               //  MinSun added back to LunStep 18 as: LEclNetSun_=LEclNetSun_+ MinSun
//===================================================================
 
typedef struct{
int year; int month; int day;
double EclipseFrac; double EclipseBegin; double MaxPhase; double EclipseEnd;
double Winter_solitice;
int SunStatus; int MoonStatus;
double observed_time; double HalfDayFan;
double ReqSun; double ReqMoon;
double half_dayoff;
 
double MoonSp;
double noon_offset;
 
double ConjuncT24;
double SunX; double MoonX; double SunMoonXS;
} EclipseTime;
 
EclipseTime Shoushi[100];
 
 EclipseTime Theo[]
= {
   {1568,  9, 21,     0.109955,  11.1367,   11.7366,   12.3428},  // from eclipse.org.uk
   {1572,  6, 10,     0.916,      6.66,      7.93,      9.93},
   {1575,  5, 10,     0.701,     13.05,     14.29,     15.45},
   {1582,  6, 20,     0.723024,  12.2055,   13.6179,   14.9154},  // from eclipse.org.uk
   {1583, 12, 14,     0.825,     11.23,     13.08,     14.77},
   {1587, 10,  2,     0.735831,  11.3223,   12.6436,   13.9343},  // from eclipse.org.uk
   {1589,  2, 15,     0.301344,   6.4499,    7.1961,    7.99914}, // from eclipse.org.uk
   {1590,  7, 31,     0.430839,  16.0803,   17.1059,   18.0372},  // from eclipse.org.uk
   {1594,  5, 20,     0.319,      9.83,     10.84,     11.9},
   {1596,  9, 22,     0.949,      9.87,     11.17,     12.5},
 
   {1603,  5, 11,     0.903,      7.4,       8.59,      9.93},
   {1604,  4, 29,     0.078077,  15.8827,   16.4079,   16.907},   // from eclipse.org.uk
   {1607,  2, 26,     999,       17.32,     18.26,     999},      // no data from eclipse.org.uk
   {1609, 12, 26,     0.044430,  16.197,    16.594,    16.9753},  // from eclipse.org.uk but no apparent noon
   {1610, 12, 15,     0.607,     14.61,     16.07,     17.32},
   {1612, 12, 99,     999,         999,      999,      999},      // no data from eclipse.org.uk
   {1615,  3, 29,     0.889903,  15.1624,   16.4118,   17.5463},  // from eclipse.org.uk
   {1617,  8,  1,     0.470019,  19.07,     19.9041,   20.6909},  // eclipse_fraction, eclipse_maximum, eclipse_end from eclipse.org.uk
   {1621,  5, 21,     0.711,     17.27,     18.23,     19.19},
   {1629,  6, 21,     0.166,     10.98,     11.71,     12.45},
 
   {1631, 10, 25,     0.108,     12.34,     12.91,     13.48},
   {1634,  3, 29,     0.81964,    8.77022,   9.94358,  11.1899 },  // from eclipse.org.uk
   {1636,  8,  1,     0.632763,   8.37958,   9.62653,  10.9643},   // from eclipse.org.uk
   {1637,  1, 26,     0.199,     12.71,     13.68,     14.58},
   {1638,  1, 15,     0.656,     11.65,     13.32,     14.84},
   {1641, 11,  3,     0.827,     12.65,     14.04,     15.33},
   {1643,  3, 20,     0.53,       7.83,      8.85,      9.94},
   {1644,  9,  1,     0.303,     10.93,     12.14,     13.34},
 
   {99999, 0,  0,     0,        0,        0,        0}
  };
 
 EclipseTime Observe[]
= {
   {1568,  9, 21,     0.109955,  11.1367,   11.7366,   12.3428},  // from eclipse.org.uk
   {1572,  6, 10,     0.916,      6.66,      7.93,      9.93},
   {1575,  5, 10,     0.701,     13.05,     14.29,     15.45},
   {1582,  6, 20,     0.723024,  12.2055,   13.6179,   14.9154},  // from eclipse.org.uk
   {1583, 12, 14,     0.825,     11.23,     13.08,     14.77},
   {1587, 10,  2,     0.735831,  11.3223,   12.6436,   13.9343},  // from eclipse.org.uk
   {1589,  2, 15,     0.301344,   6.4499,    7.1961,    7.99914}, // from eclipse.org.uk
   {1590,  7, 31,     0.430839,  16.0803,   17.1059,   18.0372},  // from eclipse.org.uk
   {1594,  5, 20,     0.319,      9.83,     10.84,     11.9},
   {1596,  9, 22,     0.949,      9.87,     11.17,     12.5},
 
   {1603,  5, 11,     0.903,      7.4,       8.59,      9.93},
   {1604,  4, 29,     0.078077,  15.8827,   16.4079,   16.907},   // from eclipse.org.uk
   {1607,  2, 26,     999,       17.32,     18.26,     999},      // no data from eclipse.org.uk
   {1609, 12, 26,     0.044430,  16.197,    16.594,    16.9753},  // from eclipse.org.uk but no apparent noon
   {1610, 12, 15,     0.607,     14.61,     16.07,     17.32},
   {1612, 12, 99,     999,         999,      999,      999},      // no data from eclipse.org.uk
   {1615,  3, 29,     0.889903,  15.1624,   16.4118,   17.5463},  // from eclipse.org.uk
   {1617,  8,  1,     0.470019,  19.07,     19.9041,   20.6909},  // eclipse_fraction, eclipse_maximum, eclipse_end from eclipse.org.uk
   {1621,  5, 21,     0.8,       17.24,     9999,      9999},
   {1629,  6, 21,     0.166,     10.98,     11.71,     12.45},
 
   {1631, 10, 25,     0.108,     12.34,     12.91,     13.48},
   {1634,  3, 29,     0.81964,    8.77022,   9.94358,  11.1899 },  // from eclipse.org.uk
   {1636,  8,  1,     0.632763,   8.37958,   9.62653,  10.9643},   // from eclipse.org.uk
   {1637,  1, 26,     0.199,     12.71,     13.68,     14.58},
   {1638,  1, 15,     0.656,     11.65,     13.32,     14.84},
   {1641, 11,  3,     0.827,     12.65,     14.04,     15.33},
   {1643,  3, 20,     0.53,       7.83,      8.85,      9.94},
   {1644,  9,  1,     0.303,     10.93,     12.14,     13.34},
 
   {99999, 0,  0,     0,        0,        0,        0}
  };
 
enum Status{LagStart, LagEnd, LeadStart, LeadEnd};
struct Status_type{char *sun, *moon;};
 
class Ref_data
{
  private:
    Status SunStatus, MoonStatus;
 
  protected:
    int calculated_year, MonthNum;
    double SumDay, Winter_solitice, LeapDay;
    double RefEclipse, RefSun, RefMoonDay, RefMoon;
    double ReqEclipse, ReqSun, ReqMoonDay, ReqMoon;
    double UseSun, SunAdd, SunSum, SunX;
    double MoonSum, SoD;
    double moon_x12_2;
    int FrontEclipseSign;
 
  public:
    Ref_data(int m,int n){calculated_year = m;  MonthNum = n;}
    ~Ref_data(){ }
 
void check_SunStatus(double half_year,double ReqSun, double& UseSun, int& SunStatus)
    {
      SunStatus=LagStart;
      UseSun = 0;
 
      if (ReqSun < half_year)
      {
        SunStatus = LagStart;
        UseSun = ReqSun;
        if (ReqSun > 93.712025){UseSun = half_year - UseSun; SunStatus = LagEnd;}
      }
//  推初末限:　視盈曆在盈初縮末限已下，縮曆在縮初盈末限已下，各為初。已上用減半歲周為末。
     if (ReqSun > half_year)
     {
        SunStatus = LeadStart;
        UseSun = ReqSun - half_year ;
        if (UseSun > 88.909225){UseSun = half_year - UseSun; SunStatus = LeadEnd;}
     }
   }//end check_SunStatus
 
void show_SunStatus(int SunStatus)
    {
    cout << "SunStatus = ";
    switch (SunStatus)
    {
      case LagStart:cout << "SunLagStart\n"; break;
      case LagEnd:cout << "SunLagEnd\n"; break;
      case LeadStart:cout << "SunLeadStart"; break;
      case LeadEnd:cout << "SunLeadEnd"; break;
      default: cout << "no SunStatus"; break;
    }
    cout << endl;
    }//end show_SunStatus
 
void get_SunX(int SunStatus,double UseSun,int IUseSun,double& SunAdd,double& SunSum,double& SunX)
    {
    SunAdd=SunSum=SunX=0;
    if (fmod(SunStatus, 3.0)==0) // SunStaus=LagStart(0) & LeadEnd(3)
    {
    SunAdd = 484.8473-(4.4362 + 0.0162*(IUseSun - 1)/2)*IUseSun ;
    SunSum = (487.06-(2.21 + 0.0027*IUseSun )*IUseSun )*IUseSun ;
    }
    if (fmod(SunStatus, 3.0) > 0) // SunStatus=LagEnd(1) & LeadStart(2)
    {
    SunAdd = 510.8569-(4.9386 + 0.0186*(IUseSun - 1)/2)*IUseSun ;
    SunSum =(513.32-(2.46 + 0.0031*IUseSun)*IUseSun)*IUseSun;
    }
    SunX =((UseSun - IUseSun)*SunAdd + SunSum)/10000.0;
    }// end get_SunX
//--------------------------------------------------------------
 void get_LSunX(int LSunStatus,double LUseSun, double& LSunX)
 {
    double LSunSum=0;
 
    if (fmod(LSunStatus, 3.0) == 0){// for LSunStaus = LagStart(0) & LeadEnd(3)
                                    LSunSum = (487.06-(2.21 + 0.0027*LUseSun)*LUseSun)*LUseSun ;
                                   }
    if (fmod(LSunStatus, 3.0) > 0){// for LSunStatus = LagEnd(1) & LeadStart(2)
                                    LSunSum =(513.32-(2.46 + 0.0031*LUseSun)*LUseSun)*LUseSun;
                                  }
    LSunX=LSunSum/10000.0;
 }// end get_LSunX
//--------------------------------------------------------------
void get_UseMoon(double ReqMoon,double& UseMoon,double& SoD,double& moon_x12_2,double& MoonSolar,
                 int& MoonStatus, double& MoonCh, double& MoonAdd)
    {
    double moon_84_check;
    if (ReqMoon > half_Zhuanzhong){UseMoon = ReqMoon - half_Zhuanzhong; MoonStatus = LagStart;}
    if (ReqMoon < half_Zhuanzhong){UseMoon = ReqMoon; MoonStatus = LeadStart;}
 
//  推遲疾限:　各置遲疾曆，以日轉限乘之，即得限數。(日轉限一十二限二十。)
    moon_x12_2=UseMoon*12.2;
    SoD=int(moon_x12_2)*0.082008;
 
    moon_84_check=moon_x12_2;
    if (moon_x12_2 > 84){
                         if (MoonStatus==LeadStart){MoonStatus=LeadEnd;}
                         if (MoonStatus==LagStart){MoonStatus=LagEnd;}
                         moon_84_check = 167 - moon_84_check;
                        }
    MoonCh = moon_84_check;
    MoonSolar = UseMoon - SoD;
    if (MoonSolar < 0){
//  求遲疾差:　置遲疾曆，以立成日率減之，(如不及減，則退一位。)餘以其下損益分乘之為實，八百二十分為法除之，得數以加其下遲疾積，即遲疾差。
                        SoD=SoD-0.082008;
                        MoonSolar=UseMoon-SoD;
                        if (fmod(MoonStatus, 2.0)==0){MoonCh=int(MoonCh);}
                        if (fmod(MoonStatus, 2.0) > 0){MoonCh=167-int(moon_x12_2);}
                       }
    MoonAdd = 11.081575-(0.05815 + 0.00195*(MoonCh - 1)/2)*MoonCh;
    }// end get_UseMoon
 
void get_MoonX(int MoonStatus, double MoonCh, double MoonAdd,double MoonSolar,
               double& MoonSum, double& MoonX, double& MoonSp)
    {
    MoonSum=0, MoonX=0, MoonSp=0;
    switch(MoonStatus)
     {
      case LagStart:
           MoonSp=0.9855+(0.05815+0.00195*(MoonCh-1)/2)*MoonCh/100;
           MoonSum=(0.1111-(0.000281+0.00000325*MoonCh)*MoonCh)*MoonCh;
           MoonX=MoonSum+MoonSolar*MoonAdd/820;
           break;
 
      case LagEnd:
           MoonSp=1.2071-(0.05815+0.00195*(MoonCh-1)/2)*MoonCh/100;
           MoonCh=MoonCh+1;
           MoonSum=(0.1111-(0.000281+0.00000325*MoonCh)*MoonCh)*MoonCh;
           MoonX=MoonSum-MoonSolar*MoonAdd/820;
           break;
 
      case LeadStart:
           MoonSp=1.2071-(0.05815+0.00195*(MoonCh-1)/2)*MoonCh/100;
           MoonSum=(0.1111-(0.000281+0.00000325*MoonCh)*MoonCh)*MoonCh;
           MoonX=MoonSum+MoonSolar*MoonAdd/820;
           break;
 
      case LeadEnd:
           MoonSp=0.9855+(0.05815+0.00195*(MoonCh-1)/2)*MoonCh/100;
           MoonCh=MoonCh+1;
           MoonSum=(0.1111-(0.000281+0.00000325*MoonCh)*MoonCh)*MoonCh;
           MoonX=MoonSum-MoonSolar*MoonAdd/820;
           break;
 
      default:
           break;
     }// end switch(MoonStatus)
    }//  end get_MoonX
 //----------------------------------------------------
 
void get_LMoonX(int LMoonStatus,double LMoonCh_more,double& LMoonX,double& LMoonSp,double& LMoonCh_)
  {
    LMoonX=0,LMoonSp=0;
    if (fmod(LMoonStatus,2) > 0){LMoonCh_more = LMoonCh_more + 1;}
 
    LMoonSp=0.9855 +(0.05815 + 0.00195*(LMoonCh_more-1)/2)*LMoonCh_more/100;
    if(fmod(LMoonStatus, 3) > 0){LMoonSp=1.2071-(0.05815 + 0.00195*(LMoonCh_more-1)/2)*LMoonCh_more/100;}
 
    LMoonX =(0.1111-(0.000281 + 0.00000325*LMoonCh_more)*LMoonCh_more)*LMoonCh_more;
    LMoonCh_=LMoonCh_more;
  } // end get_LMoonX
//-----------------------------------------------------------
 
void get_SunMoonT(int SunStatus,int MoonStatus,double MoonX,double SunX,double MoonSp,
                  double& SunXS, double& MoonXS, double& SunMoonT)
//  推加減差:　視經朔弦望下所得盈縮差、遲疾差，以盈遇遲、縮遇疾為同相併，盈遇疾、縮遇遲為異相較，各以八百二十分乘之為實，
    {
    SunXS=SunX;   if (SunStatus < LeadStart){SunXS = -SunX;}
    MoonXS=MoonX; if( MoonStatus > LagEnd){MoonXS = -MoonX;}
 
    SunMoonT =(MoonXS+SunXS)*0.082/MoonSp;
    } // end get_SunMoonT
//--------------------------------------------------------
 
void print_SolarEclipseTable(int p)
    {
        int index=p;
    printf("Reference:\n");
    printf("(1)Accuracy of Solar Eclipse Calculations, 1644-1785, with the Chinese Shoushi Calandar(Table 2)\n");
    printf("http://articles.adsabs.harvard.edu//full/1999JHA....30..161L/0000161.000.html\n\n");
 
    printf("(2)Accuracy of Solar Eclipse Observations made by Jesuit Astronomers in China\n");
    printf("http://adsabs.harvard.edu/full/1995JHA....26..227S\n");
    printf("Observation data from Table 2\n");
    printf("Observation error data from Table 3\n");
    printf("Theorectical=Observation-Observation error\n\n");
 
    printf("\t Shoushi(Observation)\t\tRemark : 999 or 9999 = no data\n");
    printf("\t-)Theorectical(Obsevation-Shoushi)\n");
    printf("\t--------------\n");
    printf("Shoushi-Theorectical(Observation-Theorectical)\n");
    printf("==================================================================================================================\n");
    printf("N  year(month)\t  EclipseFraction\tEclipseBegin\t\tMaxPhase\t\tEclipseEnd\n");
    printf("      month/day\n");
    printf("==================================================================================================================\n");
    int record;
    for (record=0; record < index; record ++)
    {
      if(Shoushi[record].year==0){Shoushi[record].year=Observe[record].year;}
      printf("%2d %4d(%1d)",record,Shoushi[record].year, Shoushi[record].month);
      printf(" \t   %5.4f", Shoushi[record].EclipseFrac);
      printf("(%5.3f)\t",Observe[record].EclipseFrac);
      printf("%5.4f", Shoushi[record].EclipseBegin);
      printf("(%5.3f)   \t", Observe[record].EclipseBegin);
      printf("%5.4f", Shoushi[record].MaxPhase);
      printf("(%5.3f)   \t", Observe[record].MaxPhase);
      printf("%5.4f", Shoushi[record].EclipseEnd);
      printf("(%5.3f)", Observe[record].EclipseEnd);
      printf("\n");
      printf("    %4d/%1d",Theo[record].month, Theo[record].day);
 
      printf("\t-) ");
      printf("%5.3f", Theo[record].EclipseFrac);
      printf("(%5.4f)   \t", Observe[record].EclipseFrac-Shoushi[record].EclipseFrac);
 
      printf("%5.3f", Theo[record].EclipseBegin);
      printf("(%5.4f)   \t",Observe[record].EclipseBegin-Shoushi[record].EclipseBegin);
 
      printf("%5.3f", Theo[record].MaxPhase);
      printf("(%5.4f)   \t",Observe[record].MaxPhase-Shoushi[record].MaxPhase);
      printf("%5.3f", Theo[record].EclipseEnd);
      printf("(%5.4f)",Observe[record].EclipseEnd-Shoushi[record].EclipseEnd);
      printf("\n");
      printf("\t\t----------------------------------------------------------------------------------------\n");
      printf("\t\t   ");
      printf("%5.4f", Shoushi[record].EclipseFrac-Theo[record].EclipseFrac);
      printf("(%5.3f)   \t", Observe[record].EclipseFrac-Theo[record].EclipseFrac);
 
      printf("%5.4f",Shoushi[record].EclipseBegin-Theo[record].EclipseBegin);
      printf("(%5.3f)   \t",Observe[record].EclipseBegin-Theo[record].EclipseBegin);
 
      printf("%5.4f",Shoushi[record].MaxPhase-Theo[record].MaxPhase);
      printf("(%5.3f)   \t",Observe[record].MaxPhase-Theo[record].MaxPhase);
 
      printf("%5.4f",Shoushi[record].EclipseEnd-Theo[record].EclipseEnd);
      printf("(%5.3f)",Observe[record].EclipseEnd-Theo[record].EclipseEnd);
      printf("\n\n");
    }; // end for loop
   }  // end print_SolarEclipseTable
  };// end class Ref_data
 
void print_CalculatedRecords_1617()
{
 printf("\t\t《大明神宗顯皇帝實錄卷之五百四十七》\n");
 printf("\t\t\t\t邢雲路獻七政真數..(望交泛分:13.845663)\n");
 printf("\t\t以此法布算今歲七月十六日戊寅夜望月食,推得是月:\n");
 printf("\t\t\t望交泛分一十三日八十四刻五十六分六十三秒(13.845663)\n");
 printf("\t\t\t陰歷交前度一度二十二分零三秒九十五微(1.220395)\n");
 printf("\t\t\t月食分一十三分五十九秒六十九微(1.35969)\n");
 printf("\t\t\t定用分六刻六十八分九十五秒二十一微(0.06689521)\n");
 printf("\t\t\t既內分二刻八十六分六十五秒\n");
 printf("\t\t\t既外分三刻八十一分七十九秒\n");
 printf("\t\t\t初虧分六刻四十三分二十秒八十九微(0.064[3]2089,應修正為:復圓分-2*定用分=0.19801131-0.06689521*2=0.064[2]2089)\n");
 printf("\t\t\t食既分一十刻二十五分零四十微\n");
 printf("\t\t\t食甚分一十三刻一十一分六十六秒一十微(0.1311[6]610,應修正為:復圓分-定用分=0.19801131-0.06689521=0.1311[1]610)\n");
 printf("\t\t\t生光分一十五刻九十八分三十一秒八十微\n");
 printf("\t\t\t復圓分一十九刻八十分一十一秒三十一微(0.19801131)\n\n");
 
 printf("\t\t\t\t   初虧丑初二刻(1.48h<=>0.064[2]2089*24=1.54130136,仍然算是丑初二刻,所以該記錄是可靠的)\n");
 printf("\t\t\t食既丑正一刻,食甚寅初二刻(3.48h? <=>0.1311[1]610*24=3.1467864h=>應修正為寅初一刻)\n");
 printf("\t\t\t生光寅初三刻,復明寅正三刻(4.72h<=>0.19801131*24=4.75227144h=寅正三刻,所以該記錄亦是可靠的)\n");
}// end  print_CalculatedRecords_1617
 
void print_CalculatedRecords_1620()
{
 printf("\t\t\t《大明熹宗哲皇帝實錄卷之七》\n");
 printf("\t\t\t\t按新法推泰昌元年庚申歲十一月十六日己丑夜望月食:\n");
 printf("\t\t\t\t\t初虧漏下二百七十三籌五十三分,計九十七刻二十四分(0.9724)\n");
 printf("\t\t\t\t\t食既漏下二百九十四籌六十分,計八十八分\n");
 printf("\t\t\t\t\t生光漏下三百三十二籌四十分,計七刻四十分\n");
 printf("\t\t\t\t\t復圓漏下三百五十三籌四十八分,計一十一刻(0.11)\n");
 printf("\t\t\t\t各以發斂求之,得:\n");
 printf("\t\t\t\t\t初虧夜子初一刻(23.24h)\n");
 printf("\t\t\t\t\t食既子正初刻\n");
 printf("\t\t\t\t\t食甚子正四刻(0.96h)\n");
 printf("\t\t\t\t\t生光丑初三刻\n");
 printf("\t\t\t\t\t復圓丑正二刻(2.48h,但見有改自三刻的痕跡)\n");
 printf("\t\t\t\t\t月食一十四分九十九秒(1.499)\n");
 printf("\t\t\t\t\t食甚月離黃道畢宿一十四度三十分一十秒(14.3010)\n");
 printf("\t\t\t\t而授時則推:\n");
 printf("\t\t\t\t\t初虧子正一刻(0.24h)\n");
 printf("\t\t\t\t\t食既丑初一刻\n");
 printf("\t\t\t\t\t食甚丑正初刻(2.00h?)\n");
 printf("\t\t\t\t\t生光丑正三刻\n");
 printf("\t\t\t\t\t復圓寅初三刻(3.72h?)\n");
 printf("\t\t\t\t\t月食一十三分三十一秒(1.331)\n");
 printf("\t\t\t\t\t月離黃道畢宿一十四度三十四分四十四秒(14.3444),與天不合\n");
}   // end print_CalculatedRecords_1620
 
void print_CalculatedRecords_1623()
{
 printf("\t\t\t《大明熹宗哲皇帝實錄卷之七》\n");
 printf("\t\t\t\t以新法推天啟三年癸亥歲九月十五壬寅夜望月食:\n");
 printf("\t\t\t\t\t初虧申正三刻\n");
 printf("\t\t\t\t\t食甚酉初三刻(17.72?,該記錄不可靠)\n");
 printf("\t\t\t\t\t復圓戍初初刻\n");
 printf("\t\t\t\t\t月食五分三十八秒\n");
 printf("\t\t\t\t\t月未出巳食五分二十五秒,月已出見食一十三秒,不及一分不救\n");
 printf("\t\t\t\t而授時則推月食:\n");
 printf("\t\t\t\t\t初虧酉初一刻\n");
 printf("\t\t\t\t\t食甚酉正二刻(18.48h,該記錄可靠)\n");
 printf("\t\t\t\t\t復圓戌正一刻\n");
 printf("\t\t\t\t\t月食七分八秒\n");
 printf("\t\t\t\t\t月未出已食一分九十八秒,月已出見食五分一十秒,與天不合\n");
}   // end print_CalculatedRecords_1623
 
void print_FindingRecords_1621()
{
 printf("\t\t《大明熹宗哲皇帝實錄卷之七》\n");
 printf("\t\t原任陝西按察使邢雲路獻七政真數.....復以此法推天啟元年辛酉歲四月壬申朔日食:\n");
 printf("\t\t\t\t初虧申正一刻(16.24h)\n");
 printf("\t\t\t\t食甚酉初刻(17.00h)\n");
 printf("\t\t\t\t復圓酉正初刻(18.00h)\n");
 printf("\t\t\t\t日食一分八十六秒(0.186)\n");
 printf("\t\t\t\t復圓日在天未入地,食不及三分不救(以上記錄缺乏其它數據可供證驗,是否可靠,甚有疑問?)\n");
 printf("\t\t而授時則推是日食:\n");
 printf("\t\t\t\t初虧申正三刻(16.72h)\n");
 printf("\t\t\t\t食甚酉正初刻(18.00h\n");
 printf("\t\t\t\t復圓戍初刻(19.00h)\n");
 printf("\t\t\t\t日食三分九十一秒(0.391)\n");
 printf("\t\t\t\t日未入已復光三分一十一秒,日已入未復光八十秒,與天不合(以上授時數據可按計算方法證驗)\n\n");
 printf("\t\t《明熹宗悊皇帝實錄卷之十六》\n");
 printf("\t\t庚申原任陝西按察使閒住邢雲路奏:\n");
 printf("\t\t\t\t臣本年春言四月朔日食歷數,臨期策驗分刻乃有不合\n");
 printf("\t\t\t\t臣始悟歷元玄法原自無差,而其差者,乃今大統歷沿郭守敬授時歷之誤也\n");
 printf("\t\t\t\t蓋日食者,月食之也,故正歷一準於月食,而月之分數與限度定法,皆其綱領最大者\n");
 printf("\t\t\t\t授時定月食分三十分及十五分,月食限十三度五分,定限八十七分\n");
 printf("\t\t\t\t此守敬率南北日官以儀器測驗于天,自丙子至庚辰又經甲午,再驗幾二十年,自謂推步之精,可永久無弊\n");
 printf("\t\t\t\t以故大統歷沿用之,臣今所步歷亦用之,及四月朔食測驗,方見其差,乃悟以異乘同\n");
 printf("\t\t\t\t除另變其法,變月食分三十分為三十五分(30->35)\n");
 printf("\t\t\t\t         變十五分為十七分五十秒\n");
 printf("\t\t\t\t         變定限十三度五分為十三度十二分五十秒(13.05->13.1250)\n");
 printf("\t\t\t\t         變定法八十七分為七十五分(8.7->7.5)\n");
 printf("\t\t\t\t諸數既變,乃依舊歷元玄法布歷,一一方合,然後知歷元法原自無差,而其差者,大統相沿授時法誤之也,故必如臣所推:\n");
 printf("\t\t\t\t如五千七百四十為日食成限,以四千九百二十為月食成限,以三十五分為月食分,而黜郭守敬之三十分,方為正法,下禮部\n\n");
 
 printf("\t\t\t\t  LunFraction=(13.05 - abs(SexLunEclA_))/8.7\n");
 printf("\t\t\t月食食分=定限(13.05->13.1250) - 視入陰陽曆/定法(8.7->7.5)\n\n");
 
 printf("\t\t\t\tLHalfEclT=(sqrt((30-LunFraction*10)*LunFraction*10)*49.20)/(LEclM_SunSp*10000)\n");
 printf("\t\t\t\t定用分=(sqrt(三十分(30->35)-月食分)*月食分)*49.20)/(定限行度*10000)\n");                   
}   // end print_FindingRecords_1621
 
void print_LunStep_13_Extra(int k)
{
    int calculated_year=k;
    if (calculated_year==1617){
                               printf("\t\t\t\t\t\t\t\t\t\t\t邢雲路食甚分 0.13111610*24=3.1467864h\n");
                               printf("\t\t\t\t\t\t\t\t\t\t\t\t\t   理論值:3.219\n");
                              }
    if (calculated_year==1620){
          printf("   邢雲路..按新法推..初虧..計九十七刻二十四分..復圓..計一十一刻=>食甚分=0.11-(0.11+(1-0.9724))/2=0.0412");
          printf("=>食甚=0.0412*24=0.9888h=子正四刻,與記錄一致\n");
          printf("\t\t\t\t\t\t\t\t\t\t\t\t\t    理論值:1.17h \n");
                             }
    if (calculated_year==1623){
                               printf("\t\t\t\t\t\t\t\t\t\t\t  而授時則推...食甚酉正二刻(18.48h),該記錄可靠\n");
                               printf("\t\t\t\t\t\t\t\t\t\t\t\t\t   理論值:19.23h \n");
                              }
}   // print_LunStep_13_Extra
 
void print_LunStep_23_Extra(int k)
{
    int calculated_year=k;
    if (calculated_year == 1617){
                                 printf("\t\t\t\t\t邢雲路獻七政真數..以此法布算..月食分一十三分五十九秒六十九微(1.35969)\n");
                                 printf("\t\t\t\t\t\t\t\t\t\t\t 理論值:1.403\n");        
                            }
    if (calculated_year == 1620){
                                 printf("\t\t\t\t\t\t\t    授時則推..月食一十三分三十一秒(1.331)\n");
                                 printf("\t\t\t\t\t\t邢雲路.....按新法推....月食一十四分九十九秒(1.499)\n");
                                 printf("\t\t\t\t\t\t\t\t\t\t   理論值:1.6\n");
                                }
    if (calculated_year == 1623){
                                 printf("\t\t\t\t以新法推天啟三年癸亥歲九月十五壬寅夜望月食...月食五分三十八秒(0.538)\n");
                                 printf("\t\t\t\t\t\t\t\t\t\t\t理論值:0.6\n");
                                }
}   // end  print_LunStep_23_Extra
 
void print_LunStep_25_Extra(int k)
{
    int calculated_year=k;
    if (calculated_year == 1617)
    {
        printf("\t\t\t\t\t(邢雲路的初虧*24=1.54130136h=>算是丑初二刻\t   4.75227144h=寅正三刻,兩項記錄都可靠)\n");
        printf("\t\t\t\t\t\t理論值: 1.49h      \t\t   4.95h\n\n");
                  
double  SexLunEclA_1617=13.05-1.35969*8.7; // 計算出原有的 SexLunEclA_1617(不分陰陽曆)               
double  NewFraction=(13.125-1.220395)/7.5;              
                                           
double  LEclM_SunSp_1617=(sqrt((30-13.5969)*13.5969)*49.20)/(0.06689521*10000);  // 計算出原有的 LEclM_SunSp_1617
double  NewHalfT=(sqrt((35-NewFraction*10)*NewFraction*10)*49.20)/(LEclM_SunSp_1617*10000);
                  
double  NewBeginT=0.13111610-NewHalfT;
double  NewBeginT24=NewBeginT*24;                  
                  
double  NewEndT=0.13111610+NewHalfT;
double  NewEndT24=NewEndT*24;
        printf("\t\t\t\t邢雲路原有的陰陽曆 SexLunEclA_1617(不分陰陽)=13.05-1.35969*8.7=%1.5f,與記錄所指的 1.220395 基本一致\n\n",SexLunEclA_1617);
 
        printf("\t\t\t\t邢雲路變法後的食分 NewFraction=(13.125-1.220395)/7.5=%1.5f,大於理論值\n",NewFraction);
        printf("\t\t\t\t\t\t\t\t\t     理論值=1.403\n");
        printf("\t\t\t\t\t\t\t\t   邢雲路原有的食分: 1.35969\n\n");
 
        printf("\t\t  邢雲路原有的定限行度 LEclM_SunSp_1617=(sqrt((30-13.5969)*13.5969)*49.20)/(0.06689521*10000)=%1.5f\n",LEclM_SunSp_1617);                    
        printf("\t\t  邢雲路變法後改大的定用分 NewHalfT=(sqrt((35-NewFraction*10)*NewFraction*10)*49.20)/(LEclM_SunSp_1617*10000)=%1.5f",NewHalfT);
        printf("\n");
        printf("\t\t  邢雲路變法後的 NewBeginT=0.13111610-NewHalfT=%1.5f\t NewEndT=0.13111610+NewHalfT=%1.5f\n\n",NewBeginT,NewEndT);    
        printf("\t\t\t變法後的初虧=NewBeginT*24=%1.5f\t\t食甚不變=0.13111610*24=%1.5f\t\t復圓=NewEndT*24=%1.5f\n",NewBeginT24,0.13111610*24, NewEndT24);
        printf("\t\t\t\t     理論值: 初虧=1.49h\t\t\t\t  食甚=3.219h\t\t\t   復圓=4.95h\n");                   
        printf("\t\t\t\t 邢雲路原有的初虧=1.54130136h\t\t\t\t\t\t\t   復圓=4.75227144h\n");                                             
    }
    if (calculated_year == 1620)
    {
        printf("\t\t (邢雲路...按新法推..初虧=0.9724*24 = 23.3376h\t\t\t\t復圓=0.11*24=2.64h)\n");
        printf("\t\t\t\t\t理論值:初虧= 23.28h\t\t\t\t\t復圓=3.05h\n\n");
//-----------
double  SexLunEclA_1620=13.05-1.499*8.7; // 計算出原有的 SexLunEclA_1620(不分陰陽曆)                
double  NewFraction=(13.125-SexLunEclA_1620)/7.5;              
 
double  HalfT_1620=(0.11+(1-0.9724))/2;
double  LEclM_SunSp_1620=(sqrt((30-14.99)*14.99)*49.20)/(HalfT_1620*10000);  // 計算出原有的 LEclM_SunSp_1620
double  NewHalfT=(sqrt((35-NewFraction*10)*NewFraction*10)*49.20)/(LEclM_SunSp_1620*10000);
                  
double  NewBeginT=1+0.0412-NewHalfT;
double  NewBeginT24=NewBeginT*24;                  
                  
double  NewEndT=0.0412+NewHalfT;
double  NewEndT24=NewEndT*24;
        printf("\t\t\t\t邢雲路原有的陰陽曆 SexLunEclA_1620(不分陰陽)=13.05-1.499*8.7=%1.5f\n",SexLunEclA_1620);
        printf("\t\t\t\t邢雲路變法後的食分 NewFraction=(13.125-SexLunEclA_1620)/7.5=%1.5f\n\n",NewFraction);
 
        printf("\t\t\t\t邢雲路原有的定用分 HalfT_1620=(0.11+(1-0.9724))/2=%1.5f\n",HalfT_1620);
        printf("\t\t\t\t邢雲路原有的定限行度 LEclM_SunSp_1620=(sqrt((30-14.99)*14.99)*49.20)/(HalfT_1620*10000)=%1.5f\n",LEclM_SunSp_1620);                    
        printf("\t\t  邢雲路變法後的定用分 NewHalfT=(sqrt((35-NewFraction*10)*NewFraction*10)*49.20)/(LEclM_SunSp_1620*10000)=%1.5f\n",NewHalfT);
        printf("\t\t\t\t邢雲路變法後的 NewBeginT=1+0.0412-NewHalfT=%1.5f\tNewEndT=0.0412+NewHalfT=%1.5f\n\n",NewBeginT,NewEndT);
 
        printf("\t\t邢雲路變法後的初虧=NewBeginT*24=%1.5f\t\t食甚不變=0.0412*24=%1.5f\t     復圓=NewEndT*24=%1.5f\n",NewBeginT24,0.0412*24,NewEndT24);
        printf("\t\t\t\t  理論值:初虧= 23.28h\t\t\t     食甚=1.17h\t\t\t\t復圓=3.05h\n");
        printf("\t\t    邢雲路原有的初虧=0.9724*24=%1.5f\t\t\t\t\t\t\t復圓=0.11*24=%1.5fh\n",0.9724*24,0.11*24); 
    }
    if (calculated_year == 1623)
    {
        printf("    以新法推天啟三年癸亥歲九月十五壬寅夜望月食[初虧申正三刻](16.72)不可靠, 因為這應該是授時的計算值\n");
        printf("\t\t\t     而授時則推月食[初虧酉初一刻](17.24)不可靠, 因為與計算值不符合,很可能是新法推的初虧\n");
        printf("\t\t\t 以新法推天啟三年..[初虧申正三刻](16.72h ?)\t ....  復圓[戍初初刻](19h ?)\n");
        printf("\t\t\t\t\t\t理論值:17.7h\t\t\t\t     20.79h\n\n");
    }
}   //  end print_LunStep_25_Extra
 //---------------------------------------------------------------------
    int main()
    {
    cout << "\tSolar eclipse time is in 24 hours\n"
         << "\tBeijing local noon time is 12 as reference\n"
         << "\tmonth means total months from winter solitice\n"
         << "\tReqMDay=ReqMoonDay" << endl;
    cout << "\tt" << "\tt" << "\tt" << "\tt" << "\tt" << "\tt" << "\tt" << "\tt" << "\tt" << "\tt" << "\tt" << "\tt" << "\tt"
         << "\tt" << "\tt" << "\tt" << endl
         << " N year(month)\n";
 
    int index = 0; int calculated_year, MonthNum;
    for (calculated_year = start_year; calculated_year - 1< end_year; calculated_year ++)
    {
    for (MonthNum = 0; MonthNum < 14; MonthNum ++)
    {
//--------------
//     Step 1:
 
       Ref_data year(calculated_year,MonthNum);
double YearDays=Tropical_year-(calculated_year-Base_year)*0.0001/100;  //0.0001/100
double half_year=YearDays/2;
double SumDay=(calculated_year-Base_year)*YearDays;    
double Winter_solitice=fmod((SumDay + Qiying_P + int(abs(SumDay))*60), 60);
 
//     推天正閏餘:　置中積，加閏應，滿朔策去之，至不滿之數，為天正閏餘。
double LeapDay = fmod((SumDay + Runying_P + int(abs(SumDay))*Shuoshi_C), Shuoshi_C);
 
//--------------
//     Step 2:
 
//     推天正盈縮:　置半歲周，內減其年閏餘全分，餘為所求天正縮曆。
double RefSun_ = half_year-LeapDay;
double MinSun=MinSun_drift*(calculated_year-Base_year);
double RefSun=RefSun_- MinSun;
        
//     推天正經朔:　置冬至，減閏餘，遇不及減，加紀法減之，為天正經朔。
double RefMoonDay = Winter_solitice-LeapDay;
       if (RefMoonDay < 0){RefMoonDay = RefMoonDay + 60;}
 
//     推天正遲疾:　置中積，加轉應，減去其年閏餘全分，餘滿轉終去之，即天正入轉。        
double RefMoon = fmod((SumDay + Zhuanying_P - LeapDay + int(abs(SumDay))*Zhuanzhong_C), Zhuanzhong_C);
 
//     推天正入交:　置中積，減閏餘，加交應，滿交終去之，即天正入交汎日。        
double RefEclipse = fmod((SumDay + Jiaoying_P - LeapDay + int(abs(SumDay))*Jiaozhong_C), Jiaozhong_C);
 
//--------------
//     Step 3:
 
//     推各月盈縮曆:　置天正縮曆，加二朔策，去半歲周，即得正月經朔下盈曆。
double ReqSun = RefSun + Shuoshi_C*MonthNum;
       ReqSun = fmod(ReqSun, YearDays);            
             
//     推各月經朔:　置天正經朔，加二朔策，滿紀法去之，即得正月經朔。
double ReqMoonDay = fmod((RefMoonDay + Shuoshi_C*MonthNum), 60);
 
//     推各月遲疾曆:　置天正經朔遲疾曆，加二轉差，得正月經朔下遲疾曆。
double ReqMoon = fmod(( RefMoon + MoonRef_C*MonthNum), Zhuanzhong_C);
 
//     推各月入交:　置天正經朔入交汎日加二交差，得正月經朔下入交汎日。
double ReqEclipse = fmod((RefEclipse + Jiao_eclipse_C*MonthNum), Jiaozhong_C);
 
//--------------
//     Step 4:
//     推交常度　置有食之朔入交汎分，以月平行度乘之，即得。
double EclAngle = ReqEclipse * MoonMeanSpeed;
 
//--------------
//     Step 5:
//     如滿半歲周去之交縮，滿半歲周又去之即復交盈。(quoted from Step 3 for 推各月盈縮曆:)
double ChSun = ReqSun;
       if (ReqSun > half_year){ChSun = ReqSun - half_year;}
 
//--------------
//     Step 6:
 
//     推初末限:　視盈曆在盈初縮末限已下，縮曆在縮初盈末限已下，各為初。已上用減半歲周為末。
       double UseSun; int SunStatus;
       year.check_SunStatus(half_year, ReqSun, UseSun, SunStatus);
 
//--------------
//     Step 7:
 
int    IUseSun = int(UseSun);
double SunAdd, SunSum, SunX;
       year.get_SunX(SunStatus,UseSun,IUseSun,SunAdd,SunSum, SunX);
             
//--------------
//     Step 8:
 
double EclA_S; int FrontEclipseSign;
 
    if (SunStatus < LeadStart){EclA_S = EclAngle - SunX;}
    if (SunStatus > LagEnd){EclA_S = EclAngle + SunX;}
 
    if (EclA_S < 7){FrontEclipseSign = 1;}
    if (EclA_S > 342){FrontEclipseSign = 1;}
    if (EclA_S > 175){if (EclA_S < 202){FrontEclipseSign = -1;}}
//     Eclipse conditions:
//     342 < EclA_S ----> FrontEclipseSign = 1
//     202 < EclA_S < 342 ----> No eclipse
//     175 < EclA_S < 202 ----> FrontEclipseSign = -1
//       7 < EclA_S < 175 ----> No eclipse
//     EclA_S < 7 ----> FrontEclipseSign = 1
 
//--------------
//     Step 9:
 
int    MoonStatus = 0;double UseMoon = 0; double MoonCh,MoonAdd, moon_x12_2,SoD, MoonSolar;
       year.get_UseMoon(ReqMoon,UseMoon,SoD,moon_x12_2,MoonSolar,MoonStatus,MoonCh,MoonAdd);
 
//--------------
//     Step 10:
 
double MoonSum, MoonX, MoonSp;
       year.get_MoonX(MoonStatus,MoonCh,MoonAdd,MoonSolar,MoonSum,MoonX,MoonSp);
 
//---------------
//     Step 11:
 
double SunXS, MoonXS, SunMoonT;
       year.get_SunMoonT(SunStatus, MoonStatus, MoonX,SunX, MoonSp,SunXS,MoonXS, SunMoonT);

double SunMoonT24 = SunMoonT * 24;
 
//---------------
//     Step 12:

double Conjunc0 = ReqMoonDay - int(ReqMoonDay);
double Conjunc0_24 = Conjunc0 * 24;
 
double ConjuncT = ReqMoonDay + SunMoonT;
       if (ConjuncT < 0){ConjuncT = 1 + ConjuncT;}
 
//double ConjuncTDay = int(ConjuncT);
double ConjuncTHour = ConjuncT - int(ConjuncT);
double ConjuncT24 = ConjuncTHour * 24;
 
//------------------
//     Step 13:
 
double time_from_noon =(ConjuncT-int(ConjuncT))-0.5;
double observed_time =(0.5-abs(time_from_noon))* time_from_noon/0.96;
double MaxPhase = ConjuncTHour + observed_time;
 
double TNoon = time_from_noon * 24;
double SeeT = observed_time * 24 * 60;
double MaxPhase24 = MaxPhase * 24;
 
//------------------
//     Step 14:
 
//     Calculate half day hours
double Sunrise_time
     = 2068.3 +(0.122928-(IUseSun -53)*0.000608)*(IUseSun - 1)/2 *IUseSun + 0.06*IUseSun ;
 
double Sunset_time = 10000 - Sunrise_time;
 
double HalfDayFan = Sunset_time - 5000;
 
       if (SunStatus == LagEnd){HalfDayFan = 5000 - HalfDayFan;}
       if (SunStatus == LeadStart){HalfDayFan = 5000 - HalfDayFan;}
 
//----------------
//     Step 15:
 
double Net_time_from_noon = time_from_noon + observed_time;
       Net_time_from_noon = abs(Net_time_from_noon*10000);
 
//----------------
//     Step 16:
 
//     推食甚入盈縮曆:　置原得盈縮曆，加入定朔大餘及食甚定分，即得。
double EclRef = ReqSun + int(ConjuncT)+ MaxPhase;
double EclSun = EclRef - ReqMoonDay;
double EclUseSun;
int    EclSunStatus;
       year.check_SunStatus(half_year,EclSun, EclUseSun, EclSunStatus);
 
//----------------
//     Step 17:
 
//     推食甚盈縮差:　依步氣朔求之。
double Ecl_SunAdd, Ecl_SunSum,Ecl_SunX;
int    IEclUseSun = int(EclUseSun);
       year.get_SunX(EclSunStatus, EclUseSun,IEclUseSun, Ecl_SunAdd, Ecl_SunSum,Ecl_SunX);
 
//-----------------
//     Step 18:
 
//     推食甚入盈縮曆行定度:　置食甚入盈縮曆，以盈縮差，盈加縮減之，即得。
double NetSun = 0;
       if (SunStatus < LeadStart){NetSun = EclSun - Ecl_SunX;}
       if (SunStatus > LagEnd){NetSun = EclSun + Ecl_SunX;}
//------------------
//     Step 19:
 
double NS_Ref;
int NetSunStatus;
//     推南北汎差:　視食甚入盈縮曆行定度，在周天象限已下為初限(NS_Ref)，已上與半歲周相減為末限(NS_Ref)。
//               以初末限自之，如一千八百七十度而一，得數，置四度四十六分減之，餘為南北汎差(NS_X)。
       year.check_SunStatus(half_year,NetSun, NS_Ref, NetSunStatus);
double NS_X = 4.46 -((NS_Ref*NS_Ref)/1870);
 
//     推南北定差:　置南北汎差，以距午定分乘之，如半晝分而一，以減汎差，餘為南北定差。若汎差數少，即反減之。
double NS_Sh = NS_X-(NS_X* Net_time_from_noon)/HalfDayFan;
double NS_ = NS_Sh;
 
//------------------
//     Step 20:
 
//     推東西汎差:　置半歲周，減去食甚入盈縮曆行定度，餘以食甚入盈縮曆行定度乘之，以一千八百七十除之為度，即東西汎差(EW_X)。
double EW_Ref = NetSun;
       if (EW_Ref > half_year){EW_Ref = EW_Ref - half_year;}
 
double EW_X =((half_year - EW_Ref)*EW_Ref)/1870;
 
//     推東西定差:　置東西汎差，以距午定分乘之，如二千五百度而一，視得數在東西汎差以下，即為東西定差。
//     若在汎差已上，倍汎差減之，餘為定差。
double EW_H = (EW_X * Net_time_from_noon)/2500;
 
double EW_Sh = EW_H;
       if( EW_H > EW_X ){EW_Sh = 2 * EW_X - EW_H;}
double EW_ = EW_Sh;
 
//------------------
//     Step 21:
 
//     盈曆中前，縮曆中後者，正交減，中交加。盈曆中後，縮曆中前者，正交加，中交減。
int    AfterNoonSign = 1;
       if (time_from_noon < 0){AfterNoonSign = -1;}
 
       switch(SunStatus)
       {
        case LagStart:
                     NS_Sh = FrontEclipseSign * NS_Sh;
                     EW_Sh = -AfterNoonSign * FrontEclipseSign * EW_Sh; break;
 
        case LagEnd:
                     NS_Sh = -FrontEclipseSign * NS_Sh;
                     EW_Sh = -AfterNoonSign * FrontEclipseSign * EW_Sh; break;
 
        case LeadStart:
                     NS_Sh = -FrontEclipseSign * NS_Sh;
                     EW_Sh = AfterNoonSign * FrontEclipseSign * EW_Sh; break;
 
        case LeadEnd:
                     NS_Sh = FrontEclipseSign * NS_Sh;
                     EW_Sh = AfterNoonSign * FrontEclipseSign * EW_Sh; break;
        default:
                     break;
       }
//----------------
//     Step 22:
 
double EclRefA = 0;
double EclipseFrac = 0;
double NetEclA_S = EclA_S;
//     推正交中交定限度:　視日食在正交者置正交度，在中交者置中交度，以南北東西二定差加減之，即得。
       if(FrontEclipseSign == -1){EclRefA = 188.05 + NS_Sh + EW_Sh;}
       if (FrontEclipseSign == 1){EclRefA = 357.64 + NS_Sh + EW_Sh;}
       if (EclA_S < 7){NetEclA_S = EclA_S + 363.7934196;}
 
//-----------------
//     Step 23:
 
double PosA = NetEclA_S - EclRefA;
double Ecl_Pos = FrontEclipseSign * PosA;
 
int    EclipseRef = 6;
       if (Ecl_Pos < 0){EclipseRef = 8;}
 
       EclipseFrac = (EclipseRef - abs(PosA))/EclipseRef;
 
//-----------------
//     Step 24:
 
double EclM = ReqMoon + SunMoonT;
int    EclMoonStatus;
double EclM_Used,EclSoD,EclM_x12_2,EclMoonSolar,EclM_Ch,EclMoonAdd;
       year.get_UseMoon(EclM,EclM_Used,EclSoD,EclM_x12_2,EclMoonSolar,EclMoonStatus,EclM_Ch,EclMoonAdd);
 
double EclMoonSum,EclMoonX,EclM_Sp;
       year.get_MoonX(EclMoonStatus,EclM_Ch,EclMoonAdd,EclMoonSolar, EclMoonSum, EclMoonX, EclM_Sp);
 
//     定限行度:以定限,取立成內行度，遲用遲，疾用疾，內減日行分八分二十秒，得之。
double EclM_SunSp = EclM_Sp - 0.082;
 
//------------------
//     Step 25:
double HalfEclT =(sqrt((20 - EclipseFrac*10)*EclipseFrac*10)*57.40)/EclM_SunSp;
 
double EclBegin_hour = MaxPhase - HalfEclT/10000;
double EclBegin24 = EclBegin_hour * 24;
 
double EclEnd_hour = MaxPhase + HalfEclT/10000;
double EclEnd24 = EclEnd_hour * 24;
 
//MMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMM
if ((calculated_year == 1605)||(calculated_year == 1617)||(calculated_year == 1620)
    ||(calculated_year == 1621)||(calculated_year == 1623))
{
//     LunStep 3:
double LunSun = fmod(ReqSun + Shuoshi_C/2, YearDays);
double LunEcl = fmod(ReqEclipse + Shuoshi_C/2, Jiaozhong_C);
double LunDay = fmod(ReqMoonDay + Shuoshi_C/2, 60);
double LunMoon = fmod(ReqMoon + Shuoshi_C/2, Zhuanzhong_C);
//-------------
//     LunStep 4:
//     推交常度:　置望下入交汎分，乘月平行，如日食法。
double LunEclA = LunEcl * MoonMeanSpeed;
 
//-------------
//     LunStep 5:
double LunChSun = LunSun;
       if (LunSun > half_year){LunChSun = LunSun - half_year;}
 
//-------------
//     LunStep 6:
double LUseSun; int LSunStatus;
year.check_SunStatus(half_year,LunSun, LUseSun, LSunStatus);
 
//-------------
//     LunStep 7:
double LSunX = 0;
year.get_LSunX(LSunStatus,LUseSun,LSunX);
 
//-----------------------------
//     LunStep 8:
//     推交定度:　置交常度，以望下盈縮差盈加縮減之即得。不及減者，加交終度減之。
double LSunXS = -LSunX;
       if (LSunStatus > LagEnd){LSunXS = LSunX;}
double LunEclA_S = LunEclA + LSunXS;
 
//-----------------------------
//     LunStep 9:
double LUseMoon,LSoD,Lmoon_x12_2,LMoonSolar,LMoonCh,LMoonAdd; int LMoonStatus;
year.get_UseMoon(LunMoon,LUseMoon,LSoD,Lmoon_x12_2,LMoonSolar,LMoonStatus, LMoonCh,LMoonAdd);
 
//-----------------------------
//     LunStep 10:
double LMoonX, LMoonSp,LMoonCh_;
double LMoonCh_more = LMoonCh+0.08306;
year.get_LMoonX(LMoonStatus,LMoonCh_more,LMoonX, LMoonSp,LMoonCh_);
 
//-----------------------------
//     LunStep 11:
 
double LMoonXS = -LMoonX;
if (LMoonStatus > LagEnd){LMoonXS = LMoonX;}
 
double LSunMoonX = LSunXS - LMoonXS;
double LSunMoonT=LSunMoonX*0.082/LMoonSp;
 
//------------------------------
//     LunStep 12:
double LConjuncT = LunDay + LSunMoonT;
if (LConjuncT < 0){LConjuncT = 1 + LConjuncT;}
 
int LConjuncTDay=int(LConjuncT);
double LConjuncTime=LConjuncT - int(LConjuncT);
 
//------------------------------
//     LunStep 13:
//     月食：视定望分在日周四分之一已下，为卯前； 已上，覆减半周，为卯后；
//         在四分之三已下，减去半周，为酉前；已上，覆减日周，为酉后。
//         以卯酉前后分自乘，退二位，如四百七十八而一，为时差；
//         子前以减，子后以加，皆加减定望分，为食甚定分；各依发敛求之，即食甚辰刻。
double LUseTime;
       LUseTime = LConjuncTime;
       if(LConjuncTime > 0.25){LUseTime = 0.5 - LConjuncTime;}
       if(LConjuncTime > 0.5){LUseTime = LConjuncTime - 0.5;}
       if(LConjuncTime > 0.75){LUseTime = 1 - LConjuncTime;}
double LunSeeT = LUseTime*LUseTime/4.78;
       if(LConjuncTime > 0.5){LunSeeT = -LunSeeT;}
double LunMaxT = LConjuncTime + LunSeeT;
double LunMaxT24 = LunMaxT * 24;
 
//--------------------------------------------------
//     LunStep 16:
//     推食甚入盈縮曆:　置原得盈縮曆，加入定朔大餘及食甚定分，即得。
double LunEclRef = LunSun + int(LConjuncT)+ LunMaxT;
double LunEclSun = LunEclRef - LunDay;
double LEcl_ChSun = LunEclSun;
       if (LunEclSun > half_year){LEcl_ChSun = LunEclSun - half_year;}
double LEclUseSun;
int    LEclSunStatus;
       year.check_SunStatus(half_year,LunEclSun, LEclUseSun, LEclSunStatus);
 
//-----------------
//     LunStep 17:
 
//     推食甚盈縮差:　依步氣朔求之。
double LEcl_SunX = 0;
year.get_LSunX(LEclSunStatus,LEclUseSun,LEcl_SunX);
 
//-----------------
//     LunStep 18:
 
//     推食甚入盈縮曆行定度:　置食甚入盈縮曆，以盈縮差，盈加縮減之，即得。
double LEcl_SunX_=-LEcl_SunX;
if (LEclSunStatus > LagEnd){LEcl_SunX_=LEcl_SunX;}
double LEclNetSun = LEcl_ChSun + LEcl_SunX_;       
 
double LEclNetSun_= LEclNetSun;
       if(LEclSunStatus > LagEnd){LEclNetSun_=LEclNetSun+182.62875;}

LEclNetSun_=LEclNetSun_+ MinSun;  
    
int total_year=calculated_year-Base_year;
double Added_drift=0.015*total_year;
double Ref_drift=10-Added_drift;
double Ref_drift_check=fmod(SumDay+315.1075,365.2575)-305.1075;
double Next_Ref_drift=Ref_drift_check-4.3445;
 
double Sun_drift=Next_Ref_drift/1.0849;
double Net_drift=4 + Sun_drift;
double LEclNetSun_All=LEclNetSun_+ Net_drift;
 
double LEclNetSun_From=0;
       if(calculated_year==1605){LEclNetSun_From=LEclNetSun_All-277.9775;}
       if(calculated_year==1617){LEclNetSun_From=LEclNetSun_All-51.08;}
       if(calculated_year==1620){LEclNetSun_From=LEclNetSun_All-160.8175;}
       if(calculated_year==1623){LEclNetSun_From=LEclNetSun_All-103.6975;}
//------------------
//     LunStep 22:
 
//求月食入陰陽曆去交前後度
// 視交定度，在交中度(181.8967)已下，為陽曆；已上，減去交中，為陰曆。
// 視入陰陽曆，在後准十五度半已下，為交後度；
// 前准一百六十六度三十九分六十八秒已上，覆減交中，余為交前度及分。
 
double SexLunEclA=LunEclA_S;
if (LunEclA_S > 181.8967){SexLunEclA = 181.8967 - LunEclA_S;}
 
int LunEclFront=-1;
    if(abs(SexLunEclA)<15.5){LunEclFront=1;}
 
double SexLunEclA_= SexLunEclA;
if (LunEclA_S > 166.3968){SexLunEclA_=181.8967-LunEclA_S;LunEclFront=-1;}
if (SexLunEclA < -166.3968){SexLunEclA_=-181.8967-SexLunEclA;LunEclFront=-1;}
 
//-------------------------------------------------
//     LunStep 23:
//求月食分秒
// 視去交前後度，
// 不用南北東西差者。用減食限，不及減者不食。
// 余如定法而一，為月食之分秒。
 
double LunFraction=(13.05-abs(SexLunEclA_))/8.7;
//------------------
//     LunStep 24:
 
double LEclM = LunMoon+LSunMoonT;
int    LEclMoonStatus;
double LEclM_Used,LEclSoD,LEclM_x12_2,LEclMoonSolar,LEclM_Ch,LEclMoonAdd;
       year.get_UseMoon(LEclM,LEclM_Used,LEclSoD,LEclM_x12_2,LEclMoonSolar,LEclMoonStatus,LEclM_Ch,LEclMoonAdd);
 
double LEclMoonX,LEclM_Sp,LEclM_Ch_;
double LEclM_Ch_more=LEclM_Ch + 0.08306;
       year.get_LMoonX(LEclMoonStatus,LEclM_Ch_more,LEclMoonX, LEclM_Sp,LEclM_Ch_);
 
//     定限行度:以定限,取立成內行度，遲用遲，疾用疾，內減日行分八分二十秒，得之。
double LEclM_SunSp = LEclM_Sp - 0.082;
 
//     LunStep 25:
double LHalfEclT=(sqrt((30-LunFraction*10)*LunFraction*10)*49.20)/(LEclM_SunSp*10000);  // Shoushi uses 57.40
double LEclBeginT=LunMaxT - LHalfEclT;
    if(LEclBeginT<0){LEclBeginT=1+LEclBeginT;}
double LEclBeginT24=LEclBeginT*24;
double LEclEndT=LunMaxT+LHalfEclT;
double LEclEndT24=LEclEndT*24;
 
if (((calculated_year==1620)&(MonthNum==13))||(calculated_year==1621))
{
printf("%2d %4d(%2d) LunEcl=%6.4f\tLunMoon=%6.4f\tLunFraction=%6.4f\tLunMaxT24=%6.4fh\tLEclNetSun_=%6.4f\n",
            index,calculated_year,MonthNum,LunEcl,LunMoon,LunFraction,LunMaxT24,LEclNetSun_);
}
//------------------------------------------------
if (((calculated_year==1605)&(MonthNum==3))
||((calculated_year==1617)&(MonthNum==8))
||((calculated_year==1620)&(MonthNum==12))
||((calculated_year==1623)&(MonthNum==10))
||((calculated_year==1621)&(MonthNum==6)))// lunar(about 15 days earlier) and solar elipse occurred in 1621(6)
{
printf("%2d %4d(%2d) ReqEcl(交泛)=%9.6f\tLunEcl(食月望交泛)=%10.7f\n", index,calculated_year,MonthNum,ReqEclipse,LunEcl);
if (calculated_year==1617){print_CalculatedRecords_1617();}
if (calculated_year==1620){print_CalculatedRecords_1620();}
if (calculated_year==1623){print_CalculatedRecords_1623();}
 
cout << endl;
printf("\tSumDay=%8.3f\t[LunStep 3]: LunDay(食月經望)=%10.7f\t\tLunSun=%5.5f\tLunMoon=%5.5f\n",SumDay,LunDay,LunSun,LunMoon);
printf("\tWinter_solitice=%4.3f\t[LunStep 4]: LunEclA(交常度)=%10.8f\n", Winter_solitice,LunEclA);
printf("\tLeapDay=%10.7f\t[LunStep 5]: LunChSun(食月盈縮曆)=%10.8f\n",LeapDay,LunChSun);
printf("\tSumDay + Q=%8.3f\t[LunStep 6]: LUseSun=%10.8f(%1d)\n", SumDay + Qiying_P,LUseSun, LSunStatus);
printf("\tRefMoonDay=%9.6f\t[LunStep 7]: LSunX(盈縮差)=%11.9f\n",RefMoonDay, LSunX);
printf("\tRefSun=%9.6f\t[LunStep 8]: LSunXS=%11.9f \tLunEclA_S(交定度)=%11.9f\n",RefSun,LSunXS,LunEclA_S);
printf("\tRefMoon=%9.6f\t[LunStep 9]: LUseMoon(食月望遲疾曆)=%10.7f(%1d) \tLmoon_x12_2(遲疾限)=%10.7f\tLMoonCh=%5.5f\n",
                                                         RefMoon,LUseMoon,LMoonStatus,Lmoon_x12_2, LMoonCh);
printf("\tRefEclipse=%9.6f\t[LunStep 10]:(+0.08306=>)LMoonCh_=%10.7f\tLMoonX(遲疾差)=%10.8f\tLMoonSp=%9.7f\n",
                                                    RefEclipse,LMoonCh_, LMoonX, LMoonSp);
if (calculated_year==1605)
{
  printf("\t\t\t\t\t\t古今律厯考卷五十(787-563)原計算推遲末限=總限(168.08306限)減遲限(138.3628779限)=29.7201821限\n");
  printf("\t\t\t\t\t\t有別大統以168限或167限減遲限計算,因此 LMoonCh_ 加回0.08306\n");
  printf("\t\t\t\t\t\t該遲行度以29限整數計算,得出1.1823。現改回29.7201821限計算,於是 LMoonSp=1.1814954\n\n");
}
printf("\t\t\t\t[LunStep 11]:LMoonXS=%10.8f\tLSunMoonX=%10.8f \tLSunMoonT(加減差)=%10.8f\n",LMoonXS,LSunMoonX,LSunMoonT);
printf("\tReqMoonDay=%9.6f\t[LunStep 12]:LConjuncT(食月定望)=%10.8f\tLConjuncTDay=%2d\t\tLConjuncTime=%10.8f\n",
          ReqSun,LConjuncT,LConjuncTDay,LConjuncTime);
printf("\tReqSun=%9.6f\t[LunStep 13]:LunSeeT(時差)=%2.7f\tLunMaxT(食甚定分)=%2.8f  =>\tLunMaxT24=%6.7fh",
          UseSun, LunSeeT,LunMaxT,LunMaxT24);
    if(calculated_year==1620){printf("[但丑正初刻略遲]");}
printf("\n");
printf("\tReqMoon=%6.6f\n",ReqMoon);  print_LunStep_13_Extra(calculated_year); 
printf("\n");
printf("\t\t\t\t[LunStep 16]:LunEclSun=%8.5f \tLEcl_ChSun(食甚入盈縮曆)=%10.8f \tLEclUseSun=%10.8f(%1d)\n",
          LunEclSun,LEcl_ChSun,LEclUseSun,LEclSunStatus);
printf("\t\t\t\t[LunStep 17]:LEcl_SunX(食甚盈縮差)=%2.8f\n",LEcl_SunX);
printf("\t\t\t\t[LunStep 18]:LEcl_SunX_=%2.8f \tLEclNetSun(食甚入盈縮曆定度)=%2.8f(%1d)\n",LEcl_SunX_,LEclNetSun,LEclSunStatus);
if(calculated_year == 1605){printf("\t\t\t\t\t     原計算=103.57495347(食甚入盈縮曆)+2.33579615(食甚盈縮差)=105.91074962\n");
                            printf("\t\t\t\t\t\t     但古今律厯考卷五十(787-565)給出錯誤食甚入盈曆定度=105.93074962\n\n");}
printf("\t\t\t\t\t\tLEclNetSun_= LEclNetSun\n");
printf("\t\t\t\t\t\tif(LEclSunStatus > LagEnd){LEclNetSun_=LEclNetSun+182.62875}\n");
printf("\t\t\t\t\t\t\tAs LEclSunStatus=%1d\tLEclNetSun_=%10.8f\n\n",LEclSunStatus,LEclNetSun_);
 
printf("\t\t\t\t\t\ttotal_year=calculated_year-Base_year=%3d => Added_drift=0.015*total_year=%7.5f\n",total_year,Added_drift);
printf("\t\t\t\t\t\tRef_drift(天正加時赤道度)=10-Added_drift=%7.5f(邢雲路計算方法)\n",10-Added_drift);
printf("\t\t\t\t\t\t  (授時曆計算方法:Ref_drift_check=fmod(SumDay+315.1075,365.2575)-305.1075=%7.5f)\n",Ref_drift_check);
printf("\t\t\t\t\t\tNext_Ref_drift=Ref_drift_check-4.3445=%7.5f\n\n",Next_Ref_drift);
 
printf("\t\t\t\t\t\tSun_drift=Next_Ref_drift/1.0849=%7.5f\n",Sun_drift);
printf("\t\t\t\t\t\tNet_drift=4+Sun_drift=%7.5f\n",Net_drift);
printf("\t\t\t\t\t\tLEclNetSun_All=LEclNetSun_ + Net_drift=%7.5f\n",LEclNetSun_All);
 
if(calculated_year==1605){printf("\t\t\t\t\t\tLEclNetSun_From(食甚月躔黃道軫宿)=LEclNetSun_All-277.9775(翼宿)=%7.5f\n",LEclNetSun_From);}
if(calculated_year==1617) {printf("\t\t\t\t\t\tLEclNetSun_From(食甚月躔黃道虚宿)=LEclNetSun_All-51.08(女宿)=%7.5f\n",LEclNetSun_From);}
if(calculated_year==1620)
                {printf("\t\t\t\t\tLEclNetSun_From(食甚月躔黃道畢宿)=LEclNetSun_All-160.8175(昴宿)=%7.5f",LEclNetSun_From);
                 printf("(幾乎=14.3444)\n");
                 printf("\t\t\t\t\t\t\t而授時則推....月離黃道畢宿一十四度三十四分四十四秒(14.3444 與計算值符合=>記錄可靠)\n");
                 printf("\t\t\t\t\t\t  邢雲路...按新法推...食甚月離黃道畢宿一十四度三十分一十秒(14.3010)\n");
                 printf("\t\t\t\t\t\t按日行每天一度,估計授時比新法遲(14.3444-14.3010)*24=1.0416h\n");
                 printf("\t\t\t\t\t\t與食甚遲的時刻比較,差距亦大致相同 LunMaxT24-0.9888h=%1.7fh,表示該些記錄亦可靠\n",LunMaxT24-0.9888);}
if(calculated_year==1623)
                {printf("\t\t\t\t\t\tLEclNetSun_From(食甚月躔黃道奎宿)=LEclNetSun_All-103.6975(壁宿)=%7.5f\n",LEclNetSun_From);}
 
printf("\t\t\t\t[LunStep 22]:\n");
      printf("\t\t\t\t\t求交常交定度:(refer to LunStep 4 and LunStep 8)\n");
      printf("\t\t\t\t\t\t置經朔望入交泛日及分秒(LunEcl)，以月平行度(%10.8f)乘之，為交常度(LunEclA)；\n",MoonMeanSpeed);
      printf("\t\t\t\t\t\t以盈縮差(LSunXS)盈加縮減之，為交定度(LunEclA_S)。\n\n");
 
      printf("\t\t\t\t\t求月食入陰陽曆去交前後度:\n");
      printf("\t\t\t\t\t\t視交定度，在交中度(181.8967)已下，為陽曆；已上，減去交中，為陰曆。\n");
      printf("\t\t\t\t\t\t視入陰陽曆，在後準十五度半已下，為交後度；\n");
      printf("\t\t\t\t\t\t前準一百六十六度三十九分六十八秒已上，覆減交中，餘為交前度及分。\n\n");
      printf("\t\t\t\t\t\t交中度=181.8967\n");
      printf("\t\t\t\t\t\tSexLunEclA=LunEclA_S=%11.9f\n",LunEclA_S);
      printf("\t\t\t\t\t\tif (LunEclA_S > 181.8967){SexLunEclA=181.8967-LunEclA_S}\n");
      printf("\t\t\t\t\t\tResult:SexLunEclA(陰陽曆,負號[-]表示陰曆;正數表示陽曆)=%11.9f\n\n",SexLunEclA);
 
      printf("\t\t\t\t\t\tLunEclFront=-1(負數[-1]表示交前度);\n");
      printf("\t\t\t\t\t\tif(abs(SexLunEclA)<15.5){LunEclFront=1(正數[1]表示交後度);}\n\n");
 
      printf("\t\t\t\t\t\tSexLunEclA_=SexLunEclA=%11.9f\n",SexLunEclA);
      printf("\t\t\t\t\t\tif (LunEclA_S > 166.3968){SexLunEclA_=181.8967-LunEclA_S; LunEclFront=-1}\n");
      printf("\t\t\t\t\t\tif (SexLunEclA < -166.3968){SexLunEclA_=-181.8967-SextLunEclA; LunEclFront=-1}\n");
      printf("\t\t\t\t\t\tResult: SexLunEclA_(LunEclFront) = 月食入陰陽曆(交前後度)= %11.9f(%1d)\n", SexLunEclA_,LunEclFront);
      if (calculated_year == 1617) { printf("\t\t\t\t\t\t邢雲路獻七政真數..以此法布算..陰歷交前度一度二十二分零三秒九十五微=>-1.220395(-1)\n");}
 
printf("\t\t\t\t[LunStep 23]:\n");
                             printf("\t\t\t\t\t求月食分秒:\n");
                             printf("\t\t\t\t\t\t視去交前後度，\n");
                             printf("\t\t\t\t\t\t不用南北東西差者。用減食限，不及減者不食。\n");
                             printf("\t\t\t\t\t\t餘如定法而一，為月食之分秒。\n");
      printf("\t\t\t\t\t\tLunFraction=(13.05-abs(SexLunEclA_))/8.7=%11.9f\n", LunFraction);
      print_LunStep_23_Extra(calculated_year);
 
printf("\t\t\t\t[LunStep 24]:\n");
      printf("\t\t\t\t\tLEclM=%2.8f(%1d)\tLEclM_Used(定入遲疾曆)=%2.8f\n",LEclM, LEclMoonStatus,LEclM_Used);
      printf("\t\t\t\t\tLEclM_x12_2(定入遲疾限)=%2.8f\tLEclM_Ch=%2.8f\tLEclM_Sp=%2.8f\n",LEclM_x12_2,LEclM_Ch,LEclM_Sp);
      printf("\t\t\t\t\tLEclM_SunSp(定限行度)=%2.5f\t    LHalfEclT(定用分)=%2.8f",LEclM_SunSp,LHalfEclT);
      printf("(邢雲路用 4920,不是 5740)\n");
      if (calculated_year == 1617)
      {
          printf("\t\t\t\t\t邢雲路獻七政真數..以此法布算..定用分六刻六十八分九十五秒二十一微(0.06689521)\n");
      }
 
printf("\t\t\t\t[LunStep 25]:\n");
      printf("\t\t\t\t\tLEclBeginT(初虧)=%2.6f\t   LEclEndT(復圓)=%2.6f\n",LEclBeginT,LEclEndT);
      if (calculated_year == 1617){printf("\t\t\t\t\t   (邢雲路的初虧: 0.064[2]2089\t\t    復圓: 0.19801131)\n\n");}
      if (calculated_year == 1620){printf("\t\t邢雲路...按新法推..初虧..計九十七刻二十四分(0.9724)\t\t複圓....計一十一刻(0.11 )\n\n");}
      printf("\t\t\t\t\tLEclBeginT24=%6.8fh",LEclBeginT24);if(calculated_year==1620){printf("[初虧子正一刻]是可靠的");}
      printf("\tLEclEndT24=%6.8fh",LEclEndT24);if(calculated_year==1620){printf("[復圓寅初三刻]略為稍遲");}
      printf("\n");                                                             
      print_LunStep_25_Extra(calculated_year);
cout << "----------------------------------------------------------\n";
 }
} //-----------------
       if (EclipseFrac > 0)
       {
         if (calculated_year >= Theo[index].year)
         {
           if (calculated_year > Theo[index].year) { Shoushi[index].EclipseFrac = -99; index = index + 1;}
 
         if (calculated_year == Theo[index].year)
         {
           if (MonthNum >= Theo[index].month)
           {
             Shoushi[index].year = calculated_year;
             Shoushi[index].month = MonthNum;
 
             Shoushi[index].EclipseFrac = EclipseFrac;
             Shoushi[index].EclipseBegin = EclBegin24;
             Shoushi[index].MaxPhase = MaxPhase24;
             Shoushi[index].EclipseEnd = EclEnd24;
 
             Shoushi[index].Winter_solitice = Winter_solitice;
             Shoushi[index].SunStatus = SunStatus;
 
             Shoushi[index].observed_time = observed_time;
             Shoushi[index].HalfDayFan = HalfDayFan;
             Shoushi[index].ReqSun = ReqSun;
             Shoushi[index].ReqMoon = ReqMoon;
 
             Shoushi[index].MoonStatus = MoonStatus;
             Shoushi[index].ConjuncT24 = ConjuncT24;
             Shoushi[index].observed_time = observed_time;
 
             index = index + 1;
           }
         }
             if (EclipseFrac > 0)
             {
                 if (calculated_year == 1621){print_FindingRecords_1621();}
 
    printf("%2d %4d(%1d)", index,calculated_year,MonthNum);
//-------------------------- Next Row for Step [1] to [6]-------------------------
    printf("\t\t\t\t[2]MinSun=%2.6f\n",MinSun);             
    printf("\t[1](365.)SumDay(WinSol)LeapDay\t[2]RefS(RefMDay)RefM(RefEcl)");
    printf("\t[3]ReqS(ReqMDay)ReqM(ReqEcl)\t[4](EclA)[5]ChSun[6]UseSun(S)\n");
    printf("\t(%5.6f)%5.2f(%5.2f)%5.2f\t  %5.2f(%5.2f)%5.2f(%5.2f)",
           YearDays-365,SumDay,Winter_solitice,LeapDay,RefSun,RefMoonDay,RefMoon,RefEclipse);
    printf("\t   %4.2f(%5.2f)%5.2f(%5.2f)\t   (%5.2f)  %5.2f  %5.2f(%1d)\n\n",
           ReqSun,ReqMoonDay,ReqMoon,ReqEclipse,EclAngle,ChSun,UseSun,SunStatus);
//-------------------------- Next Row for Step [7] to [11]-------------------------
    printf("\t[7]Sun_Add/Sum(X)[8]Ft(EclA_S)\t[9]UseM/SoD(M12_2)MSo(M)");
    printf("\t(MoonCh)MAdd[10]Sum(X)Speed\t[11]SunXS/MoonXS/SunMoonT\n");
    printf("\t%5.2f/%5.2f(%4.2f)%1d(%4.2f)",
             SunAdd,SunSum,SunX,FrontEclipseSign,EclA_S);
    printf("\t %5.2f/%5.2f(%5.2f)%1.2f(%1d)\t(%5.2f)%5.2f",
           UseMoon,SoD,moon_x12_2,MoonSolar,MoonStatus,MoonCh,MoonAdd);
    printf("  %1.2f(%1.2f)%1.4f\t   %1.3f/%1.3f/%5.3f(%5.3fh)\n\n",
           MoonSum,MoonX,MoonSp,SunXS,MoonXS,SunMoonT, SunMoonT24);
//-------------------------- Next Row for Step [12] to [16] ------------------------
    printf("\t[12]Conjunc0(T)  [13]TNoon\tSeeT(MaxPhase) [14]Srise/Sset");
    printf("\tHalfDayFan[15]NetTNoon\t\t[16]Ecl)Ref/Sun(S)UseSun\n");
    printf("\t  %1.4fh(%1.4fh)%5.2fh \t  %5.2fm(%5.4fh)  %4.0f/%4.0f",
           Conjunc0_24, ConjuncT24, TNoon,SeeT,MaxPhase24,Sunrise_time,Sunset_time);
    printf("\t  %5.2f\t%5.2f\t\t   %5.2f/%5.2f(%1d)%5.2f\n\n",
           HalfDayFan,Net_time_from_noon,EclRef,EclSun,EclSunStatus,EclUseSun);
//---------------------------- Next Row for Step [17] to [21]-----------------------
    printf("\t[17]Ecl_)SAdd/SSum(SunX)\t[18]NSun(S)[19]NS_R/NS_X(NS_)");
    printf("\t[20]EW_R/EW_X/EW_H(EW_)\t\t[21]A12_Sgn/NS_Sh/EW_Sh\n");
    printf("\t  %5.2f/%5.2f(%4.2f)\t\t  %5.2f(%1d)  %1.2f/%1.2f/(%1.2f)",
           Ecl_SunAdd,Ecl_SunSum,Ecl_SunX,NetSun,NetSunStatus,NS_Ref, NS_X, NS_);
    printf("\t %1.2f/%1.2f/%1.2f(%1.2f)",EW_Ref,EW_X,EW_H,EW_);
    printf("\t\t\t %1d/%5.2f/%5.2f\n\n",AfterNoonSign,NS_Sh,EW_Sh);
//---------------------------- Next Row for Step [22] to [25]-----------------------
    printf("\t[22]EclRefA/NetEclA_S/PosA\t[23]EclPos(Ref)Frac");
    printf("\t\t[24]EclM(M)_Ch(_Sp)(_SunSp)\t[25]HalfEclT(EclBegin)(EclEnd)\n");
    printf("\t  %5.2f/%5.2f/%5.2f\t\t   %5.2f(%1d)%5.4f",
           EclRefA, NetEclA_S,PosA,Ecl_Pos, EclipseRef,EclipseFrac);
    printf("\t\t %5.2f(%1d)%5.2f(%1.2f)(%1.2f)\t   %5.2f(%5.4fh)(%5.4fh)\n\n",
           EclM,EclMoonStatus,EclM_Ch,EclM_Sp,EclM_SunSp,HalfEclT,EclBegin24,EclEnd24);
    printf("----------------------------------------------\n");
    }
   }
  }
   }//Start_End for loop end(1)
  }//Start_End for loop end(2)
 
     Ref_data Eclipse_Comparison(calculated_year,MonthNum);
     int index_end=index;
     Eclipse_Comparison.print_SolarEclipseTable(index_end);
    return(0);
 }// main end

Show compiler warnings [ + ] Compiler args [ + ] Show input
Compilation time: 1.24 sec, absolute running time: 0.15 sec, cpu time: 0.01 sec, memory peak: 3 Mb, absolute service time: 1.42 sec

Solar eclipse time is in 24 hours Beijing local noon time is 12 as reference month means total months from winter solitice ReqMDay=ReqMoonDay t t t t t t t t t t t t t t t t N year(month) 1 1568(10) [2]MinSun=6.730724 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242213)104824.52(59.58)11.11 164.78(48.46) 8.82(17.85) 94.84(43.77) 1.03(13.82) (184.77) 94.84 87.78(1) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 11.62/23997.74(2.40)-1(182.37) 1.03/ 0.98(12.54)0.04(2) (12.54)10.21 1.34(1.34)1.1984 -2.401/-1.343/-0.256(-6.147h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 18.3859h(12.2387h) 0.24h 7.31m(12.3605h) 2456/7544 2456.06 150.22 138.35/94.59(1)88.03 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 5.06/24009.36(2.40) 92.19(0) 92.19/-0.08/(-0.08) 92.19/4.46/0.27(0.27) 1/-0.08/ 0.27 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 188.24/182.37/-5.87 5.87(6)0.0221 0.77(2) 9.41(1.20)(1.12) 107.29(12.1030h)(12.6180h) ---------------------------------------------- 1 1572(2) [2]MinSun=6.824532 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242209)106285.48(20.54)25.08 150.71(55.46)22.99(22.60) 209.77(54.52)26.94( 0.03) ( 0.34) 27.15 27.15(2) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 370.99/12005.28(1.21)1(1.55) 13.17/13.12(160.63)0.05(1) ( 6.37)10.68 0.80(0.80)1.2031 1.206/0.801/0.137(3.284h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 12.4894h(15.7731h) 3.77h 80.84m(17.1204h) 2119/7881 2118.62 2133.49 264.49/209.97(2)27.35 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 370.99/12005.28(1.21) 211.18(2) 28.56/4.02/(-0.03) 28.56/2.35/2.01(2.01) 1/ 0.03/ 2.01 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 359.68/365.34/ 5.66 5.66(6)0.0562 27.08(1) 4.70(1.20)(1.12) 169.11(16.7145h)(17.5262h) ---------------------------------------------- 2 1572(8) [2]MinSun=6.824532 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242209)106285.48(20.54)25.08 150.71(55.46)22.99(22.60) 21.72(51.70)11.25(13.94) (186.30) 21.72 21.72(0) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 388.29/9228.65(0.95)-1(185.35) 11.25/11.24(137.21)0.01(3) (29.79) 8.51 3.06(3.06)1.0112 -0.951/-3.059/-0.325(-7.804h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 16.8948h(9.0908h)-2.91h -68.87m(7.9429h) 2099/7901 2900.54 1690.45 73.05/21.34(0)21.34 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 388.29/9228.65(0.94) 20.41(0) 20.41/4.24/(1.77) 20.41/1.77/1.20(1.20) -1/-1.77/-1.20 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 185.09/185.35/ 0.27 -0.27(8)0.9668 10.92(3)33.75(1.02)(0.93) 614.28(6.4686h)(9.4172h) ---------------------------------------------- 2 1575(0) [2]MinSun=6.894888 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242206)107381.21(36.27)28.18 147.55( 8.09)13.45(26.74) 147.55( 8.09)13.45(26.74) (357.53) 147.55 35.07(1) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 326.94/14819.79(1.48)1(356.05) 13.45/13.37(164.07)0.08(3) ( 3.00)10.90 0.44(0.44)0.9873 -1.484/-0.439/-0.160(-3.833h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 2.1876h(22.3545h)10.35h 44.37m(23.0940h) 2150/7850 2150.05 4622.51 155.51/147.42(1)35.20 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 326.94/14819.79(1.49) 145.93(1) 36.69/3.74/(-4.30) 145.93/2.86/5.29(0.43) 1/ 4.30/-0.43 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 361.51/356.05/-5.46 -5.46(8)0.3174 13.29(3) 4.87(0.99)(0.91) 462.76(21.9834h)(24.2046h) ---------------------------------------------- 3 1575(6) [2]MinSun=6.894888 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242206)107381.21(36.27)28.18 147.55( 8.09)13.45(26.74) 324.73( 5.27)25.31(13.44) (179.70) 142.11 40.51(3) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 294.76/15773.60(1.59)-1(181.29) 11.53/11.48(140.65)0.05(1) (26.35) 8.90 2.76(2.76)1.1853 1.592/2.761/0.301(7.229h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 6.5929h(13.8215h) 1.82h 48.28m(14.6262h) 2173/7827 2827.25 1094.25 330.34/325.07(3)40.18 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 294.76/15773.60(1.58) 326.65(3) 38.59/3.66/(2.25) 144.03/2.97/1.30(1.30) 1/-2.25/-1.30 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 184.50/181.29/-3.21 3.21(6)0.4650 25.61(1)22.67(1.19)(1.11) 438.02(13.5750h)(15.6774h) ---------------------------------------------- 4 1582(7) [2]MinSun=7.059052 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242199)109937.90(12.96)15.71 159.85(57.25)20.05(10.74) 1.32(23.97) 6.33(26.97) (360.59) 1.32 1.32(0) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 480.41/484.85(0.06)1(360.53) 6.33/ 6.31(77.27)0.02(2) (77.27) 0.84 5.41(5.41)1.1047 -0.064/-5.408/-0.406(-9.747h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 23.2054h(13.4579h) 1.46h 40.02m(14.1250h) 2068/7932 2931.64 885.41 24.91/ 0.95(0) 0.95 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 484.85/ 0.00(0.05) 0.90(0) 0.90/4.46/(3.11) 0.90/0.09/0.03(0.03) 1/ 3.11/-0.03 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 360.72/360.53/-0.19 -0.19(8)0.9759 5.93(2)72.31(1.11)(1.03) 555.62(12.7915h)(15.4585h) ---------------------------------------------- 5 1583(13) [2]MinSun=7.082504 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242198)110303.14(18.20)26.58 148.95(51.62)16.21(11.35) 167.61(15.52)14.35(14.28) (190.89) 167.61 15.01(1) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 434.82/7135.84(0.71)-1(190.18) 0.57/ 0.49( 7.00)0.08(0) ( 7.00)10.63 0.76(0.76)0.9900 -0.714/0.764/0.004(0.099h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 12.4215h(12.5206h) 0.52h 15.56m(12.7799h) 2085/7915 2084.53 324.97 183.14/167.63(1)15.00 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 440.02/6695.81(0.71) 166.91(1) 15.71/4.33/(3.65) 166.91/1.40/0.18(0.18) 1/ 3.65/ 0.18 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 191.89/190.18/-1.71 1.71(6)0.7154 14.36(0) 7.05(0.99)(0.91) 606.00(11.3255h)(14.2343h) ---------------------------------------------- 6 1587(10) [2]MinSun=7.176312 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242194)111764.11(39.17)11.02 164.42(28.15) 4.81(18.42) 94.49(23.46)24.57(14.39) (192.42) 94.49 88.13(1) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 5.06/24009.36(2.40)-1(190.02) 10.79/10.74(131.68)0.05(1) (35.32) 7.85 3.51(3.51)1.1747 -2.401/3.508/0.077(1.855h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 10.9304h(12.7857h) 0.79h 22.95m(13.1681h) 2463/7537 2462.69 486.72 118.04/94.58(1)88.04 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 5.06/24009.36(2.40) 92.18(0) 92.18/-0.08/(-0.07) 92.18/4.46/0.87(0.87) 1/-0.07/ 0.87 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 188.85/190.02/ 1.17 -1.17(8)0.8537 24.65(1)34.38(1.18)(1.09) 519.07(11.9223h)(14.4139h) ---------------------------------------------- 7 1589(2) [2]MinSun=7.223216 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242192)112494.60(49.66) 3.24 172.16(46.41)26.66(21.96) 231.22(45.48) 3.06(26.59) (355.52) 48.60 48.60(2) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 252.82/18628.68(1.88)1(357.40) 3.06/ 3.03(37.30)0.02(2) (37.30) 7.59 3.58(3.58)1.1722 1.878/-3.585/-0.119(-2.866h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 11.4124h(8.5468h)-3.45h -76.86m(7.2658h) 2213/7787 2213.27 1972.56 276.52/231.05(2)48.42 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 252.82/18628.68(1.87) 232.92(2) 50.30/3.11/(0.34) 50.30/3.56/2.81(2.81) -1/-0.34/-2.81 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 354.49/357.40/ 2.91 2.91(6)0.5157 2.94(2)35.85(1.17)(1.09) 459.85(6.1622h)(8.3695h) ---------------------------------------------- 7 1590(2) [2]MinSun=7.246668 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242191)112859.84(54.90)14.12 161.26(40.78)22.82(22.56) 220.32(39.84)26.77(27.20) (363.65) 37.70 37.70(2) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 315.74/15468.08(1.57)1(365.22) 13.00/12.96(158.54)0.04(1) ( 8.46)10.53 1.02(1.02)1.2016 1.569/1.023/0.177(4.245h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 20.2232h(0.4677h)-11.53h -14.05m(0.2336h) 2159/7841 2158.87 4902.66 260.33/220.49(2)37.87 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 315.74/15468.08(1.57) 222.06(2) 39.44/3.63/(-4.61) 39.44/3.02/5.92(0.12) -1/ 4.61/-0.12 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 362.13/365.22/ 3.09 3.09(6)0.4853 26.95(1) 6.30(1.20)(1.12) 438.98(-0.8199h)(1.2872h) ---------------------------------------------- 8 1590(8) [2]MinSun=7.246668 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242191)112859.84(54.90)14.12 161.26(40.78)22.82(22.56) 32.26(37.03)11.08(13.90) (185.82) 32.26 32.26(0) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 334.85/13234.41(1.33)-1(184.49) 11.08/11.07(135.12)0.00(3) (31.88) 8.27 3.23(3.23)1.0136 -1.332/-3.234/-0.369(-8.865h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 0.6286h(15.7636h) 3.76h 80.73m(17.1090h) 2138/7862 2862.47 2128.76 68.97/31.95(0)31.95 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 339.79/12894.61(1.32) 30.63(0) 30.63/3.96/(1.01) 30.63/2.49/2.12(2.12) 1/-1.01/ 2.12 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 189.16/184.49/-4.67 4.67(6)0.2223 10.71(3)36.39(1.02)(0.94) 385.00(16.1850h)(18.0330h) ---------------------------------------------- 8 1594(0) [2]MinSun=7.340476 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242187)114320.80(15.86)28.08 147.20(47.78) 9.44( 0.10) 147.20(47.78) 9.44( 0.10) ( 1.39) 147.20 35.42(1) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 326.94/14819.79(1.50)1(-0.11) 9.44/ 9.43(115.14)0.01(3) (51.86) 5.49 4.61(4.61)1.0414 -1.496/-4.608/-0.481(-11.535h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 18.7321h(7.1975h)-4.80h -90.02m(5.6972h) 2150/7850 2150.05 2626.16 194.43/146.65(1)35.97 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 326.94/14819.79(1.51) 145.14(1) 37.48/3.71/(-0.82) 145.14/2.91/3.06(2.76) -1/ 0.82/ 2.76 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 361.22/363.69/ 2.46 2.46(6)0.5896 8.96(3)57.73(1.05)(0.97) 540.18(4.4008h)(6.9937h) ---------------------------------------------- 9 1594(6) [2]MinSun=7.340476 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242187)114320.80(15.86)28.08 147.20(47.78) 9.44( 0.10) 324.38(44.96)21.29(14.01) (187.35) 141.76 40.86(3) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 294.76/15773.60(1.60)-1(188.95) 7.52/ 7.46(91.71)0.05(1) (75.29) 1.25 5.40(5.40)1.1088 1.603/5.397/0.518(12.424h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 23.1375h(11.5617h)-0.44h -13.20m(11.3417h) 2173/7827 2827.25 274.28 369.85/324.89(3)40.35 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 294.76/15773.60(1.59) 326.48(3) 38.77/3.66/(3.30) 143.86/2.98/0.33(0.33) -1/-3.30/ 0.33 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 185.08/188.95/ 3.88 -3.88(8)0.5152 21.81(1)69.00(1.12)(1.04) 483.09(10.1823h)(12.5011h) ---------------------------------------------- 10 1596(10) [2]MinSun=7.387380 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242185)115051.29(26.35)20.30 154.93( 6.05) 3.73( 3.64) 84.99( 1.35)23.50(26.82) (358.58) 84.99 84.99(0) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 55.73/23718.98(2.38)1(356.20) 9.72/ 9.68(118.58)0.04(1) (48.42) 6.03 4.41(4.41)1.1566 -2.377/4.412/0.144(3.462h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 8.4302h(11.8921h)-0.11h -3.34m(11.8365h) 2436/7564 2563.84 68.14 86.49/85.14(0)85.14 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 49.94/23774.71(2.38) 82.76(0) 82.76/0.80/(0.78) 82.76/4.42/0.12(0.12) -1/ 0.78/ 0.12 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 358.54/356.20/-2.33 -2.33(8)0.7084 23.64(1)46.66(1.16)(1.08) 509.72(10.6131h)(13.0598h) ---------------------------------------------- 11 1603(5) [2]MinSun=7.551544 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242178)117607.98( 3.04) 7.83 167.24(55.21)10.34(14.85) 314.89(22.86)20.22(26.44) (353.51) 132.27 50.35(3) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 243.19/18490.50(1.86)1(355.37) 6.44/ 6.40(78.61)0.05(0) (78.61) 0.56 5.42(5.42)1.0907 1.858/5.418/0.547(13.129h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 20.6372h(9.7658h)-2.23h -56.82m(8.8188h) 2224/7776 2775.88 1325.52 338.26/315.40(3)49.85 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 248.42/18242.08(1.85) 317.24(3) 48.00/3.23/(1.69) 134.62/3.46/1.83(1.83) -1/ 1.69/-1.83 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 357.49/355.37/-2.13 -2.13(8)0.7343 20.77(1)81.72(1.10)(1.01) 546.11(7.5081h)(10.1294h) ---------------------------------------------- 12 1604(5) [2]MinSun=7.574996 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242177)117973.22( 8.28)18.71 156.34(49.57) 6.50(15.46) 303.99(17.23)16.38(27.05) (361.64) 121.37 61.25(3) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 184.59/20874.40(2.09)1(363.73) 2.60/ 2.54(31.77)0.06(0) (31.77) 8.28 3.14(3.14)1.0135 2.092/3.142/0.423(10.164h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 5.4480h(15.6119h) 3.61h 78.90m(16.9268h) 2288/7712 2711.98 2052.85 321.70/304.47(3)60.77 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 190.00/20684.40(2.08) 306.55(3) 58.69/2.62/(0.64) 123.93/3.89/3.19(3.19) 1/ 0.64/ 3.19 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 361.47/363.73/ 2.26 2.26(6)0.6227 16.80(0)36.93(1.02)(0.94) 566.77(15.5666h)(18.2871h) ---------------------------------------------- 12 1605( 3) ReqEcl(交泛)=25.341045 LunEcl(食月望交泛)=12.8941175 SumDay=118338.465 [LunStep 3]: LunDay(食月經望)=56.8288195 LunSun=278.32644 LunMoon=25.33046 Winter_solitice=13.525 [LunStep 4]: LunEclA(交常度)=172.37866026 LeapDay= 0.0532800 [LunStep 5]: LunChSun(食月盈縮曆)=95.70534750 SumDay + Q=118393.525 [LunStep 6]: LUseSun=86.91574050(3) RefMoonDay=13.471744 [LunStep 7]: LSunX(盈縮差)=2.386527878 RefSun=174.969360 [LunStep 8]: LSunXS=2.386527878 LunEclA_S(交定度)=174.765188134 RefMoon= 4.636492 [LunStep 9]: LUseMoon(食月望遲疾曆)=11.5532795(1) Lmoon_x12_2(遲疾限)=140.9500099 LMoonCh=26.04999 RefEclipse=18.385938 [LunStep 10]:(+0.08306=>)LMoonCh_=27.1330501 LMoonX(遲疾差)=2.74268889 LMoonSp=1.1844087 古今律厯考卷五十(787-563)原計算推遲末限=總限(168.08306限)減遲限(138.3628779限)=29.7201821限有別大統以168限或167限減遲限計算,因此 LMoonCh_ 加回0.08306 該遲行度以29限整數計算,得出1.1823。現改回29.7201821限計算,於是 LMoonSp=1.1814954 [LunStep 11]:LMoonXS=-2.74268889 LSunMoonX=5.12921676 LSunMoonT(加減差)=0.35511034 ReqMoonDay=263.561139 [LunStep 12]:LConjuncT(食月定望)=57.18392984 LConjuncTDay=57 LConjuncTime=0.18392984 ReqSun=80.940051 [LunStep 13]:LunSeeT(時差)=0.0070774 LunMaxT(食甚定分)=0.19100728 => LunMaxT24=4.5841748h ReqMoon=10.565167 [LunStep 16]:LunEclSun=278.68862 LEcl_ChSun(食甚入盈縮曆)=96.06753528 LEclUseSun=86.55355272(3) [LunStep 17]:LEcl_SunX(食甚盈縮差)=2.38497926 [LunStep 18]:LEcl_SunX_=2.38497926 LEclNetSun(食甚入盈縮曆定度)=98.45251454(3) 原計算=103.57495347(食甚入盈縮曆)+2.33579615(食甚盈縮差)=105.91074962 但古今律厯考卷五十(787-565)給出錯誤食甚入盈曆定度=105.93074962 LEclNetSun_= LEclNetSun if(LEclSunStatus > LagEnd){LEclNetSun_=LEclNetSun+182.62875} As LEclSunStatus=3 LEclNetSun_=288.67971254 total_year=calculated_year-Base_year=324 => Added_drift=0.015*total_year=4.86000 Ref_drift(天正加時赤道度)=10-Added_drift=5.14000(邢雲路計算方法) (授時曆計算方法:Ref_drift_check=fmod(SumDay+315.1075,365.2575)-305.1075=5.03502) Next_Ref_drift=Ref_drift_check-4.3445=0.69052 Sun_drift=Next_Ref_drift/1.0849=0.63649 Net_drift=4+Sun_drift=4.63649 LEclNetSun_All=LEclNetSun_ + Net_drift=293.31620 LEclNetSun_From(食甚月躔黃道軫宿)=LEclNetSun_All-277.9775(翼宿)=15.33870 [LunStep 22]: 求交常交定度:(refer to LunStep 4 and LunStep 8) 置經朔望入交泛日及分秒(LunEcl)，以月平行度(13.36878311)乘之，為交常度(LunEclA)；以盈縮差(LSunXS)盈加縮減之，為交定度(LunEclA_S)。求月食入陰陽曆去交前後度: 視交定度，在交中度(181.8967)已下，為陽曆；已上，減去交中，為陰曆。視入陰陽曆，在後準十五度半已下，為交後度；前準一百六十六度三十九分六十八秒已上，覆減交中，餘為交前度及分。交中度=181.8967 SexLunEclA=LunEclA_S=174.765188134 if (LunEclA_S > 181.8967){SexLunEclA=181.8967-LunEclA_S} Result:SexLunEclA(陰陽曆,負號[-]表示陰曆;正數表示陽曆)=174.765188134 LunEclFront=-1(負數[-1]表示交前度); if(abs(SexLunEclA)<15.5){LunEclFront=1(正數[1]表示交後度);} SexLunEclA_=SexLunEclA=174.765188134 if (LunEclA_S > 166.3968){SexLunEclA_=181.8967-LunEclA_S; LunEclFront=-1} if (SexLunEclA < -166.3968){SexLunEclA_=-181.8967-SextLunEclA; LunEclFront=-1} Result: SexLunEclA_(LunEclFront) = 月食入陰陽曆(交前後度)= 7.131511866(-1) [LunStep 23]: 求月食分秒: 視去交前後度，不用南北東西差者。用減食限，不及減者不食。餘如定法而一，為月食之分秒。 LunFraction=(13.05-abs(SexLunEclA_))/8.7=0.680285992 [LunStep 24]: LEclM=25.68557384(1) LEclM_Used(定入遲疾曆)=11.90838984 LEclM_x12_2(定入遲疾限)=145.28235601 LEclM_Ch=21.71764399 LEclM_Sp=1.18899494 LEclM_SunSp(定限行度)=1.10699 LHalfEclT(定用分)=0.05583191(邢雲路用 4920,不是 5740) [LunStep 25]: LEclBeginT(初虧)=0.135175 LEclEndT(復圓)=0.246839 LEclBeginT24=3.24420901h LEclEndT24=5.92414049h ---------------------------------------------------------- 13 1607(3) [2]MinSun=7.645352 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242174)119068.95(24.01)21.80 153.17( 2.21)24.51(19.60) 241.76(30.80) 2.89(26.56) (355.04) 59.14 59.14(2) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 187.65/21085.95(2.11)1(357.15) 2.89/ 2.87(35.21)0.02(2) (35.21) 7.86 3.42(3.42)1.1749 2.111/-3.422/-0.091(-2.195h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 19.1461h(16.9514h) 4.95h 90.89m(18.4662h) 2276/7724 2275.93 2694.25 272.53/241.74(2)59.12 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 187.65/21085.95(2.11) 243.85(2) 61.23/2.46/(-0.45) 61.23/3.97/4.28(3.67) 1/ 0.45/ 3.67 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 361.76/357.15/-4.60 -4.60(8)0.4244 2.79(2)34.09(1.18)(1.09) 428.93(17.4368h)(19.4956h) ---------------------------------------------- 13 1609(8) [2]MinSun=7.692256 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242172)119799.43(34.49)14.02 160.91(20.47)18.81(23.14) 31.91(16.72) 7.06(14.47) (193.47) 31.91 31.91(0) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 339.79/12894.61(1.32)-1(192.15) 7.06/ 7.05(86.18)0.01(3) (80.82) 0.09 5.43(5.43)1.0954 -1.320/-5.429/-0.505(-12.126h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 17.1731h(5.0473h)-6.95h -91.39m(3.5242h) 2134/7866 2866.46 3531.58 48.06/31.34(0)31.34 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 339.79/12894.61(1.30) 30.04(0) 30.04/3.98/(-0.92) 30.04/2.45/3.46(1.44) -1/ 0.92/-1.44 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 187.53/192.15/ 4.62 -4.62(8)0.4228 6.56(2)80.00(1.10)(1.02) 460.92(2.4180h)(4.6304h) ---------------------------------------------- 15 1610(13) [2]MinSun=7.715708 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242171)120164.67(39.73)24.90 150.01(14.84)14.97(23.75) 168.66(38.74)13.11(26.67) (356.57) 168.66 13.96(1) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 445.20/6250.61(0.67)1(355.90) 13.11/13.04(159.88)0.07(3) ( 7.12)10.63 0.88(0.88)0.9901 -0.668/-0.881/-0.128(-3.077h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 17.6550h(14.5777h) 2.58h 63.25m(15.6318h) 2081/7919 2080.57 1513.26 207.32/168.58(1)14.04 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 440.02/6695.81(0.67) 167.91(1) 14.71/4.34/(1.18) 167.91/1.32/0.80(0.80) 1/-1.18/-0.80 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 355.66/355.90/ 0.25 0.25(6)0.9588 12.98(3) 8.68(0.99)(0.91) 630.79(14.1179h)(17.1457h) ---------------------------------------------- 15 1612(6) [2]MinSun=7.762612 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242169)120895.16(50.22)17.12 157.74(33.10) 9.27( 0.07) 334.93(30.29)21.12(13.98) (186.87) 152.31 30.32(3) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 344.71/12549.90(1.27)-1(188.14) 7.35/ 7.30(89.62)0.05(1) (77.38) 0.82 5.42(5.42)1.1045 1.266/5.417/0.496(11.907h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 6.8712h(18.7785h) 6.78h 92.17m(20.3147h) 2130/7870 2870.34 3464.46 365.77/335.49(3)29.76 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 349.62/12200.28(1.25) 336.73(3) 28.51/4.03/(-0.83) 154.11/2.35/3.26(1.44) 1/ 0.83/-1.44 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 187.44/188.14/ 0.70 -0.70(8)0.9130 21.62(1)71.33(1.12)(1.03) 552.64(18.9884h)(21.6410h) ---------------------------------------------- 16 1613(6) [2]MinSun=7.786064 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242168)121260.40(55.46)27.99 146.85(27.47) 5.43( 0.68) 324.03(24.65)17.28(14.59) (195.00) 141.41 41.21(3) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 289.68/16068.36(1.61)-1(196.61) 3.51/ 3.44(42.77)0.06(0) (42.77) 6.85 3.98(3.98)1.0278 1.613/3.984/0.447(10.718h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 15.6820h(2.3996h)-9.60h -59.99m(1.3997h) 2178/7822 2822.46 4416.78 349.09/324.43(3)40.81 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 294.76/15773.60(1.60) 326.03(3) 39.21/3.64/(-2.05) 143.41/3.01/5.31(0.70) -1/ 2.05/ 0.70 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 190.81/196.61/ 5.81 -5.81(8)0.2740 17.73(0)48.22(1.04)(0.95) 413.91(0.4063h)(2.3931h) ---------------------------------------------- 17 1615(4) [2]MinSun=7.832968 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242166)121990.88( 5.94)20.21 154.58(45.73)27.28( 4.21) 272.70(43.86) 7.63(13.48) (180.27) 90.08 92.54(3) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 8.90/24001.62(2.40)-1(182.67) 7.63/ 7.63(93.06)0.00(3) (73.94) 1.52 5.38(5.38)1.0811 2.401/-5.380/-0.226(-5.424h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 20.5693h(15.1458h) 3.15h 72.54m(16.3547h) 2489/7511 2510.86 1814.47 316.38/272.53(3)92.72 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 8.90/24001.62(2.40) 274.93(3) 90.32/0.10/(0.03) 92.31/4.46/3.24(3.24) 1/-0.03/-3.24 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 184.79/182.67/-2.12 2.12(6)0.6471 7.40(3)76.70(1.09)(1.00) 534.55(15.0718h)(17.6376h) ---------------------------------------------- 17 1617( 8) ReqEcl(交泛)=26.292590 LunEcl(食月望交泛)=13.8456625 《大明神宗顯皇帝實錄卷之五百四十七》邢雲路獻七政真數..(望交泛分:13.845663) 以此法布算今歲七月十六日戊寅夜望月食,推得是月: 望交泛分一十三日八十四刻五十六分六十三秒(13.845663) 陰歷交前度一度二十二分零三秒九十五微(1.220395) 月食分一十三分五十九秒六十九微(1.35969) 定用分六刻六十八分九十五秒二十一微(0.06689521) 既內分二刻八十六分六十五秒既外分三刻八十一分七十九秒初虧分六刻四十三分二十秒八十九微(0.064[3]2089,應修正為:復圓分-2*定用分=0.19801131-0.06689521*2=0.064[2]2089) 食既分一十刻二十五分零四十微食甚分一十三刻一十一分六十六秒一十微(0.1311[6]610,應修正為:復圓分-定用分=0.19801131-0.06689521=0.1311[1]610) 生光分一十五刻九十八分三十一秒八十微復圓分一十九刻八十分一十一秒三十一微(0.19801131) 初虧丑初二刻(1.48h<=>0.064[2]2089*24=1.54130136,仍然算是丑初二刻,所以該記錄是可靠的) 食既丑正一刻,食甚寅初二刻(3.48h? <=>0.1311[1]610*24=3.1467864h=>應修正為寅初一刻) 生光寅初三刻,復明寅正三刻(4.72h<=>0.19801131*24=4.75227144h=寅正三刻,所以該記錄亦是可靠的) SumDay=122721.367 [LunStep 3]: LunDay(食月經望)=15.0095485 LunSun=48.08149 LunMoon=24.59484 Winter_solitice=16.427 [LunStep 4]: LunEclA(交常度)=185.09965897 LeapDay=12.4275960 [LunStep 5]: LunChSun(食月盈縮曆)=48.08149050 SumDay + Q=122776.427 [LunStep 6]: LUseSun=48.08149050(0) RefMoonDay= 3.999508 [LunStep 7]: LSunX(盈縮差)=1.800930527 RefSun=162.313614 [LunStep 8]: LSunXS=-1.800930527 LunEclA_S(交定度)=183.298728443 RefMoon=21.574112 [LunStep 9]: LUseMoon(食月望遲疾曆)=10.8176565(1) Lmoon_x12_2(遲疾限)=131.9754093 LMoonCh=35.02459 RefEclipse= 7.745638 [LunStep 10]:(+0.08306=>)LMoonCh_=36.1076507 LMoonX(遲疾差)=3.49220640 LMoonSp=1.1737438 [LunStep 11]:LMoonXS=-3.49220640 LSunMoonX=1.69127587 LSunMoonT(加減差)=0.11815579 ReqMoonDay=33.316194 [LunStep 12]:LConjuncT(食月定望)=15.12770429 LConjuncTDay=15 LConjuncTime=0.12770429 ReqSun=33.316194 [LunStep 13]:LunSeeT(時差)=0.0034118 LunMaxT(食甚定分)=0.13111609 => LunMaxT24=3.1467861h ReqMoon=9.829544 邢雲路食甚分 0.13111610*24=3.1467864h 理論值:3.219 [LunStep 16]:LunEclSun=48.20306 LEcl_ChSun(食甚入盈縮曆)=48.20305809 LEclUseSun=48.20305809(0) [LunStep 17]:LEcl_SunX(食甚盈縮差)=1.80403655 [LunStep 18]:LEcl_SunX_=-1.80403655 LEclNetSun(食甚入盈縮曆定度)=46.39902153(0) LEclNetSun_= LEclNetSun if(LEclSunStatus > LagEnd){LEclNetSun_=LEclNetSun+182.62875} As LEclSunStatus=0 LEclNetSun_=54.27889353 total_year=calculated_year-Base_year=336 => Added_drift=0.015*total_year=5.04000 Ref_drift(天正加時赤道度)=10-Added_drift=4.96000(邢雲路計算方法) (授時曆計算方法:Ref_drift_check=fmod(SumDay+315.1075,365.2575)-305.1075=4.84710) Next_Ref_drift=Ref_drift_check-4.3445=0.50260 Sun_drift=Next_Ref_drift/1.0849=0.46327 Net_drift=4+Sun_drift=4.46327 LEclNetSun_All=LEclNetSun_ + Net_drift=58.74217 LEclNetSun_From(食甚月躔黃道虚宿)=LEclNetSun_All-51.08(女宿)=7.66217 [LunStep 22]: 求交常交定度:(refer to LunStep 4 and LunStep 8) 置經朔望入交泛日及分秒(LunEcl)，以月平行度(13.36878311)乘之，為交常度(LunEclA)；以盈縮差(LSunXS)盈加縮減之，為交定度(LunEclA_S)。求月食入陰陽曆去交前後度: 視交定度，在交中度(181.8967)已下，為陽曆；已上，減去交中，為陰曆。視入陰陽曆，在後準十五度半已下，為交後度；前準一百六十六度三十九分六十八秒已上，覆減交中，餘為交前度及分。交中度=181.8967 SexLunEclA=LunEclA_S=183.298728443 if (LunEclA_S > 181.8967){SexLunEclA=181.8967-LunEclA_S} Result:SexLunEclA(陰陽曆,負號[-]表示陰曆;正數表示陽曆)=-1.402028443 LunEclFront=-1(負數[-1]表示交前度); if(abs(SexLunEclA)<15.5){LunEclFront=1(正數[1]表示交後度);} SexLunEclA_=SexLunEclA=-1.402028443 if (LunEclA_S > 166.3968){SexLunEclA_=181.8967-LunEclA_S; LunEclFront=-1} if (SexLunEclA < -166.3968){SexLunEclA_=-181.8967-SextLunEclA; LunEclFront=-1} Result: SexLunEclA_(LunEclFront) = 月食入陰陽曆(交前後度)= -1.402028443(-1) 邢雲路獻七政真數..以此法布算..陰歷交前度一度二十二分零三秒九十五微=>-1.220395(-1) [LunStep 23]: 求月食分秒: 視去交前後度，不用南北東西差者。用減食限，不及減者不食。餘如定法而一，為月食之分秒。 LunFraction=(13.05-abs(SexLunEclA_))/8.7=1.338847305 邢雲路獻七政真數..以此法布算..月食分一十三分五十九秒六十九微(1.35969) 理論值:1.403 [LunStep 24]: LEclM=24.71299629(1) LEclM_Used(定入遲疾曆)=10.93581229 LEclM_x12_2(定入遲疾限)=133.41690993 LEclM_Ch=33.58309007 LEclM_Sp=1.17556264 LEclM_SunSp(定限行度)=1.09356 LHalfEclT(定用分)=0.06709524(邢雲路用 4920,不是 5740) 邢雲路獻七政真數..以此法布算..定用分六刻六十八分九十五秒二十一微(0.06689521) [LunStep 25]: LEclBeginT(初虧)=0.064021 LEclEndT(復圓)=0.198211 (邢雲路的初虧: 0.064[2]2089 復圓: 0.19801131) LEclBeginT24=1.53650019h LEclEndT24=4.75707192h (邢雲路的初虧*24=1.54130136h=>算是丑初二刻 4.75227144h=寅正三刻,兩項記錄都可靠) 理論值: 1.49h 4.95h 邢雲路原有的陰陽曆 SexLunEclA_1617(不分陰陽)=13.05-1.35969*8.7=1.22070,與記錄所指的 1.220395 基本一致邢雲路變法後的食分 NewFraction=(13.125-1.220395)/7.5=1.58728,大於理論值理論值=1.403 邢雲路原有的食分: 1.35969 邢雲路原有的定限行度 LEclM_SunSp_1617=(sqrt((30-13.5969)*13.5969)*49.20)/(0.06689521*10000)=1.09838 邢雲路變法後改大的定用分 NewHalfT=(sqrt((35-NewFraction*10)*NewFraction*10)*49.20)/(LEclM_SunSp_1617*10000)=0.07805 邢雲路變法後的 NewBeginT=0.13111610-NewHalfT=0.05307 NewEndT=0.13111610+NewHalfT=0.20916 變法後的初虧=NewBeginT*24=1.27362 食甚不變=0.13111610*24=3.14679 復圓=NewEndT*24=5.01995 理論值: 初虧=1.49h 食甚=3.219h 復圓=4.95h 邢雲路原有的初虧=1.54130136h 復圓=4.75227144h ---------------------------------------------------------- 18 1617(8) [2]MinSun=7.879872 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242164)122721.37(16.43)12.43 162.31( 4.00)21.57( 7.75) 33.32( 0.24) 9.83(26.29) (351.50) 33.32 33.32(0) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 329.90/13569.26(1.37)1(350.13) 9.83/ 9.76(119.92)0.07(3) (47.08) 6.23 4.33(4.33)1.0340 -1.367/-4.330/-0.452(-10.844h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 5.8620h(19.0180h) 7.02h 91.05m(20.5355h) 2142/7858 2858.39 3556.48 34.17/33.93(0)33.93 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 329.90/13569.26(1.39) 32.54(0) 32.54/3.89/(-0.95) 32.54/2.61/3.72(1.51) 1/-0.95/-1.51 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 355.18/350.13/-5.05 -5.05(8)0.3689 9.38(3)52.59(1.04)(0.96) 463.56(19.4230h)(21.6481h) ---------------------------------------------- 18 1620(12) ReqEcl(交泛)=12.496823 LunEcl(食月望交泛)= 0.0498955 《大明熹宗哲皇帝實錄卷之七》按新法推泰昌元年庚申歲十一月十六日己丑夜望月食: 初虧漏下二百七十三籌五十三分,計九十七刻二十四分(0.9724) 食既漏下二百九十四籌六十分,計八十八分生光漏下三百三十二籌四十分,計七刻四十分復圓漏下三百五十三籌四十八分,計一十一刻(0.11) 各以發斂求之,得: 初虧夜子初一刻(23.24h) 食既子正初刻食甚子正四刻(0.96h) 生光丑初三刻復圓丑正二刻(2.48h,但見有改自三刻的痕跡) 月食一十四分九十九秒(1.499) 食甚月離黃道畢宿一十四度三十分一十秒(14.3010) 而授時則推: 初虧子正一刻(0.24h) 食既丑初一刻食甚丑正初刻(2.00h?) 生光丑正三刻復圓寅初三刻(3.72h?) 月食一十三分三十一秒(1.331) 月離黃道畢宿一十四度三十四分四十四秒(14.3444),與天不合 SumDay=123817.093 [LunStep 3]: LunDay(食月經望)=25.7638615 LunSun=163.03997 LunMoon=22.95696 Winter_solitice=32.153 [LunStep 4]: LunEclA(交常度)=0.66704211 LeapDay=15.5211300 [LunStep 5]: LunChSun(食月盈縮曆)=163.03997400 SumDay + Q=123872.153 [LunStep 6]: LUseSun=19.58110650(1) RefMoonDay=16.631449 [LunStep 7]: LSunX(盈縮差)=0.908488692 RefSun=159.149722 [LunStep 8]: LSunXS=-0.908488692 LunEclA_S(交定度)=-0.241446579 RefMoon=12.031333 [LunStep 9]: LUseMoon(食月望遲疾曆)= 9.1797775(1) Lmoon_x12_2(遲疾限)=111.9932855 LMoonCh=55.00671 RefEclipse=11.888619 [LunStep 10]:(+0.08306=>)LMoonCh_=56.0897745 LMoonX(遲疾差)=4.77402895 LMoonSp=1.1443566 [LunStep 11]:LMoonXS=-4.77402895 LSunMoonX=3.86554026 LSunMoonT(加減差)=0.27698911 ReqMoonDay=148.274677 [LunStep 12]:LConjuncT(食月定望)=26.04085061 LConjuncTDay=26 LConjuncTime=0.04085061 ReqSun=34.346403 [LunStep 13]:LunSeeT(時差)=0.0003491 LunMaxT(食甚定分)=0.04119972 => LunMaxT24=0.9887933h[但丑正初刻略遲] ReqMoon=8.191665 邢雲路..按新法推..初虧..計九十七刻二十四分..復圓..計一十一刻=>食甚分=0.11-(0.11+(1-0.9724))/2=0.0412=>食甚=0.0412*24=0.9888h=子正四刻,與記錄一致理論值:1.17h [LunStep 16]:LunEclSun=163.31731 LEcl_ChSun(食甚入盈縮曆)=163.31731222 LEclUseSun=19.30376828(1) [LunStep 17]:LEcl_SunX(食甚盈縮差)=0.89700279 [LunStep 18]:LEcl_SunX_=-0.89700279 LEclNetSun(食甚入盈縮曆定度)=162.42030943(1) LEclNetSun_= LEclNetSun if(LEclSunStatus > LagEnd){LEclNetSun_=LEclNetSun+182.62875} As LEclSunStatus=1 LEclNetSun_=170.37053743 total_year=calculated_year-Base_year=339 => Added_drift=0.015*total_year=5.08500 Ref_drift(天正加時赤道度)=10-Added_drift=4.91500(邢雲路計算方法) (授時曆計算方法:Ref_drift_check=fmod(SumDay+315.1075,365.2575)-305.1075=4.80008) Next_Ref_drift=Ref_drift_check-4.3445=0.45558 Sun_drift=Next_Ref_drift/1.0849=0.41993 Net_drift=4+Sun_drift=4.41993 LEclNetSun_All=LEclNetSun_ + Net_drift=174.79046 LEclNetSun_From(食甚月躔黃道畢宿)=LEclNetSun_All-160.8175(昴宿)=13.97296(幾乎=14.3444) 而授時則推....月離黃道畢宿一十四度三十四分四十四秒(14.3444 與計算值符合=>記錄可靠) 邢雲路...按新法推...食甚月離黃道畢宿一十四度三十分一十秒(14.3010) 按日行每天一度,估計授時比新法遲(14.3444-14.3010)*24=1.0416h 與食甚遲的時刻比較,差距亦大致相同 LunMaxT24-0.9888h=-0.0000067h,表示該些記錄亦可靠 [LunStep 22]: 求交常交定度:(refer to LunStep 4 and LunStep 8) 置經朔望入交泛日及分秒(LunEcl)，以月平行度(13.36878311)乘之，為交常度(LunEclA)；以盈縮差(LSunXS)盈加縮減之，為交定度(LunEclA_S)。求月食入陰陽曆去交前後度: 視交定度，在交中度(181.8967)已下，為陽曆；已上，減去交中，為陰曆。視入陰陽曆，在後準十五度半已下，為交後度；前準一百六十六度三十九分六十八秒已上，覆減交中，餘為交前度及分。交中度=181.8967 SexLunEclA=LunEclA_S=-0.241446579 if (LunEclA_S > 181.8967){SexLunEclA=181.8967-LunEclA_S} Result:SexLunEclA(陰陽曆,負號[-]表示陰曆;正數表示陽曆)=-0.241446579 LunEclFront=-1(負數[-1]表示交前度); if(abs(SexLunEclA)<15.5){LunEclFront=1(正數[1]表示交後度);} SexLunEclA_=SexLunEclA=-0.241446579 if (LunEclA_S > 166.3968){SexLunEclA_=181.8967-LunEclA_S; LunEclFront=-1} if (SexLunEclA < -166.3968){SexLunEclA_=-181.8967-SextLunEclA; LunEclFront=-1} Result: SexLunEclA_(LunEclFront) = 月食入陰陽曆(交前後度)= -0.241446579(1) [LunStep 23]: 求月食分秒: 視去交前後度，不用南北東西差者。用減食限，不及減者不食。餘如定法而一，為月食之分秒。 LunFraction=(13.05-abs(SexLunEclA_))/8.7=1.472247520 授時則推..月食一十三分三十一秒(1.331) 邢雲路.....按新法推....月食一十四分九十九秒(1.499) 理論值:1.6 [LunStep 24]: LEclM=23.23395061(1) LEclM_Used(定入遲疾曆)=9.45676661 LEclM_x12_2(定入遲疾限)=115.37255260 LEclM_Ch=51.62744740 LEclM_Sp=1.14987339 LEclM_SunSp(定限行度)=1.06787 LHalfEclT(定用分)=0.06909749(邢雲路用 4920,不是 5740) [LunStep 25]: LEclBeginT(初虧)=0.972102 LEclEndT(復圓)=0.110297 邢雲路...按新法推..初虧..計九十七刻二十四分(0.9724) 複圓....計一十一刻(0.11 ) LEclBeginT24=23.33045365h[初虧子正一刻]是可靠的 LEclEndT24=2.64713302h[復圓寅初三刻]略為稍遲 (邢雲路...按新法推..初虧=0.9724*24 = 23.3376h 復圓=0.11*24=2.64h) 理論值:初虧= 23.28h 復圓=3.05h 邢雲路原有的陰陽曆 SexLunEclA_1620(不分陰陽)=13.05-1.499*8.7=0.00870 邢雲路變法後的食分 NewFraction=(13.125-SexLunEclA_1620)/7.5=1.74884 邢雲路原有的定用分 HalfT_1620=(0.11+(1-0.9724))/2=0.06880 邢雲路原有的定限行度 LEclM_SunSp_1620=(sqrt((30-14.99)*14.99)*49.20)/(HalfT_1620*10000)=1.07267 邢雲路變法後的定用分 NewHalfT=(sqrt((35-NewFraction*10)*NewFraction*10)*49.20)/(LEclM_SunSp_1620*10000)=0.08027 邢雲路變法後的 NewBeginT=1+0.0412-NewHalfT=0.96093 NewEndT=0.0412+NewHalfT=0.12147 邢雲路變法後的初虧=NewBeginT*24=23.06240 食甚不變=0.0412*24=0.98880 復圓=NewEndT*24=2.91520 理論值:初虧= 23.28h 食甚=1.17h 復圓=3.05h 邢雲路原有的初虧=0.9724*24=23.33760 復圓=0.11*24=2.64000h ---------------------------------------------------------- 18 1620(13) LunEcl=2.3683 LunMoon=24.9332 LunFraction=-2.1950 LunMaxT24=13.1410h LEclNetSun_=201.2793 18 1621( 0) LunEcl=0.0499 LunMoon=22.9570 LunFraction=1.4721 LunMaxT24=0.9871h LEclNetSun_=170.3698 18 1621( 1) LunEcl=2.3683 LunMoon=24.9332 LunFraction=-2.1949 LunMaxT24=13.1392h LEclNetSun_=201.2785 18 1621( 2) LunEcl=4.6866 LunMoon=26.9094 LunFraction=-5.8883 LunMaxT24=23.8593h LEclNetSun_=231.8570 18 1621( 3) LunEcl=7.0050 LunMoon=1.3313 LunFraction=-9.5249 LunMaxT24=9.4888h LEclNetSun_=261.8989 18 1621( 4) LunEcl=9.3234 LunMoon=3.3075 LunFraction=-13.0994 LunMaxT24=18.6015h LEclNetSun_=291.3826 18 1621( 5) LunEcl=11.6417 LunMoon=5.2837 LunFraction=-16.6140 LunMaxT24=4.6552h LEclNetSun_=320.3892 18 1621( 6) LunEcl=13.9601 LunMoon=7.2599 LunFraction=0.8322 LunMaxT24=14.8351h LEclNetSun_=348.9415 18 1621( 6) ReqEcl(交泛)=26.407037 LunEcl(食月望交泛)=13.9601095 SumDay=124182.334 [LunStep 3]: LunDay(食月經望)=22.9474195 LunSun=340.20042 LunMoon=7.25994 Winter_solitice=37.394 [LunStep 4]: LunEclA(交常度)=186.62967609 LeapDay=26.3958350 [LunStep 5]: LunChSun(食月盈縮曆)=157.57933950 SumDay + Q=124237.394 [LunStep 6]: LUseSun=25.04174050(3) RefMoonDay=10.998565 [LunStep 7]: LSunX(盈縮差)=1.076856479 RefSun=148.251565 [LunStep 8]: LSunXS=1.076856479 LunEclA_S(交定度)=187.706532570 RefMoon= 8.191665 [LunStep 9]: LUseMoon(食月望遲疾曆)= 7.2599435(3) Lmoon_x12_2(遲疾限)=88.5713107 LMoonCh=78.42869 RefEclipse=12.496823 [LunStep 10]:(+0.08306=>)LMoonCh_=79.5117493 LMoonX(遲疾差)=5.42352139 LMoonSp=1.0926015 [LunStep 11]:LMoonXS=5.42352139 LSunMoonX=-4.34666491 LSunMoonT(加減差)=-0.32621823 ReqMoonDay=325.435123 [LunStep 12]:LConjuncT(食月定望)=22.62120127 LConjuncTDay=22 LConjuncTime=0.62120127 ReqSun=39.807037 [LunStep 13]:LunSeeT(時差)=-0.0030732 LunMaxT(食甚定分)=0.61812810 => LunMaxT24=14.8350745h ReqMoon=20.049015 [LunStep 16]:LunEclSun=339.87113 LEcl_ChSun(食甚入盈縮曆)=157.25004810 LEclUseSun=25.37103190(3) [LunStep 17]:LEcl_SunX(食甚盈縮差)=1.08905677 [LunStep 18]:LEcl_SunX_=1.08905677 LEclNetSun(食甚入盈縮曆定度)=158.33910487(3) LEclNetSun_= LEclNetSun if(LEclSunStatus > LagEnd){LEclNetSun_=LEclNetSun+182.62875} As LEclSunStatus=3 LEclNetSun_=348.94153487 total_year=calculated_year-Base_year=340 => Added_drift=0.015*total_year=5.10000 Ref_drift(天正加時赤道度)=10-Added_drift=4.90000(邢雲路計算方法) (授時曆計算方法:Ref_drift_check=fmod(SumDay+315.1075,365.2575)-305.1075=4.78440) Next_Ref_drift=Ref_drift_check-4.3445=0.43990 Sun_drift=Next_Ref_drift/1.0849=0.40548 Net_drift=4+Sun_drift=4.40548 LEclNetSun_All=LEclNetSun_ + Net_drift=353.34701 [LunStep 22]: 求交常交定度:(refer to LunStep 4 and LunStep 8) 置經朔望入交泛日及分秒(LunEcl)，以月平行度(13.36878311)乘之，為交常度(LunEclA)；以盈縮差(LSunXS)盈加縮減之，為交定度(LunEclA_S)。求月食入陰陽曆去交前後度: 視交定度，在交中度(181.8967)已下，為陽曆；已上，減去交中，為陰曆。視入陰陽曆，在後準十五度半已下，為交後度；前準一百六十六度三十九分六十八秒已上，覆減交中，餘為交前度及分。交中度=181.8967 SexLunEclA=LunEclA_S=187.706532570 if (LunEclA_S > 181.8967){SexLunEclA=181.8967-LunEclA_S} Result:SexLunEclA(陰陽曆,負號[-]表示陰曆;正數表示陽曆)=-5.809832570 LunEclFront=-1(負數[-1]表示交前度); if(abs(SexLunEclA)<15.5){LunEclFront=1(正數[1]表示交後度);} SexLunEclA_=SexLunEclA=-5.809832570 if (LunEclA_S > 166.3968){SexLunEclA_=181.8967-LunEclA_S; LunEclFront=-1} if (SexLunEclA < -166.3968){SexLunEclA_=-181.8967-SextLunEclA; LunEclFront=-1} Result: SexLunEclA_(LunEclFront) = 月食入陰陽曆(交前後度)= -5.809832570(-1) [LunStep 23]: 求月食分秒: 視去交前後度，不用南北東西差者。用減食限，不及減者不食。餘如定法而一，為月食之分秒。 LunFraction=(13.05-abs(SexLunEclA_))/8.7=0.832203153 [LunStep 24]: LEclM=6.93372527(3) LEclM_Used(定入遲疾曆)=6.93372527 LEclM_x12_2(定入遲疾限)=84.59144833 LEclM_Ch=82.40855167 LEclM_Sp=1.10120211 LEclM_SunSp(定限行度)=1.01920 LHalfEclT(定用分)=0.06483787(邢雲路用 4920,不是 5740) [LunStep 25]: LEclBeginT(初虧)=0.553290 LEclEndT(復圓)=0.682966 LEclBeginT24=13.27896554h LEclEndT24=16.39118343h ---------------------------------------------------------- 《大明熹宗哲皇帝實錄卷之七》原任陝西按察使邢雲路獻七政真數.....復以此法推天啟元年辛酉歲四月壬申朔日食: 初虧申正一刻(16.24h) 食甚酉初刻(17.00h) 復圓酉正初刻(18.00h) 日食一分八十六秒(0.186) 復圓日在天未入地,食不及三分不救(以上記錄缺乏其它數據可供證驗,是否可靠,甚有疑問?) 而授時則推是日食: 初虧申正三刻(16.72h) 食甚酉正初刻(18.00h 復圓戍初刻(19.00h) 日食三分九十一秒(0.391) 日未入已復光三分一十一秒,日已入未復光八十秒,與天不合(以上授時數據可按計算方法證驗) 《明熹宗悊皇帝實錄卷之十六》庚申原任陝西按察使閒住邢雲路奏: 臣本年春言四月朔日食歷數,臨期策驗分刻乃有不合臣始悟歷元玄法原自無差,而其差者,乃今大統歷沿郭守敬授時歷之誤也蓋日食者,月食之也,故正歷一準於月食,而月之分數與限度定法,皆其綱領最大者授時定月食分三十分及十五分,月食限十三度五分,定限八十七分此守敬率南北日官以儀器測驗于天,自丙子至庚辰又經甲午,再驗幾二十年,自謂推步之精,可永久無弊以故大統歷沿用之,臣今所步歷亦用之,及四月朔食測驗,方見其差,乃悟以異乘同除另變其法,變月食分三十分為三十五分(30->35) 變十五分為十七分五十秒變定限十三度五分為十三度十二分五十秒(13.05->13.1250) 變定法八十七分為七十五分(8.7->7.5) 諸數既變,乃依舊歷元玄法布歷,一一方合,然後知歷元法原自無差,而其差者,大統相沿授時法誤之也,故必如臣所推: 如五千七百四十為日食成限,以四千九百二十為月食成限,以三十五分為月食分,而黜郭守敬之三十分,方為正法,下禮部 LunFraction=(13.05 - abs(SexLunEclA_))/8.7 月食食分=定限(13.05->13.1250) - 視入陰陽曆/定法(8.7->7.5) LHalfEclT=(sqrt((30-LunFraction*10)*LunFraction*10)*49.20)/(LEclM_SunSp*10000) 定用分=(sqrt(三十分(30->35)-月食分)*月食分)*49.20)/(定限行度*10000) 19 1621(6) [2]MinSun=7.973680 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242160)124182.33(37.39)26.40 148.25(11.00) 8.19(12.50) 325.44( 8.18)20.05(26.41) (353.03) 142.81 39.81(3) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 299.83/15473.77(1.57)1(354.60) 6.27/ 6.23(76.52)0.04(0) (76.52) 1.00 5.40(5.40)1.0863 1.572/5.400/0.526(12.629h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 4.3710h(17.0004h) 5.00h 91.15m(18.5196h) 2168/7832 2831.96 2716.49 334.21/326.02(3)39.22 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 299.83/15473.77(1.55) 327.58(3) 37.66/3.70/(0.15) 144.96/2.92/3.17(2.67) 1/ 0.15/ 2.67 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 360.46/354.60/-5.86 -5.86(8)0.2680 20.58(0)82.94(1.10)(1.02) 384.15(17.5976h)(19.4415h) ---------------------------------------------- 19 1621( 7) LunEcl=16.2785 LunMoon=9.2362 LunFraction=-2.5819 LunMaxT24=2.2018h LEclNetSun_=11.8802 19 1621( 8) LunEcl=18.5968 LunMoon=11.2124 LunFraction=-6.0092 LunMaxT24=15.3730h LEclNetSun_=40.2466 19 1621( 9) LunEcl=20.9152 LunMoon=13.1886 LunFraction=-9.4863 LunMaxT24=7.3509h LEclNetSun_=69.1614 19 1621(10) LunEcl=23.2336 LunMoon=15.1648 LunFraction=-13.0180 LunMaxT24=0.4774h LEclNetSun_=98.6033 19 1621(11) LunEcl=25.5520 LunMoon=17.1411 LunFraction=-16.6120 LunMaxT24=17.4603h LEclNetSun_=128.5868 19 1621(12) LunEcl=0.6581 LunMoon=19.1173 LunFraction=0.6413 LunMaxT24=10.2391h LEclNetSun_=159.0944 19 1621(13) LunEcl=2.9765 LunMoon=21.0935 LunFraction=-3.0682 LunMaxT24=1.4936h LEclNetSun_=190.0187 19 1623(10) ReqEcl(交泛)=12.003066 LunEcl(食月望交泛)=26.7683625 《大明熹宗哲皇帝實錄卷之七》以新法推天啟三年癸亥歲九月十五壬寅夜望月食: 初虧申正三刻食甚酉初三刻(17.72?,該記錄不可靠) 復圓戍初初刻月食五分三十八秒月未出巳食五分二十五秒,月已出見食一十三秒,不及一分不救而授時則推月食: 初虧酉初一刻食甚酉正二刻(18.48h,該記錄可靠) 復圓戌正一刻月食七分八秒月未出已食一分九十八秒,月已出見食五分一十秒,與天不合 SumDay=124912.818 [LunStep 3]: LunDay(食月經望)=39.3346165 LunSun=100.81492 LunMoon=9.46173 Winter_solitice=47.878 [LunStep 4]: LunEclA(交常度)=357.86043247 LeapDay=18.6146460 [LunStep 5]: LunChSun(食月盈縮曆)=100.81491750 SumDay + Q=124967.878 [LunStep 6]: LUseSun=81.80616150(1) RefMoonDay=29.263390 [LunStep 7]: LSunX(盈縮差)=2.383266051 RefSun=155.985849 [LunStep 8]: LSunXS=-2.383266051 LunEclA_S(交定度)=355.477166424 RefMoon= 2.488554 [LunStep 9]: LUseMoon(食月望遲疾曆)= 9.4617325(3) Lmoon_x12_2(遲疾限)=115.4331365 LMoonCh=51.56686 RefEclipse=16.031600 [LunStep 10]:(+0.08306=>)LMoonCh_=52.6499235 LMoonX(遲疾差)=4.59614481 LMoonSp=1.0426297 [LunStep 11]:LMoonXS=4.59614481 LSunMoonX=-6.97941086 LSunMoonT(加減差)=-0.54891173 ReqMoonDay=86.049621 [LunStep 12]:LConjuncT(食月定望)=38.78570477 LConjuncTDay=38 LConjuncTime=0.78570477 ReqSun=86.049621 [LunStep 13]:LunSeeT(時差)=-0.0096072 LunMaxT(食甚定分)=0.77609756 => LunMaxT24=18.6263416h ReqMoon=22.250804 而授時則推...食甚酉正二刻(18.48h),該記錄可靠理論值:19.23h [LunStep 16]:LunEclSun=100.25640 LEcl_ChSun(食甚入盈縮曆)=100.25639857 LEclUseSun=82.36468043(1) [LunStep 17]:LEcl_SunX(食甚盈縮差)=2.38587969 [LunStep 18]:LEcl_SunX_=-2.38587969 LEclNetSun(食甚入盈縮曆定度)=97.87051887(1) LEclNetSun_= LEclNetSun if(LEclSunStatus > LagEnd){LEclNetSun_=LEclNetSun+182.62875} As LEclSunStatus=1 LEclNetSun_=105.89110287 total_year=calculated_year-Base_year=342 => Added_drift=0.015*total_year=5.13000 Ref_drift(天正加時赤道度)=10-Added_drift=4.87000(邢雲路計算方法) (授時曆計算方法:Ref_drift_check=fmod(SumDay+315.1075,365.2575)-305.1075=4.75304) Next_Ref_drift=Ref_drift_check-4.3445=0.40854 Sun_drift=Next_Ref_drift/1.0849=0.37657 Net_drift=4+Sun_drift=4.37657 LEclNetSun_All=LEclNetSun_ + Net_drift=110.26767 LEclNetSun_From(食甚月躔黃道奎宿)=LEclNetSun_All-103.6975(壁宿)=6.57017 [LunStep 22]: 求交常交定度:(refer to LunStep 4 and LunStep 8) 置經朔望入交泛日及分秒(LunEcl)，以月平行度(13.36878311)乘之，為交常度(LunEclA)；以盈縮差(LSunXS)盈加縮減之，為交定度(LunEclA_S)。求月食入陰陽曆去交前後度: 視交定度，在交中度(181.8967)已下，為陽曆；已上，減去交中，為陰曆。視入陰陽曆，在後準十五度半已下，為交後度；前準一百六十六度三十九分六十八秒已上，覆減交中，餘為交前度及分。交中度=181.8967 SexLunEclA=LunEclA_S=355.477166424 if (LunEclA_S > 181.8967){SexLunEclA=181.8967-LunEclA_S} Result:SexLunEclA(陰陽曆,負號[-]表示陰曆;正數表示陽曆)=-173.580466424 LunEclFront=-1(負數[-1]表示交前度); if(abs(SexLunEclA)<15.5){LunEclFront=1(正數[1]表示交後度);} SexLunEclA_=SexLunEclA=-173.580466424 if (LunEclA_S > 166.3968){SexLunEclA_=181.8967-LunEclA_S; LunEclFront=-1} if (SexLunEclA < -166.3968){SexLunEclA_=-181.8967-SextLunEclA; LunEclFront=-1} Result: SexLunEclA_(LunEclFront) = 月食入陰陽曆(交前後度)= -8.316233576(-1) [LunStep 23]: 求月食分秒: 視去交前後度，不用南北東西差者。用減食限，不及減者不食。餘如定法而一，為月食之分秒。 LunFraction=(13.05-abs(SexLunEclA_))/8.7=0.544111083 以新法推天啟三年癸亥歲九月十五壬寅夜望月食...月食五分三十八秒(0.538) 理論值:0.6 [LunStep 24]: LEclM=8.91282077(3) LEclM_Used(定入遲疾曆)=8.91282077 LEclM_x12_2(定入遲疾限)=108.73641341 LEclM_Ch=58.26358659 LEclM_Sp=1.05377118 LEclM_SunSp(定限行度)=0.97177 LHalfEclT(定用分)=0.05852605(邢雲路用 4920,不是 5740) [LunStep 25]: LEclBeginT(初虧)=0.717572 LEclEndT(復圓)=0.834624 LEclBeginT24=17.22171630h LEclEndT24=20.03096682h 以新法推天啟三年癸亥歲九月十五壬寅夜望月食[初虧申正三刻](16.72)不可靠, 因為這應該是授時的計算值而授時則推月食[初虧酉初一刻](17.24)不可靠, 因為與計算值不符合,很可能是新法推的初虧以新法推天啟三年..[初虧申正三刻](16.72h ?) .... 復圓[戍初初刻](19h ?) 理論值:17.7h 20.79h ---------------------------------------------------------- 19 1629(1) [2]MinSun=8.161296 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242152)127104.27(19.33)24.80 149.66(54.53)10.96(24.32) 179.19(24.06)12.93(26.64) (356.09) 179.19 3.43(1) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 495.99/1517.74(0.17)1(355.92) 12.93/12.88(157.79)0.06(3) ( 9.21)10.47 1.10(1.10)0.9916 -0.173/-1.101/-0.105(-2.529h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 1.3888h(22.8601h)10.86h 32.24m(23.3974h) 2069/7931 2068.94 4748.92 203.16/179.11(1) 3.52 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 495.99/1517.74(0.18) 178.93(1) 3.69/4.45/(-5.77) 178.93/0.35/0.67(0.04) 1/ 5.77/-0.04 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 363.37/355.92/-7.46 -7.46(8)0.0681 12.83(3)10.50(0.99)(0.91) 228.58(22.8488h)(23.9460h) ---------------------------------------------- 20 1629(7) [2]MinSun=8.161296 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242152)127104.27(19.33)24.80 149.66(54.53)10.96(24.32) 356.37(21.24)24.79(13.33) (178.26) 173.75 8.87(3) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 448.90/3753.66(0.41)-1(178.67) 11.01/10.99(134.37)0.02(1) (32.63) 8.18 3.30(3.29)1.1781 0.414/3.295/0.258(6.196h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 5.7942h(11.9906h)-0.01h -0.29m(11.9857h) 2073/7927 2927.01 5.98 377.87/356.63(3) 8.61 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 448.90/3753.66(0.40) 357.03(3) 8.21/4.42/(4.41) 174.41/0.77/0.00(0.00) -1/-4.41/ 0.00 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 183.64/178.67/-4.96 4.96(6)0.1726 25.05(1)29.48(1.18)(1.10) 293.15(11.2821h)(12.6892h) ---------------------------------------------- 20 1631(6) [2]MinSun=8.208200 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242150)127834.75(29.81)17.02 157.39(12.79) 5.25( 0.64) 334.58( 9.98)17.11(14.55) (194.52) 151.95 30.67(3) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 344.71/12549.90(1.28)-1(195.80) 3.33/ 3.28(40.68)0.05(0) (40.68) 7.14 3.84(3.84)1.0249 1.278/3.836/0.409(9.820h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 23.4157h(9.2359h)-2.76h -66.48m(8.1279h) 2130/7870 2870.34 1613.37 344.91/334.94(3)30.30 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 344.71/12549.90(1.27) 336.20(3) 29.04/4.01/(1.76) 153.58/2.38/1.54(1.54) -1/-1.76/ 1.54 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 187.83/195.80/ 7.96 -7.96(8)0.0044 17.52(0)45.67(1.03)(0.95) 56.68(7.9919h)(8.2639h) ---------------------------------------------- 21 1632(11) [2]MinSun=8.231652 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242149)128199.99(35.05)27.90 146.49( 7.16) 1.41( 1.25) 106.09(32.00)23.15(26.75) (357.62) 106.09 76.53(1) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 82.51/23442.53(2.35)1(355.27) 9.38/ 9.35(114.39)0.03(1) (52.61) 5.37 4.65(4.65)1.1500 -2.349/4.648/0.164(3.934h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 23.8977h(3.8321h)-8.17h -81.51m(2.4736h) 2383/7617 2383.35 3969.35 138.19/106.20(1)76.43 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 82.51/23442.53(2.35) 103.85(1) 78.77/1.14/(-0.76) 103.85/4.37/6.95(1.80) -1/ 0.76/ 1.80 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 360.20/355.27/-4.93 -4.93(8)0.3835 23.32(1)50.61(1.15)(1.07) 421.92(1.4610h)(3.4862h) ---------------------------------------------- 22 1634(4) [2]MinSun=8.278556 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242147)128930.48(45.54)20.11 154.23(25.42)23.27( 4.78) 272.35(23.55) 3.62(14.06) (187.92) 89.73 92.89(3) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 8.90/24001.62(2.40)-1(190.32) 3.62/ 3.61(44.12)0.01(2) (44.12) 6.66 4.08(4.08)1.1629 2.401/-4.076/-0.118(-2.834h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 13.1138h(10.2797h)-1.72h -46.05m(9.5121h) 2489/7511 2510.86 1036.62 295.75/272.20(3)93.04 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 2.98/24010.53(2.40) 274.60(3) 90.64/0.07/(0.04) 91.98/4.46/1.85(1.85) -1/-0.04/ 1.85 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 189.86/190.32/ 0.46 -0.46(8)0.9424 3.50(2)42.68(1.16)(1.08) 529.16(8.2421h)(10.7821h) ---------------------------------------------- 23 1636(8) [2]MinSun=8.325460 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242145)129660.96(56.02)12.33 161.96(43.69)17.56( 8.32) 32.97(39.93) 5.82(26.86) (359.15) 32.97 32.97(0) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 334.85/13234.41(1.36)1(357.79) 5.82/ 5.74(70.98)0.08(2) (70.98) 2.11 5.31(5.31)1.1174 -1.356/-5.308/-0.489(-11.737h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 22.4066h(10.6698h)-1.33h -36.96m(10.0537h) 2138/7862 2862.47 810.95 72.38/32.45(0)32.45 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 334.85/13234.41(1.34) 31.11(0) 31.11/3.94/(2.83) 31.11/2.52/0.82(0.82) -1/ 2.83/ 0.82 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 361.28/357.79/-3.49 -3.49(8)0.5639 5.33(2)65.00(1.13)(1.05) 493.47(8.8694h)(11.2380h) ---------------------------------------------- 25 1638(1) [2]MinSun=8.372364 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242143)130391.45( 6.51) 4.55 169.70( 1.95)11.86(11.85) 199.23(31.48)13.84(14.17) (189.45) 16.61 16.61(2) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 429.61/7570.66(0.78)-1(190.23) 0.06/ 0.00( 0.72)0.06(0) ( 0.72)11.04 0.08(0.08)0.9859 0.783/0.080/0.072(1.723h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 11.6227h(13.3460h) 1.35h 37.34m(13.9684h) 2087/7913 2086.71 820.16 230.81/199.33(2)16.70 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 429.61/7570.66(0.79) 200.11(2) 17.49/4.30/(2.61) 17.49/1.54/0.51(0.51) 1/ 2.61/-0.51 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 190.15/190.23/ 0.08 -0.08(8)0.9898 13.91(0) 1.59(0.99)(0.90) 634.62(12.4453h)(15.4915h) ---------------------------------------------- 26 1641(11) [2]MinSun=8.442720 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242140)131487.17(22.23) 7.64 166.53(14.59) 2.32(16.00) 126.13(39.42)24.06(14.29) (190.98) 126.13 56.49(1) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 205.65/20486.95(2.06)-1(188.92) 10.28/10.25(125.39)0.03(1) (41.61) 7.01 3.97(3.97)1.1664 -2.059/3.972/0.134(3.228h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 10.1316h(13.3594h) 1.36h 37.67m(13.9872h) 2258/7742 2258.16 827.98 165.71/126.29(1)56.33 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 205.65/20486.95(2.06) 124.23(1) 58.39/2.64/(1.67) 124.23/3.88/1.28(1.28) 1/ 1.67/ 1.28 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 191.00/188.92/-2.08 2.08(6)0.6526 24.19(1)40.00(1.17)(1.09) 495.34(12.7983h)(15.1760h) ---------------------------------------------- 27 1643(4) [2]MinSun=8.489624 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242138)132217.65(32.71)29.39 144.74( 3.32)22.19(17.21) 262.86( 1.44) 2.54(26.49) (354.08) 80.24 80.24(2) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 56.99/23734.40(2.37)1(356.45) 2.54/ 2.54(31.02)0.00(2) (31.02) 8.37 3.08(3.08)1.1800 2.375/-3.079/-0.049(-1.174h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 10.6136h(9.4396h)-2.56h -62.94m(8.3906h) 2410/7590 2409.68 1503.91 264.21/262.77(2)80.15 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 56.99/23734.40(2.37) 265.14(2) 82.52/0.82/(0.31) 82.52/4.42/2.66(2.66) -1/-0.31/-2.66 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 354.68/356.45/ 1.78 1.78(6)0.7035 2.49(2)30.42(1.18)(1.10) 498.95(7.1931h)(9.5881h) ---------------------------------------------- 27 1644(3) [2]MinSun=8.513076 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242137)132582.90(37.96)10.74 163.37(27.22)20.33(20.14) 251.96(55.81)26.26(27.09) (362.21) 69.34 69.34(2) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 126.46/22688.64(2.27)1(364.48) 12.48/12.47(152.26)0.01(1) (14.74)10.03 1.67(1.67)1.1965 2.273/1.667/0.270(6.480h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 19.4244h(1.9045h)-10.10h -50.07m(1.0700h) 2338/7662 2338.01 4554.18 308.01/252.20(2)69.58 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 126.46/22688.64(2.28) 254.47(2) 71.85/1.70/(-1.61) 71.85/4.26/7.75(0.76) -1/ 1.61/-0.76 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 358.49/364.48/ 5.99 5.99(6)0.0014 26.53(1)11.45(1.20)(1.12) 27.23(1.0046h)(1.1353h) ---------------------------------------------- 28 1644(9) [2]MinSun=8.513076 [1](365.)SumDay(WinSol)LeapDay [2]RefS(RefMDay)RefM(RefEcl) [3]ReqS(ReqMDay)ReqM(ReqEcl) [4](EclA)[5]ChSun[6]UseSun(S) (0.242137)132582.90(37.96)10.74 163.37(27.22)20.33(20.14) 63.90(52.99)10.56(13.79) (184.38) 63.90 63.90(0) [7]Sun_Add/Sum(X)[8]Ft(EclA_S) [9]UseM/SoD(M12_2)MSo(M) (MoonCh)MAdd[10]Sum(X)Speed [11]SunXS/MoonXS/SunMoonT 173.73/21238.16(2.14)-1(182.24) 10.56/10.50(128.83)0.06(3) (38.17) 7.48 3.73(3.72)1.0215 -2.140/-3.725/-0.471(-11.297h) [12]Conjunc0(T) [13]TNoon SeeT(MaxPhase) [14]Srise/Sset HalfDayFan[15]NetTNoon [16]Ecl)Ref/Sun(S)UseSun 23.8298h(12.5325h) 0.53h 15.90m(12.7976h) 2300/7700 2699.72 332.32 116.44/63.44(0)63.44 [17]Ecl_)SAdd/SSum(SunX) [18]NSun(S)[19]NS_R/NS_X(NS_) [20]EW_R/EW_X/EW_H(EW_) [21]A12_Sgn/NS_Sh/EW_Sh 173.73/21238.16(2.13) 61.31(0) 61.31/2.45/(2.15) 61.31/3.98/0.53(0.53) 1/-2.15/ 0.53 [22]EclRefA/NetEclA_S/PosA [23]EclPos(Ref)Frac [24]EclM(M)_Ch(_Sp)(_SunSp) [25]HalfEclT(EclBegin)(EclEnd) 186.43/182.24/-4.19 4.19(6)0.3014 10.09(3)43.91(1.03)(0.95) 433.50(11.7572h)(13.8380h) ---------------------------------------------- Reference: (1)Accuracy of Solar Eclipse Calculations, 1644-1785, with the Chinese Shoushi Calandar(Table 2) http://articles.adsabs.harvard.edu//full/1999JHA....30..161L/0000161.000.html (2)Accuracy of Solar Eclipse Observations made by Jesuit Astronomers in China http://adsabs.harvard.edu/full/1995JHA....26..227S Observation data from Table 2 Observation error data from Table 3 Theorectical=Observation-Observation error Shoushi(Observation) Remark : 999 or 9999 = no data -)Theorectical(Obsevation-Shoushi) -------------- Shoushi-Theorectical(Observation-Theorectical) ================================================================================================================== N year(month) EclipseFraction EclipseBegin MaxPhase EclipseEnd month/day ================================================================================================================== 0 1568(10) 0.0221(0.110) 12.1030(11.137) 12.3605(11.737) 12.6180(12.343) 9/21 -) 0.110(0.0878) 11.137(-0.9663) 11.737(-0.6239) 12.343(-0.2752) ---------------------------------------------------------------------------------------- -0.0878(0.000) 0.9663(0.000) 0.6239(0.000) 0.2752(0.000) 1 1572(8) 0.9668(0.916) 6.4686(6.660) 7.9429(7.930) 9.4172(9.930) 6/10 -) 0.916(-0.0508) 6.660(0.1914) 7.930(-0.0129) 9.930(0.5128) ---------------------------------------------------------------------------------------- 0.0508(0.000) -0.1914(0.000) 0.0129(0.000) -0.5128(0.000) 2 1575(6) 0.4650(0.701) 13.5750(13.050) 14.6262(14.290) 15.6774(15.450) 5/10 -) 0.701(0.2360) 13.050(-0.5250) 14.290(-0.3362) 15.450(-0.2274) ---------------------------------------------------------------------------------------- -0.2360(0.000) 0.5250(0.000) 0.3362(0.000) 0.2274(0.000) 3 1582(7) 0.9759(0.723) 12.7915(12.206) 14.1250(13.618) 15.4585(14.915) 6/20 -) 0.723(-0.2528) 12.206(-0.5860) 13.618(-0.5071) 14.915(-0.5431) ---------------------------------------------------------------------------------------- 0.2528(0.000) 0.5860(0.000) 0.5071(0.000) 0.5431(0.000) 4 1583(13) 0.7154(0.825) 11.3255(11.230) 12.7799(13.080) 14.2343(14.770) 12/14 -) 0.825(0.1096) 11.230(-0.0955) 13.080(0.3001) 14.770(0.5357) ---------------------------------------------------------------------------------------- -0.1096(0.000) 0.0955(0.000) -0.3001(0.000) -0.5357(0.000) 5 1587(10) 0.8537(0.736) 11.9223(11.322) 13.1681(12.644) 14.4139(13.934) 10/2 -) 0.736(-0.1179) 11.322(-0.6000) 12.644(-0.5245) 13.934(-0.4796) ---------------------------------------------------------------------------------------- 0.1179(0.000) 0.6000(0.000) 0.5245(0.000) 0.4796(0.000) 6 1589(2) 0.5157(0.301) 6.1622(6.450) 7.2658(7.196) 8.3695(7.999) 2/15 -) 0.301(-0.2144) 6.450(0.2877) 7.196(-0.0697) 7.999(-0.3703) ---------------------------------------------------------------------------------------- 0.2144(0.000) -0.2877(0.000) 0.0697(0.000) 0.3703(0.000) 7 1590(8) 0.2223(0.431) 16.1850(16.080) 17.1090(17.106) 18.0330(18.037) 7/31 -) 0.431(0.2086) 16.080(-0.1047) 17.106(-0.0031) 18.037(0.0042) ---------------------------------------------------------------------------------------- -0.2086(0.000) 0.1047(0.000) 0.0031(0.000) -0.0042(0.000) 8 1594(6) 0.5152(0.319) 10.1823(9.830) 11.3417(10.840) 12.5011(11.900) 5/20 -) 0.319(-0.1962) 9.830(-0.3523) 10.840(-0.5017) 11.900(-0.6011) ---------------------------------------------------------------------------------------- 0.1962(0.000) 0.3523(0.000) 0.5017(0.000) 0.6011(0.000) 9 1596(10) 0.7084(0.949) 10.6131(9.870) 11.8365(11.170) 13.0598(12.500) 9/22 -) 0.949(0.2406) 9.870(-0.7431) 11.170(-0.6665) 12.500(-0.5598) ---------------------------------------------------------------------------------------- -0.2406(0.000) 0.7431(0.000) 0.6665(0.000) 0.5598(0.000) 10 1603(5) 0.7343(0.903) 7.5081(7.400) 8.8188(8.590) 10.1294(9.930) 5/11 -) 0.903(0.1687) 7.400(-0.1081) 8.590(-0.2288) 9.930(-0.1994) ---------------------------------------------------------------------------------------- -0.1687(0.000) 0.1081(0.000) 0.2288(0.000) 0.1994(0.000) 11 1604(5) 0.6227(0.078) 15.5666(15.883) 16.9268(16.408) 18.2871(16.907) 4/29 -) 0.078(-0.5446) 15.883(0.3161) 16.408(-0.5189) 16.907(-1.3801) ---------------------------------------------------------------------------------------- 0.5446(0.000) -0.3161(0.000) 0.5189(0.000) 1.3801(0.000) 12 1607(3) 0.4244(999.000) 17.4368(17.320) 18.4662(18.260) 19.4956(999.000) 2/26 -) 999.000(998.5756) 17.320(-0.1168) 18.260(-0.2062) 999.000(979.5044) ---------------------------------------------------------------------------------------- -998.5756(0.000) 0.1168(0.000) 0.2062(0.000) -979.5044(0.000) 13 1609(0) -99.0000(0.044) 0.0000(16.197) 0.0000(16.594) 0.0000(16.975) 12/26 -) 0.044(99.0444) 16.197(16.1970) 16.594(16.5940) 16.975(16.9753) ---------------------------------------------------------------------------------------- -99.0444(0.000) -16.1970(0.000) -16.5940(0.000) -16.9753(0.000) 14 1610(13) 0.9588(0.607) 14.1179(14.610) 15.6318(16.070) 17.1457(17.320) 12/15 -) 0.607(-0.3518) 14.610(0.4921) 16.070(0.4382) 17.320(0.1743) ---------------------------------------------------------------------------------------- 0.3518(0.000) -0.4921(0.000) -0.4382(0.000) -0.1743(0.000) 15 1612(0) -99.0000(999.000) 0.0000(999.000) 0.0000(999.000) 0.0000(999.000) 12/99 -) 999.000(1098.0000) 999.000(999.0000) 999.000(999.0000) 999.000(999.0000) ---------------------------------------------------------------------------------------- -1098.0000(0.000) -999.0000(0.000) -999.0000(0.000) -999.0000(0.000) 16 1615(4) 0.6471(0.890) 15.0718(15.162) 16.3547(16.412) 17.6376(17.546) 3/29 -) 0.890(0.2428) 15.162(0.0906) 16.412(0.0571) 17.546(-0.0913) ---------------------------------------------------------------------------------------- -0.2428(0.000) -0.0906(0.000) -0.0571(0.000) 0.0913(0.000) 17 1617(8) 0.3689(0.470) 19.4230(19.070) 20.5355(19.904) 21.6481(20.691) 8/1 -) 0.470(0.1011) 19.070(-0.3530) 19.904(-0.6314) 20.691(-0.9572) ---------------------------------------------------------------------------------------- -0.1011(0.000) 0.3530(0.000) 0.6314(0.000) 0.9572(0.000) 18 1621(6) 0.2680(0.800) 17.5976(17.240) 18.5196(9999.000) 19.4415(9999.000) 5/21 -) 0.711(0.5320) 17.270(-0.3576) 18.230(9980.4804) 19.190(9979.5585) ---------------------------------------------------------------------------------------- -0.4430(0.089) 0.3276(-0.030) 0.2896(9980.770) 0.2515(9979.810) 19 1629(7) 0.1726(0.166) 11.2821(10.980) 11.9857(11.710) 12.6892(12.450) 6/21 -) 0.166(-0.0066) 10.980(-0.3021) 11.710(-0.2757) 12.450(-0.2392) ---------------------------------------------------------------------------------------- 0.0066(0.000) 0.3021(0.000) 0.2757(0.000) 0.2392(0.000) 20 1631(0) -99.0000(0.108) 0.0000(12.340) 0.0000(12.910) 0.0000(13.480) 10/25 -) 0.108(99.1080) 12.340(12.3400) 12.910(12.9100) 13.480(13.4800) ---------------------------------------------------------------------------------------- -99.1080(0.000) -12.3400(0.000) -12.9100(0.000) -13.4800(0.000) 21 1634(4) 0.9424(0.820) 8.2421(8.770) 9.5121(9.944) 10.7821(11.190) 3/29 -) 0.820(-0.1228) 8.770(0.5281) 9.944(0.4315) 11.190(0.4078) ---------------------------------------------------------------------------------------- 0.1228(0.000) -0.5281(0.000) -0.4315(0.000) -0.4078(0.000) 22 1636(8) 0.5639(0.633) 8.8694(8.380) 10.0537(9.627) 11.2380(10.964) 8/1 -) 0.633(0.0689) 8.380(-0.4898) 9.627(-0.4272) 10.964(-0.2737) ---------------------------------------------------------------------------------------- -0.0689(0.000) 0.4898(0.000) 0.4272(0.000) 0.2737(0.000) 23 1637(0) -99.0000(0.199) 0.0000(12.710) 0.0000(13.680) 0.0000(14.580) 1/26 -) 0.199(99.1990) 12.710(12.7100) 13.680(13.6800) 14.580(14.5800) ---------------------------------------------------------------------------------------- -99.1990(0.000) -12.7100(0.000) -13.6800(0.000) -14.5800(0.000) 24 1638(1) 0.9898(0.656) 12.4453(11.650) 13.9684(13.320) 15.4915(14.840) 1/15 -) 0.656(-0.3338) 11.650(-0.7953) 13.320(-0.6484) 14.840(-0.6515) ---------------------------------------------------------------------------------------- 0.3338(0.000) 0.7953(0.000) 0.6484(0.000) 0.6515(0.000) 25 1641(11) 0.6526(0.827) 12.7983(12.650) 13.9872(14.040) 15.1760(15.330) 11/3 -) 0.827(0.1744) 12.650(-0.1483) 14.040(0.0528) 15.330(0.1540) ---------------------------------------------------------------------------------------- -0.1744(0.000) 0.1483(0.000) -0.0528(0.000) -0.1540(0.000) 26 1643(4) 0.7035(0.530) 7.1931(7.830) 8.3906(8.850) 9.5881(9.940) 3/20 -) 0.530(-0.1735) 7.830(0.6369) 8.850(0.4594) 9.940(0.3519) ---------------------------------------------------------------------------------------- 0.1735(0.000) -0.6369(0.000) -0.4594(0.000) -0.3519(0.000) 27 1644(9) 0.3014(0.303) 11.7572(10.930) 12.7976(12.140) 13.8380(13.340) 9/1 -) 0.303(0.0016) 10.930(-0.8272) 12.140(-0.6576) 13.340(-0.4980) ---------------------------------------------------------------------------------------- -0.0016(0.000) 0.8272(0.000) 0.6576(0.000) 0.4980(0.000)