読者です 読者をやめる 読者になる 読者になる

muttenz's blog

スイス星空だより

EASolverに図を

EASolverで結果を出しても、数字だけだとピンときませんね。

そこで周期とその確かさを表すようにしてみました。

plot(d[m]/(N1+(1:(k+1))-1),1/dev[1:(k+1)],pch=16,col=3,grid(col=4),xlab="P")

こういう行を、繰り返しの後にでも入れます。

複雑に見えますが、それぞれの周期とその周期でのdevの逆数をグラフにするだけです。逆数を取ることによって小さな値の違いがはっきりグラフに出てきます。

これはOQ Gemの場合です。

f:id:muttenz:20170128044902p:plain

P: 2.76....がはっきりと他の候補より上というのがわかります。

 

これはCW Oriの場合です。

f:id:muttenz:20170128045053p:plain

上位2つの候補が比較的近いのがわかりますね。2番めの方が実際の周期でした。