SZ Psc

昼間、Mayhillで、SZ Pscを撮りました。
AVEで、まばらすぎて(?)極小時刻が決められなかったので、もう、2次関数fitでいいやと思って、Rで計算したのですが、やりかたがわからなくて苦労しました。
どうにか、計算できました。
もっと、簡潔に書けそうですが、いまはこれで精一杯。

V1 V2
1 2455827.6037 7.51
2 2455827.6597 7.62
3 2455827.7177 7.71
4 2455827.8007 7.68
5 2455827.8147 7.66
6 2455827.8277 7.64
7 2455827.8417 7.61
8 2455827.8927 7.51

V1 V2
1 0.60369999986 7.51
2 0.65970000019 7.62
3 0.71769999992 7.71
4 0.80070000002 7.68
5 0.81469999999 7.66
6 0.82769999979 7.64
7 0.84170000022 7.61
8 0.89270000020 7.51

Nonlinear regression model
model: y ~ a + b * x + c * x^2
data: parent.frame()
a b c
2.3920011 14.2138749 -9.5078879
residual sum-of-squares: 0.00048151991

Number of iterations to convergence: 1
Achieved convergence tolerance: 1.8364413e-06

[1] 2455827.7474
>

(追記)
最初、日進補正しないデータで解析してしまって、やり直したのですが、Rだと、やり直しは簡単です。
で、GCVSだと、O-C、大きすぎ

O=2455827.7474
C=2455828.57245
E=5065 (0.792)
O-C=-0.82505

Krakowだと以下のようになります。ほぼあってますが、予報より少し早いです。
このため、丁度極小部分が欠ける観測結果に。

O=2455827.7474
C=2455827.817266
E=839 (0.982)
O-C=-0.06987