muttenz's blog

スイス星空だより

2018-03-14から1日間の記事一覧

O - C

昨日、 T(N) = a*N^3 + b*N^2 + c*N + d (1) の係数を、データからRのnls関数で直接求めようとしたら、エラーで出来なかったと書きました。 そこで、P0を元期での周期、 E0を元期のHJDとして、この式を少し変形します。 T(N) = a*N^3 + b*N^2 + (c' + P0) * …