MCMCのスクリプトをEri RYで試してみた

皆既食のあるEA、Eri RYで先日のMCMC法のスクリプトを試してみた。

なんどか繰り返して得られたパラメータは、

p　4.97930685　      epo 3.00766085e+03  mag0 10.0994833  d0 5.68470602e-02

mag1 10.9285666    d2  6.22827045e-02      mag2 10.3120917  d1 2.29853353e-02

したのphase plotと見比べてください。かなり良い線を行っている感じです。

皆既が4.6%(d1*2)ということになります。ちなみに,PDMで得られた周期は4.97929でした。

paramsdを参考のために入れておきます。

[1] 6.93263376e-06 1.32795030e-03 2.24997437e-03 6.53248992e-04 1.92241962e-02
[6] 2.54509605e-03 9.81185897e-03 1.00539617e-03