!! 月齢
旧暦の暦がつくりたくなったので
（六曜をしりたくなったので)
月齢の計算式がきになったので調べてみた

まんまのGITがここに
*https://github.com/fgshun/qreki_py
 $ pip3.9 install  git+https://github.com/fgshun/qreki_py.git@v0.5.1#egg=qreki

http://koyomi8.com/reki_doc/doc_0250.htm

! ユリウス暦で調べる
 JDE= 2451550.09766+29.530588561*K + 0.00015437*T^2 - 0.000000150*T^3 + 0.00000000073*T^4
 らしい
 Ｔは、J2000.0 からのユリウス世紀
 J2000.0 = JD2451545.0 = 2000年1月1日12:00の事

 *jy は ユリウス年
 jc = jy /100;
 jd = jy * 365.25;
 moon_ege = ( jd + 2451544.5 - 
             ( 2451550.09766 + 
               0.00015437    * jc * jc - 
               0.000000150   * jc * jc * jc +
               0.00000000073 * jc * jc * jc * jc)
             ) % 29.530588561;



! Key words は 

海上保安庁水路部 略算式

*https://github.com/komasaru/CalendarPy/tree/master/ephemeris_jcg_i
*https://www1.kaiho.mlit.go.jp/KOHO/
*http://www.kotenmon.com/cal/jpl_DE.htm

日の出・日の入りの計算―天体の出没時刻の求め方
*https://amzn.to/39vzwnW

! ユリウス日
https://eco.mtk.nao.ac.jp/koyomi/wiki/A5E6A5EAA5A6A5B9C6FC.html

! Maxima で
                  ^
 ## f(x):= ( x + 2451544.5 - ( 2451550.09766 + 0.00015437*(x/36525)*(x/36525) - 0.000000150* (x/36525)*(x/36525)*(x/36525) + 0.00000000073 * (x/36525)*(x/36525)*(x/36525) *(x/36525) ));
 (%i1) f(x):= ( x + 2451544.5 - ( 2451550.09766 + 0.00015437*(x/36525)*(x/36525) - ( 0.000000150 * (x/36525)*(x/36525)*(x/36525)) + (0.00000000073*(x/36525)*(x/36525)*(x/36525)*(x/36525))));
 
 (%i11) solve(f(x)=0,x);