位相計算

Rで、位相で整理して光度曲線を描く。

> jdm <- read.table("ancamapr.txt")
> p <-20.9986
> epoch <- 2426003.4700
> jdm[,1] <- *1, xlab="phase", ylab="mag.")

txtファイルの中身は、JDと光度。
食外、変動している？

*1:jdm[,1]-epoch)%%p)/p > plot(jdm, xlim=c(0,1), ylim=rev(c(9.5,10.5