AV Pup

またしてもコンパイラのバグでした。0.52+0.01が0.52になっていました。加算してくれない。
犯人は(例によって)16bit時代のms-basicのコンパイラです。インタープリターも同じ動作をしました。
Contactの方も同じコーディングなのですが、こちらでは問題が発生しません。気まぐれです。
あちこちに呪いの呪文を入れてなんとかしました。

で、再開！

それにしても、(今は単にO-Cだけを見てますが)徐々にD,SD,Cの値が離れつつあります。妙です。
こんな感じに

inclination(deg.)
   D 83
  SD 81
   C 81

mass ratio
   D 0.90
  SD 0.83
   C 0.88

Contactだと思いますが（あの論文の）q=0.8にはなりそうにありません。なぜ、あの論文（のAbstractには）は小数点一桁までしか書いてないのだろう。
解析が終わったら著者にメール出してみようかな～