🚀 go-pugleaf

  金子＠仙台電波です。

    Followup-To: fj.sci.astro

してあります。

<01bd40e0$0a8e6cc0$LocalHost@mebius>の記事において
sugayan@interlink.or.jpさんは書きました。
>> 初めて投書します、菅原と申します。

>> １）日付より六曜の求め方

  日付をいわゆる旧暦の日付に変換できれば一意に決まります。
  で、旧暦ですが、まずは二十四節季(春分とか、雨水とか)と、新
月になる日時を求めて、それから計算することになります。
  二十四節季や月の動きは、9惑星を考慮した運動方程式を解いて
求める(正しくは推定するだけ)しかないんですが、解法自体はいく
つかの本に掲載されています(理解するの大変だけど)。

  でまじめにやりすぎると大変なので

    「暦の便利帳」(出版社等は、家に帰ればわかるんだけど)

という本をおすすめします。B5判(だと思った)の青い本です。BASI
Cでのサンプルプログラムが載っているので多少は役に立つでしょ
う。
  過去の旧暦日付については、一覧表を手に入れるしかないでしょ
う。現在の計算式とは異なる式で計算しているので狂ってくる可能
性がある。なにかいい文献知りませんかね(便乗質問)。占い系の本
に換算表が載っていたりするんだけど、過去数十年分しかないんで
すよ。

>> ２）春分の日、秋分の日の求め方

  厳密には、国立天文台だったかの計算に基づいて、前年の官報で
公表されるまでは日付は確定しません。でも、それでは大変なので、
太陽その他がいまのまま運動しつづけたと仮定して計算した結果か
ら導出された近似式が上記の文献に書いてあります。

  ま、適度に近似式を求めて、それでお茶を濁した方がいいですよ。
まじめにやろうとすると、奥が深すぎます。

double x,y,z=1.8;main(a,b){x=(time( /*       金子 裕(Yutaka Kaneko)          */
0)%935)/472;while(!0){y=x;x=2-z*(x- /* Sendai National College of Technology */
1)*(x-1);a=(x*631);a%=78;b=(y*95);b /*     kaneko@cc.sendai-ct.ac.jp         */
%=24;printf("\033[%d;%dHx\n",b,a);}}/*Evanx08,Megax05,Escax07,Fushx16,WedPx07*/