muttenz's blog

スイス星空だより

2015-01-26から1日間の記事一覧

EAの副極小のフェーズもパラメータにしてみた

下のスクリプトのようにやってみた。p[8]は副極小が0.5からどれだけずれているかというパラメータ。 fphの代わりになんとかph <- {(x-p[2])/p[1]}%%1だけで出来ないかと abs{ph-round(ph)}とやってみた。しかしこれで計算させたらどうしても数値がおかしくな…