303です
いつもの人サンクス

予想より値が高く出たな
近似の精度をもっと上げないと駄目か

立式は以下を参考にしました
https://eco.mtk.nao.ac.jp/koyomi/wiki/C3E6BDA9A4CECCBEB7EEA4C8A4CF2FCCBEB7EEC9ACA4BAA4B7A4E2CBFEB7EEA4CAA4E9A4BA.html