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の場合です。

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

これはCW Oriの場合です。

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